A curious integral

Without resorting to the Abel-Plana formula, you may simply notice that

$$ I=\int_{0}^{+\infty}\frac{\sin x}{e^x-1}\,dx = \sum_{n\geq 1}\int_{0}^{+\infty}\sin(x)e^{-nx}\,dx = \sum_{n\geq 1}\frac{1}{n^2+1} \tag{1}$$ and since $\frac{\sin(\pi x)}{\pi x}=\prod_{n\geq 1}\left(1-\frac{x^2}{n^2}\right)$, by applying $\frac{d}{dx}\log(\cdot)$ to both sides we get: $$ -\frac{1}{x}+\pi\cot(\pi x) = \sum_{n\geq 1}\frac{2x}{x^2-n^2}\tag{2} $$ as well as: $$ \sum_{n\geq 1}\frac{1}{n^2+z^2} = \frac{-1+\pi z \coth(\pi z)}{2 z^2}\tag{3}$$ from which $I=\frac{-1+\pi\coth(\pi)}{2}$ clearly follows.
Have a look at this thread for further proofs.

$\newcommand{\bbx}[1]{\,\bbox[15px,border:1px groove navy]{\displaystyle{#1}}\,} \newcommand{\braces}[1]{\left\lbrace\,{#1}\,\right\rbrace} \newcommand{\bracks}[1]{\left\lbrack\,{#1}\,\right\rbrack} \newcommand{\dd}{\mathrm{d}} \newcommand{\ds}[1]{\displaystyle{#1}} \newcommand{\expo}[1]{\,\mathrm{e}^{#1}\,} \newcommand{\ic}{\mathrm{i}} \newcommand{\mc}[1]{\mathcal{#1}} \newcommand{\mrm}[1]{\mathrm{#1}} \newcommand{\pars}[1]{\left(\,{#1}\,\right)} \newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}} \newcommand{\root}[2][]{\,\sqrt[#1]{\,{#2}\,}\,} \newcommand{\totald}[3][]{\frac{\mathrm{d}^{#1} #2}{\mathrm{d} #3^{#1}}} \newcommand{\verts}[1]{\left\vert\,{#1}\,\right\vert}$ \begin{align} \int_{0}^{\infty}{\sin\pars{x} \over \expo{x} - 1}\,\dd x & = {1 \over 2\ic}\int_{0}^{\infty}{\expo{-\ic x} - \expo{\ic x} \over 1 - \expo{-x}}\,\pars{-\expo{-x}}\,\dd x \,\,\,\stackrel{t\ =\ \expo{\large -x}}{=}\,\,\, -\,{1 \over 2\ic}\int_{0}^{1}{t^{\large\ic} - t^{\large -\ic} \over 1 - t} \,\dd t \\[5mm] & = -\,{1 \over 2\ic}\pars{H_{-\ic} - H_{\ic}}\qquad\qquad\qquad \pars{~H_{z}:\ Harmonic Number~} \\[5mm] & = -\,{1 \over 2\ic}\bracks{H_{-\ic} - \pars{H_{\ic - 1} + {1 \over \ic}}} \qquad\qquad\qquad\pars{~H_{z}\ Recurrence~} \\[5mm] & = -\,{1 \over 2} - {1 \over 2\ic}\bracks{H_{-\ic} - H_{\ic - 1}} = -\,{1 \over 2}-\,{1 \over 2\ic}\bracks{\pi\cot\pars{\pi\ic}} \pars{\substack{Euler \\[1mm] Reflection\ Formula}} \\[5mm] & = -\,{1 \over 2}-\,{\pi \over 2\ic}\bracks{-\ic\coth\pars{\pi}} = \bbx{{\pi \over 2}\,\coth\pars{\pi} - {1 \over 2}} \approx 1.077 \end{align}

Contour integration solution: We consider the integral

$$I = \oint_{C\left(\epsilon,R\right)} \frac{\exp(i z)dz}{\exp(z) - 1}$$

with $C\left(\epsilon,R\right)$ a rectangle from $\epsilon$ to $R$ on the real axis, from there to $R + 2\pi i$ parallel to the imaginary axis, from there to $\epsilon + 2\pi i$ parallel to the real axis, then a clockwise quarter turn with radius $\epsilon$ and center $2\pi i$ to the point $2\pi i - i\epsilon$, from there we move on the imaginary axis to the point $i\epsilon$ and then we take a clockwise quarter turn with radius $\epsilon$ and center the origin to move back to the starting point at $\epsilon$.

The integrand is analytic inside the contour, therefore $I = 0$. The sum of the two parts parallel to the real axis is:

$$I_r = \left[1-\exp(-2\pi)\right]\int_\epsilon^R \frac{\exp(i x)dx}{\exp(x) - 1}$$

So, the desired integral will follow from the imaginary part of $I_r$. The part of the contour integral from $R$ to $R + i$ tends to zero in the limit of $R\to\infty$ so this can be disregarded. The part of the contour integral along the imaginary axis can be written as:

$$I_i = -\int_{\epsilon}^{2\pi-\epsilon}\frac{\exp(-y)\exp\left(-i \frac{y}{2}\right)}{2\sin\left(\frac{y}{2}\right)}dy$$

We then see that

$$\lim_{\epsilon\to 0}\operatorname{Im}I_i = \frac{1-\exp(-2\pi)}{2}$$

The two quarter circles can be evaluated, we can borrow from the derivation of the residue theorem that each of them in the limit $\epsilon\to 0$ can be evaluated as $-\frac{\pi}{2} i$ times the residue at the poles that are at the centers of the quarter circles (unlike the case of a complete contour, this is only valid in the limit $\epsilon\to 0$). The residue at the pole at $z = 0$ is $1$, while the residue at $z = 2\pi i$ is $\exp(-2\pi)$. Having completed the evaluation of all parts of the contour integral, we can now equate the imaginary part of the contour integral to zero, and take the limit of $\epsilon\to 0$ and $R\to \infty$. This yields:

$$\int_0^{\infty}\frac{\sin(x)dx}{\exp(x)-1} = \frac{\pi}{2}\coth(\pi) - \frac{1}{2}$$