integration - Simple, yet evasive integral from zero to $pi/2$

Q: Evaluate$\newcommand{\dx}{\mathrm dx}\newcommand{\du}{\mathrm du}\newcommand{\dv}{\mathrm dv}\newcommand{\dtheta}{\mathrm d\theta}\newcommand{\dw}{\mathrm dw}$$$I=\int\limits_0^{\pi/2}\dx\,\frac {\left(\log\sin x\log\cos x\right)^2}{\sin x\cos x}$$

I'm out of ideas on what to do. I tried using the astute limit identity for integration, but that lead to nowhere, because it's simply the same integral i.e adding the two together, doesn't yield a simplification whatsoever$$I=\int\limits_0^{\pi/2}\dx\,\frac {\log^2\sin x\log^2\cos x}{\sin x\cos x}=\int\limits_0^{\pi/2}\du\,\frac {\log^2\cos u\log^2\sin u}{\cos u\sin u}$$I've looked into making a u-substitution, but have no idea where to begin. I plan to use a trigonometric identity, namely$$\sec x\csc x=\cot x+\tan x$$but I'm not sure what to do afterwards because we get$$I=\int\limits_0^{\pi/2}\dx\,\cot x\log^2\sin x\log^2\cos x+\int\limits_0^{\pi/2}\dx\,\tan x\log^2\sin x\log^2\cos x$$I would appreciate it if you guys gave me an idea on where to begin!

Also, as a side note, I've used \newcommand on \dx,\du,\dv, and \dw to automatically change to $\dx$, $\du$, $\dv$, and $\dw$ respectively, if you guys don't mind! Just to make it easier for the differential operator!

Answer

Change variable to $t = \sin^2 x$ and notice
$$\begin{align}\frac{dx}{\sin x\cos x} &= \frac{\sin x\cos x dx}{\sin^2 x\cos^2 x} = \frac12\frac{dt}{t(1-t)}\\
\log\sin x \log\cos x &= \frac14\log\sin^2 x \log \cos^2 x = \frac14 \log t\log(1-t)
\end{align}

$$
The integral at hand can be rewritten as
$$\mathcal{I} \stackrel{def}{=} \int_0^{\pi/2} \frac{(\log\sin x\log\cos x)^2}{\sin x\cos x} dx = \frac{1}{32}\int_0^1 \frac{\log^2 t \log^2(1-t)}{t(1-t)} dt$$
Since $\displaystyle\frac{1}{t(1-t)} = \frac1t + \frac{1}{1-t}$, by replacing $t$ by $1-t$ in part of the integral, we find

$$\begin{align}
\mathcal{I}
&= \frac{1}{16}\int_0^1 \frac{\log^2 t\log^2(1-t)}{t} dt
= \frac{1}{48}\int_0^1 \log^2(1-t) d\log^3 t\\
&\stackrel{\text{I.by.P}}{=} \frac{1}{48}\left\{

\left[\log^2(1-t)\log^3 t\right]_0^1 + 2
\int_0^1 \log^3 t\frac{\log(1-t)}{1-t}dt\right\}\\
&= -\frac{1}{24}\int_0^1 \frac{\log^3 t}{1-t}\sum_{n=1}^\infty\frac{t^n}{n} dt
= -\frac{1}{24}\int_0^1 \log^3 t\sum_{n=1}^\infty H_nt^n dt\\
&\stackrel{t = e^{-y}}{=} \frac{1}{24} \int_0^\infty y^3 \sum_{n=1}^\infty H_n e^{-(n+1)y} dy
= \frac14 \sum_{n=1}^\infty \frac{H_n}{(n+1)^4}
\end{align}
$$
where $H_n$ are the $n^{th}$ harmonic number. By rearranging its terms, the last sum should be expressible in terms of zeta functions. I'm lazy, I just ask WA to evaluate the sum. As expected, last sum equals to $2\zeta(5) - \frac{\pi^2}{6}\zeta(3)$${}^\color{blue}{[1]}$.

As a consequence, the integral at hand equals to:

$$\mathcal{I}
= \frac{12\zeta(5) - \pi^2\zeta(3)}{24}
\approx 0.024137789997360933616411382857235691008...$$

Notes

• $\color{blue}{[1]}$ It turns out we can compute this sum using an
  
  identity by Euler.
  $$2\sum_{n=1}^\infty \frac{H_n}{n^m} = (m+2)\zeta(m+1) - \sum_{n=1}^{m-2}\zeta(m-n)\zeta(n+1),\quad\text{ for } m = 2, 3, \ldots$$
  In particular, for $m = 4$, this identity becomes
  $$\sum_{n=1}^\infty \frac{\zeta(n)}{n^4}
  = 3\zeta(5) - \zeta(2)\zeta(3)$$
  and we can evaluate our sum as
  $$\begin{align}\mathcal{I}
  &= \frac14\sum_{n=1}^{\infty}\frac{H_{n}}{(n+1)^4}
  = \frac14\sum_{n=1}^{\infty}\left(\frac{H_{n+1}}{(n+1)^4}-\frac{1}{(n+1)^5}\right)
  = \frac14\left(\sum_{n=1}^\infty \frac{H_n}{n^4} - \zeta(5)\right)\\
  
  &= \frac14(2\zeta(5) - \zeta(2)\zeta(3))
  = \frac{12\zeta(5) - \pi^2\zeta(3)}{24}
  \end{align}
  $$