First, is evaluating the following integral of a gaussian - which occurs very commonly in branches of physics and engineering;
I=∫∞−∞e−x2dx
Let's now have a look at what I2
looks like;
I2=(∫∞−∞e−x2dx)2=∫∞−∞e−x2dx⋅∫∞−∞e−x2dx
Now, using a dummy variable y in the second integral , we get;
I2=∫∞−∞e−x2dx⋅∫∞−∞e−y2dy
Note: We can use dummy variables since they are exactly that; a variable which is only used to integrate over and (presumably) does not appear in the evaluated integral.
We can bring the exponents together - the integrals are independent.
I2=∫∞−∞∫∞−∞e−(x2+y2)dxdy
Now we use a change of variable from cartesian coordinates to polar with;
x=rcos(θ)
y=rsin(θ)
and
y=rsin(θ)
With the determinant of the Jacobian matrix being rdθdr
and noting the change of terminals of integration.
Hence I2=∫2π0∫∞0e−r2rdrdθ
Noting that nothing in the integrand depends on theta.
Now, letting u=r2
we get dudr=2r
so dr=du2r
We get I2=2π∫u=∞u=0e−urdu2r
Finally; I2=2π2[−e−u]u=∞u=0=π[−(e−∞−e0)]=π
Hence I2=π→I=√π
Note: This isn't really a strict or rigorous derivation - as I just dealt with the Improper Integrals without taking limits and checking their convergence. This was intended just to get the general idea across.
To be continued...
It's been a long time since I did any complex maths.
ReplyDeleteYou know what they say, use it or lose it.
@Raw: That is so true :x I remember some parts, but only barely..
ReplyDeleteIf all your posts will be like this, I might get addicted to your blog.
ReplyDeleteWow, a couple of years ago this would have made perfect sense to me.
ReplyDeleteAgreed with everyone else, I've lost this part of my brain.
ReplyDelete