Definir
I(α)=∫∞0logxlog(1−e−αx)dx.I(α)=∫∞0logxlog(1−e−αx)dx.
La integración por partes, seguido por la sustitución αx↦xαx↦x, tenemos
\begin{align*}
I(\alpha)
&= \alpha \int_{0}^{\infty} \frac{x - x\log x}{e^{\alpha x} - 1} \, dx \\
&= \frac{1}{\alpha} \int_{0}^{\infty} \frac{(1+\log \alpha) x - x \log x}{e^{x} de {- 1} \, dx\\
&= \frac{1}{\alpha} \left\{ (1+\log\alpha)\zeta(2) - \left.\frac{d \zeta(s)\Gamma(s)}{s}\right|_{s=2} \right\}\\
&= \frac{1}{\alpha} \left\{ (\gamma+\log\alpha)\zeta(2) - \zeta'(2) \right\}.
\end{align*}
Luego de ello se sigue que
∫∞0logxlog(1+12coshx)dx=∫∞0logxlog(1−e−3x1−e−x⋅1−e−21−e−4x)dx=I(2)+I(3)−I(1)I(4)=512ζ′(2)−572γπ2+118π2log(3).
Conectar algunas identidades relativas ζ′(2) y la Glaisher-Kinkelin constante de A, esto se reduce a Vladimir de la respuesta.
Anexo - Algo que quieras saber:
La siguiente identidad juega un papel clave en esta prueba.
$$ \int_{0}^{\infty} \frac{x^{m-1}}{e^{x} de {- 1} \, dx = \Gamma(s)\zeta(s), $$
que vale para ℜs>1.