scipy.special.gammaln

scipy.special.gammaln(z) returns the base e logarithm of the absolute value of the gamma function of z: ln(abs(gamma(z))) = <ufunc 'gammaln'>

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

See also

gammasgn