临风暖阳的博客

function fy=VDP(t,x) fy=[x(2);7*(1-x(1)^2)*x(2)-x(1)]; Y0=[1;0]; [t,x]=ode45(‘VDP’,[0,40],Y0); y=x(:,1); dy=x(:,2); plot(t,y,t,dy); Y0=[1;0]; [t,x]=ode23(‘VDP’,[0,40],Y0); y=x(:,1); dy=x(:,2); plot(t,y,t,dy); Y0=[1;0]; [t,x]=ode113(‘VDP’,[0,40],Y0); y=x(