Puedo ver que es correcto utilizando algunos conocimientos de la función Gamma. Tenemos
$$ \Gamma(\alpha ) = \int_0^\infty e^{-x} x^{\alpha - 1 } dx . $$
Diferenciando con respecto a $\alpha$ obtenemos
$$ \frac{d\Gamma}{d \alpha} = \int_0^\infty e^{-x} x^{\alpha - 1 }\ln x dx . $$
Configuración $\alpha = 1 $ obtenemos
$$ \int_0^\infty e^{-x} \ln x dx = \frac{d\Gamma}{d\alpha }\bigg|_{\alpha = 1 }= - \gamma. $$
Pero el conocimiento $d \Gamma/d \alpha |_{\alpha = 1}=-\gamma$ es un misterio para mí.
¿Alguien puede encontrar una prueba elemental?