ssir function

Stochastic SIR simulation