gammaln - The logarithm of gamma function.
gammaln(x) evaluates the logarithm of gamma function at all the elements of x , avoiding underflow and overflow. x must be real.
gammaln(0.5)
gamma , dlgamma ,
W. J. Cody and L. Stoltz (code from Netlib (specfun))