问题分析：函数文件和ode45语句；ode45要对这个函数进行仿真，所以不能放在函数文件中。
正确代码：
1. 保存函数文件：rigid.m
function dy=rigid(t,y)
dy=zeros(3,1);
dy(1)=0.7*y(2)*y(1)/(0.02+y(2));
dy(2)=0.7*y(2)/(0.02+y(2))/0.5*y(1);
dy(3)=0.1*dy(1)+0.02*y(1);
end
2. command window 输入：
[t,y]=ode45('rigid',[0 100],[0 0 0])
内容全部来源于网络收集，如有侵权，请联系网站删除：QQ：24596024