Explore the Gamma function Γ(α), the most important special function in analysis, its extension of the factorial, the Gaussian integral, and applications to probability and physics.
The factorial function is only defined for non-negative integers. But what if we want to compute or ?
The Gamma function provides this extension, defined for all positive real numbers (and beyond, via analytic continuation).
The Key Relationship
For positive integers n
Famous Value
Related to Gaussian integral
The Gamma function appears throughout mathematics and science:
For , the Gamma function is defined by the improper integral:
This integral combines a Type II singularity at (when ) with a Type I infinite limit at .
The integral converges for all .
Split at t = 1:
Near t = 0:
Since converges for (p = 1-α < 1), the first part converges.
Near t = ∞:
For any , exponential decay dominates polynomial growth:
Since converges, the second part converges. ∎
Problem: Evaluate
Solution:
Problem: Evaluate
Solution:
Integration by parts:
So
The choice of variable t is conventional. Alternative forms include:
For all :
Use integration by parts on :
Let , , so , .
The boundary term vanishes: at , ; at , exponential decay dominates.
∎
For positive integers :
By induction using and :
In general: ∎
The recurrence lets us define factorial for non-integers:
Problem: Evaluate
Solution:
Using the recurrence repeatedly:
By definition:
Substitute , so :
This equals the Gaussian integral from 0 to ∞, which is :
∎
Let . Consider :
Convert to polar coordinates: , :
Let :
Therefore . ∎
Using the recurrence :
For , the Beta function is defined by:
The Beta function is related to Gamma by:
Start with the product :
Substitute and where and :
This equals . Solving: . ∎
The Beta function is symmetric:
This follows from the substitution in the integral.
Problem: Evaluate
Solution:
Problem: Evaluate
Solution:
This can also be verified directly:
For and :
Substitute , so and :
∎
Problem: Evaluate
Solution:
Apply the formula with n = 3 and a = 2:
Problem: Evaluate
Solution:
Here n = 1/2 and a = 1:
For :
Problem: Evaluate
Solution:
Here m = 4, n = 2:
Problem: Evaluate
Solution:
This is by definition:
Problem: Evaluate
Solution:
Use the formula :
Problem: Evaluate
Solution:
This is :
Problem: Evaluate
Solution:
Substitute , so and :
Problem: Evaluate
Solution:
Use the formula with m = 3, n = 2:
Problem: Evaluate
Solution:
Substitute , so and :
For all :
This can also be written as:
Problem: Verify the duplication formula for
Solution:
Left side:
Right side:
Both sides match. ✓
Problem: Use the duplication formula with z = 1/2 to find
Solution:
With z = 1/2:
Since and :
The formula checks out. Now use z = 1:
So
The duplication formula is useful for:
Problem: Find the Laplace transform of
Solution:
The Laplace transform is:
This fundamental result connects Gamma to transform methods.
Problem: Evaluate
Solution:
Substitute :
Using :
Problem: Evaluate for
Solution:
This is :
For example, with n = 2:
Problem: Estimate
Solution:
Using the recurrence from :
Continuing:
This equals
Problem: Show that the volume of an n-sphere with radius R involves Gamma
Solution:
The volume formula is:
For n = 2: (area of circle)
For n = 3:
Problem: Evaluate for a, p > 0
Solution:
Substitute , so and :
, not . The exponent in the integral is , causing the off-by-one.
Memorize: , , and use recurrence for others.
For , use . For , use .
lets you reduce to known values.
where .
When the exponent of e is not just -t, account for the scaling: .
Always check your answer with known values. For n = 0, 1, 2 you can verify directly.
Wrong: Thinking
Right: , so
Wrong:
Right: since the definition uses
Wrong:
Right:
Wrong:
Right:
Wrong: or
Right:
Wrong: Computing
Right: is undefined (has a pole). Similarly for all non-positive integers.
| α | Γ(α) | Notes |
|---|---|---|
| 1 | 1 | Base case: Γ(1) = 0! |
| 2 | 1 | Γ(2) = 1! = 1 |
| 3 | 2 | Γ(3) = 2! = 2 |
| 4 | 6 | Γ(4) = 3! = 6 |
| 5 | 24 | Γ(5) = 4! = 24 |
| 1/2 | ≈ 1.7725 | |
| 3/2 | ≈ 0.8862 | |
| 5/2 | ≈ 1.3293 | |
| 7/2 | ≈ 3.3234 |
Challenge 1
Evaluate using Gamma function:
Use Γ(n+1)/a^{n+1} with n=5, a=3.
Challenge 2
Evaluate using Beta function:
This is B(4,5) = Γ(4)Γ(5)/Γ(9).
Challenge 3
Prove using Gamma functions:
Substitute u = x⁴.
Challenge 4
Evaluate the integral:
Use the trigonometric Beta formula with m=1/2, n=0.
A random variable X has a Gamma distribution with parameters α > 0 and β > 0 if:
The Gamma function in the denominator normalizes the distribution.
A random variable X has a Beta distribution with parameters α, β > 0 if:
The Beta function normalizes this probability density on [0,1].
The chi-squared distribution with k degrees of freedom is Gamma(k/2, 2):
This is fundamental in statistical hypothesis testing.
| p | q | B(p,q) | Notes |
|---|---|---|---|
| 1 | 1 | 1 | Γ(1)²/Γ(2) = 1 |
| 2 | 2 | 1/6 | Γ(2)²/Γ(4) = 1/6 |
| 3 | 3 | 1/30 | Γ(3)²/Γ(6) = 4/120 |
| 1/2 | 1/2 | π | Γ(1/2)²/Γ(1) = π |
| 1 | 2 | 1/2 | Γ(1)Γ(2)/Γ(3) = 1/2 |
| 2 | 3 | 1/12 | Γ(2)Γ(3)/Γ(5) = 2/24 |
| 3 | 4 | 1/60 | Γ(3)Γ(4)/Γ(7) = 12/720 |
| 1/2 | 3/2 | π/2 | Used in arc length integrals |
The Gamma function was introduced by Leonhard Euler in 1729 as an extension of the factorial to non-integer arguments. He discovered the integral representation while trying to find a smooth curve passing through the points (n, n!) for positive integers n.
The notation Γ(z) was introduced by Adrien-Marie Legendre in the early 19th century. The duplication formula is also named after Legendre.
Today, the Gamma function is one of the most important special functions in mathematics, appearing in number theory, complex analysis, probability, physics, and engineering applications.
Γ(α) = ∫₀^∞ t^{α-1}e^{-t}dt for α > 0. It extends the factorial to non-integers: Γ(n+1) = n! for positive integers n.
It connects the Gamma function to the Gaussian integral ∫e^{-x²}dx = √π, which appears in probability (normal distribution), physics, and many other areas.
B(p,q) = ∫₀¹ x^{p-1}(1-x)^{q-1}dx appears in probability (Beta distribution), combinatorics, and simplifies many integrals via B(p,q) = Γ(p)Γ(q)/Γ(p+q).
Transform to the form ∫₀^∞ t^{α-1}e^{-t}dt using substitution. For example, ∫₀^∞ x^n e^{-ax}dx = Γ(n+1)/a^{n+1} = n!/a^{n+1}.
Γ(n+1) = n! for positive integers. This means Γ(1) = 0! = 1, Γ(2) = 1! = 1, Γ(3) = 2! = 2, Γ(4) = 3! = 6, etc.
Use Γ(1/2) = √π and the recurrence. For example: Γ(3/2) = (1/2)Γ(1/2) = √π/2, Γ(5/2) = (3/2)(1/2)Γ(1/2) = 3√π/4.
Near t=0, t^{α-1} is integrable for α > 0 (p-test). Near t=∞, e^{-t} dominates any polynomial growth, ensuring convergence.
Γ(2z) = (2^{2z-1}/√π) Γ(z)Γ(z+1/2). This relates Gamma at 2z to its values at z and z+1/2.
Yes, via Γ(α) = Γ(α+1)/α. This gives poles at 0, -1, -2, ... but defines Gamma for all other real (and complex) numbers.
The Gamma distribution has PDF f(x) = x^{α-1}e^{-x/β}/(β^α Γ(α)). The Chi-squared, Exponential, and Erlang distributions are special cases.