SEIR function

SEIR model (2.6).