1546.sagews

#Madison Schaeffer

from random import gauss
def exponential(N,replacement):
    return N * replacement

# declare data
lam = 1.2
time = 100
N0 = 100

pop_sizes = [N0] # y data
time_list = range(time) # x data
meta_list = []

# main (orders things before)
for years in range(1,time): #[1,2...time]
    pop_sizes.append(exponential(pop_sizes[-1],lam))
    
for values in range(len(pop_sizes)):
    meta_list.append([values, pop_sizes[values]])
    
#print meta_list
line(meta_list, marker = '.') + text('Time', (50, -600), color='black')