A friend of mine and I wanted to solve the following indefinite integral but got stuck:
$$ \int \frac{x}{e^x - 1} dx. $$
My approach:
Let
$$ I = \int \frac{x}{e^x - 1} dx.\\ \implies I = x\int \frac{dx}{e^x - 1} - \int \left ( \int \frac{dx}{e^x - 1} \right ) dx. $$
Now, let $I_2 = \int \frac{dx}{e^x - 1}$. Also, let $z = e^x \implies dx = {dz \over z}$. Then,
$$ I_2 = \int \frac{dz}{z(z-1)} \\ \implies I_2 = \int {dz \over {z - 1}} - \int {dz \over z} \\ \implies I_2 = \ln (e^x - 1) - x. $$
Substituting the value of $I_2$ in $I$, we get,
$$ I = x[\ln (e^x - 1) - x] + {x^2 \over 2} - \int \ln (e^x - 1) dx. $$
I got stuck right here. Is it possible to proceed further?