请教高手下面的程序错误是啥意思

liming1985

opions=odeset('Reltol',[1e-4],'Abstol',[1e-6])
T0=903;
P0=0.1e+6;
y0=[0;0;0;T0;0;P0];
[t,y]=ode45('yifan',[0.8 1.4],y0,opions,N1,Neb01,rou_b,L1,qn0);
h=size(y);
e=h(1);
f=y(e,:);
X1=f(1)+f(2)+f(3);


执行结果是：

opions =
              AbsTol: 1.0000e-006
                 BDF: []
              Events: []
         InitialStep: []
            Jacobian: []
           JConstant: []
            JPattern: []
                Mass: []
        MassConstant: []
        MassSingular: []
            MaxOrder: []
             MaxStep: []
         NonNegative: []
         NormControl: []
           OutputFcn: []
           OutputSel: []
              Refine: []
              RelTol: 1.0000e-004
               Stats: []
          Vectorized: []
    MStateDependence: []
           MvPattern: []
        InitialSlope: []
??? Error using ==> feval
Undefined command/function 'yifan'.
Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...
Error in ==> ZUOYE22 at 30
[t,y]=ode45('yifan',[0.8 1.4],y0,opions,N1,Neb01,rou_b,L1,qn0);

[ 本帖最后由 eight 于 2007-12-6 18:39 编辑 ]

eight

原帖由 liming1985 于 2007-12-6 08:43 发表

                               
登录/注册后可看大图

下面是部分程序：
function rt=yifan(t,y,flag,n,Neb0,rou_b,L,qn0)
Kp=exp(19.76-1.537e+4/T-0.5223*log(T))
Pe=y(5)*(1-y(1)-y(2)-y(3))/(n+1+y(1)+y(2));
Ps=y(5)*(y(1))/(n+1+y(1)+y(2));
Pb=y(5)*y(2 ...

请搜索版面

yufeng

如果你自己编龙格库塔没有必要嵌入函数吧 我的qq604443022 交流一下