Consider the following integral: $$\mathcal{I}=\int_1^\infty\int_0^1\frac{\mathrm dy\,\mathrm dx}{\sqrt{x^2-1}\sqrt{1-y^2}\sqrt{1-y^2+4\,x^2y^2}}.$$ It can be represented as $$\mathcal{I}=\int_1^\infty\frac{K(\sqrt{1-4\,x^2})}{\sqrt{x^2-1}}\mathrm dx=\Re\int_0^1\frac{K\left(\sqrt\frac{1+3\,y^2}{1-y^2}\right)}{1-y^2}\mathrm dy,$$ where $K(x)$ is the complete elliptic integral of the 1st kind.
I was not able to further simplify any of these integrals, but numerical integration suggests the following conjectural closed form: $$\mathcal{I}\stackrel{?}{=}\frac{3\,\Gamma(1/3)^6}{2^{17/3}\pi^2}.$$ Can you suggest a proof of this conjecture?
Using Mathematica I obtained a closed form for $\mathcal{I}$ in terms of the Meijer G-function: $$\mathcal{I}=\frac{\pi^{3/2}}{4}G_{3,3}^{2,1}\left(\frac{1}{4}\left|\begin{array}{c}1,\ \ 1,\ \ 1\\\frac{1}{2},\frac{1}{2},\frac{1}{2}\\\end{array}\right.\right),$$ but I have no manual proof for that and do not know how to reduce it to $\frac{3\,\Gamma(1/3)^6}{2^{17/3}\pi^2}.$