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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
%微分方程组求解主程序
clc;clear all;clf;close all;
+ c' R) U) W# v%windows时钟自动计时
4 N( H  w- n( U3 _4 k# ~7 MT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
: R: S+ J( F* ]% z/ u7 x; y1 T[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
t=tt(191:end,
ys=(191:end,1)
yp=(191:end,2)
. G/ q( e. r3 L# v$ gyh=(191:end,3)
6 a: x) O/ E' t3 X0 C0 N) M) |4 F8 Xyr=(191:end,4)
, l& k" y/ X6 j  @* D+ Uyo=(191:end,5)
0 h7 @. T8 R, w* u! O+ w9 e: Z7 Wye=(191:end,6)
T2=clock;
( a* b1 Y; m3 L, wAPI_elapsed_time=T2-T1;
: O1 c) `. S" P% j8 u- K. ^if API_elapsed_time(6)<0
& I2 A9 `. X( ?( ~2 ]' B* KAPI_elapsed_time(6)=API_elapsed_time(6)+60;
9 ~* J# t0 n+ Z* l8 q3 CAPI_elapsed_time(5)=API_elapsed_time(5)-1;
5 J# @  ^# E! E# w9 O, c4 s  Uend
if API_elapsed_time(4)<0
4 F& I8 ~0 a! {# a, \, nAPI_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
8 v( z5 W4 p: ]6 j% U8 Xend
2 E- h8 k& \+ `6 Z; T. \str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
6 [# |5 `/ N: edisp(str);
# ~: H7 H& r- I1 f%子函数：微分方程组odefile文件并命名为lbwfun.m
8 c7 K/ `% K8 H/ n- \! Gfunction ydot=lbwfun(t,y);
/ T& p1 v1 I' c& Q+ e0 [ydot=[121.793*pi*cos(2*pi*t);
# d4 I0 Z. x. L# j" A: {# w! Jy(1)-4.03*y(2);
* k- X5 [2 w* a/ E$ P0.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)];
运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
; P: y7 Z) p5 a; c% w4 xUndefined 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.
1 i7 c- s( P8 ^0 }+ U9 o! A
Error in ==> ode45 at 173
3 R3 [/ H: `; v3 }4 W+ F0 h[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...