seirmod function

Gradient-function for the SEIR model