这个程序运行怎么要出问题呢，问题在哪
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编辑过]