mlnbinom function

Negative binomial distribution maximum likelihood estimation