Я играю на сверхмедленной скорости (1 ход - 1 год), поэтому и формулу пишу под такую скорость.
Основная формула имеет вид P (n+1) = B (n+1) + P (n) - D (n+1), где P - население города, n - номер хода, B - количество родившихся жителей, D - количество умерших жителей. При этом P по своему составу неоднородно: в нём имеются дети C (0-18 лет; едят, не работают, потомство не производят), взрослые R (18-45 лет; едят, работают, производят потомство), пожилые люди O (едят, в воспроизводстве населения не участвуют, работают или нет - активно над этим думаю).
B зависит от количества взрослых людей и коэффициента рождаемости m (я использую его в перерасчёте только на взрослых, т.к. все остальные все равно не репродуктивны). D, соответственно, от коэффициента смертности d.
Сейчас думаю по поводу того, от каких величин зависят m и d. Пока известно, что m находится в отрезке [0; 2/3].