2020/11/13にるめなる君が こちら で記事にした一般化です。試しに解いてみたところ方針が少し違ったので記事にすることにしました。
$$ \displaystyle \int_0^\infty \frac{x^{n-1}e^x}{\sinh x\cosh x}dx $$
[解説]
$ \begin{eqnarray*} &&\d\int_0^\infty \frac{x^{n-1}e^x}{\sinh x\cosh x}dx\\ &=&2\int_0^\infty \frac{x^{n-1}e^x}{\sinh2x}dx\\ &=&4\int_0^\infty x^{n-1}e^{-x}\sum_{k=0}^\infty e^{-4kx}dx\\ &=&4\sum_{n=0}^\infty \int_0^\infty x^{n-1}e^{-(4k+1)x}dx\\ &=&4\sum_{k=0}^\infty \frac1{(4k+1)^n}\Gamma(n)\\ &=&2\Gamma(n)((1-2^{-n})\zeta(n)+\beta(n))\\ \end{eqnarray*} $
よって、
$\d\int_0^\infty \frac{x^{n-1}e^x}{\sinh x\cosh x}dx=2\Gamma(n)((1-2^{-n})\zeta(n)+\beta(n))$
が証明されました。