|
#2
dsldsl2006-09-16 13:17
我也有同样的问题，同情
 |
format long
c0=1.;
y(1)=c0;
y(2)=0.0000;
y(3)=0.0000;
[t,y]=ode15s('diffreac',[0.0,0.4,1.0,4.0,10.0],[y(1);y(2);y(3)])
函数程序
fuction dy=diffreac(t,y)
format long;
k(1)=0.4;
k(2)=10000.;
k(3)=3.0*10^7;
dy=[(-k(1)*y(1)+k(2)*y(2)*y(3));(k(1)*y(1)-k(2)*y(2)*y(3)-k(3)*y(2)^2);(k(3)*y(2)^2)]
错误
Error in ==> C:\MATLAB6p5\work\zhu.m
On line 6 ==> [t,y]=ode15s('diffreac',[0.0,0.4,1.0,4.0,10.0],[y(1);y(2);y(3)])
[此贴子已经被作者于2006-9-16 11:43:39编辑过]