Beta and Gamma functions

Some definite integrals refuse every standard trick — reduction formulae, trigonometric identities, integration by parts — yet turn out to be old friends in disguise. The Beta and Gamma functions are the two master integrals that unmask them. The Gamma function, Γ(n)=∫0∞e−xxn−1 dx\Gamma(n) = \int_0^\infty e^{-x}x^{n-1}\,dx, stretches the familiar factorial n!n! off the integers and onto every positive real number, so "(4.5)!(4.5)!" suddenly makes sense. The Beta function, B(m,n)=∫01xm−1(1−x)n−1 dxB(m,n) = \int_0^1 x^{m-1}(1-x)^{n-1}\,dx, is its close cousin, and a single identity converts one into the other. Once you know these two integrals and the identity linking them, whole families of sin⁡pθcos⁡qθ\sin^p\theta\cos^q\theta integrals and algebraic-power integrals collapse into a couple of lines of Gamma arithmetic — no reduction formula required.

The core insight

The Gamma function Γ(n)=∫0∞e−xxn−1dx\Gamma(n)=\int_0^\infty e^{-x}x^{n-1}dx generalizes n!n! to all positive reals, the Beta function B(m,n)=∫01xm−1(1−x)n−1dxB(m,n)=\int_0^1 x^{m-1}(1-x)^{n-1}dx is symmetric and linked to it by B(m,n)=Γ(m)Γ(n)/Γ(m+n)B(m,n)=\Gamma(m)\Gamma(n)/\Gamma(m+n), and this single relation converts hard trigonometric- and algebraic-power integrals into fast Gamma arithmetic.

The core idea

The two defining integrals. The Beta function is

B(m,n)=∫01xm−1(1−x)n−1 dx,m,n>0,B(m,n) = \int_0^1 x^{m-1}(1-x)^{n-1}\,dx, \qquad m,n>0,

where: m,nm,n are the two positive parameters and the domain of integration is fixed at [0,1][0,1]. Swapping x→1−xx\to 1-x inside the integral shows it is symmetric, B(m,n)=B(n,m)B(m,n)=B(n,m). The Gamma function is

Γ(n)=∫0∞e−xxn−1 dx,n>0,\Gamma(n) = \int_0^\infty e^{-x}x^{n-1}\,dx, \qquad n>0,

where: nn is the single positive parameter and e−xe^{-x} is the ordinary exponential decay. Both are improper integrals — Γ(n)\Gamma(n) has an infinite upper limit, and for 0<n<10<n<1 also a singularity of xn−1x^{n-1} at x=0x=0 — and the restriction m,n>0m,n>0 is exactly what the earlier convergence tests in this unit certify as safe: the exponential tail always dies fast enough at infinity, and the power singularity at 00 stays integrable only while the exponent on xx exceeds −1-1.

The bridge between them. The relation that does almost all the work in this topic is

B(m,n)=Γ(m) Γ(n)Γ(m+n),B(m,n) = \frac{\Gamma(m)\,\Gamma(n)}{\Gamma(m+n)},

where: every symbol is the same m,nm,n as above. It means a Beta value never needs to be attacked as a defining integral — convert it to three Gamma values instead.

Recurrence and the factorial link. Integrating Γ(n)\Gamma(n) by parts gives

Γ(n+1)=n Γ(n),\Gamma(n+1) = n\,\Gamma(n),

where: n>0n>0 is any positive real. Since Γ(1)=∫0∞e−xdx=1\Gamma(1)=\int_0^\infty e^{-x}dx=1, repeated use of the recurrence for a positive integer nn collapses to

Γ(n+1)=n!,\Gamma(n+1) = n!,

which is the precise sense in which Gamma "is" the factorial — shifted by one and defined for every positive real in between, not only the integers. The special half-integer seed value is

Γ ⁣(12)=π.\Gamma\!\left(\tfrac12\right) = \sqrt{\pi}.

Trigonometric form. Substituting x=sin⁡2θx=\sin^2\theta (so 1−x=cos⁡2θ1-x=\cos^2\theta and dx=2sin⁡θcos⁡θ dθdx = 2\sin\theta\cos\theta\,d\theta) into the defining Beta integral turns it into

B(m,n)=2∫0π/2sin⁡2m−1θ cos⁡2n−1θ dθ,B(m,n) = 2\int_0^{\pi/2} \sin^{2m-1}\theta\,\cos^{2n-1}\theta\,d\theta,

where: m,n>0m,n>0 fix the odd exponents 2m−12m-1 and 2n−12n-1 on sine and cosine. This single substitution is why a whole class of trigonometric-power integrals — otherwise handled only by tedious reduction formulae — becomes a two-line Gamma computation.

Why the domain restrictions matter. None of the four boxed results above are arbitrary rules to memorize — they all trace back to the convergence tests for improper integrals covered earlier in this unit. In Γ(n)=∫0∞e−xxn−1dx\Gamma(n)=\int_0^\infty e^{-x}x^{n-1}dx, the exponential factor e−xe^{-x} guarantees convergence at the infinite upper limit for every real nn, no matter how large; the only thing that can go wrong is at the lower limit x=0x=0, where xn−1x^{n-1} blows up once n<1n<1, and the standard pp-test for integrals near a singular point shows that stays integrable only while n−1>−1n-1>-1, i.e. n>0n>0. The same logic bounds B(m,n)B(m,n): both endpoints x=0x=0 and x=1x=1 carry potential singularities, of orders xm−1x^{m-1} and (1−x)n−1(1-x)^{n-1} respectively, and each needs its exponent to exceed −1-1 — exactly m>0m>0 and n>0n>0. Recognizing a given integral as "really" a Beta or Gamma integral in disguise, then checking its parameters satisfy these domain conditions, is itself an exam skill worth practicing.

Figure
A smooth curve through the points one one, two one, three two, four six and five twenty four showing the Gamma function interpolating the factorial values
Gamma(n+1) interpolates n! between the integers n Gamma(n+1) 1 2 3 4 5 1! = 1 2! = 1 3! = 2 4! = 6 5! = 24 The smooth blue curve is Gamma of n plus one; the red dots are the ordinary factorials it passes through.
Derivation AThe trigonometric form of the Beta function

The goal is to show, starting from nothing but the defining integral B(m,n)=∫01xm−1(1−x)n−1dxB(m,n)=\int_0^1 x^{m-1}(1-x)^{n-1}dx and the single substitution x=sin⁡2θx=\sin^2\theta, that B(m,n)B(m,n) equals 2∫0π/2sin⁡2m−1θcos⁡2n−1θ dθ2\int_0^{\pi/2}\sin^{2m-1}\theta\cos^{2n-1}\theta\,d\theta — the identity that turns every sin⁡pθcos⁡qθ\sin^p\theta\cos^q\theta integral on [0,π/2][0,\pi/2] into a Beta-function evaluation.

Step 1. Choose the substitution x=sin⁡2θx=\sin^2\theta. As xx runs over [0,1][0,1], θ\theta runs over [0,π/2][0,\pi/2], so the limits transform cleanly with no extra pieces to patch in.

Step 2. Rewrite every factor in the integrand in terms of θ\theta. Directly, xm−1=sin⁡2m−2θx^{m-1}=\sin^{2m-2}\theta. For the second factor, 1−x=1−sin⁡2θ=cos⁡2θ1-x = 1-\sin^2\theta=\cos^2\theta, so (1−x)n−1=cos⁡2n−2θ(1-x)^{n-1}=\cos^{2n-2}\theta. Differentiating x=sin⁡2θx=\sin^2\theta gives dx=2sin⁡θcos⁡θ dθdx = 2\sin\theta\cos\theta\,d\theta.

Step 3. Substitute all three pieces into the integral. B(m,n) = \int_0^{\pi/2} \sin^{2m-2}\theta\,\cos^{2n-2}\theta\,\cdot 2\sin\theta\cos\theta\,d\theta. Combine the powers of sine (2m−22m-2 and 11) and of cosine (2n−22n-2 and 11): the sine exponent becomes 2m−12m-1 and the cosine exponent becomes 2n−12n-1.

Step 4. Read off the result. B(m,n) = 2\int_0^{\pi/2}\sin^{2m-1}\theta\,\cos^{2n-1}\theta\,d\theta, exactly the trigonometric form quoted above. Because the derivation only used a substitution and algebra on exponents, it holds for every m,n>0m,n>0, which is why matching a given sin⁡pθcos⁡qθ\sin^p\theta\cos^q\theta integral to this form is simply a matter of solving 2m−1=p2m-1=p and 2n−1=q2n-1=q for mm and nn.

Worked example

Problem. Evaluate ∫0π/2sin⁡4θ cos⁡2θ dθ\displaystyle\int_0^{\pi/2}\sin^4\theta\,\cos^2\theta\,d\theta using the Beta function, then use the result to evaluate the related integral ∫01x3(1−x)4 dx\displaystyle\int_0^1 x^3(1-x)^4\,dx.

Step 1. Match exponents. Comparing sin⁡4θcos⁡2θ\sin^4\theta\cos^2\theta against sin⁡2m−1θcos⁡2n−1θ\sin^{2m-1}\theta\cos^{2n-1}\theta gives 2m−1=4⇒m=522m-1=4\Rightarrow m=\tfrac52 and 2n−1=2⇒n=322n-1=2\Rightarrow n=\tfrac32, so the integral equals 12B(52,32)\tfrac12 B(\tfrac52,\tfrac32).

Step 2. Convert to Gamma values. B(52,32)=Γ(5/2)Γ(3/2)Γ(4)B(\tfrac52,\tfrac32)=\dfrac{\Gamma(5/2)\Gamma(3/2)}{\Gamma(4)}. From the recurrence starting at Γ(1/2)=π\Gamma(1/2)=\sqrt\pi: Γ(3/2)=12π\Gamma(3/2)=\tfrac12\sqrt\pi and Γ(5/2)=32⋅12π=34π\Gamma(5/2)=\tfrac32\cdot\tfrac12\sqrt\pi=\tfrac34\sqrt\pi; and Γ(4)=3!=6\Gamma(4)=3!=6.

Step 3. Substitute and simplify. \int_0^{\pi/2}\sin^4\theta\cos^2\theta\,d\theta = \frac12\cdot\frac{\left(\tfrac34\sqrt\pi\right)\left(\tfrac12\sqrt\pi\right)}{6} = \frac12\cdot\frac{\tfrac38\pi}{6} = \frac{3\pi}{96} = \frac{\pi}{32}. So ∫0π/2sin⁡4θcos⁡2θ dθ=π32\displaystyle\int_0^{\pi/2}\sin^4\theta\cos^2\theta\,d\theta = \dfrac{\pi}{32}.

Competitive stretch

Step 1. The algebraic-power integral ∫01x3(1−x)4 dx\int_0^1 x^3(1-x)^4\,dx matches the Beta definition ∫01xm−1(1−x)n−1dx\int_0^1 x^{m-1}(1-x)^{n-1}dx directly with m−1=3⇒m=4m-1=3\Rightarrow m=4 and n−1=4⇒n=5n-1=4\Rightarrow n=5, so the integral is B(4,5)B(4,5).

Step 2. Convert to Gamma values: B(4,5)=Γ(4)Γ(5)Γ(9)=3! 4!8!=6×2440320=14440320B(4,5)=\dfrac{\Gamma(4)\Gamma(5)}{\Gamma(9)}=\dfrac{3!\,4!}{8!}=\dfrac{6\times24}{40320}=\dfrac{144}{40320}.

Step 3. Simplify the fraction: 14440320=1280\dfrac{144}{40320}=\dfrac{1}{280}. So ∫01x3(1−x)4 dx=1280\displaystyle\int_0^1 x^3(1-x)^4\,dx = \dfrac{1}{280} — an integer-power Beta integral that would otherwise need four rounds of integration by parts, done here with two factorial lookups.

Common traps

  • Writing Γ(n)=(n−1)!\Gamma(n)=(n-1)! but then substituting nn where n+1n+1 was meant — keep the shift explicit: Γ(n+1)=n!\Gamma(n+1)=n!, so Γ(n)=(n−1)!\Gamma(n)=(n-1)! only for positive integers nn.
  • Dropping the factor of 22 in the trigonometric form B(m,n)=2∫0π/2sin⁡2m−1θcos⁡2n−1θ dθB(m,n)=2\int_0^{\pi/2}\sin^{2m-1}\theta\cos^{2n-1}\theta\,d\theta — the single most common slip, halving an otherwise correct answer.
  • Misreading exponents when converting sin⁡pθcos⁡qθ\sin^p\theta\cos^q\theta to Beta form: the correspondence is m=p+12m=\tfrac{p+1}{2}, n=q+12n=\tfrac{q+1}{2}, not m=pm=p, n=qn=q.
  • Attacking a Beta-function value by brute-force integration of its defining integral instead of reaching for B(m,n)=Γ(m)Γ(n)/Γ(m+n)B(m,n)=\Gamma(m)\Gamma(n)/\Gamma(m+n), which is almost always faster once Γ\Gamma at integers and half-integers is known.
  • Forgetting Γ(1)=1\Gamma(1)=1 as the base case, or omitting −1-1 when converting an algebraic-power integral ∫01xa(1−x)b dx\int_0^1 x^a(1-x)^b\,dx to B(a+1,b+1)B(a+1,b+1).
  • Assuming the Beta or Gamma definitions hold for all real parameters — both integrals require strictly positive arguments to converge; a negative or zero parameter is outside the domain covered here, not a special case to compute.

Ask the tutor: if matching a sin⁡pθcos⁡qθ\sin^p\theta\cos^q\theta or xa(1−x)bx^a(1-x)^b integral to Beta form still feels shaky, ask the AI tutor for a few more conversions to practice the exponent bookkeeping.

Why does Gamma of n plus one equal n factorial only for integers, when the integral itself is defined for every positive real n?

The integral Γ(n)=∫0∞e−xxn−1dx\Gamma(n)=\int_0^\infty e^{-x}x^{n-1}dx is perfectly well-defined for any real n>0n>0, and the recurrence Γ(n+1)=nΓ(n)\Gamma(n+1)=n\Gamma(n) holds for all of them too. What is special about integers is only the base case: unwinding the recurrence down to Γ(1)=1\Gamma(1)=1 lands you exactly on the ordinary factorial product n(n−1)(n−2)⋯1n(n-1)(n-2)\cdots1. For a non-integer like n=2.5n=2.5, unwinding lands on Γ(1/2)=π\Gamma(1/2)=\sqrt\pi instead, so you get a factorial-like number involving π\sqrt\pi, not a whole-number factorial.

How do I know whether to use the trigonometric form or the direct definition of the Beta function?

Look at what the integral in front of you looks like. If it is over [0,π/2][0,\pi/2] with powers of sine and cosine, match it to 2∫0π/2sin⁡2m−1θcos⁡2n−1θ dθ2\int_0^{\pi/2}\sin^{2m-1}\theta\cos^{2n-1}\theta\,d\theta and read off m,nm,n. If it is over [0,1][0,1] with powers of xx and (1−x)(1-x), match it directly to ∫01xm−1(1−x)n−1dx\int_0^1 x^{m-1}(1-x)^{n-1}dx. Either way, the very next step is the same: convert to Γ(m)Γ(n)/Γ(m+n)\Gamma(m)\Gamma(n)/\Gamma(m+n) rather than integrating further.

Is there a quick way to sanity-check a Beta or Gamma answer before trusting it?

Two quick checks catch most slips. First, symmetry: B(m,n)B(m,n) must equal B(n,m)B(n,m), so if swapping your mm and nn changes the numeric answer, you made an arithmetic error. Second, magnitude: since 0≤sin⁡θ,cos⁡θ≤10\le \sin\theta,\cos\theta\le 1 on [0,π/2][0,\pi/2], any trigonometric-power integral there must come out as a positive number well under π/2≈1.57\pi/2\approx1.57 — if your Beta-function route gives something larger or negative, recheck the Gamma values you substituted.

Before this

Where this leads