scipy.special.gammaln

scipy.special.gammaln(x[, out]) = <ufunc 'gammaln'>

y=gammaln(z) returns the base e logarithm of the absolute value of the gamma function of z: ln(|gamma(z)|)