N := 1;
X := u^2 + 1;
Xi := i;
q := 1/v;
t := u/v;
E := [1,(q^2-1)*(t^2-1)/16,(q^2-1)*(t^2-1)/16,0,0];
P := [(q+1)*(t^2-1)/8,(q+1)^2*(t-1)^2*(t+1)/32];
Q := [0,0];