This function simulates cases from an stochastic SIR model. The user may specify the initial epidemic growth rate \(r\), the rate of recovery gamma \(\gamma\), the initial number of infected cases \(I_0\), and the total population size \(N\).

## See also

Other simulate:
`simulate_exponential_cases()`

,
`simulate_secondary()`

,
`simulate_uniform_cases()`