N:=25;
R<u,v>:=PolynomialRing(Rationals(),2);
X := u*v^5 + (u^4 - 2*u^3 - u^2 + 2*u + 1)*v^4 - (2*u^6 - 2*u^4 + 4*u^3 + 2*u^2 - 2)*v^3 + (u^8 + u^7 - 2*u^6 + u^5 - u^4 - u^3 - 2*u^2 - u + 1)*v^2 + (u^8 + u^7 + 2*u^6 + u^5 - 2*u^4 + u^3 - u^2)*v + u^6;
y := (u-v-1)/u;
x := -(v+1)*(u^4 + u^3*v + u^3 - u*v^2 - v)/(u^2*(u^3 + u^2*v + u^2 - u*v - u - v^2 - v));
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];