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

度过无风

%微分方程组求解主程序
7 s( @/ T5 y. N7 j2 E9 V; _2 lclc;clear all;clf;close all;
%微分方程组求解主程序
/ S, m" M4 M" v) h, S) s* I* Z2 R7 Vclc;clear all;clf;close all;
%windows时钟自动计时
T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
% T9 ~8 c- j/ s( C& y- Jdisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
7 O  Y1 V+ }2 U( z2 r: L% G[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
' I) |% A: X' h9 u8 O; ot=tt(191:end,
ys=(191:end,1)
yp=(191:end,2)
yh=(191:end,3)
& e0 E  v; l  s' Z+ Nyr=(191:end,4)
yo=(191:end,5)
0 g) [, Q$ b! N* o8 iye=(191:end,6)
, t: h  k9 d, vT2=clock;
API_elapsed_time=T2-T1;
1 t) Z% z/ k3 U2 B8 Q) dif API_elapsed_time(6)<0
' S, `. A& X, C+ qAPI_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
end
' Z: }  z: S3 o! \' B0 Pif API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
0 x/ s$ w4 B$ a# o) e" L- fAPI_elapsed_time(4)=API_elapsed_time(4)-1;
  S; `2 |9 V& c6 P: @' Z  ^: p' Yend
/ x6 `5 `3 u( ustr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
2 @2 G) h% [& O5 r0 I6 c, iAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
. |# S8 m+ o# T7 b  Y: l) s1 M  h' `y(1)-4.03*y(2);
8 S) S. H& `/ G# R4 O& j0.48*y(2)-17.87*y(3);
5 `& W5 |5 Z. s4.85*y(3)-4.65*y(4);
6 w- F- a8 N: Z. c6 L$ q2.55*y(2)+6.12*y(3)+1.95*y(4);
7 y0 p* C6 P- t4 ^$ R1.10*y(2)+6.90*y(3)+2.70*y(4)];
3 s$ j% }! b# \2 s8 F, q% F运行结果（我是按书上写出来的程序）
4 a+ [1 r: c0 `# S8 O8 t  |计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.

Error in ==> odearguments at 110
: c; r$ D2 X9 Z- f9 \f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

8 _# J; Q8 U9 Q$ b2 M& W: S8 EError in ==> ode45 at 173
9 @2 q9 A2 n9 e  l1 v/ N+ G' c5 _[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...