N:=27; R:=PolynomialRing(Rationals(),2); X:=(u-1)^2*v^6+(u-1)^2*(u^3+2)*v^5-(u-1)^2*(u^5+2*u^4-2*u^3-u^2-2*u-1)*v^4+u*(u-1)*(u^6-3*u^5-4*u^4+u^3+u^2+3*u-2)*v^3+u*(u-1)*(u^2+u+1)*(3*u^4-4*u^3-2*u^2+u-1)*v^2+3*u^5*(u-1)*(u^2+u+1)*v+u^6*(u^2+u+1); G:=-1/u; x:=v/(u*v+u); y:=(G^4*x + G^4 + G^3*x - 2*G^2*x^2 - G*x^3 + G + x)/(G^4 + G^3*x - G^2*x^2 - G*x^2 + G + x); r:=(x^2*y - x*y + y - 1) / (x^2*y - x); s:=(x*y - y + 1) / (x*y); E:=[s-r*s+1,r*s-r^2*s,r*s-r^2*s,0,0]; P:=[0,0];