##############
# Repeated Dosages Model for Math 311
# Mark Goadrich Spring 2009
##############

# Set up file for writing output
fout = open("repeated.csv", "w")

# Initialize parameters
decayRate = 0.08
simLength = 100
population = 0
intervalSize = 0.25  # decrease to approximate continuous growth
dosage = 50
decayRatePerStep = decayRate * intervalSize

# Iterate through simulation
for i in range(int(simLength / intervalSize)):
    population = population - decayRatePerStep * population
    if (i % 30 == 0):
        population += dosage
    fout.write("%f, %f\n" % (i, population))

# Close the file
fout.close()