strogatz: Solve Strogatz's model

In PabRod/sleepR: Sleep-wake Cycle Models

Description

Solves the Strogatz's model for the given times, initial condition and parameters

Usage

strogatz(ts, y0, parms = strogatz_default_parms(), ...)

Arguments

ts	Vector of times (in h)
y0	Initial condition
parms	Model parameters (optional, see strogatz_default_parms)
...	Additional arguments passed to the ode integrator

Value

Results of the simulation, including times, states and asleep/awake status

Author(s)

Pablo Rodríguez-Sánchez (https://pabrod.github.io)

References

Strogatz, S. H. (1987). Human sleep and circadian rhythms: a simple model based on two coupled oscillators. Journal of Mathematical Biology, 25(3), 327–347. http://doi.org/10.1007/BF00276440

See Also

dStrogatz, strogatz_default_parms, ode

Examples

y0 <- c(th1 = 0.1, th2 = 0.05)
nDays <- 30
ts <- seq(0, nDays*24, length.out=nDays*24*10)
ys <- strogatz(ts, y0)

PabRod/sleepR documentation built on Dec. 21, 2020, 3:27 p.m.