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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
%微分方程组求解主程序
clc;clear all;clf;close all;
%windows时钟自动计时
; N5 F3 k2 }) b; t/ }6 E6 V  b6 `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,
ys=(191:end,1)
yp=(191:end,2)
yh=(191:end,3)
yr=(191:end,4)
7 c0 E, [/ K0 L4 S9 X  y  {4 O( xyo=(191:end,5)
6 y, k3 g# z, B: m  G4 N$ iye=(191:end,6)
T2=clock;
API_elapsed_time=T2-T1;
* }& c: L: l  O& g( z9 T$ Q, Gif API_elapsed_time(6)<0
/ T" B* h( L+ w4 A0 eAPI_elapsed_time(6)=API_elapsed_time(6)+60;
# B! {& f, A: E+ i4 lAPI_elapsed_time(5)=API_elapsed_time(5)-1;
end
if API_elapsed_time(4)<0
! {8 h/ ^8 d+ P5 J  Q, VAPI_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
1 `$ i4 l1 b& k0 d) Hend
) q# \2 k8 V7 n# r  {: Qstr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
8 E# `' A' n' i( A3 I/ bdisp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
8 z, y9 S0 I% G: ^) Q$ c1 }: G: }function ydot=lbwfun(t,y);
- [& j$ A# i! K. \ydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
7 W3 X5 u/ |+ n" d0.48*y(2)-17.87*y(3);
# N1 W% H4 Y0 a. t& f4.85*y(3)-4.65*y(4);
2.55*y(2)+6.12*y(3)+1.95*y(4);
/ E: O% N( W- |, D1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
7 D* g) R$ U- r7 e计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
  B' y0 u( V( r1 E# I  f+ I. E??? Error using ==> feval
! ?3 M; |  ?. a, I6 r7 Y% K# fUndefined function or method 'lbwfun' for input arguments of type 'double'.
7 X, `& P# ~; {: |3 @
Error in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
- P7 \& j# J+ s/ f: T
Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...