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

度过无风

%微分方程组求解主程序
3 x2 C8 C' A, uclc;clear all;clf;close all;
- H, g! J9 U. X$ p; s4 {%微分方程组求解主程序
clc;clear all;clf;close all;
0 Q# H! }; @1 Q; p8 y! M4 E%windows时钟自动计时
T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
7 p4 B+ X, M& U! A6 C8 M: v[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
5 x  X: ]. @7 O  tt=tt(191:end,
7 k; p: j. D' A6 Nys=(191:end,1)
" a; g+ K% _; a  Y( Y6 Yyp=(191:end,2)
yh=(191:end,3)
2 @2 B& q& g$ O1 G0 A6 Pyr=(191:end,4)
yo=(191:end,5)
, h8 u6 ^: }$ C7 r, L9 gye=(191:end,6)
1 g( Z- j9 t! @2 i, WT2=clock;
/ k1 t) r. }8 B; I: AAPI_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
( _# b+ J$ B' T" lAPI_elapsed_time(6)=API_elapsed_time(6)+60;
API_elapsed_time(5)=API_elapsed_time(5)-1;
end
4 U) a' y! h/ z6 Aif API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
end
# x+ B) C% u5 X( O4 o- Sstr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
' S5 u. s4 E1 `6 pfunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
+ [6 \* X* O9 K: ~; q# uy(1)-4.03*y(2);
0.48*y(2)-17.87*y(3);
1 ^9 q9 c) A) n5 ?, f4.85*y(3)-4.65*y(4);
2.55*y(2)+6.12*y(3)+1.95*y(4);
9 g' G% o) ~' l$ F$ }: y1 ]6 |1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
2 A+ g$ B5 s& d  T4 v6 r. R6 r; f计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
Undefined 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.

Error in ==> ode45 at 173
& q+ R# p  N9 s* `% X4 Z5 r- r. g[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...