N := 4;
X := (u^16 + 8*u^15 + 28*u^14 + 56*u^13 + 70*u^12 + 56*u^11 + 28*u^10 + 8*u^9 + u^8)*v^16 - (32*u^15 + 224*u^14 + 672*u^13 + 1120*u^12 + 1120*u^11 + 672*u^10 + 224*u^9 + 32*u^8)*v^15 - (16*u^16 + 96*u^15 - 224*u^14 - 2464*u^13 - 6720*u^12 - 9184*u^11 - 6944*u^10 - 2784*u^9 - 464*u^8)*v^14 + (448*u^15 + 2240*u^14 + 448*u^13 - 15680*u^12 - 38080*u^11 - 39872*u^10 - 20160*u^9 - 4032*u^8)*v^13 + (u^22 + 4*u^21 + 6*u^20 + 4*u^19 + u^18 + 112*u^16 + 448*u^15 - 4930*u^14 - 21960*u^13 - 10092*u^12 + 71224*u^11 + 134846*u^10 + 93632*u^9 + 23408*u^8 + u^6 + 4*u^5 + 6*u^4 + 4*u^3 + u^2)*v^12 - (24*u^21 + 72*u^20 + 72*u^19 + 24*u^18 + 2688*u^15 + 8064*u^14 - 33200*u^13 - 121104*u^12 - 27920*u^11 + 246352*u^10 + 287616*u^9 + 95872*u^8 + 24*u^5 + 72*u^4 + 72*u^3 + 24*u^2)*v^11 - (12*u^22 + 24*u^21 - 240*u^20 - 504*u^19 - 252*u^18 + 448*u^16 + 896*u^15 - 27800*u^14 - 56496*u^13 + 170720*u^12 + 397936*u^11 - 86408*u^10 - 570752*u^9 - 285376*u^8 + 12*u^6 + 24*u^5 - 240*u^4 - 504*u^3 - 252*u^2)*v^10 + (240*u^21 + 240*u^20 - 1520*u^19 - 1520*u^18 + 8960*u^15 + 8960*u^14 - 170720*u^13 - 170720*u^12 + 660704*u^11 + 660704*u^10 - 627456*u^9 - 627456*u^8 + 240*u^5 + 240*u^4 - 1520*u^3 - 1520*u^2)*v^9 + (56*u^22 - 2042*u^20 + 5832*u^18 + 1128*u^16 - 76288*u^14 + 655908*u^12 - 1550080*u^10 + 1027176*u^8 + 72*u^6 - 2042*u^4 + 5816*u^2)*v^8 - (896*u^21 - 896*u^20 - 8736*u^19 + 8736*u^18 + 6240*u^17 - 6240*u^16 + 24288*u^15 - 24288*u^14 - 336160*u^13 + 336160*u^12 + 1338656*u^11 - 1338656*u^10 - 1261280*u^9 + 1261280*u^8 - 6240*u^7 + 6240*u^6 - 5088*u^5 + 5088*u^4 - 14720*u^3 + 14720*u^2)*v^7 - (128*u^22 - 256*u^21 - 5456*u^20 + 11168*u^19 + 10784*u^18 - 32736*u^17 + 29920*u^16 - 27104*u^15 - 93408*u^14 + 213920*u^13 + 700224*u^12 - 1614368*u^11 - 367136*u^10 + 2348640*u^9 - 1201632*u^8 + 54624*u^7 - 49696*u^6 + 44768*u^5 - 47344*u^4 + 49920*u^3 - 24960*u^2)*v^6 + (1536*u^21 - 4608*u^20 - 10688*u^19 + 44352*u^18 - 45760*u^17 + 14912*u^16 + 36672*u^15 - 108992*u^14 - 199616*u^13 + 889152*u^12 - 95040*u^11 - 2182720*u^10 + 2538432*u^9 - 972096*u^8 + 177088*u^7 - 153408*u^6 + 134336*u^5 - 119872*u^4 + 84480*u^3 - 28160*u^2)*v^5 + (148*u^22 - 592*u^21 - 5480*u^20 + 24880*u^19 - 26108*u^18 - 22336*u^17 + 92704*u^16 - 157248*u^15 + 118632*u^14 + 120480*u^13 + 112016*u^12 - 1488224*u^11 + 2557288*u^10 - 1868352*u^9 + 683040*u^8 - 262976*u^7 + 227332*u^6 - 195280*u^5 + 151960*u^4 - 82512*u^3 + 20628*u^2)*v^4 - (1184*u^21 - 5920*u^20 + 3040*u^19 + 32160*u^18 - 94592*u^17 + 149376*u^16 - 166528*u^15 + 120960*u^14 - 184640*u^13 + 726592*u^12 - 1478080*u^11 + 1532608*u^10 - 861056*u^9 + 311680*u^8 - 171136*u^7 + 164480*u^6 - 147040*u^5 + 104416*u^4 - 46880*u^3 + 9376*u^2)*v^3 - (80*u^22 - 480*u^21 - 1344*u^20 + 13664*u^19 - 38352*u^18 + 58752*u^17 - 56576*u^16 + 28800*u^15 - 17120*u^14 + 137408*u^13 - 405888*u^12 + 565312*u^11 - 416672*u^10 + 150912*u^9 - 22784*u^8 + 24704*u^7 - 49776*u^6 + 54816*u^5 - 37696*u^4 + 14688*u^3 - 2448*u^2)*v^2 + (320*u^21 - 2240*u^20 + 6336*u^19 - 8512*u^18 + 1792*u^17 + 16128*u^16 - 36608*u^15 + 35072*u^14 + 8064*u^13 - 60544*u^12 + 60544*u^11 - 8064*u^10 - 35072*u^9 + 36608*u^8 - 16128*u^7 - 1792*u^6 + 8512*u^5 - 6336*u^4 + 2240*u^3 - 320*u^2)*v + u^24 - 8*u^23 + 52*u^22 - 248*u^21 + 706*u^20 - 1112*u^19 + 612*u^18 + 1368*u^17 - 4113*u^16 + 5168*u^15 - 1688*u^14 - 5168*u^13 + 8860*u^12 - 5168*u^11 - 1688*u^10 + 5168*u^9 - 4113*u^8 + 1368*u^7 + 612*u^6 - 1112*u^5 + 706*u^4 - 248*u^3 + 52*u^2 - 8*u + 1;
Xi := u^12 + i*u^11*v^6 - 6*i*u^11*v^4 + 10*i*u^11*v^2 + (-4*i - 4)*u^11 + 2*i*u^10*v^6 - 12*i*u^10*v^5 + 48*i*u^10*v^3 - 20*i*u^10*v^2 - 40*i*u^10*v + (16*i + 10)*u^10 + i*u^9*v^6 - 12*i*u^9*v^5 + 54*i*u^9*v^4 -    48*i*u^9*v^3 - 90*i*u^9*v^2 + 120*i*u^9*v + (-20*i - 20)*u^9 - u^8*v^8 + 8*u^8*v^6 - 24*u^8*v^4 - 64*i*u^8*v^3 + (200*i + 32)*u^8*v^2 - 144*i*u^8*v + 15*u^8 - 4*u^7*v^8 + 16*u^7*v^7 + 16*u^7*v^6 -     96*u^7*v^5 + (64*i + 192)*u^7*v^3 + (-196*i - 64)*u^7*v^2 + (112*i - 128)*u^7*v + (24*i + 24)*u^7 - 6*u^6*v^8 + 48*u^6*v^7 - 96*u^6*v^6 - 96*u^6*v^5 + 432*u^6*v^4 + (-64*i - 192)*u^6*v^3 + (192*i -     384)*u^6*v^2 + (-96*i + 384)*u^6*v + (-32*i - 52)*u^6 - 4*u^5*v^8 + 48*u^5*v^7 - 208*u^5*v^6 + 352*u^5*v^5 + (64*i - 704)*u^5*v^3 + (-188*i + 832)*u^5*v^2 + (96*i - 384)*u^5*v + (24*i + 24)*u^5 - u^4*v^8    + 16*u^4*v^7 - 104*u^4*v^6 + 352*u^4*v^5 - 664*u^4*v^4 + (-64*i + 704)*u^4*v^3 + (184*i - 416)*u^4*v^2 + (-112*i + 128)*u^4*v + 15*u^4 - i*u^3*v^6 + 6*i*u^3*v^4 + 64*i*u^3*v^3 - 198*i*u^3*v^2 +     144*i*u^3*v + (-20*i - 20)*u^3 - 2*i*u^2*v^6 + 12*i*u^2*v^5 - 112*i*u^2*v^3 + 212*i*u^2*v^2 - 120*i*u^2*v + (16*i + 10)*u^2 - i*u*v^6 + 12*i*u*v^5 - 54*i*u*v^4 + 112*i*u*v^3 - 106*i*u*v^2 + 40*i*u*v +     (-4*i - 4)*u + 1;
q := (-u*v - v + 2)/2;
t := (u*v + v - 2)/(2*u);
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];