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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
: S1 n* ?( c: Y: P! O( ^4 w, Q%微分方程组求解主程序
, d& @' F" |; d9 x6 n8 Qclc;clear all;clf;close all;
%windows时钟自动计时
T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
' @# L/ r' }, k% u. I9 pdisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
t=tt(191:end,
' ~$ S. d( D& a/ cys=(191:end,1)
yp=(191:end,2)
  }6 @8 n5 r: }$ x5 iyh=(191:end,3)
yr=(191:end,4)
' W# r+ O! u. E, G! e, q/ @yo=(191:end,5)
+ z+ |: O0 P8 y; |: u/ j6 {ye=(191:end,6)
2 v3 P- L" y" ~+ c3 ~) iT2=clock;
API_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
  \# `( X; M) tAPI_elapsed_time(5)=API_elapsed_time(5)-1;
& t" c. o/ \9 g4 Vend
. F- X! |- r" F: Nif API_elapsed_time(4)<0
* b9 A! Z* R% S) }API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
end
2 |, ~2 U: X' u6 E/ ?, f1 |str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
' A7 ^! v+ i5 V+ b. {+ J0 ]* m( j6 RAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
% w; I: y; }! ]0 e9 V+ N, Mdisp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
. I# D. T: i, K$ C2 E8 V7 Hfunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
: I! v- D. m! o# Y$ O( y0.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);
+ |( f, K9 A7 C$ f0 F; R1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
, c* L- H: I% |$ I" w. ^Undefined function or method 'lbwfun' for input arguments of type 'double'.
2 v+ g: X! }- V/ T! V
: y1 ~3 b1 L3 \  _3 E7 _: LError in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
4 n; q, w7 j, N3 ]6 ^! X+ S
! W, k, ^$ _% ?9 y" zError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...