# inst/doc/examples/upca.R In simecol: Simulation of Ecological (and Other) Dynamic Systems

```#################################################
# UPCA: Uniform Period Chaotic Amplitude Model
# see:
# Blasius, B., Huppert, A., and Stone, L. (1999)
# Complex dynamics and phase synchronization in
# spatially extended ecological systems. Nature,
# vol 399, p. 354--359.
# R-Implementaion by Th. Petzoldt
#################################################

upca <- new("odeModel",
main = function(time, init, parms) {
u      <- init[1]
v      <- init[2]
w      <- init[3]
with(as.list(parms), {
du <-  a * u           - alpha1 * f(u, v, k1)
dv <- -b * v           + alpha1 * f(u, v, k1) +
- alpha2 * f(v, w, k2)
dw <- -c * (w - wstar) + alpha2 * f(v, w, k2)
list(c(du, dv, dw))
})
},
equations  = list(
f1 = function(x, y, k){x*y},           # Lotka-Volterra
f2 = function(x, y, k){x*y / (1+k*x)}  # Holling II
),
times  = c(from=0, to=100, by=0.1),
parms  = c(a=1, b=1, c=10, alpha1=0.2, alpha2=1,
k1=0.05, k2=0, wstar=0.006),
init = c(u=10, v=5, w=0.1),
solver = "lsoda"
)

upca@equations\$f <- upca@equations\$f2
```

## Try the simecol package in your browser

Any scripts or data that you put into this service are public.

simecol documentation built on Nov. 17, 2017, 6:45 a.m.