Observe que
\begin{align*} \int_{0}^{\infty} \frac{\sin x}{x^s} \, dx &= \frac{1}{\Gamma(s)} \int_{0}^{\infty} \left( \int_{0}^{\infty} t^{s-1}e^{-xt} \, dt \right) \sin x \, dx \\ &= \frac{1}{\Gamma(s)} \int_{0}^{\infty} \left( \int_{0}^{\infty} e^{-tx} \sin x \, dx \right) t^{s-1} \, dt \\ &= \frac{1}{\Gamma(s)} \int_{0}^{\infty} \frac{t^{s-1}}{1+t^2} \, dt \\ &= \frac{1}{2\Gamma(s)} \beta\left(\frac{s}{2}, 1-\frac{s}{2}\right) \\ &= \frac{\pi}{2\Gamma(s)\sin\left(\frac{\pi s}{2}\right)}. \end{align*}
(Este cálculo heurístico es válido línea por línea en la franja $1 < \Re(s) <2$ en vista del teorema de Fubini, y luego se extiende a la franja mayor $0 < \Re(s) < 2$ por continuación analítica). Diferenciando ambos lados, obtenemos
$$ \int_{0}^{\infty} \frac{\sin x \log x}{x^s} \, dx \stackrel{(*)}{=} -\frac{d}{ds} \frac{\pi}{2\Gamma(s)\sin\left(\frac{\pi s}{2}\right)} = \frac{\pi}{2\Gamma(s)\sin\left(\frac{\pi s}{2}\right)} \left( \psi(s) + \frac{\pi}{2}\cot\left(\frac{\pi s}{2}\right) \right). $$
Ahora la respuesta se obtiene enchufando $s = 1$ :
$$ \int_{0}^{\infty} \frac{\sin x \log x}{x} \, dx = -\frac{\gamma \pi}{2}. $$
Justificación de $\text{(*)}$ . Demostremos que
$$F(s) := \int_{0}^{\infty} \frac{\sin x}{x^s} \, dx \tag{1}$$
es analítica en la franja abierta $S = \{ s \in \Bbb{C} : 0 < \Re(s) < 2 \}$ y su derivada se puede calcular por la regla integral de Leibniz. Una cuestión importante de $\text{(1)}$ es que la integral definitoria converge sólo condicionalmente para $0 < \Re(s) \leq 1$ y, por tanto, plantea algunas dificultades técnicas. Para evitarlo, mejoramos la velocidad de convergencia utilizando la integración por partes:
$$ F(s) = \underbrace{\left[ \frac{1-\cos x}{x^s} \right]_{0}^{\infty}}_{=0} + s \int_{0}^{\infty} \frac{1-\cos x}{x^{s+1}} \, dx. $$
Obsérvese que la integral resultante es absolutamente convergente en $S$ .
Ahora afirmamos que $g(s) := F(s)/s$ es diferenciable y su derivada se puede calcular por la regla integral de Leibniz. Sea $\epsilon$ sea tal que $\bar{B}(s, \epsilon) \subset S$ . Entonces, siempre que $0 < |h| < \epsilon$ ,
\begin{align*} &\frac{g(s+h) - g(s)}{h} - \int_{0}^{\infty} \frac{(1-\cos x)(-\log x)}{x^{1+s}} \, dx \\ &\hspace{9em} = \int_{0}^{\infty} \frac{1-\cos x}{x^{1+s}} \left( \frac{x^{-h} - 1}{h} + \log x \right) \, dx. \end{align*}
Obsérvese que el integrando está dominado por
$$ \left| \frac{1-\cos x}{x^{1+s}} \left( \frac{x^{-h} - 1}{h} + \log x \right) \right| \leq \frac{1-\cos x}{x^{1+\Re(s)}} (\max\{ x^{\epsilon}, x^{-\epsilon} \} + 1) |\log x| $$
Esta función dominante es integrable. Por lo tanto, por el teorema de convergencia dominante, como $h \to 0$ tenemos
$$ \lim_{h \to 0} \frac{g(s+h) - g(s)}{h} = \int_{0}^{\infty} \frac{(1-\cos x)(-\log x)}{x^{1+s}} \, dx. $$
Si lo corregimos, sabemos que $F(s)$ es diferenciable con
$$ F'(s) = \int_{0}^{\infty} \frac{1-\cos x}{x^{s+1}} \, dx + s \int_{0}^{\infty} \frac{(1-\cos x)(-\log x)}{x^{1+s}} \, dx. $$
Por último, realizando la integración por partes a esta última integral se obtiene la conclusión deseada:
\begin{align*} &s \int_{0}^{\infty} \frac{(1-\cos x)(-\log x)}{x^{1+s}} \, dx \\ &\hspace{5em} = \underbrace{\left[ \frac{(1-\cos x)\log x}{x^s} \right]_{0}^{\infty}}_{=0} - \int_{0}^{\infty} \left( \frac{\sin x \log x}{x^s} + \frac{1-\cos x}{x^{1+s}} \right) \, dx \end{align*}