vdpol.m

matlab从入门到精通第6章书后源码对初学者很有帮助

function ydot=vdpol(t,y,mu)
% Vad der Pol微分方程组
% ydot(1)=y(2);
% ydot(2)=mu*(1-Y(1)^2)*Y(2)-Y(1)
% 输入mu
if nargin<3
    mu=2;
end
ydot=[y(2);     mu*(1-y(1)^2)*y(2)-y(1)];