SIR function

Create an SIR model