N:=13;
X:=y^2 + (x^3 + x^2 + 1)*y - x^2 - x;
r:=-x*y + 1;
s:=1-x*y/(y + 1);
E:=[s-r*s+1,r*s-r^2*s,r*s-r^2*s,0,0];
P:=[0,0];