The question asks me to prove that $$\int_0^\infty J(x)e^{-ax}dx=\frac{1}{\sqrt{1+a^2}},$$ where $a>0$ and $J(x)=\frac{2}{\pi}\int_0^{\pi/2}\cos(x\cos{\theta})d\theta.$
I started off by considering the function $$f(x,\theta)=\frac{2}{\pi}e^{-ax}\cos(x\cos{\theta})\boldsymbol{1}_{\theta\in(0,\frac{\pi}{2})}$$ in the region $(x,\theta)\in(0,\infty)\times (0,\infty).$ Now $\vert f(x,\theta)\vert\leq \frac{2}{\pi}e^{-ax},$ so $$\int_0^\infty\left(\int_0^\infty\vert f(x,\theta)\vert dx\right)d\theta<\infty,$$ which means that by Fubini's we can do the integration in any order.
I don't think $\cos(x\cos{\theta})$ is integrable in terms of $\theta,$ so I'm thinking I should integrate $e^{-ax}\cos(x\cos{\theta})$ in terms of $x.$ But the integration of this looks quite horrible, so I was wondering (a) if my method is correct, and (b) is there a neater, more efficient way I could solve this problem?