Playing with Wolfram Alpha and inspired in [1] (I refers it if someone know how relates my problem with some of problems involving the Apéry constant in this reference, but the relation doesn't seem explicit), defining $$I_n:=-\int_0^1\frac{\log(1+x^{2n})\log x}{x}dx$$ for integers $n\geq 1$, I can calculate, as I am saying with Wolfram Alpha (but I don't know how get the indefinite integrals) $I_1$, $I_2$, $I_3$ and $I_4$. And as a conjecture $$I_8=\frac{6\zeta(3)}{8\cdot 16^2}.$$
Motivation. I would like to do a comparison with the sequence $I_1$, $I_2$, $I_3$, $I_4$ and $I_8$.
Question. If do you know that this problem was solved in the literature please add a comment: can you evaluate in a closed-form $I_5$? Many thanks.
References:
[1] Walther Janous , Around's Apéry's constant, J. Ineq. Pure and Appl. Math. 7(1) Art. 35 (2006).
By Taylor series expansion ,
For $0<x<1$ and $n\geq 1$,
$\displaystyle -\ln(1+x^{2n})=\sum_{k=1}^{+\infty} \dfrac{(-1)^kx^{2kn}}{k}$
For $k \geq 0$,
$\displaystyle \int_0^1 x^k\ln x dx=-\dfrac{1}{(k+1)^2}$
(integration by parts)
Therefore,
$\begin{align} I_n&=\int_0^1 \left(\sum_{k=1}^{+\infty} \dfrac{(-1)^kx^{2kn-1}\ln x}{k}\right) dx\\ &=\sum_{k=1}^{+\infty} \left(\int_0^1 \dfrac{(-1)^kx^{2kn-1}\ln x}{k} dx\right)\\ &=-\sum_{k=1}^{+\infty} \dfrac{(-1)^k}{k(2kn)^2}\\ &=-\dfrac{1}{4n^2}\sum_{k=1}^{+\infty} \dfrac{(-1)^k}{k^3}\\ &=-\dfrac{1}{4n^2}\left(\sum_{k=1}^{+\infty} \dfrac{1}{(2k)^3}-\sum_{k=0}^{+\infty}\dfrac{1}{(2k+1)^3}\right)\\ &=-\dfrac{1}{4n^2}\left(\dfrac{1}{8}\zeta(3)-\left(\zeta(3)-\dfrac{1}{8}\zeta(3)\right)\right)\tag{1}\\ &=-\dfrac{1}{4n^2}\times -\dfrac{3}{4}\zeta(3)\\ &=\boxed{\dfrac{3\zeta(3)}{16n^2}} \end{align}$
For (1) observe that,
$\begin{align}\sum_{k=0}^{+\infty}\dfrac{1}{(2k+1)^3}&=\sum_{k=1}^{+\infty}\dfrac{1}{k^3}-\sum_{k=1}^{+\infty}\dfrac{1}{(2k)^3}\\ &=\zeta(3)-\dfrac{1}{8}\zeta(3) \end{align}$