Estoy tratando de demostrar que
$$I:= \int_0^\infty \frac{x}{\sinh(ax)} dx = \frac{\pi^2}{4a^2}$$
Intento:
$$\sinh (ax) = \frac{1}{2}(e^{ax}-e^{-ax}) = \frac{1}{2}e^{-ax}(e^{2ax}-1)$$
Ahora tengo $$\int_0^\infty \frac{x}{\sinh(ax)} dx = \int_0^\infty 2x \sum_{q=0}^\infty \frac{\frac{a^qx^q}{q!}}{e^{2ax}-1} dx$$ .
Sustituyendo $y=2ax$ Lo entiendo: $$I= \frac{1}{2a^2} \sum_{q=0}^\infty \frac{1}{2^qq!} \int_0^\infty \frac{y^{q+1}}{e^y-1}dy=\frac{1}{2a^2} \sum_{q=0}^\infty \frac{1}{2^qq!} \Gamma(q+2) \zeta(q+2)$$
¿Estoy en el camino correcto? ¿Qué pasa con la serie infinita; por qué debe tener el valor $\frac{\pi^2}{2}$ ?
Cualquier sugerencia será muy apreciada.