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

度过无风

%微分方程组求解主程序
7 w6 Y/ P, C5 T* mclc;clear all;clf;close all;
%微分方程组求解主程序
clc;clear all;clf;close all;
%windows时钟自动计时
5 y/ W5 K; v" o0 b3 uT1=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)
yp=(191:end,2)
# A5 h/ U+ D+ X, Zyh=(191:end,3)
3 E- d4 ?8 d3 `yr=(191:end,4)
yo=(191:end,5)
: N! W- n: H' o* pye=(191:end,6)
T2=clock;
6 i* I' R. M, `' b. B' p1 JAPI_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
  }3 d, P. U- zAPI_elapsed_time(5)=API_elapsed_time(5)-1;
end
6 Z( E! E4 [* U' a7 Fif 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 秒',```
8 z7 j6 n, o; P' }API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
+ ]( W: c; c7 N- b) A5 V  C4 H3 Qdisp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
" ]' H: A: B7 K3 M& ?) _% Hfunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
; e% {9 t9 B% X. jy(1)-4.03*y(2);
1 h2 ]8 V+ H& [0.48*y(2)-17.87*y(3);
1 i# b. U. w- A0 U. f1 T4 Q; {4.85*y(3)-4.65*y(4);
7 ^4 P3 G* q' g; b+ k2 q2.55*y(2)+6.12*y(3)+1.95*y(4);
1.10*y(2)+6.90*y(3)+2.70*y(4)];
; ]! h: i2 X0 n$ G: m( M- s" d' g运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
5 v$ `; l3 z$ Z1 N: z6 hUndefined function or method 'lbwfun' for input arguments of type 'double'.

Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
  |% p- D6 N7 W) N8 w% e
7 C& m' b/ {5 D1 B# E& fError in ==> ode45 at 173
4 b, R2 G& U3 [. T" B9 W/ }* m[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...