Calcular $$\sf I=\int_0^1 \frac{\ln(1+x+x^2)\ln(1-x+x^2)}{x}dx$$
He pensado en la integral de arriba después de ver esta similar integral y creo que cambiar el cartel para tener $\sf 1+x+x^2$ podría conseguirnos una buena forma cerrada.
Así que empecé a utilizar la siguiente fórmula: $$\sf 2ab=(a+b)^2-a^2-b^2$$ $$\sf \Rightarrow 2I=\int_0^1\frac{\ln^2(1+x^2+x^4)}{x}dx-\int_0^1\frac{\ln^2(1+x+x^2)}{x}dx-\int_0^1\frac{\ln^2(1-x+x^2)}{x}dx$$ Utilizando en la primera integral la sustitución $\sf x^2\rightarrow x $ nos atrapa: $$\sf \int_0^1\frac{\ln^2(1+x^2+x^4)}{x}dx=\frac12\int_0^1\frac{\ln^2(1+x+x^2)}{x}dx$$ $$\sf \Rightarrow I=-\frac14\int_0^1\frac{\ln^2(1+x+x^2)}{x}dx-\frac12\int_0^1\frac{\ln^2(1-x+x^2)}{x}dx$$ Bueno, ahora sólo tenemos que encontrar: $$\sf I(a)=\int_0^1\frac{\ln^2(1+ax+x^2)}{x}dx $$ A continuación, establezca $a=1$ y $a=-1$ . Por supuesto, intenté usar el truco de Feynman: $$\sf I'(a)=2\int_0^1\frac{\ln(1+ax+x^2)}{1+ax+x^2}dx$$ Pero rápidamente se dio por vencido ya que no parece prometedor.
Otra forma podría ser dejar que $\sf x+\frac12=\frac{\sqrt 3}{2}t$ para conseguirlo: $$\sf \int_0^1\frac{\ln^2(1+x+x^2)}{x}dx=\int_\frac{1}{\sqrt 3}^\sqrt 3 \frac{\ln^2\left(\frac34(1+t^2)\right)}{t-\frac{1}{\sqrt 3}}dt$$ Pero bueno ¡apreciaría un poco de ayuda!
Actualización . Mientras tanto he encontrado algo, aunque sólo es un conjetura y aún no estaba probado, podríamos haberlo hecho: $$\int_0^1\frac{\ln^2(1+x+x^2)}{x}dx=\frac{2\pi}{9\sqrt3}\psi_1\left(\frac13\right)-\frac{4\pi^3}{27\sqrt3}-\frac23\zeta(3)$$ Y esto me da aún más esperanzas de que haya una buena forma cerrada. No digo que sea fácil de demostrar, pero apuesto a que hay mucho que aprender de una aproximación (no resultados de Mathematica).