function xdot=fdot2D(xi,t)

global e a

x=xi(1);y=xi(2);
xdot=zeros(2,1);

xdot(1)=y;
xdot(2)=-x+e*y.*(1-a*x.^2);

%xdot(1)=y;
%xdot(2)=-sin(x);
%xlabel('$h$','interpreter','latex')