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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
%微分方程组求解主程序
& P; ^0 b+ p3 J: {# rclc;clear all;clf;close all;
%windows时钟自动计时
8 {* R5 n: p4 j" d( IT1=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]);
t=tt(191:end,
ys=(191:end,1)
+ t; T1 B9 N; P0 G( ayp=(191:end,2)
2 ~+ q: ]( I8 F2 \) F& L! dyh=(191:end,3)
yr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
T2=clock;
+ z/ f3 @& Z8 R4 ^( A6 IAPI_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
; D' Y9 E9 g* }; d5 e! w; wAPI_elapsed_time(5)=API_elapsed_time(5)-1;
end
& @4 I: C6 U9 Cif API_elapsed_time(4)<0
9 {. [$ x( _2 N+ _7 FAPI_elapsed_time(4)=API_elapsed_time(4)+60;
$ r( V9 Y) ^& H0 N. Y- x2 V. fAPI_elapsed_time(4)=API_elapsed_time(4)-1;
) c, m8 z) X; C' ~end
+ f: Y  N; a, K1 S/ r* tstr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
1 ^9 O8 r3 S0 m7 z. Odisp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
9 }3 e* u! k3 J5 N- P- py(1)-4.03*y(2);
0.48*y(2)-17.87*y(3);
3 i' ^3 C+ d2 i; H: n4 `! a1 o4.85*y(3)-4.65*y(4);
: ]- w- R* @" R6 m3 h$ i2.55*y(2)+6.12*y(3)+1.95*y(4);
1.10*y(2)+6.90*y(3)+2.70*y(4)];
2 C: d7 V* L9 ]  J- y9 ?# E/ N运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
3 M2 q% o1 d! |. D8 l3 x5 `9 Q??? Error using ==> feval
; i, L' X9 K7 t+ b  C2 dUndefined function or method 'lbwfun' for input arguments of type 'double'.
) w2 j9 p$ c! u9 j
Error in ==> odearguments at 110
* ^, h9 K; j! W3 ~f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
& |) n; u8 `2 O$ ^) [4 h
  {% _, S$ J- LError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...