samplesize = 1000;
acceptable =921;
disp(acceptable/samplesize, "The maximum likelihood estimate of p is")