[From Bill Powers (950813.1300 MDT)]
Slight correction in program COLLIER.PAS. Where you see
for i := 0 to N2 - N1 + 1
eliminate the "+1".
Two places.
Also note that there is an analytical solution:
o = g*r/(1 + g*mealrate) and
p = g*r*mealrate/(1 + g*mealrate)
Best,
Bill P.