N:=42;
X := (x + 1)^15*y^12
        + (x + 1)^3*(x^17 + 15*x^16 + 108*x^15 + 494*x^14 + 1605*x^13 + 3918*x^12 + 7408*x^11 + 11052*x^10 + 13233*x^9 + 13035*x^8 + 10968*x^7 + 8190*x^6 + 5451*x^5 + 3075*x^4 + 1350*x^3 + 416*x^2 + 78*x + 6)*y^11
        - (x + 1)^3*(4*x^17 + 49*x^16 + 270*x^15 + 860*x^14 + 1610*x^13 + 1113*x^12 - 3122*x^11 - 12260*x^10 - 23205*x^9 - 29865*x^8 - 28806*x^7 - 22176*x^6 - 14346*x^5 - 7896*x^4 - 3510*x^3 - 1120*x^2 - 220*x - 15)*y^10
        + (10*x^20 + 125*x^19 + 685*x^18 + 2085*x^17 + 3360*x^16 + 345*x^15 - 11775*x^14 - 27969*x^13 - 25488*x^12 + 20845*x^11 + 105335*x^10 + 183295*x^9 + 209676*x^8 + 179547*x^7 + 122395*x^6 + 68555*x^5 + 31390*x^4 + 11171*x^3 + 2807*x^2 + 415*x + 20)*y^9
        - (20*x^20 + 195*x^19 + 795*x^18 + 1635*x^17 + 1155*x^16 - 2070*x^15 - 3690*x^14 + 7986*x^13 + 36630*x^12 + 58664*x^11 + 38700*x^10 - 24360*x^9 - 86859*x^8 - 108120*x^7 - 88776*x^6 - 55050*x^5 - 26985*x^4 - 10179*x^3 - 2695*x^2 - 411*x - 15)*y^8
        + (35*x^20 + 245*x^19 + 735*x^18 + 1155*x^17 + 1140*x^16 + 2529*x^15 + 9645*x^14 + 21105*x^13 + 22365*x^12 + 1400*x^11 - 25270*x^10 - 25950*x^9 + 2460*x^8 + 30735*x^7 + 37395*x^6 + 27537*x^5 + 14925*x^4 + 6075*x^3 + 1715*x^2 + 275*x + 6)*y^7
        + (x^24 + 12*x^23 + 84*x^22 + 400*x^21 + 1420*x^20 + 4082*x^19 + 9946*x^18 + 20610*x^17 + 35964*x^16 + 52771*x^15 + 67545*x^14 + 80709*x^13 + 91288*x^12 + 88150*x^11 + 63044*x^10 + 30490*x^9 + 12864*x^8 + 11832*x^7+ 13432*x^6 + 10452*x^5 + 5970*x^4 + 2563*x^3 + 763*x^2 + 131*x + 1)*y^6
        - x*(x^22 + 9*x^21 + 21*x^20 - 50*x^19 - 510*x^18 - 2178*x^17 - 6642*x^16 - 16008*x^15 - 30807*x^14 - 47543*x^13 - 60873*x^12 - 69432*x^11 - 72773*x^10 - 65379*x^9 - 45891*x^8 - 24816*x^7 - 12054*x^6 - 6966*x^5 - 4010*x^4 - 2010*x^3 - 825*x^2 - 245*x - 45)*y^5
        + x*(x^22 + 7*x^21 + 15*x^20 + 10*x^19 - 65*x^18 - 279*x^17 - 633*x^16 - 675*x^15 + 810*x^14 + 4870*x^13 + 10162*x^12 + 14334*x^11 + 17410*x^10 + 19030*x^9 + 16650*x^8 + 10881*x^7 + 5307*x^6 + 2565*x^5 + 1195*x^4 + 505*x^3 + 186*x^2 + 52*x + 10)*y^4
        - x*(x^22 + 5*x^21 + 19*x^20 + 55*x^19 + 140*x^18 + 316*x^17 + 660*x^16 + 1284*x^15 + 2205*x^14 + 2905*x^13 + 2717*x^12 + 2005*x^11 + 1198*x^10 - 40*x^9 - 1040*x^8 - 1284*x^7 - 660*x^6 - 316*x^5 - 140*x^4 - 55*x^3 - 19*x^2 - 5*x - 1)*y^3
        - x^9*(45*x^6 + 215*x^5 + 375*x^4 + 390*x^3 + 430*x^2 + 336*x + 210)*y^2
        + x^9*(9*x^6 + 17*x^5 - x^4 + 6*x^3 - 11*x^2 - 5*x - 21)*y
        - x^9*(x^6 - x^5 + x^4 - x^3 + x^2 - x + 1);
q := -x;
t := -(x+1)/(x*y);
E:=[(1+q)*t+(2-q),0,q*t^2+(1-q)*t,0,0];
P:=[-t,t^2];