coyne: Gradient-function for Coyne et al's rabies model
In epimdr: Functions and Data for "Epidemics: Models and Data in R"

Description Usage Arguments Value Examples

Gradient-function for Coyne et al's rabies model

1	coyne(t, logx, parms)

`t`	Implicit argument for time
`logx`	A vector with values for the log-states
`parms`	A vector with parameter values for the dynamical system

A list of gradients for the log system

require(deSolve)
times  = seq(0, 50, by=1/520)
paras  = c(gamma = 0.0397, b = 0.836, a = 1.34, sigma = 7.5, 
alpha = 66.36, beta = 33.25, c = 0, rho = 0.8)
start = log(c(X=12.69/2, H1=0.1, H2=0.1, Y = 0.1, I = 0.1))
out = as.data.frame(ode(start, times, coyne, paras))