N:=30;
X := y^6
        + (x^6 - 5*x^5 + 6*x^4 + 3*x^3 - 6*x^2 + 7*x + 3)*y^5
        + (x^7 - 3*x^6 - 13*x^5 + 44*x^4 - 18*x^3 + x^2 + 18*x + 3)*y^4
        + (x^8 - 3*x^7 - 13*x^6 + 27*x^5 + 46*x^4 - 32*x^3 + 21*x^2 + 15*x + 1)*y^3
        + 2*x*(x^7 - 8*x^6 + 9*x^5 + 20*x^4 + 6*x^3 - 6*x^2 + 9*x + 2)*y^2
        - 4*x^2*(2*x^5 - 7*x^4 - 3*x^3 - 1)*y
        + 8*x^6;
q := y+1;
t := 4*(y+1)*(x+y)/(x*y^3 - 4*x*y - 4*x - 3*y^3 - 6*y^2 - 4*y);
E:=[0,t^2-2*q*t-2,0,-(t^2-1)*(q*t+1)^2,0];
P:=[(t+1)*(q*t+1),t*(q*t+1)*(t+1)];