N := 4;
X := (u^8 + 8*u^7 + 27*u^6 + 50*u^5 + 56*u^4 + 40*u^3 + 19*u^2 + 6*u + 1)*v^6 + (u^10 + 8*u^9 + 33*u^8 + 82*u^7 + 115*u^6 + 74*u^5 - 9*u^4 - 48*u^3 - 35*u^2 - 14*u - 3)*v^5 + (u^12 + 8*u^11 + 34*u^10 + 90*u^9 + 164*u^8 + 184*u^7 + 86*u^6 - 30*u^5 - 31*u^4 + 20*u^3 + 29*u^2 + 14*u + 4)*v^4 - (u^12 + 6*u^11 + 29*u^10 + 90*u^9 + 208*u^8 + 358*u^7 + 401*u^6 + 262*u^5 + 108*u^4 + 50*u^3 + 29*u^2 + 10*u + 3)*v^3 + (u^12 + 4*u^11 + 12*u^10 + 30*u^9 + 59*u^8 + 120*u^7 + 198*u^6 + 190*u^5 + 104*u^4 + 40*u^3 + 17*u^2 + 4*u + 1)*v^2 + (u^10 + 2*u^9 + 6*u^8 + 14*u^7 + 3*u^6 - 14*u^5 - 12*u^4 - 2*u^3 - u^2)*v + u^8 - u^6 + u^4;
Xz := u^6*v^2 + z*u^6*v + 4*u^5*v^2 + 2*z*u^5*v - z*u^4*v^3 + (-2*z + 9)*u^4*v^2 + (4*z - 2)*u^4*v - u^4 - 4*z*u^3*v^3 + (-2*z + 10)*u^3*v^2 + (6*z - 6)*u^3*v + (-6*z - 1)*u^2*v^3 + (2*z + 6)*u^2*v^2 + (5*z - 4)*u^2*v - z*u^2 + (-4*z - 2)*u*v^3 + (4*z + 4)*u*v^2 - 2*u*v + (-z - 1)*v^3 + (z + 2)*v^2 - v;
q := u+1;
t := (v*u) / (u+1);
E := [(1+q)*t+(2-q),0,q*t^2+(1-q)*t,0,0];
P := [-t,t^2];
Q := [0,0];