function [xp] = ciclo(t,x) xp = zeros(2,1); mu=1; xp(1)=-x(2)+x(1)*(mu-x(1)^2-x(2)^2); xp(2)=x(1)+x(2)*(mu-x(1)^2-x(2)^2); end