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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
- A) b8 _$ \6 @4 x# Z& c! i; Z%微分方程组求解主程序
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]);
t=tt(191:end,
! Z: W. h0 g* F2 Iys=(191:end,1)
yp=(191:end,2)
. s$ l2 J: ?' G1 G. J9 dyh=(191:end,3)
  L) o0 d" H, V8 i, Vyr=(191:end,4)
. O5 B& Y0 [: dyo=(191:end,5)
ye=(191:end,6)
T2=clock;
API_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
3 a! o1 F5 }* w* a0 n1 nAPI_elapsed_time(6)=API_elapsed_time(6)+60;
1 s5 E9 H8 ?& DAPI_elapsed_time(5)=API_elapsed_time(5)-1;
2 Y3 l# Z2 H9 S: m& Jend
& M. P- V5 f( c* \0 T+ c) u! Wif API_elapsed_time(4)<0
- v; E( b7 |, I& r2 @API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
end
9 N7 ^) d7 {5 a6 o* Z+ _6 qstr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
, w% u3 c. ~7 A7 m+ mdisp(str);
5 g9 d/ ^) m, l%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
- X' `" N; a( Y2 H, Pydot=[121.793*pi*cos(2*pi*t);
' G% h3 s# ~4 Qy(1)-4.03*y(2);
6 N; z# T% o- c* e: q1 j2 u0.48*y(2)-17.87*y(3);
: Z9 c6 h7 I5 v' N$ m6 M3 |4.85*y(3)-4.65*y(4);
3 N7 k1 r. u( E7 N0 _/ S2.55*y(2)+6.12*y(3)+1.95*y(4);
1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
; y3 h. P. n/ Z% S$ t5 S! z  P  H2 W5 g??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.
, [" G9 _# o9 Y# k' `$ x; L4 A
Error in ==> odearguments at 110
( g$ `" O7 X8 m1 zf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
  ?! ]5 e( p0 t5 g8 I" {  p) r2 D
  B2 G* s# M# D2 c3 h5 H. B+ }" {& CError in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...