Comienzas (t = 0) con tu número de infectados (por ejemplo, 1e-6 por uno en un millón) y no (es decir, cero) recuperado. Eso significa que sus susceptibles son 0.999999 de la población.
Su simulación se ejecutará en el período de tiempo apropiado, por ejemplo, un día a la vez. Entonces en t = 1 su nueva cantidad de infectados aumentará en S * I * “infecciosidad”, y también disminuirá en I / “tasa de recuperación en días”, que son, por supuesto, las cantidades S disminuirán en y R aumentarán en.
Probablemente debería comparar y contrastar diferentes dT, por ejemplo, una hora o un mes; si el tiempo de recuperación es del orden de días de una semana puede ser suficiente, es casi seguro que los meses no, y las horas pueden ser excesivas.
Pero realmente es tan simple como simplemente hacer algunas multiplicaciones en cada paso del tiempo y mantener S + I + R = 1 (o lo que sea apropiado para su población). Sin embargo, no se preocupe por individuos parciales; use punto flotante independientemente
