odefun4.m

《matlab使用详解》的部分源码

function dy = odefun4(t,y)
% Vander Pol微分方程
dy = [0;0];
dy(1) = y(2);
dy(2) = -(y(1)^2 - 1)*y(2) - y(1);