湖泊有机物新陈代谢模型问题

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
, |2 I2 ^/ I6 y: k; S0 X% F- a%微分方程组求解主程序
5 b6 p8 @( J% b/ ]4 r6 Nclc;clear all;clf;close all;
! K0 T% F8 L% @! N3 v%windows时钟自动计时
T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
: h- f( f+ }+ m& l( K) idisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
t=tt(191:end,
# T- ~) A% o( Dys=(191:end,1)
7 G- L8 I2 U* z1 C4 \' [, `yp=(191:end,2)
yh=(191:end,3)
yr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
5 R  A8 D3 B# Z( S& o4 A- b$ nT2=clock;
API_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
" s( n6 K' k2 T! t5 C* UAPI_elapsed_time(5)=API_elapsed_time(5)-1;
end
if API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
) N2 K' x( F( e, j0 ?7 \% uend
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
" Q9 o  Q: [: Z6 d2 zAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
% k4 Y# y4 X& Z( y9 v9 {%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
0 h+ y+ Z) R$ U+ k/ u+ o( W, H* ]0.48*y(2)-17.87*y(3);
0 Y4 l$ ?  K3 V* `8 J4.85*y(3)-4.65*y(4);
8 d5 a# f+ T/ h6 X9 ]2.55*y(2)+6.12*y(3)+1.95*y(4);
/ c' V* i7 Q+ a+ x/ I1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.

0 S0 X4 m: L+ v9 b+ \7 _+ zError in ==> odearguments at 110
+ H; X! [' a! q# ]! J9 vf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
1 t$ g5 b. X' l
, f, `& t" E8 ^7 {9 d3 }5 gError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...