Fig. 1

Example of an SIR model using the ODE solver (solid lines) and the stochastic solver (dashed lines) within Broadwick. Both models start with 99 susceptible (S), 1 infected (I) and 0 recovered (R) individuals with β=0.004 and ρ=0.04 with both solvers adding 50 new individuals at t=20 using a θ-event