ModelSIRD function

SIRD model