💾 Archived View for jsreed5.org › math › 20211118-integrating-the-standard-normal-distribution › ind… captured on 2021-11-30 at 19:37:34. Gemini links have been rewritten to link to archived content
-=-=-=-=-=-=-
Download a PostScript version of this document
(PS, 230 KB)
---
One of the most prevalent probability distributions is the standard normal distribution, given by the function
f(x)=(1/sqrt(2π))e^(-x^2/2)
Like all probability distributions, the area under this function is 1, but the integral of the base function e^(-x^2) cannot be expressed in terms of elementary functions. Instead, these kinds of integrals are expressed in terms of special functions called "error functions". How, then, can we show that the area under the curve is unit?
I found the following solution online many years ago while I was still a student, but I didn't bother to go through the logic of each step. This note goes through the justification for the solution method and what can be learned from it.
We wish to determine the value of the definite integral
∫((1/sqrt(2π))e^(-x^2/2),x,-infinity,infinity)
1/sqrt(2π) is a constant factor, so it can be separated from the integrand.
Also, because x^2 is an even function, e^(-x^2) is an even function. We can thus consider the integral only from 0 to infinity.
∫((1/sqrt(2π))e^(-x^2/2),x,-infinity,infinity)=(2/sqrt(2π))∫(e^(-x^2/2),x,0,infinity)
To aid in evaluating the integral, we will introduce another term into the integrand. We do this by mutiplying the integral by itself and taking its square root.
(2/sqrt(2π))∫(e^(-x^2/2),x,0,infinity)=(2/sqrt(2π))sqrt(∫(e^(-x^2/2),x,0,infinity)∫(e^(-x^2/2),x,0,infinity))
Definite integrals are not affected by choice of variable, so we substitute x with y in the second integral. This allows us to join the two integrals into one double integral.
(2/sqrt(2π))sqrt(∫(e^(-x^2/2),x,0,infinity)∫(e^(-x^2/2),x,0,infinity))=(2/sqrt(2π))sqrt(∫(e^(-x^2/2),x,0,infinity)∫(e^(-y^2/2),y,0,infinity)) =(2/sqrt(2π))sqrt(∫(∫(e^(-(x^2+y^2)/2),x,0,infinity),y,0,infinity))
We now perform another change of variable:
y=sx, dy=xds
Substituting:
(2/sqrt(2π))sqrt(∫(∫(e^(-(x^2+y^2)/2),x,0,infinity),y,0,infinity)=(2/sqrt(2π))sqrt(∫(∫(e^(-(x^2+(s^2)(x^2))/2),x,0,infinity),s,0,infinity)) =(2/sqrt(2π))sqrt(∫(∫(xe^((-x^2/2)(1+s^2)),x,0,infinity),s,0,infinity))
Having x in front of the exponent allows us to integrate with respect to x and produce elementary functions.
(2/sqrt(2π))sqrt(∫(∫(xe^((-x^2/2)(1+s^2)),x,0,infinity),s,0,infinity))=(2/sqrt(2π))sqrt(∫([(-1/(1+s^2))e^((-x^2/2)(1+s^2)),x,0,infinity],s,0,infinity)) =(2/sqrt(2π))sqrt(∫((-1/(1+s^2))(e^(-infinity)-e^0),s,0,infinity)) =(2/sqrt(2π))sqrt(∫(-1/(1+s^2),s,0,infinity))
The remaining integrand is the derivative of the arctangent function. After evaluating the integral, the final arithmetic is simple.
(2/sqrt(2π))sqrt(∫(-1/(1+s^2),s,0,infinity))=(2/sqrt(2π))sqrt([arctan(s),s,0,infinity]) =(2/sqrt(2π))sqrt(arctan(infinity)-arctan(0)) =(2/sqrt(2π))sqrt(π/2) =2sqrt(π)/(sqrt(2)sqrt(2)sqrt(π)) =1
The core of this solution is to create a double integral and apply a clever change of variable to modify the integrand. Often the difference between an elementary integral and an nonelementary integral is no more than a leading factor--in this case the presence of x in front of the term e^(-x^2). Introducing a second integral and converting to a double integral allowed us to put in the factor we needed.
I tend to struggle with integration problems. I seem to lack the creativity needed to find an extra term or change of variable needed to simplify the problem. It might be that I simply don't remember my basic integrals well enough to realize what might be useful. At any rate, solutions like this fascinate me.
---
[Last updated: 2021-11-19]