Podemos utilizar una integral de contorno en el plano complejo, tal como he mostrado aquí para el caso de $n=3$. Ahora, sin embargo, utilizamos
$$\oint_C dz \frac{\log^2{z}}{z^n-1}$$
donde $C$ es la modificación de la cerradura de contorno sobre el eje real positivo de radio exterior $R$ y un radio interior de $\epsilon$.
Vamos a evaluar esta integral sobre los contornos. Como antes, hay $8$ piezas para evaluar, de la siguiente manera:
$$\int_{\epsilon}^{1-\epsilon} dx \frac{\log^2{x}}{x^n-1} + i \epsilon \int_{\pi}^0 d\phi \, e^{i \phi} \frac{\log^2{\left (1+\epsilon e^{i \phi}\right )}}{(1+\epsilon e^{i \phi})^n-1} \\ + \int_{1+\epsilon}^R dx \frac{\log^2{x}}{x^n-1} + i R \int_0^{2 \pi} d\theta \, e^{i \theta} \frac{\log^2{\left (R e^{i \theta}\right )}}{R^n e^{i n \theta}-1} \\ + \int_R^{1+\epsilon} dx \frac{(\log{x}+i 2 \pi)^2}{x^n-1} + i \epsilon \int_{2 \pi}^{\pi} d\phi \, e^{i \phi} \frac{(\log{\left (1+\epsilon e^{i \phi}\right )}+i 2 \pi)^2}{(1+\epsilon e^{i \phi})^n-1} \\ + \int_{1-\epsilon}^{\epsilon} dx \frac{(\log{x}+i 2 \pi)^2}{x^n-1} + i \epsilon \int_{2 \pi}^0 d\phi \, e^{i \phi} \frac{\log^2{\left (\epsilon e^{i \phi}\right )}}{\epsilon^n e^{i n \phi}-1} $$
Como $R \to \infty$, el cuarto integral desvanece como $\log^2{R}/R^{n-1}$. Como $\epsilon \to 0$, la segunda integral se desvanece como es $O(\epsilon^n)$, mientras que el octavo integral desvanece como $\epsilon \log^2{\epsilon}$. Esto deja a la primera, tercera, quinta, sexta y séptima de las integrales, que en los límites anteriores, se convierten en
$$PV \int_0^{\infty} dx \frac{\log^2{x} - (\log{x}+i 2 \pi)^2}{x^n-1} + i \frac{4 \pi^3}{n}$$
El residuo de cálculo es un poco más complicado, porque ahora tenemos $n-1$ polos en los que necesitamos para evaluar los residuos. El contorno de la integral es así
$$i 2 \pi \sum_{k=1}^{n-1} \frac{-4 \pi^2 k^2/n^2}{n e^{i 2 (n-1) \pi k/n}} = -i \frac{8 \pi^3}{n^3} \sum_{k=1}^{n-1} k^2 \, e^{-i 2 (n-1) \pi k/n} $$
La suma es factible, pero el álgebra es un poco horrible. El resultado es
$$\sum_{k=1}^{n-1} k^2 \, e^{-i 2 (n-1) \pi k/n} = \frac12 \left ( \frac{n}{\sin^2{\frac{\pi}{n}}} - n^2\right ) - i \frac12 n^2 \cot{\frac{\pi}{n}} $$
Igualando las partes reales e imaginarias de ambas ecuaciones para el contorno integral de los rendimientos
$$\int_0^{\infty} dx \frac{\log{x}}{x^n-1} = \frac{\pi^2}{n^2 \sin^2{\frac{\pi}{n}}} $$
$$PV \int_0^{\infty} dx \frac{1}{x^n-1} = -\frac{\pi}{n} \cot{\frac{\pi}{n}} $$