Моделирование эпидемий (Python)

Простейшая модель развития эпидемии описывается формулой:

$B\left(n+1\right)=\left(1+k\frac{L-B\left(n\right)-V\left(n\right)}{L}\right)B\left(n\right),$

где B (n) – количество больных в n-й день эпидемии, $V\left(n\right)$ – количество людей, имеющих иммунитет (переболевших или вакцинированных), L – общая численность населения, k – коэффициент распространения заболевания.

Графики прогноза начала эпидемии Covid-19 в трёх городах.