SIR function

Simple SIR model (P 2.1).