mlgamma function

Gamma distribution maximum likelihood estimation