【求助】关于ode45的问题

[baggio198310]

请各位大侠帮帮我

在chap1_2.m中的程序：
………………
………………
………………
[tt,xx]=ode45('chapl_2f',tSpan,xk,[],para);
………………
在chap1_2f.m中的程序是：
function dy =chap1_2f(t,y,flag,para)
u=para;
J=0.0067;
B=0.1;
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=-(B/J)*y(2)+(1/J)*u;

可是运行后总是报错：
??? Error using ==> feval
Undefined command/function 'chapl_2f'.

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 ==> chapl_2 at 14
[tt,xx]=ode45('chapl_2f',tSpan,xk,option,para);


敢问各位高手是什么问题？

[fanxing39]

你看下这个附件,是"清华大学大学数学实验"的讲稿,里面详细介绍了"库塔公式"的用法
我还把这个课件里要用到的matlab程序也上传了,希望对你有用!

[elunxu]

这两个语句有点问题吧
dy(1)=y(2);
dy(2)=-(B/J)*y(2)+(1/J)*u;