Assuming that the incubation period is a random variable with exponential distribution with parameter a (i.e. the average incubation period is a^{-1}), and also assuming the presence of vital dynamics with birth rate equal to death rate, we have the model: