How to prove that:
$$\int_0^1\frac{\ln x\ln(1+x^2)}{1-x^2}dx=\frac74\zeta(3)-\frac34\ln2 \zeta(2)-\frac{\pi}{2}G$$
where $\zeta$ is the Riemann zeta function and $G$ is Catalan constant.
I came across this integral while working on evaluating some harmonic series.
I am tagging "harmonic series" as its pretty related to logarithmic integrals.
We will start by using the following substitution: $$\frac{1-x}{1+x}=t\Rightarrow x=\frac{1-t}{1+t}\Rightarrow dx=\frac{2}{(1+t)^2}dt$$ $$\Rightarrow I=\int_0^1 \frac{\ln x\ln(1+x^2)}{1-x^2}dx=\frac12\int_0^1 \frac{[\ln(1-t)-\ln(1+t)][\ln2+\ln(1+t^2)-2\ln(1+t)]}{t}dt$$ Now we are going to use the following result to evaluate a part from above: $$\small \int_0^1 \frac{[m\ln(1+x)+n\ln(1-x)][q\ln(1+x)+p\ln(1-x)]}{x}dx=\left(\frac{mq}{4}-\frac{5}{8}(mp+nq)+2np\right)\zeta(3)$$ $$\Rightarrow I=\frac{7}{8}\zeta(3)-\frac34\zeta(2)\ln 2+\frac12{\int_0^1 \frac{[\ln(1-t)-\ln(1+t)]\ln(1+t^2)}{t}dt}$$ The last integral is $I-J=\frac74\zeta(3)-\pi G$ which appears in the following post. $$\Rightarrow I =\frac{7}{8}\zeta(3)-\frac34\zeta(2)\ln 2+\frac78\zeta(3)-\frac{\pi}{2}G=\frac74\zeta(3)-\frac34\zeta(2)\ln 2-\frac{\pi}{2}G$$