Gammafit-func
gamfit(vX, a, b) uses Newton-Raphson method to calculate Maximum Likelihood Estimates (MLE) of a (shape) and b (scale) for Gamma distributed dataset vX.
The probability density function of Gamma distribution is:
where a, b > 0, with and .