from modsim import *
def run_simulation(system, growth_func):
results = TimeSeries()
results[system.t_0] = system.p_0
for t in range(system.t_0, system.t_end):
growth = growth_func(t, results[t], system)
results[t+1] = results[t] + growth
return results
