生态系统动力学仿真matlab程序

disp（‘fox1代表幼年狐狸数量，fox2代表成年狐狸数量。rabbit1代表幼年兔子数量，rabbit2代表成年兔子数量‘）;
%f1f2rab1rab2分别代表fox1、fox2、rabbit1、rabbit2的初值。
f1=input（‘请输入幼年狐狸的初值:\n‘）;
f2=input（‘请输入成年狐狸的初值:\n‘）;
rab1=input（‘请输入幼年兔子的初值:\n‘）;
rab2=input（‘请输入成年兔子的初值:\n‘）;
fox1（1）=f1;fox2（1）=f2;rabbit1（1）=rab1;rabbit2（1）=rab2;
%设定各个流率及常数:R1代表狐狸出生率，C1为狐狸出生率常数，R2代表成年狐狸的增加率，C2为幼狐成长延滞时间
%R3为成年兔子被猎杀率R4为兔子的出生率，C3为兔子出生率常数，R5为成年兔子的增加率，C4代表幼兔成长延滞
%P3代表因为被吃导致成年兔子的死亡率C5为幼狐的食量，C6为成狐的食量R6为幼年兔子缺乏抚养和保护的死亡率并假设初值为零
%设定各个流率及常数的初值。
C1=1/24;
R1（1）=fox2（1）*C1;
C2=8;
R2（1）=fox1（1）/C2;
C3=0.4;
R4（1）=rabbit2（1）*C3;
C4=3;
R5（1）=rabbit1（1）/C4;
C5=10;
C6=60;
P3（1）=fox1（1）*C5+fox2（1）*C6;
R6（1）=0;
for i=1:130
    P5（i）=0;
end
%进入循环控制语句。
for i=2:150
%流位L1、L2、T1、T2的迭代语句。
  if i>120
    R3=7;
  else
    R3=0;
  end
fox1（i