vdpoljac.m

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

function jac=vdpoljac(t,y,mu);
% 根据Van der Pol的Jacobian矩阵求解
% 输入mu
if nargin<3
    mu=2;
end
jac=[    0                       1
         (-2*mu*y(1)*y(2)-1)  (mu*(1-y(1)^2))];
% ydot=[y(2);     mu*(1-y(1)^2)*y(2)-y(1)];