3-40.m

来自「《MATLAB 7.0编程基础》第3章 （源码实例）主要讲解MATLAB的数值计」· M 代码 · 共 8 行

tspan = [0 12];
% 在同一目录下，计算方程数值解。给出时间区间和初始值
y0 = [-1 1 0];
[t,Y] = ode45('odet3',tspan,y0);
% 采用ode45算法
plot(t,Y(:,1),'-',t,Y(:,2),'-.',t,Y(:,3),'.')
% 绘制计算结果并标注，如图3-15
legend('Y1','Y2','Y3')