[MATLAB数学相关] ODE45求解微分方程 - MATLAB爱好者论坛-LabFans.com

w43787545o · 2013-04-27, 10:35

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)，我觉得自己编的没有问题，可为什么运行不了？求指教

2013-04-27, 10:35	#1
w43787545o 初级会员注册日期: 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+E1E2/(E3k2); C=E1E2/(k1k2); dy=zeros(2,1); dy(1)=y(2); dy(2)=-By(2)/A-Cy(1)/A; y0=E1E3p0/(E1+E3); [t,y]=ode45(@myfun,[0,100],y0)，我觉得自己编的没有问题，可为什么运行不了？求指教