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