Integrate $\int_0^\pi \theta^2 \ln^2\big(2\cosh\frac{\theta}{2}\big)d \theta$ Hello I am trying to integrate
$$
I=\int_0^\pi \theta^2 \ln^2\big(2\cosh\frac{\theta}{2}\big)d \theta
$$
which is similar to Integral $\int_0^\pi \theta^2 \ln^2\big(2\cos\frac{\theta}{2}\big)d \theta$., however we have hyperbolic functions in this example.
We can possibly write the integrand as 
$$
\theta^2 \ln^2\big(2\cosh\frac{\theta}{2}\big)=\theta^2\big(\ln 2+\ln \cosh \frac{\theta}{2}\big)^2=\theta^2\left(\ln^2(2)+\ln^2\cosh\frac{\theta}{2}+2\ln (2)\ln \cosh\big(\frac{\theta}{2}\big)\right).
$$
We can now write
$$
I=\ln^2(2)\int_0^\pi \theta^2d\theta  +\int_0^\pi\theta^2 \ln^2 \cosh \frac{\theta}{2}d\theta+2\ln 2 \int_0^\pi\theta^2 \ln \cosh{\frac{\theta}{2}}d\theta.
$$
Change of variables $y=\theta/2$ we obtain
$$
I=\frac{\pi^3\ln^2(2)}{3}+16\ln(2)\int_0^{\pi/2} y^2 \ln \cosh (y) dy+8\int_0^{\pi/2} y^2 \ln^2 \cosh (y) dy.
$$
I am stuck here, thanks.
 A: $$\int_0^\pi\theta^2\ln^2\left(2\cosh\dfrac{\theta}{2}\right)d\theta$$
$$=\int_0^\pi\theta^2(\ln(e^\frac{\theta}{2}+e^{-\frac{\theta}{2}}))^2~d\theta$$
$$=\int_0^\pi\theta^2(\ln(e^\frac{\theta}{2}~(1+e^{-\theta})))^2~d\theta$$
$$=\int_0^\pi\theta^2\left(\dfrac{\theta}{2}+\ln(1+e^{-\theta})\right)^2~d\theta$$
$$=\int_0^\pi\left(\dfrac{\theta^4}{4}+\theta^3\ln(1+e^{-\theta})+\theta^2(\ln(1+e^{-\theta}))^2\right)d\theta$$
$$=\int_0^\pi\left(\dfrac{\theta^4}{4}+\sum\limits_{n=1}^\infty\dfrac{(-1)^n\theta^3e^{-n\theta}}{n}+\sum\limits_{n=1}^\infty\sum\limits_{k=1}^\infty\dfrac{(-1)^{n+k}\theta^2e^{-(n+k)\theta}}{nk}\right)d\theta$$
$$=\left[\dfrac{\theta^5}{20}-\sum\limits_{n=1}^\infty\dfrac{(-1)^n(n^3\theta^3+3n^2\theta^2+6n\theta+6)e^{-n\theta}}{n^5}-\sum\limits_{n=1}^\infty\sum\limits_{k=1}^\infty\dfrac{(-1)^{n+k}((n+k)^2\theta^2+2(n+k)\theta+2)e^{-(n+k)\theta}}{nk(n+k)^3}\right]_0^\pi$$ (according to http://en.wikipedia.org/wiki/List_of_integrals_of_exponential_functions)
$$=\dfrac{\pi^5}{20}-\sum\limits_{n=1}^\infty\dfrac{(-1)^n((\pi^3n^3+3\pi^2n^2+6\pi n+6)e^{-\pi n}-6)}{n^5}-\sum\limits_{n=1}^\infty\sum\limits_{k=1}^\infty\dfrac{(-1)^{n+k}((\pi^2(n+k)^2+2\pi(n+k)+2)e^{-\pi(n+k)}-2)}{nk(n+k)^3}$$
