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

度过无风

%微分方程组求解主程序
) _! ]+ W1 c9 W5 R- D  r% g* Oclc;clear all;clf;close all;
$ g0 _( X& f+ Q0 N1 X%微分方程组求解主程序
3 @, k* X0 H* A: V% c( h6 _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]);
5 Q* G: I8 h. e% }% Ht=tt(191:end,
& m$ `" L( k; ^: F- j5 y/ M, G: {ys=(191:end,1)
. L: }2 U# s3 ~; z3 u4 ^yp=(191:end,2)
, f+ P3 _1 ~& c& c0 I  z+ Z& {yh=(191:end,3)
yr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
T2=clock;
API_elapsed_time=T2-T1;
) A( k; |3 M6 M/ N: H8 E* uif API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
# E; d6 T# ^# s" @, Hend
, ^9 E$ i- O2 h4 |if API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
API_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);
%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
/ `- v! |) {8 r, ny(1)-4.03*y(2);
- u4 n% G5 }* |- c0.48*y(2)-17.87*y(3);
/ a7 J1 H  V) ~8 ?( E4.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)];
. a- a; M! _' }  i* o& M4 c1 @运行结果（我是按书上写出来的程序）
0 v( J# B) j, ?$ e; t' j# a计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
5 E. S/ l7 _2 J5 ?  f% @Undefined function or method 'lbwfun' for input arguments of type 'double'.
9 ^: z7 _7 S: H5 v6 {! B
" U7 K3 F" L& TError in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
" L) T' l$ N* Q0 u( Q
( p6 `6 R7 p5 |' O. I9 I3 GError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...