[交流] 【求助】ode45数值求解常微分方程总是错，找不出原因

程序如下：S1=0.4; S2=0.8; y=zeros(42,3); y(1,1)=6e10; y(2,1)=4e10; y(3,1)=0; for i=0:5; for j=1:7         if(j==6)|(j==7)         [t,y]=ode45('tumor',[i*7+j  i*7+j+1],[y(i*7+j,1) y(i*7+j,2) y(i*7+j,3)]);         else            a=y(i*7+j,1);            y(i*7+j,1)=S1*a;            b=y(i*7+j,2);            y(i*7+j,2)=S2*b;            y(i*7+j,3)=(1-S1)*a+(1-S2)*b;           [t,y]=ode45('tumor',[i*7+j i*7+j+1],[y(i*7+j,1) y(i*7+j,2) y(i*7+j,3)]);         end    end   end plot(t,y(1,,'*',t,y(2,,'.',t,y(3,,'+'); 错误为：Index exceeds matrix dimensions. Error in ==> number1 at 12         [t,y]=ode45('tumor',[i*7+j i*7+j+1],[y(i*7+j,1) y(i*7+j,2) y(i*7+j,3)]); 我估计是向量传递出错了 发现ode45返回的是t和y(:,1) y(:,2) y(:,3)的列向量， 如果想得到向量的最后一个值应该怎么办？ [ Last edited by liangbao on 2010-11-25 at 09:51 ]

回复此楼