N := 2;
X := (u^4 + 6*u^2 + 1)*v^4 - 16*u^2*v^3 - (4*u^4 - 20*u^2)*v^2 + (8*u^4 - 8*u^2)*v + u^6 - 2*u^4 + u^2;
Xi := u^3 + i*u^2*v^2 - 2*u*v^2 + 4*u*v - u + i*v^2;
q := (u^2+2*v-1) / (2*v);
t := (u^2+2*v-1) / (2*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];