avec I = int( exp(-x^2/2) , x=0..+oo)
on a : int( exp(-x^2/2) , x=-oo .. +oo) = 2*I par parité
Ensuite I^2 = int( exp(-x^2/2) , x=0..+oo) * int( exp(-y^2/2) , x=0..+oo) = int( int (exp(-(x^2+y^2)/2) , x=0..+oo), y=-oo..+oo) par Fubini puisque les fonctions sont Riemann intégrables
Du coup, tu peux faire le changement polaire des physiciens x^2+y^2 -> r^2 et dx*dy -> r dr*d(theta) (justification à partir du jacobien classique)
donc ça donne I^2=int ( int ( r*exp(-r^2/2) , r=0..+oo), theta = 0.. 2*Pi)=-2*Pi * [exp(-r^2/2)][0..+oo] = 2*Pi
donc I = sqrt(2*Pi) (tiens bizarre de tête j'aurai dit sqrt(Pi), j'ai pas la force de chercher l'erreur tant pis

)
alors l'intégrale que tu cherches vaut 2*I=sqrt(8*Pi)
(je sais que ça sert à rien mais c'est juste pour refaire des maths, ça rappelle des bons souvenirs

)