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 _ 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 _ 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.