SIRTauLeap function

SIR model with tau leap method (P 6.5).