add a function for calculate the gamma used by Hill (1987)