matlab解微分方程ode程序

function dy=vdp(t,y) %#ok<FNDEF>
dy=[y(2);1000(1-y(1)^2)*y(2)-y(1)];
end

[T,Y]=ode15s('vdp',[0 1000],[2 0]);

plot(T,Y(:,1));