total_people = 857;
days= 20;
disp(total_people/days, "The maximum likelihood estimate of lambda")