SIRBirthDeath function

SIR model with births and deaths (P 2.2).