1550.sagews

def expo(N, r, t):
    return N * exp(r*t)
    
time = 10
n0 = 100.
pop_sizes = [n0]

for years in range(1, time):
    pop_sizes.append(expo(pop_sizes[-1], 0.002, years))
    
series = []
for p in range(len(pop_sizes)):
    series.append([p, pop_sizes[p]])

line(series, marker='.')