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