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