sirmod function

Gradient-function for the SIR model