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

度过无风

%微分方程组求解主程序
2 x! g' o4 Z7 a, ~& oclc;clear all;clf;close all;
%微分方程组求解主程序
& t; W# T+ B4 L8 e, T0 Vclc;clear all;clf;close all;
+ {. v" ?- `: G& o: `%windows时钟自动计时
# |7 f) {' `7 x1 `% ?T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
$ S) E, E) V4 w' q" J[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)
/ c8 M, c$ J1 Y# yyp=(191:end,2)
yh=(191:end,3)
6 E9 f) N+ }: a; Wyr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
' G) N( b1 Y( C, `: ]% WT2=clock;
API_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
' o- D! e2 v& e& w+ C7 F& tAPI_elapsed_time(6)=API_elapsed_time(6)+60;
- u/ N" s" ?& H; d" ~8 hAPI_elapsed_time(5)=API_elapsed_time(5)-1;
3 k0 I' `* @9 X2 Tend
5 N! ]% p( ?2 a3 l. [. tif API_elapsed_time(4)<0
0 d, k% `+ C, O( PAPI_elapsed_time(4)=API_elapsed_time(4)+60;
API_elapsed_time(4)=API_elapsed_time(4)-1;
! j0 P% A  S; c. e8 J- rend
str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
! ]% u; t* r6 h7 k3 N4 x: l  ^API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
disp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
9 h7 A$ j$ G+ G$ @) A+ n3 Qfunction ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
& F9 U& h4 U' r$ T- q% Py(1)-4.03*y(2);
0.48*y(2)-17.87*y(3);
4.85*y(3)-4.65*y(4);
/ o4 S3 {: |# \6 i& y2.55*y(2)+6.12*y(3)+1.95*y(4);
1.10*y(2)+6.90*y(3)+2.70*y(4)];
运行结果（我是按书上写出来的程序）
计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
??? Error using ==> feval
0 U* t, `/ Y2 C3 n, G7 u4 wUndefined function or method 'lbwfun' for input arguments of type 'double'.

& V7 v- k2 y& I, JError in ==> odearguments at 110
' f" v: g% N! {: u( p2 |5 c. |f0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.

* r3 G" p( d5 FError in ==> ode45 at 173
; K) Q$ j/ L8 ~0 w3 N+ g6 r: U7 ^; y3 |[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...