主题: [MATLAB数学相关] ODE45求解微分方程
查看单个帖子
旧 2013-04-27, 10:35   #1
w43787545o
初级会员
 
注册日期: 2012-03-02
年龄: 43
帖子: 1
声望力: 0
w43787545o 正向着好的方向发展
帖子 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),我觉得自己编的没有问题,可为什么运行不了?求指教
w43787545o 当前离线   回复时引用此帖