There are several ways of proving $I=\int_{-\infty}^{\infty} \mathrm{e}^{-x^2}\, dx = \sqrt \pi$, of which one of the ways is to consider $$I^2=\bigg(\int_{-\infty}^{+\infty} e^{-x^2} dx\bigg)\bigg(\int_{-\infty}^{+\infty} e^{-y^2} dy\bigg)=\int_{-\infty}^{+\infty}\int_{-\infty}^{+\infty}e^{-(x^2+y^2)} dxdy$$ and then convert to polar coordinates.
The step changing the product of two single variable integrals into a double integral can be justified using Fubini's theorem, the proof of which is beyond the standard undergraduate curriculum of many universities. Is there a more elementary way to justify this step?
Clarification: This is not a duplicate of the question requiring a proof of the integral. I am asking for justification of a step used in the accepted answer without using Fubini's theorem which is not something addressed there.