%Dinâmica do Oscilador de Van der Pol function [xponto]=equacao(t,x) global mi % Para mi=0 ==> Centro mi=-2; % Para mi>0 ==> Ciclo limite estável % Para mi<0 ==> Ciclo limite instável xponto=[x(2);-x(1)+mi*(1-x(1)^2)*x(2)];