Group 6

Nadine Galler

Kirk M. Soodhalter

Jennifer S. Novak

Project Description:

Develop and simulate a Hastings model for two competing species on a 1-dimensional flat landscape. Then develop the model for the same situation on a curved landscape for 1 and 2 species. Try to incorporate altitude in the probability density, instead of only distance.

The Hastings Model

The Hasting's model assumes the following:

- Discrete time (reproductive cycle, fixed life span, etc)

- Adults are stationary (plants, crabs, etc)

- Juveniles are motile

The Hasting's model is as follows:

N_Next := proc (x, N) options operator, arrow; a*N+...

where a is the survival rate of the adults per cycle,

the probability density of the juveniles is given by k := proc (x) options operator, arrow; c*exp(-beta*...

and m is the number of juveniles reaching adulthood per adult.

Basic Assumptions:

a: = 0 ( adults do not survive from one cycle to the next)

beta := 1 ( dispersal rate)

c := 1/2 ( conservative constant)

m:=1 ( one offspring per adult)

1. Motivation:

1 species on a flat landscape:

N_NEXT := proc (x, N) options operator, arrow; c*in...

2. Two competing species on a flat landscape:

Our original model:

N_1_Next := proc (x, N_1, N_2) options operator, ar...
N_1_Next := proc (x, N_1, N_2) options operator, ar...

N_2_Next := proc (x, N_1, N_2) options operator, ar...
N_2_Next := proc (x, N_1, N_2) options operator, ar...

N_1 := proc (x) options operator, arrow; piecewise(...

N_2 := proc (x) options operator, arrow; piecewise(... [Maple Plot]

These are the 2 species after the 2nd cycle:

[Maple Plot]

These are the 2 species after the 3rd cycle:

[Maple Plot]

These are the 2 species after the 4th cycle:

[Maple Plot]

Here we changed the competing factor: instead of ( 0.1 , 0.9) we used ( 0.3, 0.7):

[Maple Plot] [Maple Plot] [Maple Plot]

Possible modifications to the model:

N_1_Next := proc (x, N_1, N_2) options operator, ar...
N_1_Next := proc (x, N_1, N_2) options operator, ar...

3. One species on a curved landscape:

[Maple Plot]

N_NEXT := proc (x, g, N) options operator, arrow; a...
N_NEXT := proc (x, g, N) options operator, arrow; a...

[Maple Plot] [Maple Plot]

[Maple Plot] [Maple Plot]

[Maple Plot] [Maple Plot]

These models are not conservative. We conjecture that the conservative model would have:

c := 1/2*beta/((1+diff(g(t),t)(x))^(1/2))

4. Two species on a curved landscape:

N_1_Next := proc (x, g, N_1, N_2) options operator,...

N_2_Next := proc (x, g, N_1, N_2) options operator,...

[Maple Plot]

5. One species on a curved landscape with altitude:

N_NEXT := proc (x, g, N) options operator, arrow; a...
N_NEXT := proc (x, g, N) options operator, arrow; a...

[Maple Plot] [Maple Plot]

N_NEXT := proc (x, g, N) options operator, arrow; a...
N_NEXT := proc (x, g, N) options operator, arrow; a...

[Maple Plot]