## Create an SIR model with scheduled events.
model <- SIR(
u0 = u0_SIR(),
tspan = 1:(4 * 365),
events = events_SIR(),
beta = 0.16,
gamma = 0.077
)
## Extract the events and display a summary.
ev <- events(model)
summary(ev)
## Plot the event schedule over time.
plot(ev)
Run the code above in your browser using DataLab