Una forma de resolver esto es reescribir la integral como
$$ \int_0 ^{ \infty } dx \, \frac {x \, e^{-x}}{1+e^{-x}} $$
y luego expandir el denominador en una serie de Taylor. La representación de esta integral se convierte
$$ \int_0 ^{ \infty } dx \, x \sum_ {k=0}^{ \infty } (-1)^k e^{-(k+1) x} $$
Debido a que la suma y la integral convergen, podemos cambiar el orden de la suma y la integración:
$$ \sum_ {k=0}^{ \infty } (-1)^k \int_0 ^{ \infty } dx \, x \, e^{-(k+1) x} = \sum_ {k=0}^{ \infty } \frac {(-1)^k}{(k+1)^2}$$
Otra forma de trabajar esto, de la que puede o no ser consciente, es usar el teorema de Cauchy. Considere el contorno integral
$$ \oint_C dz \frac {z^2}{1+e^z} $$
donde $C$ es el rectángulo con vértices $0$ , $R$ , $R+i 2 \pi $ y $i 2 \pi $ con un segmento semicircular de radio $ \epsilon $ centrado en $i \pi $ en el rectángulo. El contorno integral se convierte entonces
$$ \int_0 ^R dx \frac {x^2}{1+e^x} + i \int_0 ^{2 \pi } dy \frac {(R+i y)^2}{1+e^R e^{i y}} + \int_R ^0 dx \frac {(x+i 2 \pi )^2}{1+e^x} \\ +i \int_ {2 \pi }^{ \pi + \epsilon } dy \frac {(i y)^2}{1+e^{i y}}+ i \epsilon \int_ { \pi /2}^{- \pi /2} d \phi \, e^{i \phi } \frac {(i \pi + \epsilon e^{i \phi })^2}{1-e^{ \epsilon e^{i \phi }}} + i \int_ { \pi - \epsilon }^0 dy \frac {(i y)^2}{1+e^{i y}}$$
Como $R \to\infty $ la segunda integral se desvanece. Como $ \epsilon \to 0$ la quinta integral se convierte en
$$i \epsilon \int_ { \pi /2}^{- \pi /2} d \phi \, e^{i \phi } \frac {- \pi ^2}{- \epsilon e^{i \phi }} = -i \pi ^3$$
La integral del contorno es entonces
$$-i 4 \pi \int_0 ^{ \infty } dx \frac {x}{1+e^x} + 4 \pi ^2 \int_0 ^{ \infty } \frac {dx}{1+e^x} + i PV \int_0 ^{2 \pi } dy \frac {y^2}{1+e^{i y}}-i \pi ^3$$
Tengan en cuenta que $PV$ denota el valor principal de Cauchy. También hay que tener en cuenta que
$$ \begin {align}i PV \int_0 ^{2 \pi } dy \frac {y^2}{1+e^{i y}} &= i \frac12 \int_0 ^{2 \pi } dy \, y^2 + \frac12 PV \int_0 ^{2 \pi } dy \, y^2 \tan { \left ( \frac {y}{2} \right )} \\ &= i \frac {4 \pi ^3}{3} + \frac12 PV \int_0 ^{2 \pi } dy \, y^2 \tan { \left ( \frac {y}{2} \right )} \end {align}$$
Según el teorema de Cauchy, la integral del contorno es cero, lo que significa que tanto la parte real como la imaginaria son cero. Si la parte imaginaria de la integral de contorno es cero, tenemos
$$ \int_0 ^{ \infty } dx \frac {x}{1+e^x} = \frac { \pi ^2}{12} $$