SIR model with Gillespie (Tau-leap) algorithm