Si lo he entendido, no necesitas complicarte demasiado, simplemente debes desplazar el cálculo de la población:
pop = 100000
t = 5 # tiempo
t_m = 0.008
t_n = 0.01
población = []
nac = []
mort = []
for i_ in range(t + 1):
población.append(pop)
nacimientos = int(round(pop * t_n))
muertes = int(round(pop * t_m))
pop = pop + nacimientos - muertes
mort.append(muertes)
nac.append(nacimientos)
print(población)
print(nac)
print(mort)
[100000, 100200, 100400, 100601, 100802, 101004] [1000, 1002, 1004, 1006, 1008, 1010] [800, 802, 803, 805, 806, 808]
El cálculo en este caso es:
población[0] = pop
poblacion[1] = poblacion[0] + nac[0] - mort[0]
poblacion[2] = poblacion[1] + nac[1] - mort[1]
...
Otra opción sería que fuera:
población[0] = pop
poblacion[1] = poblacion[0] + nac[1] - mort[1]
poblacion[2] = poblacion[1] + nac[2] - mort[2]
...
pero en este caso lo apropiado es que mort[0]
y nac[0]
sean ambos 0:
pop = 100000
t = 5 # tiempo
t_m = 0.008
t_n = 0.01
población = [pop]
nac = [0]
mot = [0]
for i_ in range(t):
nacimientos = int(round(pop * t_n))
muertes = int(round(pop * t_m))
pop = pop + nacimientos - muertes
población.append(pop)
mot.append(muertes)
nac.append(nacimientos)
print(población)
print(nac)
print(mot)
[100000, 100200, 100400, 100601, 100802, 101004] [0, 1000, 1002, 1004, 1006, 1008] [0, 800, 802, 803, 805, 806]
Asumo que estás trabajando con una población, tasa de natalidad y tasa de mortalidad, por eso he redondeado al entero más cercano dado que en principio no existen individuos medio nacidos o medio muertos ni poblaciones con una fracción de un individuo :). Además resto mort
en vez de sumarlo a la población.