Prove $\int_0^{\infty} \frac{x^2}{\cosh^2 (x^2)} dx=\frac{\sqrt{2}-2}{4} \sqrt{\pi}~ \zeta \left( \frac{1}{2} \right)$ Wolfram Alpha evaluates this integral numerically as 
$$\int_0^{\infty} \frac{x^2}{\cosh^2 (x^2)} dx=0.379064 \dots$$
Its value is apparently
$$\frac{\sqrt{2}-2}{4} \sqrt{\pi}~ \zeta \left( \frac{1}{2} \right)=0.37906401072\dots$$

How would you solve this integral?

Obviously, we can make a substitution $t=x^2$
\begin{align}
\int_0^{\infty} \frac{x^2}{\cosh^2 (x^2)} dx&=\frac{1}{2} \int_0^{\infty} \frac{\sqrt{t}}{\cosh^2 (t)} dt\\[10pt]
&=\int_0^{\infty} \frac{\sqrt{t}}{\cosh (2t)+1} dt\\[10pt]
&=\frac{1}{2 \sqrt{2}}\int_0^{\infty} \frac{\sqrt{u}}{\cosh (u)+1} du
\end{align}
We could use geometric series since $\cosh (u) \geq 1$, but I don't know how it will help.

Hint:
Consider the parametric integral
\begin{equation}
I(a) =  \int_0^\infty \frac{t^{a-1}}{\cosh^{2} t}\ dt=4 \int_{0}^{\infty} \frac{t^{a-1}}{(e^{t}+e^{-t})^{2}}\ dt = 4 \int_{0}^{\infty}\frac{t^{a-1} e^{-2t}}{(1+e^{-2t})^{2}}\ dt
\end{equation}
Hence, your integral is simply
\begin{equation}
\int_0^{\infty}  \frac{x^2}{\cosh^2x^2}\ dx = -2\ \frac{\partial}{\partial b}\left[\int_0^{ \infty}\frac{t^{a-1}}{1+e^{bt}} \ dt\right]_{a=\frac{1}{2}\ ,\ b=2}
\end{equation}
I believe you can evaluate the last expression by your own.
