N:=40;
X := x^4*y^12
        - x^4*(x^9 - x^8 + 2*x^7 - 2*x^6 + 3*x^5 - 3*x^4 + 4*x^3 - 4*x^2 - 10*x - 2)*y^11
        - x^4*(6*x^10 - 2*x^9 + x^8 + 8*x^7 - 7*x^6 + 12*x^5 - 8*x^3 - 65*x^2 - 10*x - 1)*y^10
        - (15*x^15 + 3*x^14 - 20*x^13 + 44*x^12 - 21*x^11 + 15*x^10 + 58*x^9 - 114*x^8 - 123*x^7 - 119*x^6 + 69*x^5 - 39*x^4 + 19*x^3 - 9*x^2 + 3*x - 1)*y^9
        - (20*x^16 + 10*x^15 - 53*x^14 + 60*x^13 + 23*x^12 - 34*x^11 + 123*x^10 - 188*x^9 - 416*x^8 - 42*x^7 - 31*x^6 + 116*x^5 - 116*x^4 + 66*x^3 - 39*x^2 + 12*x - 6)*y^8
        - (15*x^17 + 5*x^16 - 55*x^15 + 7*x^14 + 140*x^13 - 128*x^12 + 125*x^11 - 245*x^10 - 604*x^9 - 128*x^8 + 251*x^7 - 83*x^6 + 38*x^5 - 130*x^4 + 75*x^3 - 75*x^2 + 15*x - 15)*y^7
        - (5*x^18 - 40*x^16 + 105*x^14 - 182*x^12 - 725*x^10 + 348*x^8 - 219*x^6 - 106*x^4 - 90*x^2 - 20)*y^6
        + (15*x^17 - 5*x^16 - 55*x^15 - 7*x^14 + 140*x^13 + 128*x^12 + 125*x^11 + 245*x^10 - 604*x^9 + 128*x^8 + 251*x^7 + 83*x^6 + 38*x^5 + 130*x^4 + 75*x^3 + 75*x^2 + 15*x + 15)*y^5
        - (20*x^16 - 10*x^15 - 53*x^14 - 60*x^13 + 23*x^12 + 34*x^11 + 123*x^10 + 188*x^9 - 416*x^8 + 42*x^7 - 31*x^6 - 116*x^5 - 116*x^4 - 66*x^3 - 39*x^2 - 12*x - 6)*y^4
        + (15*x^15 - 3*x^14 - 20*x^13 - 44*x^12 - 21*x^11 - 15*x^10 + 58*x^9 + 114*x^8 - 123*x^7 + 119*x^6 + 69*x^5 + 39*x^4 + 19*x^3 + 9*x^2 + 3*x + 1)*y^3
        - x^4*(6*x^10 + 2*x^9 + x^8 - 8*x^7 - 7*x^6 - 12*x^5 + 8*x^3 - 65*x^2 + 10*x - 1)*y^2
        + x^4*(x^9 + x^8 + 2*x^7 + 2*x^6 + 3*x^5 + 3*x^4 + 4*x^3 + 4*x^2 - 10*x + 2)*y
        + x^4;
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];