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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
' L8 N/ d/ K9 G7 ]6 J! Z9 e- H%微分方程组求解主程序
clc;clear all;clf;close all;
%windows时钟自动计时
T1=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)
& ~% \) z! D5 O/ yyp=(191:end,2)
% h$ P1 ?9 h* l3 Oyh=(191:end,3)
yr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
) }9 s( g1 Q' `: g) OT2=clock;
API_elapsed_time=T2-T1;
) E$ S5 b& @+ T0 \( {if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
$ s8 y1 V4 R, Eend
1 F$ P, e) v- L; r- y2 Vif API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
; ~% n8 w1 j  a8 t# [' sAPI_elapsed_time(4)=API_elapsed_time(4)-1;
end
+ f: R" K. F) H' ?str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
) y7 @8 E# a8 Y5 m" Ddisp(str);
9 V. Z4 D& ]+ U- a+ g2 c! J$ z%子函数：微分方程组odefile文件并命名为lbwfun.m
) Y, R: ]8 @4 s7 {: a# M6 ^6 d7 Ifunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
9 T2 l5 E" a- e' `* M/ v  ?y(1)-4.03*y(2);
" W/ v5 ]. ~* ~4 o1 B0.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)];
; _( `4 D; b1 E! I运行结果（我是按书上写出来的程序）
# R  \! W0 @; h1 ]计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
, }+ X/ c, X' M??? Error using ==> feval
% L$ ]9 n# K% I% V7 X( cUndefined function or method 'lbwfun' for input arguments of type 'double'.

* F3 W3 D6 t+ HError in ==> odearguments at 110
' `' O6 O  s3 T" L2 {) uf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
, l  u/ T4 D+ n5 m; E% {( D5 e
Error in ==> ode45 at 173
: D1 @) Z7 q7 K[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...