| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
 |
| Labfans是一个针对大学生、工程师和科研工作者的技术社区。 | 论坛首页 | 联系我们(Contact Us) |
|
2013-04-27, 10:35
|
#1 |
|
初级会员
注册日期: 2012-03-02
年龄: 43
帖子: 1
声望力: 0  |
ODE45求解微分方程
function dy=myfun(t,y)
E1=28.9; E2=45.36; E3=1.95e5; k1=0.45; k2=871.7; p0=0.034; A=1+E1/E3; B=E1/k2+E1/k1+E2/k2+E1*E2/(E3*k2); C=E1*E2/(k1*k2); dy=zeros(2,1); dy(1)=y(2); dy(2)=-B*y(2)/A-C*y(1)/A; y0=E1*E3*p0/(E1+E3); [t,y]=ode45(@myfun,[0,100],y0)，我觉得自己编的没有问题，可为什么运行不了？求指教 |
|
|
树形模式