function dx=ecdif1(t,x,m)
dx=zeros(2,1);
m=3;	%coeficiente de amortiguamiento
dx(1)=x(2);
dx(2)=-x(1)-(m.*x(2).*((x(1).^2)-1));
%fin
