Empecemos con el$\;\displaystyle\int \frac {dx}{\cosh(x)^2}=\frac 2{1+e^{-2x}}+C\;$ (con la opción$C=-2$) y usamos la integración por partes:
\begin{align}
I(\alpha) &= \int_0^\infty \frac{x^\alpha}{\cosh^2x}dx\\
&=\left.x^\alpha\left(\frac 2{1+e^{-2x}}-2\right)\right|_0^\infty-\int_0^\infty \alpha\,x^{\alpha-1}\left(\frac 2{1+e^{-2x}}-2\right)\,dx\\
&=2\alpha \int_0^\infty \frac {x^{\alpha-1}\,e^{-2x}}{1+e^{-2x}}\,dx\\
&=2\alpha \int_0^\infty \frac {x^{\alpha-1}}{e^{2x}+1}\,dx\\
&=2\alpha\, 2^{1-\alpha-1}\int_0^\infty \frac {(2x)^{\alpha-1}}{e^{2x}+1}\,d(2x)\\
&= 2^{1-\alpha}\alpha\,\int_0^\infty \frac {t^{\alpha-1}}{e^t+1}\,dt
\end{align}
Ahora una variante cerrada de$\zeta$ es la función de Dirichlet eta que puede ser definida por:$$\eta(s) := \frac{1}{\Gamma(s)} \int_{0}^{\infty} \frac{t^{s-1}}{e ^ t + 1}\, dt=\left(1-2^{1-s}\right)\zeta(s) $ $
Todo esto nos da:$$I(\alpha) ={2^{1-\alpha}\left(1-2^{1-\alpha}\right)}{\,\Gamma(\alpha+1)}\,\zeta(\alpha)$ $