SciPy

scipy.special.psi

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

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