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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
) u: ~# [& P+ @' @%微分方程组求解主程序
clc;clear all;clf;close all;
! `0 F1 f& u, J%windows时钟自动计时
/ ^/ {6 k& a: f8 a; M  c6 rT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
t=tt(191:end,
3 r* ~+ v) r, w% S! fys=(191:end,1)
yp=(191:end,2)
yh=(191:end,3)
yr=(191:end,4)
yo=(191:end,5)
! D7 i2 O5 U% J2 W9 `  V# T; oye=(191:end,6)
0 P$ |% A3 m  D6 H4 PT2=clock;
API_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
. |: v5 n# D* q0 p* O& ~4 EAPI_elapsed_time(6)=API_elapsed_time(6)+60;
2 l/ `& c0 v3 _9 ]  Z6 jAPI_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;
7 C! p$ T* s! p. B/ ]end
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
# W* z9 {/ W3 ?API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
6 W0 G) A! `* j8 w. ]%子函数：微分方程组odefile文件并命名为lbwfun.m
8 B. y/ O+ @2 J0 Y$ dfunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
6 W* N" n5 J# ]y(1)-4.03*y(2);
0.48*y(2)-17.87*y(3);
4.85*y(3)-4.65*y(4);
2.55*y(2)+6.12*y(3)+1.95*y(4);
1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
: y" x6 T6 T1 o) Y: @- n% ^??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.
$ k. u4 ^/ v, N' P
Error in ==> odearguments at 110
+ J! S) U8 h( jf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
0 C! \- Y* m. }. F
; C$ d% w- N" p0 wError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...