##Model Input
S_0 <- 989
E_0 <- 10
I_0 <- 1
R_0 <- 0
beta <- 3
gamma <- 1/2
tau <- 2
chi <- 0.5
parameters <- c(beta = beta, gamma = gamma,
chi = chi, tau = tau)
inits <- c(S = S_0, E = E_0, I = I_0, R_0 = R_0)
SEIRS_ode(1, inits, parameters)
Run the code above in your browser using DataLab