##Model Input
S_u_0 <- 989
I_u_0 <- 1
R_u_0 <- 0
S_v_0 <- 0
I_v_0 <- 0
R_v_0 <- 0
beta <- 3
tau <- 2
mu <- 1/81
alpha <- 0.8
lambda <- 0.7
parameters <- c(beta = beta, tau = tau, mu = mu,
alpha = 0.8, lambda = 0.7)
inits <- c(S_u = S_u_0, I_u = I_u_0, R_u_0 = R_u_0,
S_v = S_v_0, I_v = I_v_0, R_v_0 = R_v_0)
SIR_vaccination_demographics_ode(1, inits, parameters)
Run the code above in your browser using DataLab