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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
- Q. D) Y; F( x& L  c%微分方程组求解主程序
clc;clear all;clf;close all;
%windows时钟自动计时
& l' p! z. P1 R3 l' ]  MT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
7 x( w8 H+ w3 m. A& J3 R( Kdisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
0 C6 |6 {& A) b1 m- [[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
! _$ Q" V; U& r1 D9 [t=tt(191:end,
1 U( O& a9 \4 m( Z6 fys=(191:end,1)
; q* y/ H' Z& Z4 f# o5 typ=(191:end,2)
yh=(191:end,3)
1 i+ M/ ]3 ?0 @4 u2 _yr=(191:end,4)
+ l  `7 t; c9 v+ ^8 byo=(191:end,5)
ye=(191:end,6)
T2=clock;
API_elapsed_time=T2-T1;
/ b# m8 z7 F& ?/ q8 _if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
% l5 i0 O0 P- x9 o! JAPI_elapsed_time(5)=API_elapsed_time(5)-1;
end
$ B- c  M  F" V8 Y- rif API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
end
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
  t( Q1 G: {0 ^4 b, E$ iAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
- v# b" z% F! g/ ~/ C; Ndisp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
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)];
运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
( g$ U' i  p+ Z. o% s1 |! z) _7 iUndefined function or method 'lbwfun' for input arguments of type 'double'.
" L7 p8 `% Y( Q( w
1 o, [$ Y+ l( C) n! XError in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
3 b4 ]) \/ H. G1 \  \' }
Error in ==> ode45 at 173
( t0 j1 l6 F. y! s" r8 `( m[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...