May: Harvesting model with a bifurcation from May (1977)...

In cboettig/pdg_control: Pretty Darn Good Control

Description

Harvesting model with a bifurcation from May (1977) increasing parameter "a" (p[3]) will drive the bifurcation

Usage

May(x, h, p)

Arguments

x	the current population level
h	total harvest level
p	vector of parameters c(r, K, a, H, Q)

Details

try pars <- c(r = .75, k = 10, a=1, H=1, Q = 3) See the position of the bifurcation at around a = 1.9: curve(.75*(1-x/10), 0, 10) curve(1.9*x^2/(x^3+1), 0, 10, add=T, col="red") goes to nonzero alternate stable state.

Value

the population level in the next timestep

cboettig/pdg_control documentation built on May 13, 2019, 2:10 p.m.