function  v=xdot2G(psi,t);

global epsilon

v=zeros(2,1);
x=psi(1);y=psi(2);


v(1)=y;
v(2)=-x+epsilon*(y-y^3/3);
