Prove that $\int_0^{\pi/2}\ln^2(\cos x)\,dx=\frac{\pi}{2}\ln^2 2+\frac{\pi^3}{24}$ 
Prove that
  \begin{equation}
\int_0^{\pi/2}\ln^2(\cos x)\,dx=\frac{\pi}{2}\ln^2 2+\frac{\pi^3}{24}
\end{equation}


I tried to use by parts method and ended with
\begin{equation}
\int \ln^2(\cos x)\,dx=x\ln^2(\cos x)+2\int x\ln(\cos x)\tan x\,dx
\end{equation}
The latter integral seems hard to evaluate. Could anyone here please help me to prove it preferably with elementary ways (high school methods)? Any help would be greatly appreciated. Thank you.
Addendum:
I also found this nice closed-form
\begin{equation}
-\int_0^{\pi/2}\ln^3(\cos x)\,dx=\frac{\pi}{2}\ln^3 2+\frac{\pi^3}{8}\ln 2 +\frac{3\pi}{4}\zeta(3)
\end{equation}
I hope someone here also help me to prove it. (>‿◠)✌

We are all familiar with the famous Wallis integrals, $W_n=\displaystyle\int_0^\frac\pi2\sin^{n-1}x~dx=\int_0^\frac\pi2\cos^{n-1}x~dx$ 
$=\dfrac12B\bigg(\dfrac12,~\dfrac n2\bigg)$, see beta function for more details. It follows then that our integral is nothing 
other than $W''(1)$, which can be evaluated in terms of the digamma and trigamma functions, for 
arguments $\dfrac12$ and $1$. The former can be expressed as $\psi_{_0}(k+1)=H_k-\gamma$, where $\gamma\approx\dfrac1{\sqrt3}~$ is the 
Euler-Mascheroni constant, and $H_k=\displaystyle\int_0^1\frac{1-x^k}{1-x~~}dx=-\int_0^1\ln\Big(1-\sqrt[^k]x\Big)~dx$ is the harmonicnumber. Thus, that $H_1=1$ is self-evident, and that $H_\frac12=2~(1-\ln2)$ can be shown by a simple 
substitution. The values are also found here. For the latter we have $~\psi_{_1}(x)=\displaystyle\sum_{k=0}^\infty\frac1{(k+x)^2}$, which 
implies $\psi_{_1}(1)=\zeta(2)=\dfrac{\pi^2}6$, and $\psi_{_1}\bigg(\dfrac12\bigg)=4\bigg(1-\dfrac14\bigg)\zeta(2)=\dfrac{\pi^2}2$. See Basel problem for more 
information.
