Voy a seguir adelante y publicar mi método. Es similar a xpauls excepto yo solía digamma, que está relacionada con la serie armónica de todos modos.
Romper integral:
$$\int_{0}^{1}\frac{\log^{2}(x)\log(1+x)}{(1-x)(x^{2}+1)}dx+\int_{1}^{\infty}\frac{\log^{2}(x)\log(1+x)}{(1-x)(x^{2}+1)}dx$$
En el derecho integral, hacer el sub $x=1/t$. Esto nos da:
$$\int_{0}^{1}\frac{\log^{2}(x)\log(1+x)}{(x^{2}+1)}dx+\int_{0}^{1}\frac{x\log^{3}(x)}{(1-x)(x^{2}+1)}dx$$
El derecho integral:
Se rompen en $$1/2\int_{0}^{1}\frac{x\log^{3}(x)}{x^{2}+1}dx-1/2\int_{0}^{1}\frac{\log^{3}(x)}{x^{2}+1}dx+1/2\int_{0}^{1}\frac{\log^{3}(x)}{1-x}dx$$
Yo no soy de ir a trabajar a través de cada uno de estos. Pero, baste decir, que se puede hacer sin demasiado esfuerzo mediante el uso de series geométricas. Por ejemplo, la de en medio:
$$1/2\int_{0}^{1}\log^{3}(x)\sum_{k=0}^{\infty}(-1)^{k}x^{2k}dx=3\sum_{k=0}^{\infty}\frac{(-1)^{k}}{(2k+1)^{4}}$$
Haciendo así que a los tres nos lleva a la serie que evaluar en términos de $\zeta(4)$ y $\psi_{3}$. Sumando a ellos los resultados en:
$$ \boxed{\displaystyle \int_{0}^{1}\frac{x\log^{3}(x)}{(1-x)(x^{2}+1)}dx=\frac{-9\pi^{4}}{256}+\frac{1}{512}\left[\psi_{3}(1/4)-\psi_{3}(3/4)\right]}$$
La izquierda integral hasta la parte superior es un poco más difícil. Al menos así lo creo.
$$\int_{0}^{1}\frac{\log^{2}(x)\log(1+x)}{x^{2}+1}dx$$
Utilizar la serie de Taylor para $\log(1+x)$:
$$\int_{0}^{1}\frac{\log^{2}(x)}{x^{2}+1}\sum_{n=1}^{\infty}\frac{(-1)^{n}x^{n}}{n}$$
Nota la función Beta incompleta se define como: $\displaystyle \int_{0}^{1}\frac{x^{a}}{x^{2}+1}dx=1/4\left[\psi \left(\frac{a+3}{4}\right)-\psi\left(\frac{a+1}{4}\right)\right]$.
Distinción esto dos veces w.r.t 'a' introduce el registro-cuadrado plazo y da:
$$\int_{0}^{1}\frac{x^{a+n}\log^{2}(x)}{x^{2}+1}dx=1/64\left[\psi_{2} \left(\frac{a+n+3}{4} \right)-\psi_{2} \left(\frac{a+n+1}{4} \right) \right]$$.
Por lo tanto, dejar que $a=0$, $$\int_{0}^{1}\frac{\log^{2}(x)\log(1+x)}{x^{2}+1}dx=1/64\sum_{n=1}^{\infty}\frac{(-1)^{n}}{n}\left[\psi_{2}\left(\frac{n+3}{4}\right)-\psi_{2}\left(\frac{n+1}{4}\right)\right]$$
$$=\boxed{\displaystyle \frac{\pi^{2}}{6}G+\frac{\pi^{3}}{32}\log(2)-\frac{1}{768}\left[\psi_{3}\left(1/4\right)-\psi_{3}\left(3/4\right)\right]}$$
Esta serie de resultados, cuando se combina con la otra caja resultado, se obtiene la solución a la integral original.
El único problema menor que tengo es la evaluación de esta tetragamma de la serie. Como he dicho, La Flajolet-Salvy residuo método puede funcionar, pero encontrar el núcleo correcto es la primera tarea importante. Ya que los suplentes, quiero suponer algo con $\pi \csc(\pi z)$
Por supuesto, uno podría simplemente decir que el diablo con él y utilizar esto como un lexema. Pero, me gustaría evaluar.