Notes on Complex Analysis
Export PDFGitHub

§12.1 The Γ-Function

12.1 The Γ-Function

Definition 12.1.1.

The Gamma function is defined by

Γ(𝑧)=∫0∞e−𝑡𝑡𝑧−1d𝑡,(12.1.1)
where 𝑧∈ℂ.

By letting 𝑧=𝑥+i𝑦 where 𝑥,𝑦∈ℝ, we have |e−𝑡𝑡𝑧−1|=e−𝑡𝑡𝑥−1. Notice that for 𝑥>0,

|Γ(𝑥)|=∫01e−𝑡𝑡𝑥−1d𝑡+∫1∞e−𝑡𝑡𝑥−1d𝑡≤∫01𝑡𝑥−1d𝑡+∫1∞e−𝑡𝑡𝑥−1d𝑡=1𝑥+∫1∞e−𝑡𝑡𝑥−1d𝑡.

Since ∫1∞d𝑡𝑡2 is convergent and lim𝑡→∞e−𝑡𝑡𝑧−1𝑡−2=0, then by comparison, the second integral is convergent.

Therefore, Γ(𝑥) is convergent on ℝ>0. It follows that Γ(𝑧) is absolutely convergent on the right half-plane {𝑧∈ℂ:ℜ𝔢(𝑧)>0}.

Theorem 12.1.2.

The Γ-function is holomorphic on {𝑧∈ℂ:ℜ𝔢(𝑧)>0}.
Proof.

Let 𝛾⊂{𝑧∈ℂ:ℜ𝔢(𝑧)>0} be an arbitrary simple, closed, rectifiable curve. By Morera’s Theorem (Theorem 3.2.6Theorem 3.2.6Theorem 3.2.6), it suffices to show that ∮𝛾Γ(𝑧)d𝑧=0. We have

∮𝛾Γ(𝑧)d𝑧=∮𝛾(∫0∞e−𝑡𝑡𝑧−1d𝑡)d𝑧=∫0∞e−𝑡(∮𝛾𝑡𝑧−1d𝑧)d𝑡=0,
where the integral exchange is justified by uniform convergence of Γ on compact subsets of {𝑧∈ℂ:ℜ𝔢(𝑧)>0} (by the Weierstrass 𝑀–Test, Theorem 2.3.3Theorem 2.3.3Theorem 2.3.3). ⁠□

From integration by parts, we obtain

Γ(𝑧+1)=∫0∞e−𝑡𝑡𝑧d𝑡=−e−𝑡𝑡𝑧|0∞+𝑧∫0∞e−𝑡𝑡𝑧−1d𝑡=𝑧Γ(𝑧).

Additionally,

Γ(1)=∫0∞e−𝑡d𝑡=−e−𝑡|0∞=1.

Hence, we have Γ(𝑧+1)=𝑧! for nonnegative integers, and the Γ-function generalizes the factorial. We also have

Γ(𝑧+𝑛)=Γ(𝑧)∏𝑘=0𝑛−1(𝑧+𝑘),ℜ𝔢(𝑧)>0,𝑛∈ℕ.

Therefore, can derive its analytic continuation via

Γ(𝑧)=Γ(𝑧+𝑛)∏𝑘=0𝑛−1(𝑧+𝑘),ℜ𝔢(𝑧)>−𝑛.

Since the numerator is holomorphic on ℜ𝔢(𝑧)>−𝑛 and 𝑛 was arbitrary, the analytic continuation of Γ has simple poles at each of ℤ≤0. Hence, Γ(𝑧) is meromorphic on ℂ.

By (4.5.2)(4.5.2)(4.5.2), the residue at each pole is equal to

Res𝑧=−𝑛Γ(𝑧)=lim𝑧→−𝑛Γ(𝑧+𝑛+1)∏𝑘=0𝑛−1(𝑧+𝑘)=1∏𝑘=1𝑛(−𝑘)=(−1)𝑛𝑛!.

We will now study two representations for the Gamma function.

Theorem 12.1.3 (Gauss).

The Gamma function satisfies

Γ(𝑧)=lim𝑛→∞𝑛𝑧𝑛!∏𝑘=0𝑛(𝑧+𝑘),ℜ𝔢𝑧>0.(12.1.2)
Proof.

Define the sequence of functions

𝑓𝑛(𝑧)=∫0𝑛(1−𝑡𝑛)𝑛𝑡𝑧−1d𝑡=𝑛𝑧∫01(1−𝑡)𝑛𝑡𝑧−1d𝑡,ℜ𝔢(𝑧)>0.

By integration by parts, we have

𝑓𝑛(𝑧)=𝑛𝑧[𝑡𝑧𝑧(1−𝑡)𝑛|01+𝑛𝑧∫01(1−𝑡)𝑛−1𝑡𝑧d𝑡]=(𝑛𝑛−1)𝑧+1𝑓𝑛−1𝑧+1𝑧=[𝑛𝑧+1𝑛−1(𝑛−2)𝑧+2]𝑓𝑛−2𝑧+2𝑧(𝑧+1)=𝑛𝑧+1(𝑛−1)!𝑓1𝑧+𝑛−1∏𝑘=0𝑛−2(𝑧+𝑘)=𝑛𝑧𝑛!∏𝑘=0𝑛(𝑧+𝑘).(12.1.3)

Let us now analyze the difference

lim𝑛→∞[∫0𝑛e−𝑡𝑡𝑧−1d𝑡−𝑓𝑛(𝑧)]=lim𝑛→∞∫0𝑛e−𝑡𝑡𝑧−1[1−e𝑡(1−𝑡𝑛)𝑛]d𝑡.(12.1.4)

Since

de𝑡(1−𝑡𝑛)𝑛d𝑡=e𝑡(1−𝑡𝑛)𝑛−e𝑡(1−𝑡𝑛)𝑛−1=−e𝑡𝑡𝑛(1−𝑡𝑛)𝑛−1,

we have

1−e𝑡(1−𝑡𝑛)𝑛=1𝑛∫0𝑡𝑢e𝑢(1−𝑢𝑛)𝑛−1d𝑢.(12.1.5)

Additionally, since

de𝑢(1−𝑢𝑛)𝑛−1d𝑢=e𝑢(1−𝑢𝑛)𝑛−1−𝑛−1𝑛e𝑢(1−𝑢𝑛)𝑛−2=e𝑢𝑛(1−𝑢𝑛)𝑛−2(1−𝑢)

has zeros at 𝑢=1 and at 𝑢=𝑛, and

d2e𝑢(1−𝑢𝑛)𝑛−1d𝑢2=e𝑢𝑛2(1−𝑢𝑛)𝑛−3(𝑢2−2𝑢−𝑛+2),

evaluates to −e(𝑛−1)𝑛−2𝑛𝑛−1<0 at 𝑢=1 and evaluates to

e𝑛𝑛𝑛−1(𝑛−𝑢)𝑛−3(𝑛−2)(𝑛−1)→0+

as 𝑢→𝑛−, e𝑢(1−𝑢𝑛)𝑛−1 attains its maximum of e(𝑛−1𝑛)𝑛−1 at 𝑢=1. For 𝑛>1, e(𝑛−1𝑛)𝑛−1≤e. From (12.1.5)(12.1.5)(12.1.5), we have

1−e𝑡(1−𝑡𝑛)𝑛≤e𝑡22𝑛.

Moreover, since for 0≤𝑢≤𝑡≤𝑛,

𝑢e𝑢(1−𝑢𝑛)𝑛−1>0,

it follows that 1−e𝑡(1−𝑡𝑛)𝑛 is positive. By (12.1.4)(12.1.4)(12.1.4), we have

|∫0𝑛e−𝑡𝑡𝑧−1[1−e𝑡(1−𝑡𝑛)𝑛]d𝑡|≤e2𝑛|∫0𝑛e−𝑡𝑡𝑧+1d𝑡|<12𝑛|Γ(𝑧+2)|→0
as 𝑛→∞. From (12.1.3)(12.1.3)(12.1.3), we have Γ(𝑧)=lim𝑛→∞𝑛𝑧𝑛!∏𝑘=0𝑛(𝑧+𝑘), or (12.1.2)(12.1.2)(12.1.2) ⁠□

The Weierstrass formula is a direct consequence of the Gauss formula.

Theorem 12.1.4 (Weierstrass).

The reciprocal Γ-function has the entire Weierstrass factorization of

1Γ(𝑧)=𝑧∏𝑘=1∞[(1+𝑧𝑘)e−𝑧𝑘]e𝑧𝛾,(12.1.6)
where 𝛾=∫1∞(1⌊𝑥⌋−1𝑥)d𝑥.
Proof.

Since the Gauss formula agrees with (12.1.1)(12.1.1)(12.1.1) on the right half-plane, the analytic continuation of Γ(𝑧) is unique on the entire complex plane except for the poles at ℤ≤0 by the Identity Theorem (Theorem 3.3.3Theorem 3.3.3Theorem 3.3.3). Since

𝑛𝑧𝑛!∏𝑘=0𝑛(𝑧+𝑘)=exp(𝑧log(𝑛))𝑧∏𝑘=1𝑛(1+𝑧𝑘)=exp(𝑧∫1𝑛d𝑥𝑥)𝑧∏𝑘=1𝑛(1+𝑧𝑘)exp(−𝑧∑𝑘=1𝑛1𝑘)∏𝑘=1𝑛exp(−𝑧𝑘)=1𝑧∏𝑘=1𝑛[(1+𝑧𝑘)e−𝑧𝑘]exp(−𝑧(∫1𝑛(1⌊𝑥⌋−1𝑥)d𝑥)).

Therefore,

1Γ(𝑧)=lim𝑛→∞∏𝑘=0𝑛(𝑧+𝑘)𝑛𝑧𝑛!=𝑧∏𝑘=1𝑛[(1+𝑧𝑘)e−𝑧𝑘]lim𝑛→∞exp(𝑧(∫1𝑛(1⌊𝑥⌋−1𝑥)d𝑥))=𝑧∏𝑘=1𝑛[(1+𝑧𝑘)e−𝑧𝑘]exp(𝑧𝛾).

The constant 𝛾=∫1∞(1⌊𝑥⌋−1𝑥)d𝑥 is known as the Euler–Mascheroni constant. For each 𝑥, 0≤1⌊𝑥⌋−1𝑥<1⌊𝑥⌋−1⌊𝑥⌋+1, and the integral is bounded by ∑𝑛=1∞(1𝑛−1𝑛+1)<1, and hence the integral is convergent. The properties of 𝛾 still largely remain unknown.

By the Weierstrass Factorization Theorem (Theorem 4.4.3.5Theorem 4.4.3.5Theorem 4.4.3.5), if we let 𝑎𝑛=−𝑛 and 𝑝𝑛=1, it follows that

∑𝑛=1∞|𝑅𝑎𝑛|𝑝𝑛+1=𝑅2∑𝑛=1∞1𝑛2=𝑅2π26

is convergent. Thus, the Weierstrass formula defines an entire function with zeros at each of ℤ≤0. ⁠□

We have two famous identities on the Γ-function:

Theorem 12.1.5 (Euler’s Reflection Formula).

The Γ-function can be analytically continued to the left half-plane with the functional equation of

Γ(𝑧)Γ(1−𝑧)=πcsc(π𝑧)(12.1.7)
for 𝑧∈ℂ∖ℤ.
Proof.

By the Weierstrass Formula (Theorem 12.1.4Theorem 12.1.4Theorem 12.1.4), we have

1Γ(𝑧)=𝑧∏𝑘=1𝑛[(1+𝑧𝑘)e−𝑧𝑘]e𝑧𝛾,1Γ(−𝑧)=−𝑧∏𝑘=1𝑛[(1−𝑧𝑘)e𝑧𝑘]e−𝑧𝛾.

Since the Weierstrass elementary factors form an absolutely convergent infinite product, we may rearrange its terms. Hence, by Theorem 4.4.5.4Theorem 4.4.5.4Theorem 4.4.5.4, we have

1Γ(𝑧)Γ(1−𝑧)=−1𝑧Γ(𝑧)Γ(−𝑧)=𝑧∏𝑘=1𝑛(1−𝑧2𝑘2)=sin(π𝑧)π,
which confirms (12.1.7)(12.1.7)(12.1.7). ⁠□

Example 12.1.6.

Evaluate Γ(12).
Solution to Example 12.1.6Example 12.1.6Example 12.1.6.

By the Reflection Formula (Theorem 12.1.5Theorem 12.1.5Theorem 12.1.5), we have that

Γ(12)2=πcsc(π2)=π,
and it follows that Γ(12)=π as it is positive. ⁠□

Theorem 12.1.7 (Legendre’s Duplication Formula).

For any 𝑧∈ℂ∖(−ℕ2), we have

Γ(𝑧)Γ(𝑧+12)=21−2𝑧πΓ(2𝑧).(12.1.8)
Proof.

From Theorem 12.1.3Theorem 12.1.3Theorem 12.1.3, we have

Γ(𝑧)Γ(𝑧+12)=lim𝑛→∞𝑛2𝑧+12𝑛!2∏𝑘=0𝑛(𝑧+𝑘)(𝑧+𝑘+12)=lim𝑛→∞22𝑛+2𝑛2𝑧+12𝑛!2∏𝑘=02𝑛+1(2𝑧+𝑘)

where the left-hand side is defined since 𝑧∈ℂ∖(−ℕ2). By expansion of the value, we have

Γ(𝑧)Γ(𝑧+12)=lim𝑛→∞(2𝑛)2𝑧(2𝑛)!∏𝑘=02𝑛(2𝑧+𝑘)⋅𝑛12𝑛!222𝑛+2−2𝑧(2𝑧+2𝑛+1)(2𝑛)!=Γ(2𝑧)lim𝑛→∞𝑛12𝑛!222−2𝑧(2𝑧+2𝑛+1)∏𝑘=0𝑛−1(𝑘+12)∏𝑘=1𝑛𝑘=22−2𝑧Γ(2𝑧)lim𝑛→∞𝑛12𝑛!∏𝑘=0𝑛(𝑘+12)⋅𝑛+122𝑧+2𝑛+1=21−2𝑧Γ(2𝑧)Γ(12)=21−2𝑧Γ(2𝑧)π,
where the last step is derived from Example 12.1.6Example 12.1.6Example 12.1.6. ⁠□

The identity above is a special case of the following result:

Theorem 12.1.8 (Gauss Multiplication Theorem).

Suppose 𝑚∈ℕ≥2. Let 𝑧∈ℂ∖(−ℕ𝑚). Then we have

Γ(𝑧)Γ(𝑧+1𝑚)⋯Γ(𝑧+𝑚−1𝑚)=(2π)𝑚−12𝑚12−𝑚𝑧Γ(𝑚𝑧).(12.1.9)

The Gamma function as in (12.1.1)(12.1.1)(12.1.1) is commonly referred to as the Euler Integral of the Second Kind.

Definition 12.1.9 (Beta Function).

The Euler Integral of the First Kind is also known as the Beta function, and is defined by

Β(𝑧1,𝑧2)=∫01𝑡𝑧1−1(1−𝑡)𝑧2−1d𝑡.

By a change of variables (by letting 𝜏=1−𝑡), we derive the symmetry of the Beta function:

Β(𝑧1,𝑧2)=∫01𝜏𝑧2−1(1−𝜏)𝑧1−1d𝜏=Β(𝑧2,𝑧1).

The Beta function is commonly treated as an auxiliary function in many cases of integral evaluation due to its connection with the Gamma function:

Theorem 12.1.10.

For any ℜ𝔢(𝑧1),ℜ𝔢(𝑧2)>0, we have

Β(𝑧1,𝑧2)=Γ(𝑧1)Γ(𝑧2)Γ(𝑧1+𝑧2).
Proof.

Consider the product Γ(𝑧1)Γ(𝑧2). By letting 𝑠=𝑢𝑡 and 𝑣=𝑡(𝑢+1), we have

Γ(𝑧1)Γ(𝑧2)=∫0∞e−𝑠𝑠𝑧2−1[∫0∞e−𝑡𝑡𝑧1−1d𝑡]d𝑠=∫0∞𝑢𝑧2−1[∫0∞e−𝑣(𝑣𝑢+1)𝑧1+𝑧2−1d(𝑣𝑢+1)]d𝑢=∫0∞𝑢𝑧2−1(𝑢+1)𝑧1+𝑧2[∫0∞e−𝑣𝑣𝑧1+𝑧2−1d𝑣]d𝑢.(12.1.10)

Let 𝑟=𝑢𝑢+1, 𝑢=𝑟1−𝑟, and d𝑢=1(1−𝑟)2d𝑟. Then we have

Γ(𝑧1)Γ(𝑧2)=Γ(𝑧1+𝑧2)∫01𝑟𝑧2−1(1−𝑟)𝑧1−1d𝑟=Γ(𝑧1+𝑧2)Β(𝑧1,𝑧2).□

Example 12.1.11 (MIT Integration Bee 2023 Finals #1).

Evaluate

∫0π2tan(𝑥)13(sin(𝑥)+cos(𝑥))2d𝑥.
Solution to Example 12.1.11Example 12.1.11Example 12.1.11.

By rewriting the integral, and letting 𝑢=tan(𝑥), we have

𝐼=∫0π2tan(𝑥)13(sin(𝑥)+cos(𝑥))2d𝑥=∫0∞𝑢13sec2(𝑥)d𝑥(𝑢+1)2=∫0∞𝑢13d𝑢(𝑢+1)2=Β(43,23),

where the last step recognizes the form of (12.1.10)(12.1.10)(12.1.10). Theorem 12.1.10Theorem 12.1.10Theorem 12.1.10 then gives

𝐼=Γ(23)Γ(43)Γ(2)=13Γ(13)Γ(23).

Lastly, the Reflection Formula (Theorem 12.1.5Theorem 12.1.5Theorem 12.1.5) gives that

𝐼=π3sin(π3)=2π39.□

Example 12.1.12.

Evaluate

𝐼=∫0π2𝑥tan𝑥d𝑥.
Proof of Example 12.1.12Example 12.1.12Example 12.1.12.

Apply the substitution 𝑢=tan𝑥, 𝑠=𝑢2 and the integral representation arctan𝑢=∫0𝑢d𝑣1+𝑣2=∫01𝑢d𝑣1+𝑢2𝑣2 to get

𝐼=∫0∞arctan𝑢𝑢d𝑢1+𝑢2=∫0∞(∫01𝑢1+𝑢2𝑣2d𝑣)𝑢d𝑢1+𝑢2=∫01∫0∞𝑢32(1+𝑢2𝑣2)(1+𝑢2)d𝑢d𝑣=∫01∫0∞𝑠142(1+𝑠𝑣2)(1+𝑠)d𝑠d𝑣.

The integral switch is justified since the integrand ≥0 throughout. Consider the integral for

∫0∞𝑠14d𝑠(1+𝑠𝑣2)(𝑠+1)=11−𝑣2∫0∞𝑠𝛽(11+𝑠−1𝑐+𝑠)d𝑠,

where 𝑐=𝑣−2 and 𝛽=14. Let the integral in the last expression be a complex function of 𝛽, more explicitly 𝐽(𝛽), defined for −1<ℜ𝔢𝛽<1.

Notice that the integral maybe written as the difference of two convergent integrals for −1<ℜ𝔢𝛽<0:

𝐽(𝛽)=∫0∞𝑠𝛽d𝑠1+𝑠−∫0∞𝑠𝛽𝑐𝛽d𝑠1+𝑠=(1−𝑐𝛽)Β(−𝛽,𝛽+1)=(𝑐𝛽−1)πcsc(π𝛽)

by the reflection formula (Theorem 12.1.5Theorem 12.1.5Theorem 12.1.5), which derives an analytic continuation of 𝐽|−1<ℜ𝔢𝛽<1 to everywhere in ℂ∖ℤ.

One can show by differentiation under the integral sign, that the integral 𝐽(𝛽) is holomorphic for −1<ℜ𝔢𝛽<1. The identity theorem (Theorem 3.3.3Theorem 3.3.3Theorem 3.3.3) then allows the computed cosecant expression to be assigned to 𝐽(𝛽):

𝐼=∫01(𝑣−12−1)πcsc(π4)2−2𝑣2d𝑣=π22∫01𝑣−12−11−𝑣2d𝑣=π2∫011−𝑢1−𝑢4d𝑢=π2∫01d𝑢(1+𝑢)(1+𝑢2)=π22∫01(11+𝑢+1−𝑢1+𝑢2)d𝑢=π22(log2+arctan1−12log2)=π228+π24log2.□