function dx=rigid(t,x)
dx =zeros(5,1);
dx(1)=-x(4)*x(1)*x(3)/12000000;
dx(2)=x(4)*x(1)*x(3)/12000000-1/5.3*x(2);
dx(3)=1/5.3*x(2)-1/5.61*x(3);
dx(4)=1/5.3/5.3*x(3);
dx(2)=1/5.3*x(2);
x(4)=m(t);
下面是定义的函数m.m
function y=m(x);
if x<150
y=0.27;
else
y=0.27*exp(-0.048*(1-x));
end
我是画图画不出来