Question

matlab新手编程，请教．．．

function dy=rigid1(t,y)
dy=zeros(3,1);
dy(1)=5-0.2*y(1)-0.0014453*y(1)*y(3);
dy(2)=0.0014453*y(1)*y(3)-0.32*y(2);
dy(3)=-3*y(3)+100*0.32*y(2);

>>options=odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
>> [T,Y]=ode23(@rigid1,[0 80],[80 20 20],options);
>> plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
以M文件运行后，出现这个错误～不懂
??? Error: File: C:\MATLAB6p5\work\rigid1.m Line: 7 Column: 1
"End of Input" expected, ">" found.

Answer 1

%你是要解决微分方程，程序我运行了都对！
%你的问题是函数
function dy=rigid1(t,y)
dy=zeros(3,1);
dy(1)=5-0.2*y(1)-0.0014453*y(1)*y(3);
dy(2)=0.0014453*y(1)*y(3)-0.32*y(2);
dy(3)=-3*y(3)+100*0.32*y(2);
%要保存到一个M文件中命名为rigid1.m
options=odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y]=ode23(@rigid1,[0 80],[80 20 20],options);
plot(T,Y(:,1),'含磨余-',T,Y(:,2),'-.',T,Y(:,3),'.')
%以上三条语句不能放入rigid1.m中，要么再建一个M文件或者直接在命令窗口输入即可！要是建立了M文件，谈滚直接在命令窗口输入文件名即可游丛！