Notes on Complex Analysis
Export PDFGitHub

§4.5 The Residue Theorem

4.5 The Residue Theorem

After Riemann and Weierstrass refined the understanding of analytic functions and the formal characterization of Jordan curves, the Cauchy Residue Theorem was consequently formalized. Cauchy had the informal notion of a residue, which we will now formally introduce.

Definition 4.5.1 (Residue).

For some 𝑟∈ℝ>0, 𝑎∈𝑈, suppose 𝑓:𝐷∗(𝑎,𝑟)→ℂ is holomorphic. Then the residue of 𝑓 at 𝑎, denoted by Res𝑧=𝑎𝑓(𝑧) or Res(𝑓,𝑎), is equal to

Res𝑧=𝑎𝑓(𝑧)=12πi∮𝜕𝐷(𝑎,𝜌)𝑓(𝑧)d𝑧,(4.5.1)

where 0<𝜌<𝑟 is arbitrary. Since 𝑓 has a Laurent expansion at 𝑎, being

∑𝑛=−∞∞𝑐𝑛(𝑧−𝑎)𝑛,𝑐𝑛=12πi∮𝜕𝐷(𝑎,𝜌)𝑓(𝑧)d𝑧(𝑧−𝑎)𝑛+1,
we get that the residue of 𝑓 at 𝑎 is equal to the first term 𝑐−1 of the principal part of its Laurent expansion.

It then follows that the residue at a removable singularity is 0. As a direct consequence of (4.5.1)(4.5.1)(4.5.1), we can derive explicit formulas for the calculation of residues at poles. If 𝑈⊆ℂ is open, 𝑎∈𝑈 is an isolated singularity (a pole of order 𝑚≠∞) of 𝑓:𝑈∖{𝑎}→ℂ that is holomorphic, then locally:

𝑓(𝑧)=𝑐−𝑚(𝑧−𝑎)−𝑚+𝑐1−𝑚(𝑧−𝑎)1−𝑚+⋯+𝑐−1(𝑧−𝑎)−1+….

Multiplying by (𝑧−𝑎)𝑚, we obtain that

(𝑧−𝑎)𝑚𝑓(𝑧)=𝑐−𝑚+𝑐1−𝑚(𝑧−𝑎)+⋯+𝑐−1(𝑧−𝑎)𝑚−1+….

By the definition of a Taylor series, we find that

𝑐−1=Res𝑧=𝑎𝑓(𝑧)=1(𝑚−1)!lim𝑧→𝑎d𝑚−1d𝑧𝑚−1[(𝑧−𝑎)𝑚𝑓(𝑧)].(4.5.2)

Let 𝑧=∞ be an isolated singularity of 𝑓(𝑧), which is holomorphic in ℂ∖𝐷(0,𝑅)‾, for sufficiently large finite 𝑅. Then for finite 𝜌>𝑅, the residue at 𝑧=∞ is defined as (notice the orientation)

Res𝑧=∞𝑓(𝑧)=12πi∲𝜕𝐷(0,𝜌)𝑓(𝑧)d𝑧.

Let 𝜁=1𝑧. Then we get that

Res𝑧=∞𝑓(𝑧)=−12πi∳𝜕𝐷(0,𝜌)𝑓(1𝜁)d(1𝜁)=12πi∳𝜕𝐷(0,1/𝜌)𝑓(1/𝜁)𝜁2d𝜁=−Res𝜁=0𝑓(1/𝜁)𝜁2.

In this definition, if

𝑓(𝑧)=∑𝑛=−∞∞𝑐𝑛𝑧𝑛⟺𝑓(1/𝜁)𝜁2=∑𝑛=−∞∞𝑐𝑛𝜁−𝑛−2,

the residue at 𝑧=∞ is equal to −𝑐−1. We will later explain the reasoning behind this definition.

Theorem 4.5.2 (Residue Theorem).

Let 𝑈⊂ℂ be an open set with a simple closed boundary curve 𝜕𝑈. Suppose {𝑧𝑛}⊂𝑈 is a finite set and 𝑓(𝑧) is holomorphic on 𝑈∖{𝑧𝑛} and continuous on 𝑈‾∖{𝑧𝑛}. Then,

∮𝜕𝑈𝑓(𝑧)d𝑧=2πi∑𝑘=1𝑛Res𝑧=𝑧𝑘𝑓(𝑧)
Proof.

Since 𝑈 is open, there exists a small disk centered at each isolated singularity 𝑧𝑘 of radii 𝛿𝑘. By the Cauchy–Goursat Theorem (Theorem 3.1.9Theorem 3.1.9Theorem 3.1.9), we get that

∫⋃𝑘=1𝑛𝐷(𝑧𝑘,𝛿𝑘)−∪𝜕𝑈+𝑓(𝑧)d𝑧=0.
From rearrangement, ∮𝜕𝑈𝑓(𝑧)d𝑧=∑𝑘=1𝑛∮𝜕𝐷(𝑧𝑘,𝛿𝑘)𝑓(𝑧)d𝑧, and the conclusion follows. ⁠□

This result itself is fairly trivial. Now we will explain the significance of the residue at infinity.

Theorem 4.5.3 (Global Residue Theorem).

If {𝑧1,…𝑧𝑛,∞} is discrete and finite, and 𝑓:ℂ̂∖{𝑧1,…𝑧𝑛,∞}→ℂ is holomorphic, and these points are the isolated singularities of 𝑓, then the sum of the residues at each of these isolated singularities is zero, or

∑𝑘=1𝑛Res𝑧=𝑧𝑘𝑓(𝑧)+Res𝑧=∞𝑓(𝑧)=0.
Proof.

Let 𝑅>max𝑗∈ℕ≤𝑛|𝑧𝑛| be arbitrary. By the Residue Theorem (Theorem 4.5.2Theorem 4.5.2Theorem 4.5.2),

−Res𝑧=∞𝑓(𝑧)=12πi∮𝜕𝐷(0,𝑅)𝑓(𝑧)d𝑧=∑𝑘=1𝑛Res𝑧=𝑧𝑘𝑓(𝑧)
as desired. This is merely a restatement of Theorem 4.5.2Theorem 4.5.2Theorem 4.5.2. ⁠□

There is not a directly trivial reason for the definition of the residue at ∞, except for the fact that it seemingly “unifies” the Riemann sphere.

Figure 4.2: The orientation of a neighborhood that does not enclose ∞ after projection.

However, if we take a neighborhood of an arbitrary point in ℂ on the Riemann sphere and traverse its boundary clockwise (from the perspective of outside the sphere), its projection onto ℂ will be counterclockwise (Figure 4.2Figure 4.2Figure 4.2). However, the boundary of a neighborhood of ∞ in 𝑆2 will have a clockwise projection (hence the difference in orientation). We define its equality with the residue of −𝑓(1/𝜁)𝜁2 at 𝜁=0, rather than 𝑓(1/𝜁), because we compose the differential form 𝑓(𝑧)d𝑧 with the inversion, as opposed to 𝑓(𝑧).

For any closed rectifiable curve 𝛾⊂𝑈 (here we are not bound under the assumption of simpleness), the Residue Theorem can be generalized into:

∮𝛾𝑓(𝑧)d𝑧=2πi∑𝑘Ind𝛾(𝑧𝑘)Res𝑧=𝑧𝑘𝑓(𝑧)

where 𝑧𝑘 are the singularities of 𝑓 in 𝑈 and Ind𝛾 is the winding index.

Residues are extremely important as they allow for simple evaluation of definite (most commonly improper) real-valued integrals. This is because oftentimes, residues at poles are generally easy to calculate and have an integral representation. We can integrate over a contour (a smooth closed curve) that contains the important part of the real interval. Oftentimes this is the most non-trivial step.

Example 4.5.4.

Evaluate the improper integral 𝐼=∫−∞∞1(𝑥2+1)𝑛+1d𝑥, where 𝑛∈ℕ.
Figure 4.3: A semicircular contour with orientation marked.
Solution to Example 4.5.4Example 4.5.4Example 4.5.4.

Consider 𝛾 to be a closed semicircle with radius 𝑅≥2 as in Figure 4.3Figure 4.3Figure 4.3. Notice that the function 𝑧↦1(𝑧2+1)𝑛+1 has singularities at only 𝑧=i and 𝑧=−i, both of which are poles of order 𝑛+1. By (4.5.2)(4.5.2)(4.5.2), the residue at 𝑧=i is

Res𝑧=i1(𝑧2+1)𝑛+1=1𝑛!d𝑛d𝑧𝑛(𝑧+i)−𝑛−1|𝑧=i=1𝑛!(−1)𝑛∏𝑘=1𝑛(𝑛+𝑘)(2i)2𝑛+1=(−1)𝑛(2𝑛)!(𝑛!)2(2i)2𝑛+1=(2𝑛)!22𝑛+1i(𝑛!)2.

The singularity at 𝑧=−i is not relevant, as it is not enclosed by the contour. By the Residue Theorem (Theorem 4.5.2Theorem 4.5.2Theorem 4.5.2), we have

∮𝛾1(𝑧2+1)𝑛+1d𝑧=∫−𝑅𝑅1(𝑥2+1)𝑛+1d𝑥+∫0π𝑅i(𝑅2e2i𝜃+1)𝑛+1ei𝜃d𝜃=2πiRes𝑧=i1(𝑧2+1)𝑛+1=(2𝑛)!π22𝑛(𝑛!)2.

We will now show that the integral over the semicircle vanishes as 𝑅→∞. Under the assumption that 𝑅≥2, since

|𝑅iei𝜃(𝑅2e2i𝜃+1)𝑛+1|=𝑅|𝑅2e2i𝜃+1|𝑛+1≤𝑅|𝑅2−1|≤23,

which is integrable over [0,π], and we can commute the limit with the integral. Therefore, we have

∫−∞∞1(𝑥2+1)𝑛+1d𝑥=lim𝑅→∞∮𝛾1(𝑧2+1)𝑛+1d𝑧−∫0πlim𝑅→∞𝑅i(𝑅2e2i𝜃+1)𝑛+1ei𝜃d𝜃=(2𝑛)!π22𝑛(𝑛!)2.□

Example 4.5.5 (Dirichlet Integral).

Evaluate the integral ∫0∞sin𝑥𝑥d𝑥.
Figure 4.4: An indented semicircular contour with orientation marked.
Solution to Example 4.5.5Example 4.5.5Example 4.5.5.

It is common to use integration with parameters to approach this integral. However, we will now provide a solution via contour integration.

Let 𝑓(𝑧)=ei𝑧𝑧. Consider a closed contour 𝛾 in the form of Figure 4.4Figure 4.4Figure 4.4, consisting of a semicircle of radius 𝑅 in ℍ+‾ (𝐶𝑅), a line segment from −𝑅 to −𝜀, a smaller semicircle of radius 𝜀 in the upper half-plane (𝐶𝜀), and a line segment from 𝜀 to 𝑅.

By the Cauchy–Goursat Theorem (Theorem 3.1.9Theorem 3.1.9Theorem 3.1.9), we have that

∮𝛾𝑓(𝑧)d𝑧=∫𝐶𝑅𝑓(𝑧)d𝑧+∫−𝑅−𝜀𝑓(𝑧)d𝑧+∫𝐶𝜀𝑓(𝑧)d𝑧+∫𝜀𝑅𝑓(𝑧)d𝑧=0.

We will now analyze each integral. The first integral is

∫𝐶𝑅𝑓(𝑧)d𝑧=∫0πexp(i𝑅ei𝜃)𝑅ei𝜃𝑅iei𝜃d𝜃=i∫0πei𝑅cos𝜃e−𝑅sin𝜃d𝜃.

Notice that 2π𝜃≤sin(𝜃)≤𝜃 over the integration range. We want to observe the behavior as 𝑅→∞:

|i∫0πei𝑅cos𝜃e−𝑅sin𝜃d𝜃|≤∫0πe−𝑅sin𝜃d𝜃=2∫0π/2e−𝑅sin𝜃d𝜃<2∫0π/2e−𝑅2π𝜃d𝜃=−π𝑅e−𝑅2π𝜃|0π/2=π𝑅(1−e−𝑅)→0.

Let us evaluate the integral on 𝛾𝜀 as 𝜀→0:

∫𝐶𝜀𝑓(𝑧)d𝑧=i∫π0exp(𝜀(icos𝜃−sin𝜃))d𝜃=i∫π0e−𝜀sin𝜃ei𝜀cos𝜃d𝜃.

Obviously,

|e−𝜀sin𝜃ei𝜀cos𝜃|≤1,

and therefore, the integral and the limit may commute:

lim𝜀→0∫𝐶𝜀𝑓(𝑧)d𝑧=i∫π0lim𝜀→0+e−𝜀sin𝜃ei𝜀cos𝜃d𝜃=i∫π0d𝜃=−iπ.

Evaluating the integral over the line segments, we have

∫−𝑅−𝜀𝑓(𝑧)d𝑧+∫𝜀𝑅𝑓(𝑧)d𝑧=∫−𝑅−𝜀ei𝑧𝑧d𝑧+∫𝜀𝑅ei𝑧𝑧d𝑧=∫−𝑅−𝜀ei𝑧𝑧d𝑧−∫−𝑅−𝜀e−i𝑧𝑧d𝑧→2i∫−∞0sin𝑧𝑧d𝑧=2i∫0∞sin𝑧𝑧d𝑧.

Hence,

−iπ+2i∫0∞sin(𝑧)𝑧d𝑧=0⟺∫0∞sin𝑧𝑧d𝑧=π2.□

Example 4.5.6 (Fresnel Integrals).

Evaluate the improper integrals

𝐼1=∫0∞cos(𝑥2)d𝑥,𝐼2=∫0∞sin(𝑥2)d𝑥.
Figure 4.5: A wedge contour with orientation marked.
Solution to Example 4.5.6Example 4.5.6Example 4.5.6.

Let 𝑓(𝑧)=ei𝑧2. Choose the wedge contour composed of

Γ1={𝑥∈ℝ:0≤𝑥≤𝑅},Γ2={𝑟eiπ/4:0≤𝑟≤𝑅},𝐶𝑅={𝑅ei𝜃:0≤𝜃≤π4}

as in Figure 4.5Figure 4.5Figure 4.5. By the Cauchy–Goursat Theorem (Theorem 3.1.9Theorem 3.1.9Theorem 3.1.9), we have that

∫Γ1𝑓(𝑧)d𝑧+∫Γ2𝑓(𝑧)d𝑧+∫𝐶𝑅𝑓(𝑧)d𝑧=0.(4.5.3)

The third integral can be written as

∫𝐶𝑅𝑓(𝑧)d𝑧=𝑅i∫0π/4exp[i(𝑅ei𝜃)2]ei𝜃d𝜃.

Using the fact that 4π𝜃<sin(2𝜃) on the integration range, it can be bounded as

|∫𝐶𝑅𝑓(𝑧)d𝑧|≤𝑅∫0π/4e−𝑅2sin(2𝜃)d𝜃<𝑅∫0π/4e−4π𝑅2𝜃d𝜃=−π4𝑅e−4π𝑅2𝜃|0π/4=π4𝑅(1−e−𝑅2).

As 𝑅→∞, this integral tends to 0. Let 𝑧=𝑟eiπ/4 on Γ2. Then, we have

lim𝑅→∞∫Γ2𝑓(𝑧)d𝑧=∫∞0exp[i(𝑟eiπ/4)2]eiπ/4d𝑟=eiπ/4∫∞0exp(−𝑟2)d𝑟.

From (4.5.3)(4.5.3)(4.5.3), we have that

∫0∞ei𝑟2d𝑟=eiπ/4∫0∞e−𝑟2d𝑟.

Since ∫0∞e−𝑟2d𝑟=π2, we have

∫0∞ei𝑟2d𝑟=(22+i22)π2.

Since ei𝑟2=cos(𝑟2)+isin(𝑟2), we have

∫0∞cos(𝑟2)d𝑟=ℜ𝔢[∫0∞ei𝑟2d𝑟]=2π4,∫0∞sin(𝑟2)d𝑟=ℑ𝔪[∫0∞ei𝑟2d𝑟]=2π4,
as desired. ⁠□

Example 4.5.7.

Evaluate the integrals ∫02πΦ(cos𝜃,sin𝜃)d𝜃, where Φ(𝜉,𝜂) is a rational function of 𝜉 and 𝜂 that is continuous on 𝜃∈[0,2π].
Solution to Example 4.5.7Example 4.5.7Example 4.5.7.

Let 𝑧=ei𝜃. Consequently, we have cos𝜃=𝑧+𝑧−12, sin𝜃=𝑧−𝑧−12i, and d𝑧=iei𝜃d𝜃, implying that d𝜃=d𝑧i𝑧. Therefore, by the Residue Theorem (Theorem 4.5.2Theorem 4.5.2Theorem 4.5.2), letting 𝑓(𝑧)=1i𝑧Φ(𝑧+𝑧−12,𝑧−𝑧−12i), we have

∫02πΦ(cos𝜃,sin𝜃)d𝜃=∮𝜕𝔻𝑓(𝑧)d𝑧=2πi∑𝑘=1𝑛Res𝑧=𝑧𝑘𝑓(𝑧),
where 𝑧𝑘 where 𝑘=1,…,𝑛 are the isolated singularities of 𝑓 in 𝔻. ⁠□

Example 4.5.8.

Evaluate 𝐼=∫0∞𝑥𝛼1+𝑥𝛽d𝑥, where 0<𝛼+1<𝛽.
Figure 4.6: An indented wedge contour with orientation marked.
Solution to Example 4.5.8Example 4.5.8Example 4.5.8.

Let 𝑓(𝑧)=𝑧𝛼1+𝑧𝛽 and let −π<Arg(𝑧)≤π in the principal branches of 𝑧𝛼=e𝛼Log(𝑧) and 𝑧𝛽=e𝛽Log(𝑧). Then except for at the zeros of 1+𝑧𝛽, 𝑓 is holomorphic.

The solutions to 𝑧𝛽=−1 are 𝑧=exp(iπ𝛽+2i𝑘π𝛽). Choose an indented wedge contour (as there is a logarithmic branch point singularity at the origin) with an angle of 2π𝛽 (as in Figure 4.5Figure 4.5Figure 4.5). The only singularity it encloses is exp(iπ𝛽). Since it is a simple zero of 1𝑓, this singularity is a simple pole.

The contour is the union of the following curves:

Γ1={𝑥∈ℝ:𝜀≤𝑥≤𝑅},Γ2={𝑟exp(i2π/𝛽):𝜀≤𝑟≤𝑅},𝐶𝑅={𝑅ei𝜃:0≤𝜃≤2π/𝛽},𝐶𝜀={𝜀ei𝜃:0≤𝜃≤2π/𝛽}

where 𝑅>1 and 0<𝜀<1. By the Residue Theorem (Theorem 4.5.2Theorem 4.5.2Theorem 4.5.2), we get that

lim𝜀→0lim𝑅→∞(∫Γ1+∫Γ2+∫𝐶𝑅+∫𝐶𝜀)𝑓(𝑧)d𝑧=2πiRes[𝑓,exp(iπ𝛽)].

By (4.5.2)(4.5.2)(4.5.2), it follows that

Res[𝑓,exp(iπ/𝛽)]=lim𝑧→exp(iπ/𝛽)𝑧−exp(iπ/𝛽)(1+𝑧𝛽)/𝑧𝛼=lim𝑧→exp(iπ/𝛽)(dd𝑧(𝑧−𝛼+𝑧𝛽−𝛼))−1=lim𝑧→exp(iπ/𝛽)𝑧𝛼+1(𝛽−𝛼)𝑧𝛽−𝛼=−1𝛽exp(iπ𝛽(𝛼+1)).

We can write the integral on Γ2 in terms of 𝐼:

lim𝑅→∞𝜀→0∫Γ2𝑓(𝑧)d𝑧=lim𝑅→∞𝜀→0∫𝑅0𝑓[𝑟exp(i2π/𝛽)]exp(i2π/𝛽)d𝑟=−exp[i2π𝛽(1+𝛼)]∫0∞𝑟𝛼1+𝑟𝛽d𝑟.

We also have

∫𝐶𝑅𝑓(𝑧)d𝑧=𝑅i∫02π/𝛽𝑓(𝑅ei𝜃)ei𝜃d𝜃=i∫02π/𝛽𝑅𝛼+11+𝑅𝛽ei𝛽𝜃exp[i𝜃(1+𝛼)]d𝜃.

It can also be shown that the integral is bounded by a vanishing function as 𝑅→∞:

|∫02π/𝛽𝑅𝛼+11+𝑅𝛽ei𝛽𝜃exp[i𝜃(1+𝛼)]d𝜃|≤∫02π/𝛽𝑅𝛼+1𝑅𝛽−1d𝜃=2π𝛽𝑅𝛼+1𝑅𝛽−1→0.

Similarly, as 𝜀→0,

|∫𝐶𝜀𝑓(𝑧)d𝑧|≤𝜀∫02π/𝛽|𝑓(𝜀ei𝜃)|d𝜃=∫02π/𝛽𝜀𝛼+11−𝜀𝛽d𝜃=2π𝛽𝜀𝛼+11−𝜀𝛽→0.

By letting 𝑅→∞ and 𝜀→0, we have

[1−exp(i2π𝛽(1+𝛼))]𝐼=−2πi𝛽exp(iπ𝛽(𝛼+1)).

It follows that

𝐼=2πi𝛽[exp(iπ𝛽(𝛼+1))−exp(−iπ𝛽(𝛼+1))]−1=π𝛽csc(π𝛽(𝛼+1)).□

Example 4.5.9.

Prove that the Fourier transform of sech(π𝑥) is itself, or that

𝐼(𝜉)=∫−∞∞exp(−2πi𝑥𝜉)sech(π𝑥)d𝑥=sech(π𝜉).
Figure 4.7: A rectangular contour with orientation marked.
Solution to Example 4.5.9Example 4.5.9Example 4.5.9.

Fix 𝜉∈ℝ and let 𝑓(𝑧)=exp(−2πi𝑧𝜉)cosh(π𝑧). Its poles in ℂ occur when eπ𝑧+e−π𝑧=0, or equivalently, when 𝑧=i(𝑛+12), where 𝑛∈ℤ.

Since

cosh(π(𝑧+i))=−cosh(π𝑧),exp(−2πi(𝑧+i)𝜉)=exp(2π𝜉)exp(−2πi𝑧𝜉),

we have that 𝑓(𝑧) is a constant multiple of 𝑓(𝑧+i). In particular, 𝑓(𝑧+i)=−exp(2π𝜉)𝑓(𝑧). Therefore, we can use a rectangular contour as shown in Figure 4.7Figure 4.7Figure 4.7. Let the sides be denoted by

Γ←={𝑥+i:−𝑅≤𝑥≤𝑅,𝑥∈ℝ},Γ→={𝑥∈ℝ:−𝑅≤𝑥≤𝑅}Γ↓={−𝑅+i𝑦:𝑦∈[0,1]},Γ↑={𝑅+i𝑦:𝑦∈[0,1]}.

The only enclosed singularity is a simple pole at 𝑧=i2 (simple by evaluation of the Taylor expansion of the denominator). By the Residue Theorem (Theorem 4.5.2Theorem 4.5.2Theorem 4.5.2), we get that

(∫Γ→+∫Γ↑+∫Γ←+∫Γ↓)𝑓(𝑧)d𝑧=2πiRes(𝑓,i2).(4.5.4)

By (4.5.2)(4.5.2)(4.5.2), we have

Res(𝑓,i2)=lim𝑧→i/2(𝑧−i2)exp(−2πi𝑧𝜉)cosh(π𝑧)=lim𝑧→i/2dd𝑧(cosh(π𝑧)exp(−2πi𝑧𝜉))−1=lim𝑧→i/2exp(−2πi𝑧𝜉)πsinh(π𝑧)+2πi𝜉cosh(π𝑧)=exp(π𝜉)πi.

The sum of the horizontal line integrals is equal to

∫−𝑅𝑅𝑓(𝑧)d𝑧+∫𝑅−𝑅𝑓(𝑧+i)d𝑧=∫−𝑅𝑅𝑓(𝑧)d𝑧−∫𝑅−𝑅e2π𝜉𝑓(𝑧)d𝑧=(1+e2π𝜉)∫−𝑅𝑅𝑓(𝑧)d𝑧.

As 𝑅→∞, we have ∫Γ→𝑓(𝑧)d𝑧+∫Γ←𝑓(𝑧)d𝑧→(1+e2π𝜉)𝐼(𝜉). The remaining two integrals can be written as

∫Γ↑𝑓(𝑧)d𝑧=∫01exp(−2πi(𝑅+i𝑧)𝜉)cosh(π(𝑅+i𝑧))d𝑧∫Γ↓𝑓(𝑧)d𝑧=∫10exp(2πi(𝑅−i𝑧)𝜉)cosh(π(−𝑅+i𝑧))d𝑧.

They can be bounded with

|∫01exp(−2πi(𝑅+i𝑧)𝜉)cosh(π(𝑅+i𝑧))d𝑧|≤2∫01exp(2π𝑧𝜉)|eπ𝑅eπi𝑧+e−π𝑅e−πi𝑧|d𝑧≤2∫01exp(2π𝑧𝜉)|eπ𝑅−e−π𝑅|d𝑧

and

|∫10exp(2πi(𝑅−i𝑧)𝜉)cosh(π(−𝑅+i𝑧))d𝑧|≤2∫01exp(2π𝑧𝜉)|e−π𝑅eπi𝑧+eπ𝑅e−πi𝑧|d𝑧≤2∫01exp(2π𝑧𝜉)|eπ𝑅−e−π𝑅|d𝑧.

Since the integrands are continuous and uniformly convergent to 0 with respect to 𝑧, we have

∫Γ↑𝑓(𝑧)d𝑧+∫Γ↓𝑓(𝑧)d𝑧→0

as 𝑅→∞. By rearrangement of (4.5.4)(4.5.4)(4.5.4),

𝐼(𝜉)(1+e2π𝜉)=2exp(π𝜉),

or that

𝐼(𝜉)=2e−π𝜉+eπ𝜉=sech(π𝜉),

which proves the result. ⁠□

Contour integration provides a powerful method for evaluating real improper integrals by leveraging the Residue Theorem (Theorem 4.5.2Theorem 4.5.2Theorem 4.5.2). The primary challenge often lies in constructing a suitable contour in the complex plane that encloses the relevant singularities of the integrand 𝑓 while ensuring that the contribution from the contributions from the remaining segments of the contour either vanishes or can be calculated with ease.

If the function 𝑓 is even and integrated on a domain such as ℝ≥0, then the integral can be extended to the entire real axis. If 𝑓 decays sufficiently rapidly in the upper half plane ℍ+, a semicircular contour is generally preferable, as illustrated in Figure 4.3Figure 4.3Figure 4.3. In the presence of singularities on the contour itself, we can insert arc indentations around them, as shown in Figure 4.4Figure 4.4Figure 4.4.

If 𝑓(𝑧) is a constant multiple of 𝑓(𝑧+i𝑦) (a type of quasi-periodicity) for some 𝑦∈ℝ, it is a strong indication to use a rectangular contour. If 𝑓(𝑧) is a constant multiple of 𝑓(𝑧ei𝜏) for some 𝜏∈ℝ, a wedge-shaped contour is an appropriate choice.

In the case that there are indentations along the contour, we have

Theorem 4.5.10.

Let 𝜆>0 and let 𝑎∈ℂ. Suppose 𝑓(𝑧) is a holomorphic function on 𝐷∗(𝑎,𝜆) with a simple pole at 𝑧=𝑎∈𝑈. Let 0<𝜀<𝜆 and define 𝛾𝜀⊆𝜕𝐷(𝑎,𝜀) be a counterclockwise-oriented, connected arc subtending an angle 𝜗. Then,

lim𝜀→0∫𝛾𝜀𝑓(𝑧)d𝑧=i𝜗⋅Res𝑧=𝑎𝑓(𝑧).
Proof.

Parameterize 𝛾𝜀 with 𝑧=𝑎+𝜀ei𝜃, where 𝜃∈[𝛼,𝛽] and 𝛽−𝛼=𝜗. Then,

∫𝛾𝜀𝑓(𝑧)d𝑧=∫𝛼𝛽𝑓(𝑎+𝜀ei𝜃)d𝑧d𝜃d𝜃=𝜀i∫𝛼𝛽𝑓(𝑎+𝜀ei𝜃)ei𝜃d𝜃.

Since 𝑓 has a simple pole at 𝑧=𝑎, we can write a Laurent expansion around 𝑎 as

𝑓(𝑧)=𝑐−1𝑧−𝑎+𝜑(𝑧),

where 𝜑(𝑧) is holomorphic in a neighborhood of 𝑎 and 𝑐−1=Res𝑧=𝑎𝑓(𝑧).

Then for 𝑧=𝑎+𝜀ei𝜃,

𝑓(𝑎+𝜀ei𝜃)=𝑐−1𝜀ei𝜃+𝜑(𝑎+𝜀ei𝜃).

So,

∫𝛾𝜀𝑓(𝑧)d𝑧=𝜀i∫𝛼𝛽(𝑐−1𝜀ei𝜃+𝜑(𝑎+𝜀ei𝜃))ei𝜃d𝜃=i𝑐−1∫𝛼𝛽d𝜃+𝜀i∫𝛼𝛽𝜑(𝑎+𝜀ei𝜃)ei𝜃d𝜃=i𝑐−1𝜗+𝜀i∫𝛼𝛽𝜑(𝑎+𝜀ei𝜃)ei𝜃d𝜃.

Let 𝜀<𝜆2. Since 𝜑 is continuous on the disk 𝐷(𝑎,𝜆2)‾, it is bounded. Therefore, letting 𝜀→0, we have

lim𝜀→0𝜀i∫𝛼𝛽𝜑(𝑎+𝜀ei𝜃)ei𝜃d𝜃=lim𝜀→0𝜀i∫𝛼𝛽𝜑(𝑎)ei𝜃d𝜃=0.

Therefore,

lim𝜀→0∫𝛾𝜀𝑓(𝑧)d𝑧=i𝜗Res𝑧=𝑎𝑓(𝑧).□

In the case that a branch point singularity is present on the contour, we may attempt to rewrite the function in a way such that the branch point is irrelevant. Otherwise, there are two types of “keyhole contours” that can be used to avoid the branch cut.

Example 4.5.11.

Evaluate 𝐼=∫0∞log(𝑥2+1)𝑥2+1d𝑥.
Solution to Example 4.5.11Example 4.5.11Example 4.5.11.

Notice that the integrand itself has branch points at 𝑧=±i coinciding with the poles from the denominator. We can rewrite the integral as

𝐼=12∫−∞∞log(𝑥2+1)𝑥2+1=∫−∞∞log(𝑥+i)(𝑥−i)𝑥2+1=∫−∞∞log|𝑥±i|𝑥2+1=ℜ𝔢∫−∞∞log(𝑥+i)𝑥2+1.(4.5.5)

Let 𝛾=Γ∪𝐶𝑅, where concretely,

Γ={𝑥∈ℝ:−𝑅≤𝑥≤𝑅},𝐶𝑅={𝑅ei𝜃:0≤𝜃≤π}

and 𝑅>2, and let 𝑓(𝑧)=Log(𝑧+i)𝑧2+1, where the branch for Log is chosen to satisfy [0,π]⊂ℑ𝔪log(ℂ∗), such as the principal branch. The only singularity of 𝑓 in the upper half plane is a simple pole at 𝑧=i. By the Residue Theorem (Theorem 4.5.2Theorem 4.5.2Theorem 4.5.2), we have

lim𝑅→∞∮𝛾𝑓(𝑧)d𝑧=lim𝑅→∞(∫Γ+∫𝐶𝑅)𝑓(𝑧)d𝑧=2πiRes𝑧=i𝑓(𝑧).

By (4.5.2)(4.5.2)(4.5.2), we have

Res𝑧=i𝑓(𝑧)=lim𝑧→i(𝑧−i)log(𝑧+i)𝑧2+1=lim𝑧→ilog(𝑧+i)𝑧+i=log(2i)2i=π4−ilog(2)2.

Additionally, for 𝑧∈𝐶𝑅, since as 𝑅→∞, |𝑓(𝑧)|=|Log(𝑧+i)𝑧2+1|≤|log|𝑧+i||+π𝑅2−1≤log|𝑅+1|+π𝑅2−1<𝑅+1+π𝑅2−1→0 by virtue of 𝑅>2, it follows that ∫𝐶𝑅𝑓(𝑧)d𝑧→0.

Since lim𝑅→∞∫Γ𝑓(𝑧)d𝑧=∫−∞∞𝑓(𝑧)d𝑧 and

∫−∞∞𝑓(𝑧)d𝑧=π2i2+πlog(2),

by (4.5.5)(4.5.5)(4.5.5), we have 𝐼=ℜ𝔢∫−∞∞𝑓(𝑧)d𝑧=πlog(2). ⁠□