sirwmod | R Documentation |
Gradient-function for the SIRWS model
sirwmod(t, logy, parameters)
t |
Implicit argument for time |
logy |
A vector with values for the log(states) |
parameters |
A vector with parameter values for the SIRWS system |
A list of gradients (in log-coordinates)
require(deSolve) times = seq(0, 26, by=1/10) paras = c(mu = 1/70, p=0.2, N = 1, beta = 200, omega = 1/10, gamma = 17, kappa=30) start = log(c(S=0.06, I=0.01, R=0.92, W = 0.01)) out = as.data.frame(ode(start, times, sirwmod, paras))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.