www.gusucode.com > 关于海航matlab和lingo的训练题 > vdp1000.m

    %[t,x]=solver('f',ts,x0,options)
%solver:ode45,ode15s,ode23s
%t为自变量，x为函数值，
%ts=[t0,tf],t0、tf为自变量的初值和终值
%options用于设定误差限

function dy=vdp1000(t,y)
% t不可少
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=1000*(1-y(1)^2)*y(2)-y(1);