%你是要解決微分方程，程式我運行了都對！

%你的問題是函式

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檔案，直接在命令視窗輸入檔名即可！