INFLOWS:
BIRTH_PREY_1 ¼ (PREY_1ÀDEATH_PREY_1) * BR_PREY_1 * (1ÀPREY_1/
CC_PREY_1)
IM_PREY_1 ¼ PREY_LF_2+PREY_UP_3
OUTFLOWS:
DEATH_PREY_1 ¼ CONSUME__1
MIG_PREY_1 ¼ MIG_RATE_PREY * (PREY_1ÀDEATH_PREY_1)
PREY_2(t) ¼ PREY_2(t À dt) + (BIRTH_PREY_2 + IM_PREY_2 À
DEATH_PREY_2 À MIG_PREY_2) * dt
INIT PREY_2 ¼ 0
INFLOWS:
BIRTH_PREY_2 ¼ (PREY_2ÀDEATH_PREY_2) * BR_PREY_1 * (1ÀPREY_2/
CC_PREY_2)
IM_PREY_2 ¼ PREY_RT_1+PREY_UP_4
OUTFLOWS:
DEATH_PREY_2 ¼ CONSUME_2
MIG_PREY_2 ¼ MIG_RATE_PREY * (PREY_2ÀDEATH_PREY_2)
PREY_3(t) ¼ PREY_3(t À dt) + (BIRTH_PREY_3 + IM_PREY_3 À
DEATH_PREY_3 À MIG_PREY_3) * dt
INIT PREY_3 ¼ 0
INFLOWS:
BIRTH_PREY_3 ¼ (PREY_3ÀDEATH_PREY_3) * BR_PREY_1 * (1ÀPREY_3/
CC_PREY_3)
IM_PREY_3 ¼ PREY_DN_1+PREY_LF_4
OUTFLOWS:
DEATH_PREY_3 ¼ CONSUME_3
MIG_PREY_3 ¼ MIG_RATE_PREY * (PREY_3ÀDEATH_PREY_3)
PREY_4(t) ¼ PREY_4(t À dt) + (IM_PREY_4 + BIRTH_PREY_4 À
DEATH_PREY_4 À MIG_PREY_4) * dt
INIT PREY_4 ¼ 0
INFLOWS:
IM_PREY_4 ¼ PREY_DN_2+PREY_RT_3
BIRTH_PREY_4 ¼ (PREY_4ÀDEATH_PREY_4) * BR_PREY_1 * (1ÀPREY_4/
CC_PREY_4)
OUTFLOWS:
DEATH_PREY_4 ¼ CONSUME_4
MIG_PREY_4 ¼ MIG_RATE_PREY * (PREY_4ÀDEATH_PREY_4)
BR_PRED ¼ 0.2
BR_PREY_1 ¼ 2
CC_PREY_1 ¼ 2250 * SIN(2 * PI * TIME/CYCLE_TIME)+2500
CC_PREY_2 ¼ 2250 * SIN(2 * PI * TIME/CYCLE_TIMEÀPI/2)+2500
CC_PREY_3 ¼ 2250 * SIN(2 * PI * TIME/CYCLE_TIMEÀ3 * PI/2)+2500
CC_PREY_4 ¼ 2250 * SIN(2 * PI * TIME/CYCLE_TIMEÀPI)+2500
CONSUME_2 ¼ MIN(PREY_2,CONSUME_RATE * PRED_2)
372
40 Modeling Spatial Dynamics of Predator–Prey Interactions. . .
BIRTH_PREY_1 ¼ (PREY_1ÀDEATH_PREY_1) * BR_PREY_1 * (1ÀPREY_1/
CC_PREY_1)
IM_PREY_1 ¼ PREY_LF_2+PREY_UP_3
OUTFLOWS:
DEATH_PREY_1 ¼ CONSUME__1
MIG_PREY_1 ¼ MIG_RATE_PREY * (PREY_1ÀDEATH_PREY_1)
PREY_2(t) ¼ PREY_2(t À dt) + (BIRTH_PREY_2 + IM_PREY_2 À
DEATH_PREY_2 À MIG_PREY_2) * dt
INIT PREY_2 ¼ 0
INFLOWS:
BIRTH_PREY_2 ¼ (PREY_2ÀDEATH_PREY_2) * BR_PREY_1 * (1ÀPREY_2/
CC_PREY_2)
IM_PREY_2 ¼ PREY_RT_1+PREY_UP_4
OUTFLOWS:
DEATH_PREY_2 ¼ CONSUME_2
MIG_PREY_2 ¼ MIG_RATE_PREY * (PREY_2ÀDEATH_PREY_2)
PREY_3(t) ¼ PREY_3(t À dt) + (BIRTH_PREY_3 + IM_PREY_3 À
DEATH_PREY_3 À MIG_PREY_3) * dt
INIT PREY_3 ¼ 0
INFLOWS:
BIRTH_PREY_3 ¼ (PREY_3ÀDEATH_PREY_3) * BR_PREY_1 * (1ÀPREY_3/
CC_PREY_3)
IM_PREY_3 ¼ PREY_DN_1+PREY_LF_4
OUTFLOWS:
DEATH_PREY_3 ¼ CONSUME_3
MIG_PREY_3 ¼ MIG_RATE_PREY * (PREY_3ÀDEATH_PREY_3)
PREY_4(t) ¼ PREY_4(t À dt) + (IM_PREY_4 + BIRTH_PREY_4 À
DEATH_PREY_4 À MIG_PREY_4) * dt
INIT PREY_4 ¼ 0
INFLOWS:
IM_PREY_4 ¼ PREY_DN_2+PREY_RT_3
BIRTH_PREY_4 ¼ (PREY_4ÀDEATH_PREY_4) * BR_PREY_1 * (1ÀPREY_4/
CC_PREY_4)
OUTFLOWS:
DEATH_PREY_4 ¼ CONSUME_4
MIG_PREY_4 ¼ MIG_RATE_PREY * (PREY_4ÀDEATH_PREY_4)
BR_PRED ¼ 0.2
BR_PREY_1 ¼ 2
CC_PREY_1 ¼ 2250 * SIN(2 * PI * TIME/CYCLE_TIME)+2500
CC_PREY_2 ¼ 2250 * SIN(2 * PI * TIME/CYCLE_TIMEÀPI/2)+2500
CC_PREY_3 ¼ 2250 * SIN(2 * PI * TIME/CYCLE_TIMEÀ3 * PI/2)+2500
CC_PREY_4 ¼ 2250 * SIN(2 * PI * TIME/CYCLE_TIMEÀPI)+2500
CONSUME_2 ¼ MIN(PREY_2,CONSUME_RATE * PRED_2)
372
40 Modeling Spatial Dynamics of Predator–Prey Interactions. . .
