mlinvgamma function

Inverse Gamma distribution maximum likelihood estimation