start = 1800 end = 2000 lifespan = end—start ci = lifespan * .25 end_i = end—ci end_f = end + ci
p_max = 0 max_correct = 0 for i in range(10000): p = i * .001 correct = 0 for year in range(lifespan): end_prediction = year * p + start if end_prediction >= end_i and end_prediction ⇐ end_f: correct += 1
start = 1800
end = 2000
lifespan = end—start
ci = lifespan * .25
end_i = end—ci
end_f = end + ci
p_max = 0
max_correct = 0
for i in range(10000):
p = i * .001
correct = 0
for year in range(lifespan):
end_prediction = year * p + start
if end_prediction >= end_i and end_prediction ⇐ end_f:
correct += 1
if correct > max_correct:
max_correct = correct
p_max = p
print(“correct %:”, max_correct / lifespan)
print(“p: ”, p_max)
works_i = (end_i—start) / p_max
works_f = (end_f—start) / p_max
print(“Works from”, start + works_i, “to”, start + works_f)