ex5_23.m

来自「源代码集合」· M 代码 · 共 5 行

h_opt=odeset; x0=[-0.2; -0.7]; t_final=20;
mu=1; [t1,y1]=ode45('vdp_eq',[0,t_final],x0,h_opt,mu);
mu=2; [t2,y2]=ode45('vdp_eq',[0,t_final],x0,h_opt,mu);
plot(y1(:,1),y1(:,2),y2(:,1),y2(:,2))