## Extract code from an mparse-generated SIR model.
## Using writeLines formats the character vector output for
## readability.
model <- mparse(
transitions = c("S -> beta * S * I/(S + I + R) -> I",
"I -> gamma * I -> R"),
compartments = c("S", "I", "R"),
gdata = c(beta = 0.16, gamma = 0.077),
u0 = data.frame(S = 99, I = 1, R = 0),
tspan = 1:10,
use_enum = TRUE
)
writeLines(C_code(model))
Run the code above in your browser using DataLab