$$\int_{0}^{\infty}{x^s-1 \over x-1}\cdot{1-e^{-ax} \over 1-e^{ax}}dx$$
$$\int_{0}^{\infty}{x^s-1 \over x-1}\cdot{1-e^{-ax} \over 1-e^{ax}}\cdot{e^{ax}\over e^{ax}}dx$$
$$\int_{0}^{\infty}-{x^s-1 \over x-1}\cdot{e^{-ax}}dx$$
Suppose $s=2$
$$\int_{0}^{\infty}-{x^2-1 \over x-1}\cdot{e^{-ax}}dx$$
$(x^2-1)=(x+1)(x-1)$
$$-\int_{0}^{\infty}{e^{-ax}}(1+x)dx$$
Integrating using by parts
$u=1+x$ $\rightarrow$ $du=dv$
$dv=e^{-ax}$ $\rightarrow$ $v=-{{e^{-ax}} \over {a}}$
$=\left.{-(1+x)e^{-ax}\over a}\right|_{0}^{\infty}-{1\over a}\int_{0}^{\infty}e^{-ax}dx$
$=\left.{-(1+x)e^{-ax}\over a}\right|_{0}^{\infty}+{1\over a^2}\left.e^{-ax}\right|_{0}^{\infty}$
$=-{1\over a}-{1\over a^2}$
Let $s=k$ for k is an integer $k\ge1$ ; restriction to integers for sake of easiness
The factorisation of $x^k-1$ is quite lengthy, so integrating via this method seem tedious. How can one evaluate this integral using another less effort strategy?