La función es medible porque es continua y, también, positiva. Por el Fubini-Tonelli Teorema,
$$
\begin{aligned}\int_{\mathbb{R}^n}\exp\left({-\sum_{i=1}^na_i^2x_i^2}\right)\,dx&=\prod_{i=1}^n\int_{\mathbb{R}}e^{-a_i^2x_i^2}\,dx_i=\prod_{i=1}^n\left(\int_{\mathbb{R}}e^{-a_i^2x^2}\,dx\right)\\&=\left(\frac{1}{|a_1\cdots a_n|}\int_{\mathbb{R}}e^{-x^2}\,dx\right)^n=\frac{(\sqrt{\pi})^n}{|a_1\cdots a_n|}.
\end{aligned}$$
En caso de que usted no está familiarizado con la Integral de Gauss, aquí es un resultado rápido.
Aplicando el teorema anterior y un cambio de variable a $(\int_{\mathbb{R}}e^{-x^2}\,dx)^2\,$ obtenemos
$$\left(\int_{\mathbb{R}}e^{-x^2}\,dx\right)^2=\int_{\mathbb{R}}e^{-x^2}\,dx\int_{\mathbb{R}}e^{-y^2}\,dy=\int_{\mathbb{R}^2}e^{-(x^2+y^2)}\,dxdy=\int_{]0,+\infty[\times]0,2\pi[}re^{-r^2}\,drd\theta=\pi.$$
Tomar la raíz cuadrada de ambos lados, y el resultado de la siguiente manera.