I need to solve the integral
$$\int_{0}^{1}\int_{0}^{1} \sqrt{1 + 4(x^2 + y^2)}\,dx\,dy$$
I am using polar coordinates here to get :
$$ \int_{0} ^{\pi/4}\int_{0}^{\sec \theta} \sqrt{(1 + 4r^2)} r \,dr \,d\theta + \int_{\pi/4} ^{\pi/2}\int_{0}^{\operatorname{cosec}\theta} \sqrt{(1 + 4r^2)} r \,dr \,d\theta$$
After this integral becomes too complex to solve further . For eg : the first integral gives :
$$\int_{0}^{\pi/4}\frac1{12}{((1 + 4\sec^2\theta)^{3/2} - 1)}\, d\theta$$
After this I am stuck how to proceed further, Please help.
Thank You.