options_jac.m

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

% options_jac.m
% 根据雅可比矩阵来求解
tspan=[0 20];
yo=[2;0];
mu=10;
options=odeset(options,'Jacobian',@vdpoljac);
[t,y]=ode15s(@vdpol,tspan,yo,options,mu);
% 绘制数据结果
plot(t,y(:,1),t,y(:,2),'--');
xlabel('t');
legend('x','dx/dt');
title(['Van der Pol solution num(t)=' num2str(length(t))])