function x=F1(y)
global ep
x(1)=y(2);
x(2)=-y(1)+ep*(1-y(1)^2)*y(2);