Model of mutualism with saturating functional response

Used internally by run_mutualism().

Usage

mutualism(time, init, params)

Arguments

time

Vector of time units over which to run model, starting at 0.

init

Named vector of initial population sizes.

• N1: Initial population size of species 1

• N2: Initial population size of species 2

params

Named vector of model parameters.

• r1: Per capita growth rate of species 1

• r2: Per capita growth rate of species 1

• a12: Effect of species 2 on species 1

• a21: Effect of species 1 on species 2

• b1: Half-saturation constant for species 1

• b2: Half-saturation constant for species 2

• d1: Rate of self-limitation for species 1

• d2: Rate of self-limitation for species 2

Value

A list where the first and only element is a vector of derivatives of the state variables (N1, N2) with respect to time, given params.

See also

run_mutualism()