N:=28;
X := x^2*(x + 1)^3*y^6 
        + x^2*(x + 1)^2*(x^4 + 11)*y^5
        + 3*x^2*(x + 1)*(3*x^4 - 2*x^2 + 15)*y^4
        - (5*x^8 - 28*x^6 + 42*x^4 - 84*x^2 + 1)*y^3
        - (x - 1)*(25*x^6 - 23*x^4 + 67*x^2 - 5)*y^2
        + 6*(x - 1)^5*(x + 1)^3*y
        - (x - 1)^6*(x + 1)^3;
q := (x*y + y + 2)/(2*x);
t := (-x*y - y - 2)/2;
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];