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

度过无风

%微分方程组求解主程序
0 L. n5 Y5 Q, Q) Q/ v; L) p; ]clc;clear all;clf;close all;
%微分方程组求解主程序
clc;clear all;clf;close all;
%windows时钟自动计时
T1=clock;      %函数clock返回的值是clock=[year month day hour minute seconds]
1 e4 P  d  f( ]6 u1 H# H% U8 ~disp('计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······');
[tt,y]=ode45('lbwfun',[0:10:2020],[95.9,0.83,0.003,0.0001,0.0,0.0]);
) e# F% z) D8 Q( q3 pt=tt(191:end,
2 {, P* X% Y# f" n1 ?5 @7 y: @ys=(191:end,1)
yp=(191:end,2)
3 z  X. f8 E; s8 I1 @yh=(191:end,3)
% ?' s) F; Y3 Xyr=(191:end,4)
yo=(191:end,5)
ye=(191:end,6)
T2=clock;
; C  j) u; j, e! Z$ ^- p, G$ |5 yAPI_elapsed_time=T2-T1;
if API_elapsed_time(6)<0
. \* W& n* ~# k' gAPI_elapsed_time(6)=API_elapsed_time(6)+60;
+ P! L  `2 c1 E. U( S/ f% I+ s5 xAPI_elapsed_time(5)=API_elapsed_time(5)-1;
  U' K7 S1 Z+ y. ^: cend
6 F8 S+ c0 z7 R1 Bif API_elapsed_time(4)<0
API_elapsed_time(4)=API_elapsed_time(4)+60;
! Q0 F# z3 w3 j* E+ {API_elapsed_time(4)=API_elapsed_time(4)-1;
/ d) {* T4 D8 Dend
; s2 ~0 {2 h, _str=sprintf('湖泊新陈代谢模拟程序共运行 %d 小时 %d 分钟 %.4f 秒',```
API_elapsed_time(4),API_elapsed_time(5),API_elapsed_time(6));
2 [. ^- F2 t4 adisp(str);
%子函数：微分方程组odefile文件并命名为lbwfun.m
function ydot=lbwfun(t,y);
ydot=[121.793*pi*cos(2*pi*t);
7 J0 o2 J% |) Q# G/ @" Q/ v3 ly(1)-4.03*y(2);
0.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);
1.10*y(2)+6.90*y(3)+2.70*y(4)];
+ u+ o/ F1 B: z运行结果（我是按书上写出来的程序）
! e! V9 W9 F& ?% @. l计算机正在准备输出湖泊有机物新陈代谢结果，请耐心等待······
3 \7 z" |1 n( g  D* q??? Error using ==> feval
Undefined function or method 'lbwfun' for input arguments of type 'double'.

Error in ==> odearguments at 110
- d0 Z. D4 H( a! t: ~4 Hf0 = feval(ode,t0,y0,args{:});   % ODE15I sets args{1} to yp0.
% |: V: I  u1 J! R3 s
0 d8 ]0 E) _  D3 H5 R! _Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, odeFcn, ...