population = [7333253 3448613 2731743 1702086 1524249 1151977 1048949 1022830 998905 992038 816884 752279 734676 702979 665070 635913 617044 614289 579307 567094 547727 520947 514013 504505 493559];
disp(perctl(population, 10), "The sample 10 percentile is")
disp(perctl(population, 80), "The sample 80 percentile is")
disp(perctl(population, 50), "The sample 80 percentile is")
disp(median(population), "The median is")