Using the user-specified numbers of people N a in the age classes and the fractions F a of people under high risk within each age class (Table 2) , we obtain the initial population sizes according to age and risk class: Using these initial values, the set of differential equations is solved numerically with a Runge-Kutta method with step-size control. 