Is there an easy way to find the value of e^(-x^2) using numerical methods?

  • Context: Undergrad 
  • Thread starter Thread starter C0nfused
  • Start date Start date
  • Tags Tags
    Integration
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
10 replies · 7K views
C0nfused
Messages
139
Reaction score
0
Hi everybody,
Do you have any idea how this is solved?
[tex]\int_{-\infty}^{+\infty} e^{-x^2} dx =?[/tex]
Thanks
 
Physics news on Phys.org
yes. first note the integral exists, let its value be T.

By some theorem

T^2 = [tex]\int_{-\infty}^{+\infty} e^{-x^2} dx\int_{-\infty}^{+\infty} e^{-y^2} dy[/tex]

which equals

[tex]\int_{-\infty}^{+\infty} e^{-x^2-y^2} dxdy[/tex]

now put it into polars
 
I=[tex]\int_{0}^{+\infty} e^{-x^2} dx[/tex], is solved using a trick, and integrating in polar coordinates, in which the Jacobian helps solving it. It is equal to [tex]\frac{\sqrt{\pi}}{2}[/tex]. Your integral is just 2I, I think.
 
Oggy,the "trick" applies for his integral.Yours,if u apply Fubini's theorem,will not lead to an integral over [itex]\mathbb{R}^{2}[/itex],but over the semiplane [itex]x\geq 0[/itex] and it wouldn't be the same...

To evaluate that integral (and to get the value u wrote),u need another "trick":the integrand is even,therefore

[tex]\int_{0}^{+\infty} e^{-x^{2}} \ dx =\frac{1}{2}\int_{\mathbb{R}} e^{-x^{2}} \ dx[/tex]



Daniel.
 
elaborating with what dexter had stated:

[tex]\int_{0}^{\infty} e^{-x^2} dx[/tex]

Convert to DI problem:

[tex]{(\int_{0}^{\infty} e^{-x^2} dx)}^2[/tex]

[tex]=(\int_{0}^{\infty} e^{-x^2} dx) (\int_{0}^{\infty} e^{-y^2} dy)[/tex]

[tex]=\int_{0}^{\infty}\int_{0}^{\infty} e^{-x^2}e^{-y^2} dxdy[/tex]

[tex]=\int_{0}^{\infty}\int_{0}^{\infty} e^{-{(x^2+y^2)}} dxdy[/tex]

From here using a Change of Variables step or noticing the region we are integrating is simply a semicircular arc, we can use polar coordinates:

[tex]=\int_{0}^{\frac{\pi}{2}}\int_{0}^{R} e^{-r^2}r drd\theta[/tex]

Where we take the limit of this as R --> infinity and the extra r in the integrand is from the Jacobian when switching to polar coordinates (proven in the change of variables theorem)

And NOW we can use a U-substitution to easily solve this iterated double integral:

[tex]{(\int_{0}^{\infty} e^{-x^2} dx)}^2 = \frac{\pi}{4}[/tex]

Thus:

[tex]\int_{0}^{\infty} e^{-x^2} dx = \frac{\sqrt \pi}{2}[/tex]

It seems like we can solve quite a few "definite" or "improper" integrals of non elementary functions using multivar. calculus. Does anyone know any other methods?
 
That's exactly what I referred to.
 
Well... you didn't explicitly state residue calculus now did you :wink:
 
Thanks for your help. I am not very familiar with integration theorems but i wanted to calculate this integral because the normal distibution function has a term like this
(e^(-x^2)) and the integral is 1. I checked it in Mathematica and the indefinite integral of the function g(x)=e^(-x^2) is (1/2)*[(pi)^(1/2)]*erf(x). Do u know what this erf(x) is? (a simple explanation as I am not very advanced in analysis)
 
The way i know it (it coincides with the definition swallowed by Mathematica and my Maple)

[tex]\mbox{erf} (x)=:\frac{2}{\sqrt{\pi}}\int_{0}^{x} e^{-t^{2}} \ dt[/tex]



Daniel.
 
Numerical values of erf(x) or the "error function" can be readily looked up in any book of math tables.

http://www.uni-mainz.de/FB/Geo/Geologie/Geophysik/Lithosphere/erfTable.htm
 
Last edited by a moderator: