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

度过无风

%微分方程组求解主程序
+ e& _( i0 j1 `+ ]clc;clear all;clf;close all;
%微分方程组求解主程序
( g8 x- x3 n/ a# A0 [4 M6 Tclc;clear all;clf;close all;
% x  m* o% ]! Q- j; ~1 m%windows时钟自动计时
T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
, B$ w$ C7 ~1 Bdisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
( ^% R8 @& k1 C" \3 B% et=tt(191:end,
6 G3 c5 P) @' o  A  y3 Dys=(191:end,1)
yp=(191:end,2)
yh=(191:end,3)
yr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
+ ^# m) N7 z0 X" g: r5 r; E* }T2=clock;
API_elapsed_time=T2-T1;
$ r" r' u% J" s; ~  o, Jif API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
end
( X5 I* y& m7 K: {7 K" g- Yif API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
0 x) f5 L1 [$ r" z+ |# o9 z( t, X. HAPI_elapsed_time(4)=API_elapsed_time(4)-1;
end
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
; F% Z. B" v7 K1 l3 QAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
: D# c$ v9 x5 V2 G%子函数：微分方程组odefile文件并命名为lbwfun.m
, k2 M( j- P& {" f) X1 j2 Wfunction ydot=lbwfun(t,y);
7 y2 ~( {$ N7 U# X& u" c" Oydot=[121.793*pi*cos(2*pi*t);
9 w3 y) B4 C* l. x* uy(1)-4.03*y(2);
9 q8 ~4 s: b6 L2 O& \! I+ _0.48*y(2)-17.87*y(3);
* J. o: ?. n: X# C: y; j# q. Q4.85*y(3)-4.65*y(4);
: o! j$ a7 I# }, \2.55*y(2)+6.12*y(3)+1.95*y(4);
: d1 r+ g* y4 H+ ]3 g7 C/ A6 a1.10*y(2)+6.90*y(3)+2.70*y(4)];
: \, H) Y* C* R% Q' K5 t运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.
( V5 R1 B# g" V( L* a/ ~
Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
! d' H1 ~0 w# e! k- F1 ~[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...