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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
%微分方程组求解主程序
& t9 N: V* b' H: p% R% G1 Rclc;clear all;clf;close all;
1 q, ]! `& G6 d# j7 E0 M%windows时钟自动计时
9 c3 h# \* j) U  V! ]! q+ L8 NT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
6 o8 C$ c; g$ y' m5 l) z% zdisp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
/ J: ^6 l4 y! H% W- It=tt(191:end,
ys=(191:end,1)
! ~1 s8 h# a' E0 P5 L$ `: J: jyp=(191:end,2)
yh=(191:end,3)
yr=(191:end,4)
6 m. }" i" n. [: L9 @# C! v6 myo=(191:end,5)
ye=(191:end,6)
T2=clock;
, `7 h( L. d* {% {" v. d# {API_elapsed_time=T2-T1;
7 m5 N1 m$ N- A0 `7 Oif API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
1 E2 A* h; o0 S0 FAPI_elapsed_time(5)=API_elapsed_time(5)-1;
end
if API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
' M; q  g0 Y; P& E1 Bend
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
- L# I# C8 [7 R0 P% |disp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
4 ?& Z* z- F6 i) u8 p) Pydot=[121.793*pi*cos(2*pi*t);
, N' D. v, k; p5 T" B9 K- `8 c, xy(1)-4.03*y(2);
8 ^- D9 V% W9 s; y8 I1 _& S4 L0.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);
6 y( o" A* Y" j8 L8 ]! u" _  V1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
) a- F9 [" j+ B4 ~+ i* k  f( b计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
9 K% K+ v% \4 B8 Y3 p. H. N: s: \$ C??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.
6 d$ U8 e7 X# n! U4 v
7 r6 e: n. R! p* a4 ~Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

: }# O' \& h; I. u- G, r5 WError in ==> ode45 at 173
6 @9 d6 P) Q8 M2 L! Q% P[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...