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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
' v- n, l; h; Y8 l%微分方程组求解主程序
clc;clear all;clf;close all;
0 A" Q% M- l8 k: z%windows时钟自动计时
& F6 e+ Y% }6 o$ V! C2 i$ s1 uT1=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]);
" f2 ]9 n+ f% L- s5 \t=tt(191:end,
ys=(191:end,1)
) j- z6 g& x# Y+ s; t! U4 Byp=(191:end,2)
yh=(191:end,3)
yr=(191:end,4)
# K% G3 f7 I* Q9 V8 lyo=(191:end,5)
ye=(191:end,6)
T2=clock;
3 L$ C) n5 K! J$ H( V5 U) F, k. LAPI_elapsed_time=T2-T1;
2 E5 L6 B- U8 f" Qif API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
$ I: c5 K# q% v5 CAPI_elapsed_time(5)=API_elapsed_time(5)-1;
end
! i# Y! n* B: i$ y7 |$ G' A  h, gif API_elapsed_time(4)<0
" Y/ q. q* v. l, \2 ^( dAPI_elapsed_time(4)=API_elapsed_time(4)+60;
7 F! x5 v5 j/ _API_elapsed_time(4)=API_elapsed_time(4)-1;
end
, e0 m6 ?* P" j6 O, q2 B) r; l" \str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
* e, k/ S6 i+ j) ^" |3 g/ rdisp(str);
% \& b8 b: A: W%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
# ?* v$ E8 T$ z! w" M- u/ Vydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
0.48*y(2)-17.87*y(3);
  {$ d9 B: w, }) I5 G0 ?4.85*y(3)-4.65*y(4);
0 L* S. r: Z; k  Q2.55*y(2)+6.12*y(3)+1.95*y(4);
! S  M3 W- `5 O- c: t0 W3 L1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
3 l9 M+ Y" a" G1 e计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.

Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

# D6 b( y. \8 d5 oError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...