EstMLEBetaBin function

Estimating the shape parameters a and b for Beta-Binomial Distribution