Starting with:
$$\zeta (s)={\frac {1}{\Gamma (s)}}\int _{0}^{\infty }{\frac {x^{s-1}}{e^{x}-1}}\,\mathrm {d} x$$
How can we prove for $s > 2$ the following conjecture:
$$\zeta (s-1)-\zeta(s)={\frac {2}{\Gamma (s+1)}}\int _{0}^{\infty }{\frac {x^s e^x}{(e^{x}-1)^3}}\,\mathrm {d} x$$
Integration by parts may be a way to go with:
$$u=x^s \qquad du=s x^{s-1} \\ dv=\frac {e^x dx}{(e^{x}-1)^3} \qquad v=- \frac{1}{2} \frac {1}{(e^{x}-1)^2}$$
Which gives us:
$$\zeta (s-1)-\zeta(s)={\frac {1}{\Gamma (s)}}\int _{0}^{\infty }{\frac {x^{s-1} }{(e^{x}-1)^2}}\,\mathrm {d} x$$
It's a little nicer, but I'm not sure how to prove it either.
It just came to me: we should probably do another integration by parts, that's it.