3-39.m

《MATLAB 7.0编程基础》第3章 （源码实例）主要讲解MATLAB的数值计算

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