Take the 2-minute tour ×
Mathematics Stack Exchange is a question and answer site for people studying math at any level and professionals in related fields. It's 100% free, no registration required.

It seems as if no one has asked this here before, unless I don't know how to search.

The Gamma function is $$ \Gamma(\alpha)=\int_0^\infty x^{\alpha-1} e^{-x}\,dx. $$ Why is $$ \Gamma\left(\frac{1}{2}\right)=\sqrt{\pi}\ ? $$ (I'll post my own answer, but I know there are many ways to show this, so post your own!)

share|improve this question
    
Possible duplicate: math.stackexchange.com/questions/9286/… –  Argon Oct 20 '12 at 17:35
1  
@Argon Did you look at the other question? It's not even close to the same thing. And, how on earth did 2 other people vote to close without even looking at it? –  Graphth Oct 20 '12 at 23:49
4  
@Graphth With a simply change of variables ($x^2 = t$), the integrals are virtually identical, except for a constant coefficient. –  Argon Oct 20 '12 at 23:52
    
Maybe not that question, but I am quite sure this has been asked before... –  J. M. Oct 21 '12 at 6:09
1  
math.uconn.edu/~kconrad/blurbs/.../gaussianintegral.pdf is a good source of methods to solve this. –  Ishan Banerjee Mar 17 '13 at 8:07

9 Answers 9

We only need Euler's formula:

$$\Gamma(1-z)\Gamma(z)=\frac{\pi}{\sin \pi z}\Longrightarrow \Gamma^2\left(\frac{1}{2}\right)=\pi$$

share|improve this answer
    
You missed a square-root. –  user17762 Oct 17 '12 at 3:33
    
Indeed I did. Thanks. –  DonAntonio Oct 17 '12 at 3:39
    
@DonAntonio , Have you any source which prove Euler's formula which you have used ? –  Maths Lover Aug 20 '13 at 3:06
    
@MathsLover, there must be thousands as that is a rather well known formula (Euler's Reflection Formula). Check the following: en.wikipedia.org/wiki/Gamma_function , mathworld.wolfram.com/GammaFunction.html (41) , en.wikipedia.org/wiki/Reflection_formula , planetmath.org/eulerreflectionformula , aw.twi.tudelft.nl/~koekoek/documents/wi4006/gammabeta.pdf (formula 18, page 7) , etc. –  DonAntonio Aug 20 '13 at 10:46

One can use the trick with spherical coordinates: $$ \Gamma\left(\frac{1}{2}\right) = \int_0^\infty \mathrm{e}^{-x} \frac{\mathrm{d} x}{\sqrt{x}} = \int_0^\infty \mathrm{e}^{-x} \mathrm{d} (2 \sqrt{x}) = \int_{-\infty}^\infty \mathrm{e}^{-u^2} \, \mathrm{d} u $$ Then: $$ \Gamma\left(\frac{1}{2}\right) \Gamma\left(\frac{1}{2}\right) = \int_{-\infty}^\infty \int_{-\infty}^\infty \mathrm{e}^{-x^2} \mathrm{e}^{-y^2} \, \mathrm{d}x \, \mathrm{d} y = \underbrace{\int_0^\infty r \mathrm{e}^{-r^2} \, \mathrm{d} r}_{\frac{1}{2}} \cdot \underbrace{\int_{0}^{2\pi} \, \mathrm{d} \phi}_{2 \pi} = \pi $$

Alternatively one can use the second Euler's integral: $$ \frac{\Gamma\left(\frac{1}{2}\right) \Gamma\left(\frac{1}{2}\right) }{\Gamma(1)} = \int_0^1 t^{-1/2} (1-t)^{-1/2} \mathrm{d}t = \int_0^1 \mathrm{d} \left(2\arcsin\left(\sqrt{t}\right)\right) = \pi $$ Now, using $\Gamma(1) = 1$ the result follows.

Yet another method is to use the duplication identity: $$ \Gamma(2s) = \frac{2^{2s-1}}{\sqrt{\pi}} \Gamma(s) \Gamma\left(s+\frac{1}{2}\right) $$ at $ s= \frac{1}{2}$.

share|improve this answer

Method 1: $$\Gamma(1/2) = \int_0^{\infty} x^{-1/2} \exp(-x) dx$$ Set $x = t^2$, to get $$\Gamma(1/2) = \int_0^{\infty} 2\exp(-t^2) dt = \sqrt{\pi}$$

Method 2: $$\beta(x,y) = \dfrac{\Gamma(x) \Gamma(y)}{\Gamma(x+y)}$$ Take $y=x=1/2$, to get $$\beta(1/2,1/2) = \dfrac{\Gamma^2(1/2)}{\Gamma(1)} \implies \Gamma(1/2) = \sqrt{\beta(1/2,1/2)}$$

$$\beta(x,y) = 2 \int_0^{\pi/2} \sin^{2x-1}(\theta) \cos^{2y-1}(\theta) d \theta$$ Hence, $$\beta(1/2,1/2) = 2 \int_0^{\pi/2} d \theta = \pi$$ $$\Gamma(1/2) = \sqrt{\pi}$$

share|improve this answer

This is a "proof". We know that the surface area of the $n-1$ dimensional unit sphere is $$ |S^{n-1}| = \frac{2\pi^{\frac{n}2}}{\Gamma(\frac{n}2)}. $$ On the other hand, we know that $|S^2|=4\pi$, which gives $$ 4\pi = \frac{2\pi^{\frac32}}{\Gamma(\frac32)} = \frac{2\pi^{\frac32}}{\frac12\Gamma(\frac12)}. $$

share|improve this answer
2  
This formula also holds for $n = 1$ if you take counting measure as your $0$-dimensional measure. –  JavaMan Dec 23 '12 at 18:09

It also follows from the famous Riemann functional equation with $s=1/2$: $$ \zeta(s) = 2^s\pi^{s-1}\ \sin\left(\frac{\pi s}{2}\right)\ \Gamma(1-s)\ \zeta(1-s) $$ but Euler's reflection formula is probably more basic.

share|improve this answer

If there's any justice in the universe, someone must have asked here how to show that $$ \int_{-\infty}^\infty e^{-x^2/2}\,dx = \sqrt{2\pi}. $$ Let's suppose that has been answered here. Let (capital) $X$ be a random variable whose probability distribution is $$ \frac{e^{-x^2/2}}{\sqrt{2\pi}}\,dx. $$ Consider the problem of finding $\mathbb E(X^2)$. It is $$ \frac{1}{\sqrt{2\pi}}\int_{-\infty}^\infty x^2 e^{-x^2/2}\,dx = \text{(by symmetry)} \frac{2}{\sqrt{2\pi}} \int_0^\infty x^2 e^{-x^2/2}\,dx $$ $$ \sqrt{\frac2\pi}\int_0^\infty xe^{-x^2/2}\Big(x\,dx\Big) = \sqrt{\frac2\pi}\int_0^\infty \sqrt{u}\ e^{-u}\,du = \sqrt{\frac2\pi}\ \Gamma\left(\frac32\right) = \frac12\sqrt{\frac2\pi} \Gamma\left(\frac12\right). $$ So it is enough to show that this expected value is $1$. That is true if the sum of two independent copies of it has expected value $2$. So: $$ \Pr\left(X^2+Y^2<w\right) = \frac{1}{2\pi}\iint\limits_\mathrm{disk} e^{-(x^2+y^2)/2}\,dx\,dy $$ where the disk has radius $\sqrt{w}$. This equals $$ \frac{1}{2\pi}\int_0^{2\pi}\int_0^{\sqrt{w}} e^{-\rho^2/2} \,\rho\,d\rho\,d\theta = \int_0^{\sqrt{w}} e^{-\rho^2/2} \,\rho\,d\rho. $$ This last equality holds because we are integrating with respect to $\theta$ something not depending on $\theta$. Differentiating this with respect to $w$ gives the probability density function of the random variable $X^2+Y^2$: $$ e^{-w/2}\sqrt{w}\frac{1}{2\sqrt{w}} = \frac{e^{-w/2}}{2}\text{ for }w>0. $$ So $$ \mathbb E(X^2+Y^2) = \int_0^\infty w \frac{e^{-w/2}}{2}\,dw =2. $$

share|improve this answer
3  
If $\int_{-\infty}^\infty e^{-u^2/2}\,du = \sqrt{2\pi}$ is given, then it seems much simpler to take $\Gamma(1/2) = \int_0^{\infty} x^{-1/2} \exp(-x) dx$, substitute $u^2=2x$, and we are done. –  leonbloy Jun 17 '13 at 1:13

$$B\left(\frac{1}{2},\frac{1}{2}\right)=\frac{\left[\Gamma\left(\frac{1}{2}\right)\right]^{2}}{\Gamma{(1)}}=\left[\Gamma\left(\frac{1}{2}\right)\right]^{2}$$ $$B\left(\frac{1}{2},\frac{1}{2}\right)=\frac{\pi}{\sin{\frac{\pi}{2}}}=\pi$$

Chris.

share|improve this answer
1  
For those who find the second line not obvious, here is an elementary explanation: by definition, $B(\frac{3}{2},\frac{3}{2}) = \intop_0^1 x^{1/2} (1-x)^{1/2} dx$, and the graph of $x^{1/2} (1-x)^{1/2}$ is a semicircle, so this directly relates to $\pi$ being defined via the area of the circle. This allows to calculate $\Gamma(\frac{3}{2})$, which is just $\frac{1}{2} \Gamma(\frac{1}{2})$ –  Alexander Shamov Jun 16 '13 at 22:27

I'm surprised that no one has mentioned that this $\sqrt{\pi}$ is also "the same one" as the $\sqrt{\pi}$ in Stirling's formula.

It is a known fact that $\Gamma$ is uniquely characterized as the function that satisfies the equation $\Gamma(z+1) = z \Gamma(z)$ and whose logarithm "has nice asymptotics at $+\infty$". The conventional way to explain what "nice asymptotics" should mean in this context involves logarithmic convexity, but I prefer to leave it like this: any function that satisfies that functional equation must coincide with $\Gamma(z)$ up to a periodic factor, which, if nontrivial, would prevent it from having an asymptotic expansion at $+\infty$ in terms of powers and logs.

Now let's turn to Stirling's formula and assume that it holds both for integer and half-integer values of $\Gamma$. Clearly,

$$\Gamma(n+1/2) = \Gamma(1/2) \cdot \frac{1}{2} \cdot \frac{3}{2} \cdot \dots \cdot \frac{2n-1}{2} = \Gamma(1/2) \cdot \frac{(2n)!}{2^{2n} n!},$$

which is basically the duplication formula. Now if we plug it into the Stirling's formula, we will find out that $\Gamma(1/2) = \sqrt{\pi}$.

share|improve this answer

There are a few details missing (showing the integral converges, justifying the differentiation under the integral sign, and so forth), but the main idea should be clear.

Consider the Fourier transform of $f(x)=e^{-x^2}$.

$$\hat f(\xi)= \int_{-\infty}^{\infty} e^{-x^2} e^{-2\pi i x\xi} \ dx =\int_{-\infty}^{\infty} e^{-x^2} (\cos(2\pi x \xi) - i \sin(2\pi x \xi)) dx =\int_{-\infty}^{\infty} e^{-x^2} \cos(2\pi x \xi) \ dx$$

Differentiate under the integral sign and integrate by parts.

$$\hat f'(\xi) =-2\pi \int e^{-x^2}\sin(2\pi \xi x) \ dx = -2\pi^2 \xi \int e^{-x^2} \cos(2\pi x \xi) \ dx = -2\pi^2\xi \hat f(\xi)$$

Separate variables.

$$\frac{\hat f'(\xi)}{\hat f(\xi) }= -2\pi^2 \xi \Rightarrow \hat f(\xi) =Ce^{-\pi^2 \xi^2}.$$

Use the inversion formula and change variables .

$$e^{-x^2} = C\int e^{-\pi^2 \xi^2} \cos(2\pi x \xi) \ d\xi = \frac{C}{\pi} \int e^{- \xi^2} \cos(2 x \xi) \ d\xi $$

$$e^{-\pi^2 x^2}\frac{C}{\pi} =\int e^{- \xi^2} \cos(2 x\pi \xi) \ d\xi =\frac{C^2}{\pi}e^{-\pi^2x^2}\Rightarrow C=\sqrt{\pi}$$

$$\hat f (0) = \sqrt{\pi} = \int e^{-x^2} \cos(0) \ dx = \int e^{-x^2} \ dx.$$

share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.