El problema es mostrar que la integral impropia $I = \int_0^{\infty} \frac{x^2}{e^x-1} dx$ converge a $2\sum_1^{\infty} \frac1{n^3}$.
Anteriormente, he calcula la siguiente integral:
$$f(x) = -\int_0^x\frac{\log(1-t)}{t}dt = \sum_1^{\infty} \frac{x^n}{n^2}$$
que tiene por $-1 \le x \le 1$. A partir de esto, deduzco lo siguiente:
$$g(x) = \int_0^x \frac{f(t)}{t}dt = \sum_1^{\infty} \frac{x^n}{n^3}$$
que tiene también para $-1 \le x \le 1$. Estoy tratando de demostrar que $I = g(1)$.
La primera cosa que voy a hacer con $I$ es hacer la sustitución $t = 1 - e^{-x}$, lo que se obtiene:
$$I = \int_0^1 \frac{\log^2(1-t)}{t} dt$$
Entonces yo intento de integrar por partes, distinguiendo $\log(1-t)$ y la integración de $\frac{\log(1-t)}{t}$, lo que se obtiene:
$$I = -\log(1-t)f(t)|_0^1-\int_0^1\frac{f(t)}{1-t}dt $$
Que no funciona. El primer término no convergen. El límite superior se evalúa a $-\log(0)\frac{pi^2}6$ y el segundo término no es g(1). El denominador debe ser $t$, no $t-1$. Así que no estoy seguro de qué hacer a continuación.