simulate_SIR_model_discrete function

SIR model