N:=37;
X:=(x + 1)^4*y^18
        + (x + 1)^3*(4*x^3 + 15*x^2 + 4*x - 8)*y^17
        + (x + 1)^3*(6*x^5 + 46*x^4 + 71*x^3 - 39*x^2 - 56*x + 28)*y^16
        + (x + 1)^3*(4*x^7 + 58*x^6 + 204*x^5 + 109*x^4 - 300*x^3 - 91*x^2 + 186*x - 56)*y^15
        + (x + 1)^2*(x^10 + 34*x^9 + 265*x^8 + 692*x^7 + 295*x^6 - 1055*x^5 - 745*x^4 + 570*x^3 + 135*x^2 - 205*x + 70)*y^14
        + (x + 1)*(7*x^12 + 130*x^11 + 755*x^10 + 1659*x^9 + 426*x^8 - 3357*x^7 - 3594*x^6 + 612*x^5 + 1760*x^4 + 274*x^3 + 32*x^2 + 6*x - 56)*y^13
        + (21*x^14 + 294*x^13 + 1420*x^12 + 2579*x^11 - 870*x^10 - 9749*x^9 - 11481*x^8 + 664*x^7 + 10483*x^6 + 6962*x^5 - 527*x^4 - 2730*x^3 - 733*x^2 + 322*x + 28)*y^12
        + (35*x^15 + 363*x^14 + 1210*x^13 + 414*x^12 - 6515*x^11 - 14529*x^10 - 6010*x^9 + 16913*x^8 + 23552*x^7 + 4413*x^6 - 13117*x^5 - 8675*x^4 + 2488*x^3 + 2270*x^2 - 494*x - 8)*y^11
        + (26*x^16 + 184*x^15 + 129*x^14 - 2320*x^13 - 8129*x^12 - 7355*x^11 + 12096*x^10 + 31592*x^9 + 16424*x^8 - 20987*x^7 - 29776*x^6 - 1252*x^5 + 14573*x^4 + 1991*x^3 - 3145*x^2 + 412*x + 1)*y^10
        - x*(x^17 + 5*x^16 + 52*x^15 + 501*x^14 + 2023*x^13 + 2681*x^12 - 4715*x^11 - 19581*x^10 - 18001*x^9 + 14660*x^8 + 40593*x^7 + 16842*x^6 - 25005*x^5 - 20440*x^4 + 8210*x^3 + 5975*x^2 - 2525*x + 210)*y^9
        - x*(x^18 + 10*x^17 + 52*x^16 + 172*x^15 + 162*x^14 - 1311*x^13 - 5622*x^12 - 7317*x^11 + 5435*x^10 + 25859*x^9 + 20839*x^8 - 16281*x^7 - 34475*x^6 - 2999*x^5 + 20366*x^4 + 2065*x^3 - 5725*x^2 + 1264*x - 66)*y^8
        + x*(3*x^17 + 32*x^16 + 180*x^15 + 595*x^14 + 789*x^13 - 1585*x^12 - 7639*x^11 - 9008*x^10 + 5183*x^9 + 22272*x^8 + 12561*x^7 - 16256*x^6 - 16253*x^5 + 7183*x^4 + 5741*x^3 - 3053*x^2 + 395*x - 12)*y^7
        - x*(3*x^16 + 39*x^15 + 243*x^14 + 807*x^13 + 1079*x^12 - 1278*x^11 - 6267*x^10 - 6155*x^9 + 4834*x^8 + 13100*x^7 + 1905*x^6 - 9941*x^5 - 1372*x^4 + 3726*x^3 - 978*x^2 + 72*x - 1)*y^6
        + x^2*(x^14 + 22*x^13 + 166*x^12 + 564*x^11 + 678*x^10 - 891*x^9 - 3469*x^8 - 2293*x^7 + 3656*x^6 + 4495*x^5 - 2053*x^4 - 2249*x^3 + 1275*x^2 - 180*x + 6)*y^5
        - 5*x^3*(x^11 + 12*x^10 + 42*x^9 + 36*x^8 - 111*x^7 - 270*x^6 - 37*x^5 + 331*x^4 + 80*x^3 - 182*x^2 + 48*x - 3)*y^4
        + x^4*(10*x^8 + 35*x^7 - 15*x^6 - 265*x^5 - 325*x^4 + 201*x^3 + 321*x^2 - 180*x + 20)*y^3
        + x^5*(26*x^4 + 77*x^3 + 30*x^2 - 72*x + 15)*y^2
        - x^6*(x^3 + 7*x^2 + 12*x - 6)*y
        + x^7;
r:=(x^2*y + x*y + y - 1)/(x^2 + x + y - 1);
s:=(x*y + y - 1)/(x + y - 1);
E:=[s-r*s+1,r*s-r^2*s,r*s-r^2*s,0,0];
P:=[0,0];