vdp.m

作业题 就是说matlab中的常微分方程问题

function fy=vdp(t,x)

fy=[x(2);7*(1-x(1)^2)*x(2)-x(1)];

%再编写m文件sy12_6.m如下：

function hehe

y0=[1;0]

[t,x]=ode45(@vdp,[0,40],y0);

y=x(:,1);dy=x(:,2);

plot(t,y,t,dy)