匿名用户
1级
2015-02-07 回答
因为t<51时你没给a赋值, 所以我把t==51改成了t<=51了, 这样的代码如下:
f=@(t,x)[-3*10^(-8)*(1-((t<=51)*0.4+(((t<=60)&(t>51))*0.7)+(t> 60)*0.9))*x(3)*x(1)-1*10^(-11)*(1-((t<=51)*0.4+(((t<=60)&(t>51))*0.7)+(t> 60)*0.9))*x(5)*x(1)-1*10^(-11)*((t<=51)*0.4+(((t<=60)&(t>51))*0.7)+(t> 60)*0.9)*0.05*x(1)*x(5);3*10^(-8)*x(3)*x(1)*(1-((t<=51)*0.4+(((t<=60)&(t>51))*0.7)+(t> 60)*0.9))+3*10^(-8)*(1-((t<=51)*0.4+(((t<=60)&(t>51))*0.7)+(t> 60)*0.9))*(1-((t<=51)*0.4+(((t<=60)&(t>51))*0.7)+(t> 60)*0.9))*x(5)*x(3)+3*10^(-8)*(1-((t<=51)*0.4+(((t<=60)&(t>51))*0.7)+(t> 60)*0.9))*x(3)*((t<=51)*0.4+(((t<=60)&(t>51))*0.7)+(t> 60)*0.9)*0.05*x(5)-0.2*x(2);0.2*x(