function xdot=limcyc(t,x)
E=x(1)^2+x(2)^2;
xdot1=x(2)+x(1)*(1-E);
xdot2=-x(1)+x(2)*(1-E);
xdot=[xdot1;xdot2];
