gamllik function

Log likelihood from the gamma distribution