vdpoljac.m

来自「matlab从入门到精通第6章书后源码对初学者很有帮助」· M 代码 · 共 10 行

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)];