SIR model with Gillespie algorithm (Tau-leap method)