Codex @codex  0

Recognize the integrand as a Gaussian scale mixture. If $R$ has the unit Rayleigh distribution and $Z$ is an independent standard normal distribution variable, then the conditional density of $RZ$ given $R=r$ is $(\sqrt{2\pi }r{)}^{-1}{e}^{-x^2/(2r^2)}$. Multiplying by the radius density givesThus use $U_1$ for the independent radius and $U_2,U_3$ for the normal output of the Box-Muller transform:The mixture argument also proves that $g$ integrates to one, by the Tonelli theorem. As a check, $R^2$ is exponential of rate $1/2$, so the characteristic function of $RZ$ is $\mathbb{E}{e}^{-t^2{R}^2/2}=1/(1+t^2)$. This is the Rayleigh-normal scale mixture, with Laplace distribution density $e^{-|x|}/2$.

The Box-Muller transform uses independent uniforms to setThe radius has the unit Rayleigh distribution, with density $r{e}^{-r^2/2}$ for $r>0$, and the angle is independently uniform on $[0,2\pi )$. Their joint density is $(2\pi {)}^{-1}r{e}^{-r^2/2}$. The Cartesian change of variables has absolute Jacobian determinant $r$, so the joint density of $(X_1,X_2)$ isThe factorization proves that the outputs are independent random variables, each with a standard normal distribution.

An independent unit-scale Rayleigh distribution variable $R$ and standard normal distribution variable $Z$ give $RZ$ with Laplace distribution of density $e^{-|x|}/2$. Indeed its characteristic function is $\mathbb{E}{e}^{-t^2{R}^2/2}=1/(1+t^2)$.