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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
%微分方程组求解主程序
, P. @7 C/ j1 @9 bclc;clear all;clf;close all;
%windows时钟自动计时
1 P; R9 k/ k( W  x- I4 }3 @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,
6 E8 J2 d, I0 V) m0 R$ cys=(191:end,1)
yp=(191:end,2)
/ j( `$ j# ?5 i) ]+ m) Fyh=(191:end,3)
% B, O( d. ~) z& ?yr=(191:end,4)
+ @9 {) d- m; l  \4 Z8 ]+ Wyo=(191:end,5)
ye=(191:end,6)
T2=clock;
8 _# [3 ~$ h$ HAPI_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
2 @0 Z/ E" r8 zend
3 E5 V  @# ?8 z/ `* \if API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
" q6 d  f& o) g( A& qAPI_elapsed_time(4)=API_elapsed_time(4)-1;
end
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
9 E% Q0 t/ I* r# b$ C" G% J* EAPI_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
. r/ X8 y1 K! s" ?5 vdisp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
. S( a* S2 O' L5 f/ C- U. {function ydot=lbwfun(t,y);
) b0 Y" C. _0 Y! W% y6 eydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
5 g2 L# R8 U7 X0 a1 r0.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);
) u5 D5 D; v9 N1.10*y(2)+6.90*y(3)+2.70*y(4)];
3 g" t9 E% d9 R9 E. @% I运行结果（我是按书上写出来的程序）
% R) |  @' ?! G6 X计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
) r0 p3 @0 k( i" R9 `Undefined function or method 'lbwfun' for input arguments of type 'double'.
: ~' h# S3 A! Z; n' r0 M, @& M5 O
; b+ I; a) N/ ~3 J( W$ U6 HError in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...