Matlab实现常微分计算

y10=2;
y20=0;
tspan=[0 20];
f=@(t,y)[y(2);(1-y(1)^2)*y(2)-y(1)];
y0=[2,0];
[t,y]=ode45(f,tspan,y0);
%绘图
plot(t,y)
xlabel \itt;ylable '{\ity}_1,{\ity}_2'
legend('{\ity}_1(\itt)','{\ity}_2(\itt)')