emv_gamma function

compute the MLE in the Gamma model using the Newton-Raphson method