examp7_14.m

来自「《高等应用数学问题MATLAB求解》的PDF版及习题解答」· M 代码 · 共 5 行

h_opt=odeset; h_opt.RelTol=1e-6; x0=[2;0]; t_final=3000;
tic, mu=1000; [t,y]=ode15s('vdp_eq',[0,t_final],x0,h_opt,mu); toc

plot(t,y(:,1)); figure; plot(t,y(:,2))