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

度过无风

%微分方程组求解主程序
clc;clear all;clf;close all;
# j* H- b0 e" k5 t%微分方程组求解主程序
clc;clear all;clf;close all;
( _+ Y0 \7 ^* }3 {* N1 ?%windows时钟自动计时
& W0 `( E/ o9 {9 a0 Q7 E1 z$ X; gT1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
& ^2 ]" S7 X. T, ^[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
) h- A2 `. ?+ ct=tt(191:end,
ys=(191:end,1)
yp=(191:end,2)
yh=(191:end,3)
, Q0 z8 c$ c) t% R2 L/ Ayr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
& i5 O, `' z" l% d  p/ V/ d" z5 k  YT2=clock;
2 n; `8 p) W. T+ F( [- ~API_elapsed_time=T2-T1;
$ w* c7 Y1 x; rif API_elapsed_time(6)<0
API_elapsed_time(6)=API_elapsed_time(6)+60;
# i  \8 r7 _7 @0 p) M5 R' @1 sAPI_elapsed_time(5)=API_elapsed_time(5)-1;
3 {1 m; {2 S4 r* M. ]6 z6 Z' rend
if API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
2 U  x, w" x3 U$ N" s; F& Lend
, d4 T  C9 o7 s* Z" Vstr=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
/ q/ B& F( C1 X, f  g0 S) x( c3 fdisp(str);
- J" y/ T! g  b$ t) w%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
y(1)-4.03*y(2);
0.48*y(2)-17.87*y(3);
4.85*y(3)-4.65*y(4);
+ n: r- E/ {/ u$ H( B% S# L2.55*y(2)+6.12*y(3)+1.95*y(4);
1.10*y(2)+6.90*y(3)+2.70*y(4)];
( ?* p- I) E- S0 b运行结果（我是按书上写出来的程序）
( P* J6 k+ Q) B' q' X计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
; n8 e: ?: ?% H: k" V1 l4 b??? Error using ==> feval
2 a, k" \# d/ q  LUndefined function or method 'lbwfun' for input arguments of type 'double'.
; L! R! x2 q3 D9 p7 O
" q4 K, o1 F# Y2 A) X/ }# l% j: G. qError in ==> odearguments at 110
f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

! O. y: K8 I) U# J; jError in ==> ode45 at 173
  c3 Y7 {5 i) I7 q8 d" K[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...