scipy.special.psi

scipy.special.psi(x[, out])

y=psi(z) is the derivative of the logarithm of the gamma function evaluated at z (also called the digamma function).