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