N:=52;
X := x^10*y^21
        + 21*x^11*y^20
        + x*(x^20 - 2*x^19 + 5*x^18 - 8*x^17 + 14*x^16 - 20*x^15 + 30*x^14 - 40*x^13 + 55*x^12 + 140*x^11 + 70*x^10 - 70*x^9 + 55*x^8 - 40*x^7 + 30*x^6 - 20*x^5 + 14*x^4 - 8*x^3 + 5*x^2 - 2*x + 1)*y^19
        + x*(10*x^21 - 13*x^20 + 24*x^19 - 13*x^18 + 12*x^17 + 18*x^16 - 16*x^15 + 26*x^14 + 54*x^13 + 1151*x^12 + 84*x^11 - 84*x^10 - 106*x^9 + 231*x^8 - 216*x^7 + 206*x^6 - 132*x^5 + 102*x^4 - 44*x^3 + 33*x^2 - 6*x + 9)*y^18
        + 3*x*(15*x^22 - 10*x^21 + 4*x^20 + 43*x^19 - 61*x^18 + 103*x^17 - 42*x^16 - 37*x^15 + 361*x^14 + 1134*x^13 + 692*x^12 - 834*x^11 + 407*x^10 - 111*x^9 + 166*x^8 - 97*x^7 + 174*x^6 - 77*x^5 + 116*x^4 - 23*x^3 + 52*x^2 +8*x + 12)*y^17
        + x*(120*x^23 - 15*x^22 - 180*x^21 + 587*x^20 - 478*x^19 + 336*x^18 + 760*x^17 - 1735*x^16 + 5106*x^15 + 9035*x^14 + 5832*x^13 - 6766*x^12 - 476*x^11 + 4674*x^10 - 2456*x^9 + 2493*x^8 - 254*x^7 + 1331*x^6 + 276*x^5 +937*x^4 + 382*x^3 + 588*x^2 + 168*x + 84)*y^16
        + x*(210*x^24 + 60*x^23 - 630*x^22 + 1045*x^21 + 19*x^20 - 1131*x^19 + 2734*x^18 - 2631*x^17 + 7968*x^16 + 26913*x^15 + 2540*x^14 - 8880*x^13 - 11870*x^12 + 19936*x^11 - 5220*x^10 + 2209*x^9 + 5644*x^8 + 729*x^7 + 4578*x^6 + 2365*x^5 + 3375*x^4 + 2285*x^3 + 1470*x^2 + 420*x + 126)*y^15
        + 3*x*(84*x^25 + 42*x^24 - 360*x^23 + 243*x^22 + 716*x^21 - 1428*x^20 + 1228*x^19 - 69*x^18 + 3434*x^17 + 14547*x^16 + 3536*x^15 - 15168*x^14 + 1764*x^13 + 2566*x^12 + 10268*x^11 - 10191*x^10 + 11402*x^9 - 1257*x^8 + 4584*x^7 + 3375*x^6 + 3444*x^5 + 3208*x^4 + 1784*x^3 + 770*x^2 + 196*x + 42)*y^14
        + (x^28 + 202*x^27 + 128*x^26 - 1276*x^25 + 48*x^24 + 3729*x^23 - 4357*x^22 - 3210*x^21 + 11876*x^20 + 5430*x^19 + 55679*x^18 + 7572*x^17 - 72839*x^16 - 21414*x^15 + 70000*x^14 - 14736*x^13 + 22830*x^12 - 28912*x^11 + 76805*x^10 - 28884*x^9 + 57264*x^8 + 12225*x^7 + 31343*x^6 + 13802*x^5 + 7475*x^4 + 2084*x^3 + 548*x^2 + 76*x + 1)*y^13
        + (3*x^29 + 98*x^28 + 90*x^27 - 1092*x^26 - 264*x^25 + 3588*x^24 - 2034*x^23 - 9424*x^22 + 15290*x^21 + 17220*x^20 + 25369*x^19 + 23804*x^18 - 116961*x^17 - 18214*x^16 + 82306*x^15 + 30144*x^14 - 67482*x^13 + 62192*x^12 - 10323*x^11 + 92060*x^10 - 27672*x^9 + 119676*x^8 + 4822*x^7 + 57544*x^6 + 3471*x^5 + 9078*x^4 + 52*x^3 + 612*x^2 - 33*x + 10)*y^12
        + 3*(x^30 + 8*x^29 + 11*x^28 - 192*x^27 - 110*x^26 + 926*x^25 - 310*x^24 - 2688*x^23 + 2710*x^22 + 9515*x^21 - 530*x^20 - 112*x^19 - 25997*x^18 - 22938*x^17 + 48365*x^16 - 564*x^15 - 20010*x^14 - 3018*x^13 + 32980*x^12 - 6224*x^11 + 36828*x^10 + 6160*x^9 + 48016*x^8 - 4332*x^7 + 21437*x^6 - 5857*x^5 + 4313*x^4 - 1176*x^3 + 429*x^2 - 84*x + 15)*y^11
        + (x^31 + 2*x^30 - 3*x^29 - 148*x^28 - 198*x^27 + 1460*x^26 - 116*x^25 - 5176*x^24 + 4534*x^23 + 19678*x^22 - 1473*x^21 - 34644*x^20 - 25349*x^19 - 51614*x^18 + 93917*x^17 + 64524*x^16 - 146738*x^15 + 34112*x^14 + 52802*x^13 + 54972*x^12 + 8640*x^11 + 142604*x^10 + 40404*x^9 + 99284*x^8 - 29683*x^7 + 46844*x^6 - 27772*x^5 + 16588*x^4 - 6727*x^3 + 2434*x^2 - 563*x + 120)*y^10
        - (4*x^30 + 12*x^29 + 30*x^28 - 362*x^27 - 270*x^26 + 2476*x^25 - 2073*x^24 - 9507*x^23 + 2964*x^22 + 26808*x^21 + 15921*x^20 - 9710*x^19 - 20736*x^18 - 68292*x^17 + 102845*x^16 + 38466*x^15 - 122637*x^14 + 7380*x^13 - 48738*x^12 - 69492*x^11 - 107872*x^10 - 5372*x^9 - 36972*x^8 + 27513*x^7 - 26155*x^6 + 20112*x^5 - 14683*x^4 + 6972*x^3 - 3108*x^2 + 756*x - 210)*y^9
        + 3*(2*x^29 + 8*x^28 + 28*x^27 - 180*x^26 + 53*x^25 + 1222*x^24 - 945*x^23 - 4016*x^22 - 1637*x^21 + 5410*x^20 + 5130*x^19 + 1120*x^18 - 5149*x^17 - 21258*x^16 + 31237*x^15 + 11000*x^14 - 5998*x^13 + 29280*x^12 + 15888*x^11 + 8124*x^10 - 5916*x^9 + 3900*x^8 - 3660*x^7 + 4608*x^6 - 2697*x^5 + 2906*x^4 - 1400*x^3 + 896*x^2 - 210*x + 84)*y^8
        - (4*x^28 + 20*x^27 + 81*x^26 - 672*x^25 + 617*x^24 + 4106*x^23 - 771*x^22 - 8089*x^21 - 6375*x^20 + 6206*x^19 + 10143*x^18 + 6224*x^17 - 13049*x^16 - 59820*x^15 + 30578*x^14 - 39010*x^13 - 47242*x^12 + 12104*x^11 + 4924*x^10 + 3654*x^9 - 9616*x^8 - 102*x^7 - 6965*x^6 + 1247*x^5 - 4057*x^4 + 1134*x^3 - 1638*x^2 + 294*x - 210)*y^7
        + (x^27 + 6*x^26 + 27*x^25 - 700*x^24 + 471*x^23 + 2760*x^22 + 604*x^21 - 3468*x^20 - 3851*x^19 + 1798*x^18 + 8001*x^17 + 12336*x^16 + 8938*x^15 - 11808*x^14 + 35512*x^13 - 12244*x^12 - 11616*x^11 + 6776*x^10 + 5398*x^9 + 5700*x^8 + 2097*x^7 + 3420*x^6 + 1062*x^5 + 1952*x^4 + 264*x^3 + 744*x^2 - 36*x + 120)*y^6
        + 3*(165*x^23 - 55*x^22 - 374*x^21 - 194*x^20 + 446*x^19 + 861*x^18 + 828*x^17 + 321*x^16 - 114*x^15 - 74*x^14 + 2136*x^13 - 4693*x^12 + 2418*x^11 + 1776*x^10 + 616*x^9 + 666*x^8 + 553*x^7 + 603*x^6 + 366*x^5 + 301*x^4 + 116*x^3 + 88*x^2 + 12*x + 15)*y^5
        - (220*x^22 - 44*x^21 - 352*x^20 - 485*x^19 - 108*x^18 + 179*x^17 + 530*x^16 + 364*x^15 + 352*x^14 + 225*x^13 + 1618*x^12 - 5418*x^11 + 872*x^10 - 458*x^9 - 948*x^8 - 878*x^7 - 708*x^6 - 443*x^5 - 278*x^4 - 128*x^3- 68*x^2 - 19*x - 10)*y^4
        + (66*x^21 + 6*x^20 - 45*x^19 - 137*x^18 - 120*x^17 - 131*x^16 - 26*x^15 + 15*x^14 + 187*x^13 + 311*x^12 + 690*x^11 - 834*x^10 + 623*x^9 + 276*x^8 + 190*x^7 + 119*x^6 + 69*x^5 + 39*x^4 + 19*x^3 + 9*x^2 + 3*x + 1)*y^3
        - 3*x^7*(4*x^13 + 2*x^12 + 2*x^11 - 3*x^10 - 4*x^9 - 9*x^8 - 8*x^7 - 10*x^6 - 4*x^5 + 14*x^3 - 65*x^2 + 12*x - 1)*y^2
        + x^7*(x^12 + x^11 + 2*x^10 + 2*x^9 + 3*x^8 + 3*x^7 + 4*x^6 + 4*x^5 + 5*x^4 + 5*x^3 + 6*x^2 - 18*x + 3)*y
        + x^7;
q :=(x*y - x - y - 1)/(x*y + x + y + 1);
t := (-x*y + x + y + 1)/(x*y + x - y - 1);
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];