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))