Para utilizar un contorno rectangular, considere la integral
$$\oint_C dz \frac{z^2}{e^z+1}$$
donde $C$ es el contorno rectangular que tiene vértices en $0$, $R$, $R+i 2 \pi$, y $i 2 \pi$, con una forma semicircular, desvío en el rectángulo de radio $\epsilon$$z=i \pi$. El contorno de la integral es entonces igual a
$$\int_0^R dx \frac{x^2}{e^x+1} + i \int_0^{2 \pi} dy \frac{(R+i y)^2}{e^{R+i y}+1} \\ + \int_R^0 dx \frac{(x+i 2 \pi)^2}{e^{x+i 2 \pi}+1} + i \int_{2 \pi}^{\pi+\epsilon} dy \frac{(i y)^2}{e^{i y}+1}\\ + i \epsilon \int_{\pi/2}^{-\pi/2} d\phi \, e^{i \phi} \frac{(i \pi+\epsilon e^{i \phi})^2}{e^{i \pi+\epsilon e^{i \phi}}+1}+i \int_{\pi-\epsilon}^0 dy \frac{(i y)^2}{e^{i y}+1}$$
Consideramos que el límite de $R\to\infty$$\epsilon \to 0$. Como $R\to\infty$, la segunda integral se desvanece. Como $\epsilon\to 0$, el quinto enfoques integrales
$$i \epsilon \int_{\pi/2}^{-\pi/2} d\phi \, e^{i \phi} \frac{-\pi^2}{-\epsilon e^{i \phi}}=-i \pi^3$$
Por Cauchy teorema, el contorno de la integral es cero. Tenemos, entonces, la expansión de la primera y la tercera de las integrales y la combinación de la cuarta y la sexta de las integrales en un valor principal de Cauchy:
$$-i 4 \pi \int_0^{\infty} dx \frac{x}{e^x+1} +4 \pi^2 \int_0^{\infty} \frac{dx}{e^x+1}\\ +i \,PV \int_0^{2 \pi} dy \frac{y^2}{e^{i y}+1}-i \pi^3=0$$
donde $PV$ denota el valor principal de Cauchy. Tenga en cuenta que
$$PV \int_0^{2 \pi} dy \frac{y^2}{e^{i y}+1} = \frac12 \int_0^{2 \pi} dy \: y^2 -i \frac12 \, PV \int_0^{2 \pi} dy \frac{y^2 \sin{y}}{1+\cos{y}}$$
La equiparación de la imaginaria, vemos que
$$-4 \pi \int_0^{\infty} dx \frac{x}{e^x+1} + \frac12 \frac{(2 \pi)^3}{3} - \pi^3=0$$
o
$$\int_0^{\infty} dx \frac{x}{e^x+1} = \frac13 \pi^2-\frac14\pi^2 = \frac{\pi^2}{12}$$