ModelSIR function

SIR model