N:=44; X:=x^7*y^15 + 15*x^8*y^14 + (x^16 - 2*x^15 + 5*x^14 - 8*x^13 + 14*x^12 - 20*x^11 + 30*x^10 + 65*x^9 + 40*x^8 - 40*x^7 + 30*x^6 - 20*x^5 + 14*x^4 - 8*x^3 + 5*x^2 - 2*x + 1)*y^13 + (6*x^17 - 7*x^16 + 12*x^15 - 3*x^14 + 4*x^13 + 6*x^12 + 24*x^11 + 381*x^10 + 24*x^9 - 24*x^8 - 56*x^7 + 106*x^6 - 84*x^5 + 74*x^4 - 36*x^3 + 27*x^2 - 6*x + 7)*y^12 + 3*(5*x^18 - 2*x^17 - 3*x^16 + 21*x^15 - 20*x^14 + 15*x^13 + 64*x^12 + 246*x^11 + 150*x^10 - 202*x^9 + 90*x^8 - 9*x^7 + 48*x^6 - 25*x^5 + 54*x^4 - 15*x^3 + 29*x^2 + 2*x + 7)*y^11 + (20*x^19 + 5*x^18 - 62*x^17 + 130*x^16 - 44*x^15 - 119*x^14 + 722*x^13 + 974*x^12 + 900*x^11 - 1176*x^10 - 100*x^9 + 893*x^8 - 292*x^7 + 389*x^6 + 86*x^5 + 275*x^4 + 112*x^3 + 205*x^2 + 50*x + 35)*y^10 + (15*x^20 + 10*x^19 - 90*x^18 + 93*x^17 + 86*x^16 - 304*x^15 + 882*x^14 + 1645*x^13 - 61*x^12 - 714*x^11 - 1544*x^10 + 2431*x^9 - 321*x^8 + 140*x^7 + 847*x^6 + 336*x^5 + 656*x^4 + 458*x^3 + 315*x^2 + 90*x + 35)*y^9 + 3*(2*x^21 + x^20 - 20*x^19 - x^18 + 64*x^17 - 106*x^16 + 208*x^15 + 527*x^14 - 138*x^13 - 699*x^12 + 20*x^11 + 151*x^10 + 870*x^9 - 777*x^8 + 784*x^7 + 173*x^6 + 376*x^5 + 374*x^4 + 208*x^3 + 95*x^2 + 26*x + 7)*y^8 + (x^22 - x^21 - 21*x^20 - 19*x^19 + 110*x^18 - 53*x^17 + 46*x^16 + 1185*x^15 - 777*x^14 - 2004*x^13 - 211*x^12 + 1844*x^11 - 19*x^10 + 654*x^9 - 1437*x^8 + 3882*x^7 - 332*x^6 + 2257*x^5 + 728*x^4 + 431*x^3 + 129*x^2 + 35*x + 7)*y^7 - (4*x^21 - 24*x^19 - 45*x^18 + 66*x^17 - 341*x^16 + 204*x^15 + 1677*x^14 - 378*x^13 - 1183*x^12 - 288*x^11 + 1105*x^10 - 774*x^9 - 543*x^8 - 2868*x^7 + 499*x^6 - 2154*x^5 + 311*x^4 - 272*x^3 + 12*x^2 - 12*x - 1)*y^6 + 3*x*(2*x^19 + 2*x^18 - 8*x^16 + 5*x^15 - 192*x^14 + 15*x^13 + 360*x^12 - 115*x^11 - 172*x^10 + 12*x^9 + 430*x^8 + 225*x^7 + 366*x^6 - 100*x^5 + 312*x^4 - 215*x^3 + 94*x^2 - 25*x + 5)*y^5 - x*(4*x^18 + 8*x^17 + 18*x^16 + 23*x^15 + 48*x^14 - 397*x^13 + 58*x^12 + 463*x^11 - 336*x^10 - 530*x^9 - 550*x^8 - 75*x^7 - 188*x^6 + 164*x^5 - 186*x^4 + 269*x^3 - 208*x^2 + 70*x - 20)*y^4 + x*(x^17 + 3*x^16 + 9*x^15 + 19*x^14 + 39*x^13 - 151*x^12 + 119*x^11 + 279*x^10 + 73*x^9 + 16*x^8 - 51*x^7 + 25*x^6 - 23*x^5 + 57*x^4 - 17*x^3 + 72*x^2 - 30*x + 15)*y^3 + 3*x*(22*x^11 - 4*x^10 - 13*x^9 - 2*x^8 + 2*x^7 + 8*x^6 + 6*x^5 + 8*x^4 + 3*x^3 + 4*x^2 - x + 2)*y^2 - x*(12*x^10 - 2*x^9 - 5*x^8 - 4*x^7 - 4*x^6 - 3*x^5 - 3*x^4 - 2*x^3 - 2*x^2 - x - 1)*y + x^10; 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];