如果你正在用 Matlab 做微分方程求解，这个小技巧可以让你事半功倍。通过创建一个 m 文件，比如vdp1000.m，你就能一些典型的微分方程问题。代码如下：

function dy=vdp1000(t,y)
  dy=zeros(2,1);
  dy(1)=y(2);
  dy(2)=1000*(1-y(1)^2)*y(2)-y(1);
。之后，使用ode15s来求解，就能得到你想要的结果：
[T,Y]=ode15s('vdp1000',[0 3000],[2 0]);
plot(T,Y(:,1),'-')
。图像生成的效果蛮直观，能你理解微分方程的解。
如果你对相关的内容有兴趣，也可以参考以下链接，获取更多的技术资源：