aufgabe5(1).m

来自「作业题 就是说matlab中的常微分方程问题」· M 代码 · 共 15 行

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)