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