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(:,1);
dy=x(:,2);
plot(t,y,t,dy);
Y0=[1;0];
[t,x]=ode23t(‘VDP’,[0,40],Y0);
y=x(:,1);
dy=x(:,2);
plot(t,y,t,dy);
Y0=[1;0];
[t,x]=ode15s(‘VDP’,[0,40],Y0);
y=x(:,1);
dy=x(:,2);
plot(t,y,t,dy);
Y0=[1;0];
[t,x]=ode23s(‘VDP’,[0,40],Y0);
y=x(:,1);
dy=x(:,2);
plot(t,y,t,dy);
Y0=[1;0];
[t,x]=ode23tb(‘VDP’,[0,40],Y0);
y=x(:,1);
dy=x(:,2);
plot(t,y,t,dy);
MATLAB概率与数据统计分析—张德丰等编著–第12页–[M]北京:机械工业出版社