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

度过无风

%微分方程组求解主程序
9 g  o* |2 U4 N) H2 Y# z. u0 oclc;clear all;clf;close all;
%微分方程组求解主程序
5 O8 \# M. s$ M$ hclc;clear all;clf;close all;
%windows时钟自动计时
- _8 t0 X' f6 k* `% V$ Z8 q. AT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
& W/ ?9 y/ b  xdisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
: |' |+ x* P6 o8 n* w[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
) S" }" T6 L  q  c' N+ T& Y+ ft=tt(191:end,
) Q) J0 B( I$ s# y; D9 ]; K! @ys=(191:end,1)
yp=(191:end,2)
yh=(191:end,3)
* O! m" V% R' O4 |" O$ d+ s7 Zyr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
T2=clock;
API_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
6 h0 q& }0 M. a1 B6 @$ L* MAPI_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
end
8 X7 X( z- n0 O& a' @6 r4 l* Q4 H9 A7 bif API_elapsed_time(4)<0
3 V3 `6 \1 q4 MAPI_elapsed_time(4)=API_elapsed_time(4)+60;
  w. R8 F1 n. oAPI_elapsed_time(4)=API_elapsed_time(4)-1;
end
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
4 v) K; C5 f# h; Z1 r%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
: f) @& `: Q4 V, H1 ^- D0 W* ^ydot=[121.793*pi*cos(2*pi*t);
5 O& O5 c$ A# P# dy(1)-4.03*y(2);
, |* }+ g1 P+ d/ q2 `0.48*y(2)-17.87*y(3);
4.85*y(3)-4.65*y(4);
. \8 n9 O, e- K! O2 v! q, F2.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
( u8 S8 j" U5 ^3 ]2 m  y6 w5 QUndefined function or method 'lbwfun' for input arguments of type 'double'.

Error in ==> odearguments at 110
) V% g; L# R) }% M! q( Zf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
$ y/ H! `, H- {/ Q0 t$ u
) ]& I  r9 u+ }4 s  IError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...