I'm struggling with this integral $$I=\int_0^1\frac{1-x^2+\left(1+x^2\right)\ln x}{\left(x+x^2\right)\ln^3x}dx.\tag1$$ Mathematica could not evaluate it in a closed form. Its numeric value is approximately $I\approx0.7804287418294087023386965471512328112...$$^\text{[more]}$, but I could not find a plausible closed form for this number using inverse symbolic calculators available online.
Update: Based on Raymond Manzoni's comment below, there is actually a conjectural closed form, numerically matching up to at least $10^3$ decimal digits: $$I\stackrel?=6\ln A-\frac{\ln4}3-\frac14,\tag2$$ where $A$ is the Glaisher-Kinkelin constant: $$A=\exp\left(\frac1{12}-\zeta'(-1)\right).\tag3$$ Could you suggest a proof of the conjecture $(2)$?
$$\begin{align*} I&=\int^1_0\frac{1-x^2+(1+x^2)\log x}{x+1}\frac{dx}{x\log^3x}\\ &=\left.\frac{-1}{2\log^2x}\frac{1-x^2+(1+x^2)\log x}{x+1}\right|^1_0-\int^1_0\frac{-1}{2\log^2x}\frac{\partial}{\partial x}\left(\frac{1-x^2+(1+x^2)\log x}{x+1}\right)dx\\ &=\int^1_0\frac{1}{2\log^2x}\frac{\partial}{\partial x}\left(1-x+\frac{(1+x^2)\log x}{x+1}\right)dx\\ &=\frac12\int^1_0\frac{1}{\log^2x}\left(-1+\frac{(1+x^2)}{x(x+1)}+\frac{(x^2+2x-1)\log x}{(x+1)^2}\right)dx\\ &=\frac12\int^1_0\left(\frac{1-x}{x(x+1)\log^2x}+\frac{1}{\log x}-\frac{2}{(x+1)^2\log x}\right)dx\\ &=\frac12\int^1_0\left(\frac{1-x}{x(x+1)\log^2x}+\frac{2}{(x+1)^2\log x}\right)dx+\frac12\int^1_0\left(\frac{1}{\log x}-\frac{4}{(x+1)^2\log x}\right)dx\\ &=\frac12(I_1+I_2). \end{align*}$$
$$\begin{align*} I_1&=\int^1_0\left(\frac{1-x}{x(x+1)\log^2x}+\frac{2}{(x+1)^2\log x}\right)dx\\ &=\int^1_0\left(\frac{1-x^2+2x\log x}{(x+1)^2}\right)\frac{dx}{x\log^2x}\\ &=\left.\frac{-1}{\log x}\frac{1-x^2+2x\log x}{(x+1)^2}\right|^1_0-\int^1_0\frac{-1}{\log x}\frac{\partial}{\partial x}\left(\frac{1-x^2+2x\log x}{(x+1)^2}\right)dx\\ &=\int^1_0\frac{-1}{\log x}\frac{-2(1-x)\log x}{(x+1)^3}dx\\ &=2\int^1_0\frac{1-x}{(x+1)^3}dx=\frac12. \end{align*}$$
Mathematica is able to evaluate $I_2=12\log A-1-\frac43\log2$, so the conjectured closed form is proved.
Edit: to evaluate $I_2$, we write $F(a)=\int^1_0\frac{x^a-1}{(x+1)^2\log x}dx$, so that $F(0)=0$, and $F'(a)=\int^1_0\frac{x^a}{(x+1)^2}dx$. We have $$\begin{align*} I_2&=\int^1_0\left(\frac{1}{\log x}-\frac{4}{(x+1)^2\log x}\right)dx\\ &=\int^1_0\frac{x^2+2x-3}{(x+1)^2\log x}dx=F(2)+2F(1). \end{align*}$$
Then we need to give closed form of $F(a)$. First we have $$\begin{align*} \int^1_0\frac{x^adx}{(1-zx)^2}&=\int^1_0\sum_{n=0}^{\infty}(n+1)z^nx^{n+a}dx\\ &=\sum_{n=0}^{\infty}\frac{(n+1)z^n}{n+a+1}\\ &=\sum_{n=0}^{\infty}z^n-a\sum_{n=0}^{\infty}\frac{z^n}{n+a+1}\\ &=\frac{1}{1-z}-\frac{a}{a+1}~_2F_1(1,a+1;a+2\mid z) \end{align*}$$ for all $|z|<1$. Taking limits, we have $$\begin{align*} F'(a)&=\frac12-\frac{a}{a+1}~_2F_1(1,a+1;a+2\mid-1)\\ &=\frac12-\frac{a}{2}\left(\psi\left(\frac{a+2}{2}\right)-\psi\left(\frac{a+1}{2}\right)\right) \end{align*}$$ and therefore $$\begin{align*} F(a)&=\int^a_0F'(b)~db\\ &=\frac{a}{2}-\int^a_0\frac{b}{2}\left(\psi\left(\frac{b+2}{2}\right)-\psi\left(\frac{b+1}{2}\right)\right)db\\ &=\frac{a}{2}-2\int^{a/2}_{0}c\left(\psi(c+1)-\psi(c+\frac12)\right)dc\\ &=\frac{a}{2}+\int^{a/2}_{0}(2\psi(c+1)-\psi(c+\frac12))dc-2\int^{a/2}_{0}\left((c+1)\psi(c+1)-(c+\frac12)\psi(c+\frac12)\right)dc\\ &=\frac{a}{2}+(2\log\Gamma(a/2+1)-2\log\Gamma(1)-\log\Gamma(a/2+1/2)+\log\Gamma(1/2))\\ &-2\left(\int^{a/2+1}_{1}-\int^{a/2+1/2}_{1/2}\right)c\psi(c)dc\\ &=\frac{a}{2}+\log\Gamma(1/2)+2\log\Gamma(a/2+1)-\log\Gamma(a/2+1/2)-2\left.(c\log\Gamma(c)-\psi^{(-2)}(c))\right|^{a/2+1,1/2}_{a/2+1/2,1}\\ &=\frac{a}{2}+\log\Gamma(1/2)+2\log\Gamma(a/2+1)-\log\Gamma(a/2+1/2)-(a+2)\log\Gamma(a/2+1)\\&+(a+1)\log\Gamma(a/2+1/2)-\log\Gamma(1/2)+2\left.(\psi^{(-2)}(c))\right|^{a/2+1,1/2}_{a/2+1/2,1}\\ &=\frac{a}{2}-a\log\Gamma(a/2+1)+a\log\Gamma(a/2+1/2)+2\left.(\psi^{(-2)}(c))\right|^{a/2+1,1/2}_{a/2+1/2,1}\\ &=\frac{a}{2}-a\log\Gamma(a/2+1)+a\log\Gamma(a/2+1/2)\\ &+2\left(\psi^{(-2)}(a/2+1)-\psi^{(-2)}(a/2+1/2)+\psi^{(-2)}(1/2)-\psi^{(-2)}(1)\right).\\ \end{align*}$$ Here $\psi^{(-2)}(a)=\int^a_0\log\Gamma(x)dx$. We can now give a closed form for $I_2$: $$\begin{align*} I_2&=F(2)+2F(1)\\ &=(1-2\log\Gamma(2)+2\log\Gamma(3/2))+2\left(1/2-\log\Gamma(3/2)+\log\Gamma(1)\right)\\ &+2(\psi^{(-2)}(2)-\psi^{(-2)}(3/2)+\psi^{(-2)}(1/2)-\psi^{(-2)}(1))\\ &+4(\psi^{(-2)}(3/2)-\psi^{(-2)}(1)+\psi^{(-2)}(1/2)-\psi^{(-2)}(1))\\ &=2+2\psi^{(-2)}(2)+2\psi^{(-2)}(3/2)+6\psi^{(-2)}(1/2)-10\psi^{(-2)}(1). \end{align*}$$ Using the known values $\psi^{(-2)}(2)=\log2+\log\pi-1$, $\psi^{(-2)}(1)=\frac12\log2+\frac12\log\pi$, $\psi^{(-2)}(3/2)=\frac34\log\pi+\frac{5}{24}\log2+\frac32\log A-\frac12$ and $\psi^{(-2)}(1/2)=\frac14\log\pi+\frac{5}{24}\log 2+\frac32\log A$, we finally conclude that $I_2=12\log A-1-\frac43\log2$, and therefore $$I=\frac12(I_1+I_2)=6\log A-\frac23\log2-\frac14.$$