data(upca)
plot(sim(upca))
# omit stabilizing parameter wstar
parms(upca) <- list(wstar=0)
plot(sim(upca))
# change functional response from
# Holling II (default) to Lotka-Volterra
equations(upca)$f <- function(x, y, k) x * y
plot(sim(upca))Run the code above in your browser using DataLab