CA
CA
F REDERICK T SZ -H O F ONG
Preface ix
Chapter 1. Preliminaries 1
1.1. Complex Numbers 1
1.2. Sequences and Series 12
1.3. Point-Set Topology of C 25
vii
Preface
「虛則實之,實則虛之。」
《孫子兵法》
This is the lecture note written for the course MATH 4023 - Complex Analysis, taught
by the author at the Hong Kong University of Science and Technology (HKUST) in Spring
2017 and Spring 2020.
The purpose of this lecture note and the course is to introduce both theory and
applications of complex-valued functions of one variable. It begins with basic notions of
complex differentiability (i.e. holomorphic) functions. The central part of the course is
the Cauchy’s integral formula, which is a fundamental theorem leading many important
and exciting results in the later half of the course. The last chapter of the course explains
the statement of the Riemann Hypothesis, a famous problem that worths US$1,000,000
in Pure Mathematics.
The prerequisites of the course include Multivariable Calculus in which students
should be familiar with line integrals and basic analysis of multivariable functions, and
Analysis I (as basic ε − δ languages are needed). It is also recommended that students
have taken Analysis II (MATH 3033) or Honors Analysis I (MATH 2043) before taking
this course, as some standard toolkits such as uniform convergence, Weierstass’s M-
test, Lebesgue Dominated Convergence Theorem, Fubini’s Theorem, etc. will appear
very frequently in the later half of the course. Students without MATH 2043/3033 are
recommended to first go through some examples in these topics. A list of theorems in
MATH 2043/3033 which are essential for this course can be found in the appendix of
this lecture note.
The author welcomes any students and/or readers to point out typographical errors
and mistakes of this lecture note.
Frederick Tsz-Ho Fong
15 January, 2020
HKUST, Clear Water Bay, Hong Kong
ix
Chapter 1
Preliminaries
√
i= −1 so that i2 = −1.
While complex numbers make their appearance for purely algebraic purposes, their uses
branch out to many scientific fields beyond Mathematics, including Quantum Mechanics,
String Theory, Electrical Engineering, Fluid Mechanics, etc.
Remark 1.2. Note that a real number is also considered as a complex number, since
a = a + 0i. In other words, we have R ⊂ C.
1
2 1. Preliminaries
Im
z = x + yi
y
x Re
−y
z = x − yi
z1 + z2 = (a + c) + (b + d)i
z1 − z2 = (a − c) + (b − d)i
z1 z2 = (a + bi)(c + di)
= (ac − bd) + (ad + bc)i
z1 a + bi c − di
= · (where z2 6= 0)
z2 c + di c − di
(ac + bd) (bc − ad)i
= 2 + 2
c + d2 c + d2
Remark 1.4. It is important to note that complex numbers are un-ordered. It does not
make sense to say z1 < z2 or z1 > z2 . However, since |z| is a real number, it makes sense
to make comparison of |z1 | and |z2 |.
Remark 1.5. Geometrically, z is obtained by reflecting z across the Re-axis (see Figure
1.1), and |z| is the magnitude of the position vector representing z.
Listed below are some very useful properties of complex numbers. Given any z, z1 , z2 ∈ C,
we have:
2
zz = |z| z=z |z| = |z|
z+z z−z
Re(z) = Im(z) =
2 2i
z1 z1
z1 ± z2 = z1 ± z2 z1 z2 = z1 z2 =
z2 z2
1.1. Complex Numbers 3
The proofs are all straight-forward and hence omitted. Simply let z = x + yi and
verify LHS and RHS are equal in each property. Let’s look at some examples on how to
make good use of these properties:
Example 1.1. Show that for any z1 , z2 ∈ C, we have:
2 2 2
|z1 + z2 | = |z1 | + |z2 | + 2Re(z1 z2 ).
Solution
2
The key step is to use the property that |z| = zz for any z ∈ C.
2
|z1 + z2 | = (z1 + z2 )(z1 + z2 ) = (z1 + z2 )(z1 + z2 )
= z1 z1 + z1 z2 + z1 z2 + z2 z2
2 2
= |z1 | + z1 z2 + z1 z2 + |z2 |
2 2
= |z1 | + |z2 | + 2Re(z1 z2 )
β
Example 1.2. Let α, β ∈ C\{0}. Show that αβ ∈ R if and only if α ∈ R.
Solution
1.1.3. Polar Form. There are two common types of coordinates in R2 , namely
rectangular and polar. Apart from the standard (rectangular) form x + yi for representing
a complex number, we can also represent a complex number by a polar form. The
conversion rule between rectangular and polar coordaintes is given by:
x = r cos θ
y = r sin θ
Therefore, a complex number z = x + yi can be written as:
z = (r cos θ) + i(r sin θ) = r(cos θ + i sin θ).
The form z = r(cos θ + i sin θ) is commonly called the polar form of z.
p
Note that |cos θ + i sin θ| = cos2 θ + sin2 θ = 1. When z = r(cos θ + i sin θ), it is
easy to see that r = |z|. However, the value of θ is not unique as both sin and cos are
periodic functions of period 2π. We define the principal argument of a complex number
to be the angle θ with a specified range described below:
Definition 1.7 (Principal Argument). Given a complex number z, the principal argument
of z, denoted by Arg(z), is defined to be the angle θ0 ∈ (−π, π] such that:
z = |z| (cos θ0 + i sin θ0 ).
√
For example, −1 − 3i has modulus 2 and so the r-coordinate is 2:
√ !
√ 1 3
−1 − 3i = 2 − − i .
2 2
√
To find the θ-coordinate, we solve cos θ = − 21 and sin θ = − 23 . From standard
trigonometry, we get θ = 4π3 + 2kπ for any integer k. The only θ that falls into the range
(−π, π] is − 2π
3 = 4π
3 − 2π. Therefore, we have:
√ 2π 2π
−1 − 3i = 2 cos − + i sin −
3 3
√
and Arg(−1 − 3i) = − 2π 3 .
Im
z = r(cos θ + i sin θ)
r
θ
Re
2 − 2π
3
√
z = −1 − 3i
1.1. Complex Numbers 5
In general, Arg(x + yi) can be found using tan−1 xy since if x = r cos θ and y = r sin θ,
then tan θ = xy . However, it is important to note that Arg(x + yi) is NOT simply equal to
tan−1 xy because by definition of the inverse tangent, tan−1 xy takes the value in (− π2 , π2 )
only. Precisely, we have (when x 6= 0):
−1 y
tan x if (x, y) is in 1st and 4th quadrants;
−1 y
Arg(x + yi) = tan x + π if (x, y) is in 2nd quadrant;
−1 y
tan x − π if (x, y) is in 3rd quadrant;
π
Furthermore, Arg(0 + yi) = 2 when y > 0; and Arg(0 + yi) = − π2 when y < 0. Note that
Arg(0 + 0i) is undefined.
Im
y y
tan−1 x +π tan−1 x
Re
y y
tan−1 x −π tan−1 x
Exercise 1.5. Express the following complex numbers in polar form, and find their
principal arguments Arg:
(a) 1 + 2i
(b) 1 − 2i
(c) cos(−π) + i sin(−π)
(d) −i
Using (1.1), we can see that given z1 = r1 (cos θ1 + i sin θ1 ) and z2 = r2 (cos θ2 +
i sin θ2 ), then we have:
Im z1 z2
z1
θ2
θ1
Re
Proof. We prove by induction for positive n’s. Clearly (1.2) is true when n = 1. Assume
that (1.2) is true when n = k for some positive integer k. Then, for n = k + 1, we have:
Hence (1.2) is true when n = k + 1. By induction, (1.2) is true for all positive integer n.
When n = 0, (1.2) also holds because (cos θ + i sin θ)0 = 1.
1.1. Complex Numbers 7
1.1.5. Roots of Complex Numbers. In the real number system, the root equation
xn = a where a 6= 0 and n ∈ N, has at most two solutions. When √ n is odd (no matter
whether a is positive or negative), the only
√ real solution
√ is x = n
a. When n is even and
a > 0, there are two real solutions x = n a or − n a. The equation has no solution when
n is even and a < 0.
However, in the complex number system, the root equation z n = a, where a ∈ C\{0}
and n ∈ N, always has n solutions! Let’s first look at the simplest equation z n = 1:
Certainly, 1 is a solution to the equation. Furthermore, using De Moivre’s Theorem,
we get:
n
2π 2π 2π 2π
cos + i sin = cos · n + i sin · n = cos(2π) + i sin(2π) = 1.
n n n n
Clearly, this shows the complex number cos 2π 2π n
n + i sin n satisfies the equation z = 1. In
2kπ 2kπ
fact, any number which can be expressed in form of cos n + i sin n , where k is an
8 1. Preliminaries
Definition 1.9 (Roots of a Complex Number). Given any a ∈ C\{0} and n ∈ N, the
n-th roots of a is a set given by:
1 p Arg(a) + 2kπ Arg(a) + 2kπ
a n = n |a| cos + i sin : k ∈ {0, 1, . . . , n − 1}
n n
1 √ 1
Example 1.3. Find i 3 and (1 − 3i) 2 .
Solution
First express i into polar form i = cos π2 + i sin π2 . Hence by Definition 1.9, we have:
π π
1 + 2kπ + 2kπ
i 3 = cos 2 + i sin 2 : k = 0, 1, 2
3 3
π π 5π 5π 3π 3π
= cos + i sin , cos + i sin , cos + i sin
| 6 {z 6} | 6 {z 6} | 2 {z 2}
k=0 k=1 k=2
(√ √ )
3+i 3−i
= , , −i
2 2
n √ 1o
Similarly, to find (1 − 3i) 2 , we first express:
√ π π
1 − 3i = 2 cos − + i sin −
3 3
Hence, by Definition 1.9, we have:
π π
√ 1 √ − 3 + 2kπ − 3 + 2kπ
(1 − 3i) = 2 2 cos + i sin : k = 0, 1
2 2
( √ ! √ !)
√ 3−i √ − 3+i
= 2 , 2
2 2
(√ √ )
3−i − 3+i
= √ , √
2 2
where n is an integer with n ≥ 2, then one can show the following identity holds:
1 + ω + ω 2 + . . . + ω n−1 = 0.
(1 − ω)(1 + ω + ω 2 + . . . + ω n−1 )
= (1 + ω + ω 2 + . . . + ω n−1 ) − ω(1 + ω + ω 2 + . . . + ω n−1 )
= (1 + ω + ω 2 + . . . + ω n−1 ) − (ω + ω 2 + . . . + ω n−1 + ω n )
= 1 − ωn
n
2π 2π
= 1 − cos + i sin
n n
= 1 − (cos(2π) + i sin(2π)) = 1 − 1 = 0.
Since ω 6= 1 as n ≥ 1, we conclude that:
1 + ω + ω 2 + . . . + ω n−1 = 0.
Using this result, one can derive some trigonometric identities. Express ω in terms of its
real and imaginary parts:
2 n−1
2π 2π 2π 2π 2π 2π
1 + cos + i sin + cos + i sin + . . . + cos + i sin =0
n n n n n n
| {z }
ω
2π 2π 4π 4π
1 + cos + i sin + cos + i sin + ...
n n n n
2(n − 1)π 2(n − 1)π
+ cos + i sin =0
n n
By equating the real and imaginary parts, we obtain two trigonometric identities:
2π 4π 2(n − 1)π
cos + cos + . . . + cos = −1
n n n
2π 4π 2(n − 1)π
sin + sin + . . . + sin =0
n n n
Exercise 1.14. Let Pk (xk , yk ), where k = 1, 2, 3, be three distinct points in C and let
zk := xk + yk i be the complex number representing Pk . Denote ω = cos 2π 2π
3 + i sin 3 .
Show that 4P1 P2 P3 is equilateral if and only if
z1 + ωz2 + ω 2 z3 = 0.
Using this, show that it is impossible for 4P1 P2 P3 being equilateral if xk , yk ∈ Q for
all k = 1, 2, 3.
12 1. Preliminaries
Remark 1.13. The definition of convergence of complex sequences is almost the same
as the that of real sequences. The only difference is now |·| represents the modulus
while for real sequence it represents the absolute value. Therefore, many computational
rules about limits carry over to complex sequences. For instance, if lim zn = L and
n→∞
lim wn = M , then we have
n→∞
lim (zn ± wn ) = L ± M
n→∞
lim (zn wn ) = LM
n→∞
zn L
lim = (whenever M 6= 0)
n→∞ wn M
Solution
First consider the case |z| < 1: if z = 0, then zn = 0 for any n and the desired result
clearly holds. From now on we assume z 6= 0. For any ε > 0, we pick a positive
log ε
integer N > log|z| . Whenever, n ≥ N , we have:
n N
|zn − 0| = |z n | = |z| ≤ |z| .
Here we have used the fact that |z| < 1 and n ≥ N . By our choice of N , we have:
log ε
N
|z| < |z| log|z|
log|z| ε
= |z| = ε.
This shows lim zn = 0 in case of |z| < 1. The case of z = 1 is trivial.
n→0
1.2. Sequences and Series 13
When |z| ≥ 1 and z 6= 1, the sequence zn = z n can be shown to diverge using the
squeezing principle (see Exercise 1.16). It can also be proved using the following useful
fact:
Now given a complex number z expressed in polar form as z = r(cos θ + i sin θ), and
suppose |z| ≥ 1 (i.e. r ≥ 1) and z 6= 1. Consider again the sequence zn = z n . By De
Moivre’s Theorem, we have:
zn = rn (cos nθ + i sin nθ).
It is well known in real analysis that when θ 6= 2kπ (where k ∈ Z), at least one of the real
sequences {cos nθ} and {sin nθ} diverges as n → ∞. Hence, when r ≥ 1 and θ 6= 2kπ
(k ∈ Z), at least one of the real sequences {rn cos nθ} and {rn sin nθ} diverges. This
shows zn diverges.
Since z n → 0 when |z| < 1, we also have the geometric series formula for complex
numbers:
a
a + az + az 2 + az 3 + · · · = .
1−z
when |z| < 1 just like the real case.
Exercise 1.15. Show that if lim zn = L, then lim zn = L and lim |zn | = |L|.
n→∞ n→∞ n→∞
Exercise 1.16. Show (without using Proposition 1.14) that if |z| ≥ 1 and z 6= 1,
then the sequence {z n } must diverge. [Hint: First prove the following inequality:
|z − 1| ≤ z n+1 − w + |z n − w|
for any z ∈ C such that |z| ≥ 1, and any w ∈ C.]
In Real Analysis, there is a notion of Cauchy sequences which describe sequences that
are closer and closer to each other. It is a priori different from convergent sequences, which
are sequences that are closer and closer to a certain limit. However, it is well-known that
for sequences in R, the Cauchy condition will guarantee convergence. This important
fact is known as completeness of real numbers.
In Complex Analysis, we have a similar notion of Cauchy sequences and completeness,
to be discussed below.
Definition 1.15 (Cauchy Sequence). A sequence {zn }∞ n=1 of complex numbers is said
to be a Cauchy sequence if and only if for any ε > 0, there exists an integer N ∈ N such
that whenever m, n ≥ N , we have |zn − zm | < ε.
14 1. Preliminaries
Theorem 1.16 (Completeness of C). Every Cauchy sequence of complex numbers con-
verges to a certain complex number. In other words, C is complete.
Proof. Let {zn } be a Cauchy sequence of complex numbers. We need to show it converges.
Write zn = xn + iyn , where xn , yn ∈ R. Since we have:
|xn − xm | ≤ |zn − zm |
|yn − ym | ≤ |zn − zm |
and given that {zn } is a Cauchy sequence, the real sequences {xn } and {yn } are also
Cauchy sequences. By Completeness of R, both {xn } and {yn } converge to some real
numbers x∞ and y∞ respectively. By Proposition 1.14, the complex sequence {zn }
converges to x∞ + iy∞ .
Exercise 1.17. Suppose {zn }∞ n=0 is a complex sequence. Suppose there exists a real
constant α ∈ [0, 1) such that:
|zn+1 − zn | ≤ α |zn − zn−1 | for any n ∈ N.
Show that the complex sequence {zn }∞
n=0 converges.
∞
X
1.2.2. Series of Complex Numbers. An (infinite) series zn of complex num-
n=1
N
X
bers zn ∈ C is the limit (if exists) of the N -th partial sums zn as N → ∞. In Real
n=1
Analysis, we learned that many series convergence tests rely on the fact that R is com-
plete. Now that we know C is also complete (Theorem 1.16), we can generalize many
series convergence tests for C.
∞
X
Proposition 1.18 (Absolute Convergence Test). If the series |zn | converges, then the
n=1
∞
X
complex series zn also converges.
n=1
Proof. The proof is almost identical to the real case, mutatis mutandis. Given that
X∞ N
X
|zn | converges, its N -th partial sum |zn | is a Cauchy sequence. Now consider
n=1 n=1
N
X
the sequence of N -th partial sums zn . We want to show the later is also a Cauchy
n=1
sequence.
For any ε > 0, there exists an integer K > 0 such that whenever M > N ≥ K, we
have
XM N
X
|zn | − |zn | < ε.
n=1 n=1
1.2. Sequences and Series 15
It implies:
M
X N
X M
X M
X M
X N
X
zn − zn = zn ≤ |zn | = |zn | − |zn | < ε.
n=1 n=1 n=N +1 n=N +1 n=1 n=1
N
X
Therefore, zn is also a Cauchy sequence. By completeness of C (Theorem 1.16), the
n=1
N
X ∞
X
N -th partial sum zn (and hence the infinite series zn ) converges.
n=1 n=1
X∞ n
i
Example 1.5. Does the series converge absolutely, conditionally, or does not
n=1
n
X∞
in
converge? How about the series ?
n=1
n2
Solution
X ∞ X∞
in 1
The series = diverges by p-test. The N -th partial sum can be decom-
n=1
n n=1
n
posed into:
k
k−1
XN
in − 1 + 1 − . . . + (−1) + 1 − 1 + 1 + . . . + (−1) i if N = 2k
2 4 2k 3 5 2k−1
= k
k+1
n − 1 + 1 − . . . + (−1) + 1 − 1 + 1 + . . . + (−1) i if N = 2k + 1
n=1 2 4 2k 3 5 2k+1
In either case, the real and imaginary parts converge by alternating series test. By
X∞ n
i
Proposition 1.14, the series converges, and so it converges conditionally.
n=1
n
X∞ X∞ X∞
in in 1
Now consider 2
. The series 2
= converges by p-test. There-
n=1
n n=1
n n=1
n2
X∞
in
fore, the series converges absolutely.
n=1
n2
One good property of an absolute convergent series is that we can rearrange the
terms as we wish without changing the value of the series. Precisely, given an absolute
∞
X
convergent series zn =: L and a bijection σ : N → N, then the rearranged series
n=1
∞
X
zσ(n) also converges absolutely to the limit L. The proof in the complex case is
n=1
∞
X
based on the analogous result in the real case: if the complex series zn converges
n=1
∞
X ∞
X
absolutely, then the real series Re(zn ) and Im(zn ) both converge absolutely by
n=1 n=1
the comparison test (as |Re(zn )| ≤ |zn | and |Im(zn )| ≤ |zn |). It then follows from the
16 1. Preliminaries
∞
X ∞
X
rearrangement of real absolute convergent series that Re(zσ(n) ) = Re(zn ) and
n=1 n=1
∞
X ∞
X X∞ ∞
X
Im(zσ(n) ) = Im(zn ). By Proposition 1.14, it proves zσ(n) = zn .
n=1 n=1 n=1 n=1
Recall from Real Analysis that the ratio test and root test follow from the absolute
convergence test and geometric series test. Now we learned that both hold on C, hence
the ratio test and root test can be extended to complex series:
∞
X
Proposition 1.19 (Ratio Test). Consider the complex series zn :
n=1
X∞
zn+1
• If lim < 1, then zn converges absolutely.
n→∞ zn n=1
X∞
zn+1
• If lim > 1, then zn diverges.
n→∞ zn n=1
zn+1
• If lim = 1, then no conclusion can be drawn.
n→∞ zn
∞
X
Proposition 1.20 (Root Test). Consider the complex series zn :
n=1
p ∞
X
n
• If lim |zn | < 1, then zn converges absolutely.
n→∞
n=1
p ∞
X
• If lim n |zn | > 1, then zn diverges.
n→∞
n=1
p
n
• If lim |zn | = 1, then no conclusion can be drawn.
n→∞
Remark 1.21. The proofs of the ratio and root tests are the same as in the real case. We
omit their proofs but we encourage readers to write down their proofs as an exercise.
X∞
zn
Example 1.6. Show that for any z ∈ C, the complex series converges
n=0
n!
absolutely.
Solution
We use the ratio test. Consider:
z n+1 /(n + 1)! z n+1 n!
lim n
= lim n
n→∞ z /n! n→∞ z (n + 1)!
z |z|
= lim = lim
n→∞ n + 1 n→∞ n + 1
Alternatively, we can also use the root test (Proposition 1.20) by showing that:
s
zn |z|
lim n = lim √ =0<1
n→∞ n! n→∞ n n!
√n
for any z ∈ C. Here we have used the fact that lim n! = ∞.
n→∞
∞
X
Example 1.7. Determine all complex numbers z such that the series nz n con-
n=0
verges.
Solution
(n + 1)z n+1 (n + 1)
Consider the limit lim n
= lim |z| = |z|. Therefore, by ratio
n→∞ nz n→∞ n
test (Proposition 1.19), the series converges absolutely when |z| < 1; and diverges
when |z| > 1.
When |z| = 1, the ratio test fails to conclude anything. In this case, we let
∞
X
z = cos θ + i sin θ where θ ∈ R. Then, the series is given by (n cos nθ + in sin nθ),
n=0
and the real and imaginary parts are:
∞
! ∞ ∞
! ∞
X X X X
n n
Re nz = n cos(nθ) and Im nz = n sin(nθ).
n=0 n=0 n=0 n=0
By Proposition 1.14, if the complex series converges, then both their real and
imaginary parts converge, and in particular we have:
lim n cos(nθ) = 0 and lim n sin(nθ) = 0.
n→∞ n→∞
By Squeeze Theorem, it will imply:
lim cos(nθ) = lim sin(nθ) = 0.
n→∞ n→∞
However, it would contradict the fact that cos2 (nθ) + sin2 (nθ) = 1; and so the series
∞
X
nz n does not converge when |z| = 1.
n=0
∞
X
Conclusion: the series nz n converges if and only if |z| < 1.
n=0
Below is a very nice example that demonstrates complex numbers could come into
play to solve a problem (about integers) which seems to have nothing to do with complex
numbers.
Solution
What is needed to show is that for any integer k ≥ 2, if ai , di ∈ N where 1 ≤
i ≤ k are positive integers such that whenever i 6= j, we have di 6= dj . Denote
S1 , S2 , . . . , Sk are the sets of infinite arithmetic sequences:
S1 = {a1 , a1 + d1 , a1 + 2d1 , . . .}
..
.
Sk = {ak , ak + dk , ak + 2dk , . . .}
We need to show Si ∩ Sj = ∅ for any i 6= j, then:
k
[
N 6= Si .
i=1
X
For each 1 ≤ i ≤ k, we first express the sum z n where |z| < 1 in terms of ai , di
n∈Si
and z with no summation sign.
X
z n = z ai + z ai +di + z ai +2di + · · ·
n∈Si
= z ai 1 + z di + (z di )2 + · · ·
z ai
=
1 − z di
di
using the geometric series formula (note that the common ratio z di = |z| < 1).
Now suppose on the contrary N = {1, 2, 3, · · · } can be decomposed into disjoint
union of Si ’s, then X X X
zn = zn + · · · + zn
n∈N n∈S1 n∈Sk
for any |z| < 1. Here we have used the absolute convergence of geometric series,
and the rearrangement theorem of absolute convergent series. However, it will
imply:
z z a1 z ak
= + · · · + .
1−z 1 − z d1 1 − z dk
Assume without loss of generality that d1 < d2 < · · · < dk so that dk is the largest
among all d’s. Letting z → cos 2π 2π
dk + i sin dk which is an dk -th root of 1, then for any
j = 1, 2, · · · , k − 1:
z cos 2π 2π
dk + i sin dk
lim =
z→cos 2π
dk +i sin 2π
dk
1−z 1 − (cos 2π 2π
dk + i sin dk )
2a π 2a π
z aj cos djk + i sin djk
lim = 2a π 2a π .
z→cos 2π
dk +i sin 2π
dk
1 − z dj 1 − (cos djk + i sin djk )
Note that the both limits exist because the denominators do not go to zero.
However, for the last term, we will get:
z ak cos 2adkkπ + i sin 2adkkπ
lim = which does not exist.
z→cos 2π
dk +i sin 2π
dk
1 − z dk 1−1
This gives a contradiction. Hence N cannot be decomposed into a finite disjoint
union of arithmetic sequences of distinct common differences.
1.2. Sequences and Series 19
Exercise 1.18. Determine whether each of the following complex series converges
absolutely, conditionally, or diverge:
X∞
(1 − 3i)n
(a)
n=0
(4 + i)2n
∞
X n2
(b)
n=1
n + n3 i
∞
X
(c) (cos n − i sin n)
n=1
Exercise 1.19. In each of the following complex series: (i) determine all complex
numbers z such that the series converges, (ii) sketch the range of these z’s on the
complex plane C.
∞
X
(a) zn
n=1
X∞ n
z
(b)
n=1
z+1
X∞
(−1)n z 5n
(c)
n=1
n!
X∞
z n!
(d)
n=1
n2
∞
zn
(b) Using (a), or otherwise, find all z ∈ C such that the sequence
1 + z 2n n=1
converges.
∞
X zn
(c) Find all z ∈ C such that the series converges.
n=1
1 + z 2n
X∞
zn
1.2.3. Euler’s Identity. The series considered in Example 1.6 is an important
n=0
n!
one – it defines the complex exponential function. When z = x is a real number, the value
of the series is given by ex . Given that the series converges for any z ∈ C, we define ez to
be the limit of this series:
Remark 1.23. Some books (especially those about physics and engineering) define
ex+yi = ex (cos y + i sin y), but here we will first define ez as an infinite series first, and
justify that ex+yi = ex (cos y + i sin y) from the series definition of ez .
Here is the famous Euler’s identity that relates complex exponentials with the polar
form of a complex number:
Proof. The key idea is to split the defining series into real and imaginary parts.
X∞ 2N n n
X
(iθ)n i θ
eiθ = = lim
n=0
n! N →∞
n=0
n!
N N −1 2k+1 2k+1
!
X i2k θ2k X i θ
= lim + [by rearrangement]
N →∞ (2k)! (2k + 1)!
k=0 k=0
N
! N −1
X (−1)k θ2k (−1)k θ2k+1 X
= lim +i lim [using i2k = (i2 )k = (−1)k ]
N →∞ (2k)! N →∞ (2k + 1)!
k=0 k=0
| {z } | {z }
=cos θ =sin θ
= cos θ + i sin θ
Remark 1.25. From (1.3), it is evident that we have:
eiπ + 1 = 0
which is a single identity involving 5 most important constants in mathematics, namely 1,
0, e, π and i.
Remark 1.26. From the Euler’s identity, we can now write down the polar form of a
complex number in a simpler way: if z = r(cos θ + i sin θ), then we can write:
z = reiθ .
In particular, any z ∈ C such that |z| = 1 can be expressed as z = eiθ for some θ ∈ R.
We are going to show that the complex exponential has the property that ez ew = ez+w
just like the real case. Informally, we express both ez and ew into two infinite series. After
multiplying the two series, we express the double sum diagonally:
∞
! ∞ ! ∞ X ∞
X zn X wm X z n wm
z w
e e = =
n=0
n! m=0
m! n=0 m=0
n!m!
∞
X X ∞
X X z n wk−n k
z n wm
= = [since m = k − n]
n!m! n=0
n!(k − n)!
k=0 m+n=k k=0
∞ X
X k
Cnk z n wk−n k!
= [since Cnk = ]
n=0
k! (n − k)!n!
k=0
X∞
(z + w)k
= = ez+w [Binomial Theorem]
k!
k=0
Although this “proof” above seems convincing and neat, there isPa little
P step we need
to justify, namely why we can rearrange the infinite double sum n m in a diagonal
1.2. Sequences and Series 21
P P P P
way: k m+n=k ? We have seen in Real Analysis that even switching n and m may
sometimes result in a different sum. Below we give a rigorous (and more refined) proof
of this fact:
XN XN
zn wm
Proof. Consider the N -th partial sums and , then:
n=0
n! m=0
m!
N
! N
! N X N
X zn X wm X z n wm
=
n=0
n! m=0
m! n=0 m=0
n!m!
| {z }
Region I in Fig. 1.2
2N
X X N
X −1 2N
X −m N −1 2N −n
z w n m
z n wm X X z n wm
= − −
n!m! m=0
n!m! n=0
n!m!
k=0 m+n=k n=N +1 m=N +1
| {z } | {z } | {z }
Region I+II+III in Fig. 1.2 Region II in Fig. 1.2 Region III in Fig. 1.2
PN PN
Here we break down the finite double sum n m into three triangular sums. See
Figure 1.2 for illustration. For the sum corresponding to the large triangle (Region
I+II+III in Figure 1.2), we can rewrite it as:
2N
X X 2N k
X X z n wk−n X X C k z n wk−n 2N k
X (z + w)k 2N
z n wm n
= = = → ez+w
n!m! n=0
n!(k − n)! n=0
k! k!
k=0 m+n=k k=0 k=0 k=0
as N → ∞.
For Region II in Figure 1.2, we can show that it converges to 0 as N → ∞:
N
X −1 2N
X −m N
X −1 2N
X −m
z n wm z n wm
≤
m=0 n=N +1
n!m! m=0
n!m!
n=N +1
N
X −1
X∞ n m
|z| |w|
≤
m=0 n=N +1
n! m!
N −1
! ∞
!
X |w|
m X |z|
n
=
m=0
m! n!
n=N +1
∞
!
X |z|
n
≤ e|w| →0
n!
n=N +1
∞
X n
|z|
as N → ∞ since the infinite sum converges (to e|z| ). The sum corresponding to
n=0
n!
Region III in Figure 1.2 can be shown to converge to 0 by switching m and n, and z and
w in the above argument.
Overall, we have shown:
N
! N !
X zn X wm
lim
N →∞
n=0
n! m=0
m!
2N
X X N
X −1 2N
X −m N
X −1 2N
X −n
z n wm z n wm z n wm
= lim − lim − lim
N →∞ n!m! N →∞
m=0
n!m! N →∞
n=0
n!m!
k=0 m+n=k n=N +1 m=N +1
z+w
=e − 0 − 0,
22 1. Preliminaries
2N
II
I III
m
O N 2N
Figure 1.2
∞
X ∞
X
Exercise 1.21. Given two series zn and wn which converge absolutely to A
n=0 n=0
and B respectively, show that the series below converges absolutely to AB:
∞ k
!
X X
zn wk−n
k=0 n=0
X∞
zn
Exercise 1.23. Using the result from Exercise 1.22, show that the series
n=1
n
converges for any z such that |z| = 1 and z 6= 1.
Using the multiplicative property ez ew = ez+w , one can show the following properties
about the complex exponential function. We leave the proofs for readers.
Remark 1.28. For any z = x + yi ∈ C where x, y ∈ R, we have:
n
• (ez ) = enz for any integer n.
1.2. Sequences and Series 23
The complex exponential az with other real base a > 0 is defined via the natural
exponential ez . Recall that a = eln a , and we define:
az := e(ln a)z .
Remark 1.30. For any positive integer n, the rational number n1 is no doubt also a
1
complex number. Therefore, now e n could mean two different things, namely the value
∞ k
X 1
n
of the series , or the n-th roots of e. It is a confusing ambiguity but fortunately
k!
k=0
we seldom deal with both of them in the same context. One way to avoid such a confusion
1
is to represent the n-th roots of e by e n , and use exp( n1 ) to represent the value of the
aforesaid series.
1.2.4. Riemann ζ Function: the first encounter. The Riemann zeta function, de-
noted by ζ(z), is of central importance in Complex Analysis and Number Theory. It is an
infinite series defined as:
X∞
1
ζ(z) := z
n=1
n
for Re(z) > 1. This complex series motivates the discussions of the famous Riemann
Hypothesis, which is a conjecture purposed by Riemann in 1859 and remains unsolved
as of today (February 29, 2020). The statement of the Riemann Hypothesis will be
explained after we learn about analytic continuation of holomorphic functions. The
Riemann zeta function has deep connections with Number Theory, in particular on the
study of distribution of prime numbers. It is used to show the renowned Prime Number
Theorem, which asserts that:
π(x)
lim =1
x→∞ x/ ln x
where π(x) is the number of positive prime numbers less than or equal to x. The
deep connection between ζ(z) and prime numbers is beyond the scope of this course.
Interested readers may consult Stein-Shakarchi’s book.
Meanwhile, we show that this series converges absolutely when Re(z) > 1 by the
(real) p-test. The main reason is as follows. Write z = x + yi where x, y ∈ R, then we
have:
1 1 1 1 1
= z log n = x log n iy log n = x = Re(z) .
nz e |e e | n n
24 1. Preliminaries
∞
X1
By (real) p-test, the series Re(z)
converges if and only if Re(z) > 1. Therefore, by
n=1
n
X∞
1
the (complex) absolute convergence test, the series converges absolutely when
n=1
nz
Re(z) > 1.
1.3. Point-Set Topology of C 25
1.3.1. Open and Closed Subsets. Intuitively, an open subset Ω in C is one that
does not have a boundary. However, this “definition” is not rigorous enough since the
term “boundary” has not been defined so far. We are going to give a rigorous definition
of open and closed subsets here. We first define:
In the figure below, the yellow set is the subset U ⊂ C. The point z ∈ U is an interior
ball because by drawing a ball with a small enough radius (i.e. the blue ball), the ball is
completely inside U . In layman terms, an interior point of U is a point z whose “nearby”
points are contained in U .
On the other hand, the point w in the figure below is a boundary point. No matter
how small the ball you draw around w, that ball must contain some points in U , and
some points not in U . In layman terms, a boundary point of U is a point w at which if
you look around it, you can see “nearby” some points in U and some point not in U .
The point η in the figure is an exterior point of U . In layman terms, it is a point
whose “nearby” are outside U .
η ε
δ z
w
26 1. Preliminaries
Remark 1.32. Since z ∈ Bε (z) for any z ∈ C and ε > 0, if z is an interior point of U , it
is automatically that z ∈ U . In other words, an interior point of a set must belong to that
set. However, a boundary point of a set can be contained or not contained in the set.
Furthermore, according to the definitions, interior points, boundary points and exterior
points are mutually exclusive.
Example 1.9. Find all interior, boundary and exterior points of the set:
U = {z ∈ C : Re(z) > 1}.
Solution
We claim that the set of interior points is U itself. For any z ∈ U , we have Re(z) > 1.
Write z = x + yi, then we have x > 1. We need to find a small ε > 0 such that
Bε (x + yi) ⊂ U . According to the figure below, an appropriate choice of ε should
be ε = x−1
2 . We next verify that it is indeed Bε (z) ⊂ U for this choice of ε.
x−1
For any α ∈ Bε (z), we have |α − z| < ε = . Then, by Re(z − α) ≤ |z − α|,
2
we know that:
x−1 x−1
Re(z − α) < =⇒ x − Re(α) < .
2 2
x−1 x+1 1+1
By rearrangement, we get Re(α) > x − = > = 1, which is equiv-
2 2 2
alently to saying that α ∈ U . This shows Bε (z) ⊂ U , and hence z is an interior
point.
η z
1 x
Exercise 1.24. Find all the interior, boundary and exterior points of each set below:
(a) U1 = {z ∈ C : Re(z) ≥ 0 and Im(z) < 0}.
(b) U2 = Br (z0 ) where z0 ∈ C is a fixed number and r > 0
(c) U3 = Br (z0 ) where z0 ∈ C is a fixed number and r > 0.
(d) U4 = ∂Br (z0 ) where z0 ∈ C is a fixed number and r > 0.
(e) U5 = C.
From now on, given any set U ⊂ C, we denote and define:
U c := C\U = the complement of U in C
U ◦ := set of all interior points of U
∂U := set of all boundary points of U
U := U ∪ ∂U = the closure of U
There is no standard notation for the set of all exterior points though. According to the
definition of interior points, we must have U ◦ ⊂ U .
We are now ready to define what are open sets and closed sets. The way we define
open sets is very common in many other textbooks, while the way we define closed sets
may sound different from some textbooks but it is more intuitive and is nonetheless
equivalent to the definition found in other textbooks.
Definition 1.33 (Open and Closed Sets). A set Ω ⊂ C is said to be open if every point
z ∈ Ω is an interior point of Ω (i.e. Ω = Ω◦ ). A set E ⊂ C is said to be closed if all
boundary points of E belong to E (i.e. ∂E ⊂ E).
Let’s look at some examples. Consider the set Ω = {z ∈ C : Re(z) > 1}:
Ω◦ = {z ∈ C : Re(z) > 1} = Ω
∂Ω = {z ∈ C : Re(z) = 1} 6⊂ Ω
Therefore, Ω is an open set, but is not closed.
Let’s look at another example: E = {z ∈ C : Re(z) ≥ 1}. By inspection (we left the
detail for readers), we can see:
E ◦ = {z ∈ C : Re(z) > 1} =
6 E
∂E = {z ∈ C : Re(z) = 1} ⊂ E
Therefore, E is not open, but is closed.
There are sets which are not open and not closed! For instance, consider the unit
circle W = {z ∈ C : Re(z) ≥ 0 and Im(z) > 0}. We can see from Figure 1.3 that:
W ◦ = {z ∈ C : Re(z) > 0 and Im(z) > 0} =
6 W
∂W = {x + 0i ∈ C : x ≥ 0} ∪ {0 + yi ∈ C : y ≥ 0} 6⊂ W.
W is neither open nor closed.
28 1. Preliminaries
Figure 1.3. The set W = {z ∈ C : Re(z) ≥ 0 and Im(z) > 0} and its interior and boundary sets
Surprisingly, there are sets which are both open and closed (so “open” and “closed”
are not exactly opposite, which is a linguistic nightmare)! For subsets of C, there are not
many though. They are the empty set ∅ and the whole C. It is easy to see that C◦ = C
and ∂C = ∅ ⊂ C (the empty-set is a subset of every set). This shows C is both open and
closed.
The argument that shows ∅ is both open and closed has a bit of philosophical favor.
We claim that ∅◦ = ∅. Suppose otherwise, then we must have ∅◦ 6⊂ ∅ (since ∅ is a subset
of every set). This means there exists z ∈ ∅◦ such that z 6∈ ∅. Then, z being an interior
point of ∅ implies there exists ε > 0 such that Bε (z) ⊂ ∅, which is clearly impossible!
This shows ∅◦ = ∅ and so the empty set is open. To show ∅ is closed as well, we claim
∂∅ = ∅. Suppose ∂∅ is non-empty, then we can pick w ∈ ∂∅, then for any δ > 0, both
Bδ (w) ∩ ∅ and Bδ (w) ∩ (C\∅) are non-empty. However, the former cannot happen! This
concludes ∂∅ = ∅, and so ∅ is closed as well!
Remark 1.35. There is an interesting YouTube video titled “Hitler learns Topology”.
Exercise 1.25. Determine whether each set U1 to U5 in Exercise 1.24 is open, closed,
neither or both.
Readers who have learned a bit point-set topology may have seen another definition
of closed sets, namely a set E is closed if its complement E c is open. We are going to
show that this is equivalent to our definition:
Therefore, from now on we can say a set is closed if and only if its complement is
open, which is more convenient sometimes. For instance, this fact can be used to show
an important and nice property about a closed set E: if there is a convergent sequence in
E, then the limit must be inside E.
Proof. We prove by contradiction. The key idea is that if w 6∈ E, then one can draw a
small ball around w such that the ball is completely outside E. However, then zn which
approaches w must go within the ball, and hence outside E, when n is large (see Figure
1.4).
Here we present the detail: suppose w 6∈ E, then w ∈ E c . By Proposition 1.36, E c
is open and so there exists ε > 0 such that Bε (w) ⊂ E c . By the fact that zn → w, there
exists N ∈ N such that whenever n ≥ N , we have |zn − w| < ε. However, it implies:
zn ∈ Bε (w) ⊂ E c =⇒ zn 6∈ E
which is clearly a contradiction. It proves w ∈ E.
w
zn
z4
z3 E
z2 z1
Below is a list of useful facts about open and closed sets. We will prove some of them
and leave the others as exercises for readers.
Proposition 1.38. Open and closed sets in C have the following properties:
[
• The union Uα of any family (finite or infinite) of open sets {Uα } in C is open.
α
N
\
• The intersection Uk of a finite family of open sets U1 , . . . , UN in C is open.
k=1
N
[
• The union Ek of a finite family of closed sets E1 , . . . , EN in C is closed.
k=1
\
• The intersection Eα of any family (finite or infinite) of closed sets {Eα } in C is
α
closed.
Proof. Let’s prove the second statement only, that if U1 , . . . , UN are open, then their
N
\
intersection is also open. Let z ∈ Uk , then z ∈ Uk for any k = 1, . . . , N . For each
k=1
k, since Uk is open, z is an interior point of Uk and so there exists εk > 0 such that
Bεk (z) ⊂ Uk . Let ε = min{ε1 , . . . , εN }, which is positive, then ε ≤ εk for any k, and so
we have:
Bε (z) ⊂ Bεk (z) ⊂ Uk for any k = 1, . . . , N .
N
\ N
\ N
\
Therefore, Bε (z) ⊂ Uk . This shows z is an interior point of Uk . As a result, Uk
k=1 k=1 k=1
is an open set.
30 1. Preliminaries
We leave the proof of the first statement as an exercise for readers. Once the
first two statements are established, the third and fourth statements about
!c closed sets
[ \
easily follow from Proposition 1.36 and De Morgan’s Rule: Ek = Ekc and
!c k k
\ [
Eα = Eαc .
α α
Exercise 1.26. Prove all the other three statements in Proposition 1.38. Give an
example of a family of open sets whose intersection is not open. Also give an
example of a family of closed sets whose union is not closed.
Exercise 1.28. Use the Bolzano-Weierstrass’s Theorem for R to show the Bolzano-
Weierstrass’s Theorem for C, which asserts that if {zn }∞
n=1 is a complex sequence in
a bounded set Ω, then there exists a convergent subsequence {znk }∞ ∞
k=1 of {zn }n=1 .
1.3.2. Connected Sets. Intuitively, a connected set is one that is in one “piece”.
However, such a definition is not rigorous as the word “piece” is quite vague. To define
connectedness, we first need to understand what it means by a disconnected set:
Remark 1.40. The condition Ω ⊂ U ∪ V means that U and V together cover the whole
set Ω. The condition Ω ∩ U and Ω ∩ V being non-empty means that Ω has something
inside U and something inside V . Since the definition requires U and V are disjoint (i.e.
separated in some sense), these sets U and V create a separation for the set Ω, and hence
we say Ω is disconnected (see Figure 1.5).
Figure 1.5. Ω is the yellow set. It is disconnected with disjoint open sets U and V that
separate Ω.
an annulus 1 < |z| < 2. However, thanks for a proposition that we will state, one can
verify that they are all connected easily. Before we state the proposition, we need to
define:
For instance, any convex set is polygonally path-connected since every pair of points
can be joined by a single line segment contained inside the set. The annulus 1 < |z| < 2
is also polygonally path-connected (see the figure below):
z1
z2
The following proposition asserts that for any open set Ω, connectedness and polygonal-
path-connectedness are equivalent:
We omit the proof in this lecture note. Interested readers may consult Stein-
Shakarchi’s book (Exercise 5 in P.25) for an outline of the proof and try to complete the
detail as an exercise. Using this proposition, it is easy to see that any convex open sets
(and many other non-convex open sets) are connected.
The last notion about sets in C to be introduced is simply-connectedness. Readers
should have encountered this concept in Multivariable Calculus (typically in the chapter
about conservative vector field).
On C, simply-connected sets have one nice property concerning simple closed curves
(“simple closed” means closed without self-intersections). If γ is a simple closed curve
contained inside a simply-connected set Ω, then the region enclosed by γ will be a subset
of Ω. Some textbooks put this as the definition of simply-connected sets in C.
Exercise 1.30. For each set described below, sketch the region on C, and determine
whether it is (i) open, (ii) closed, (iii) bounded, (iv) compact, (v) connected and
(vi) simply-connected or not.
(a) Ω1 = {z ∈ C : |z + 1| ≥ 4 |z − 1|}
(b) Ω2 = {z ∈ C : |z + 1| < 4 |z − 1|}
(c) Ω3 = {z ∈ C : |z| ≤ Re(z) + 1}
(d) Ω4 = {ez ∈ C : 1 ≤ Re(z) ≤ 2}
(e) Ω5 = {z ∈ C : z 2 − 1 ≤ 1}
Chapter 2
Holomorphic Functions
33
34 2. Holomorphic Functions
Exercise 2.1. For each function below, state its domain and find its real and imagi-
nary parts:
1
(a) f (z) =
z
2
(b) f (z) = |z|
(c) f (z) = e2z
iz + 1
(d) f (z) =
z−i
The orange graph represents the real part u(x, y) = x2 − y 2 , whereas the blue graph
represents the imaginary part v(x, y) = 2xy.
Similarly, the exponential function ez = ex cos y + iex sin y can be visualized as two
separate graphs
Again the orange graph is the real part, and the blue graph is the imaginary part.
1 x yi
Some functions such as f (z) = = 2 − 2 are not defined everywhere
z x + y2 x + y2
on C. Let’s see how its graphs look:
2.1. Complex-Valued Functions 35
From the graph, one can see easily that both real and imaginary parts tend to ±∞ as
(x, y) approaches (0, 0).
5
0.5 0.5
0.0 0.0
3
-0.5 -0.5
-1.0 0 -1.0
-1.0 -0.5 0.0 0.5 1.0 0 1 2 3 4 5 6 -1.0 -0.5 0.0 0.5 1.0
One can see the orange and blue level curves intersect each other orthogonally at
almost all points. It is in fact not coincident! This orthogonality phenomenon is related
to complex differentiability as we will see in the next section.
1.0 2.0
0.8
1.5
0.6
1.0
0.4
0.2 0.5
0.0
0.0 0.2 0.4 0.6 0.8 1.0 - 1.0 - 0.5 0.5 1.0
1
Example 2.4. Consider another example f (z) = z + . We want to find the image
z
of the circle |z| = r0 under this map. Write z = r0 eiθ where 0 ≤ θ ≤ 2π, then we
have:
1
f (z) = r0 eiθ +
r0 eiθ
1
= r0 eiθ + e−iθ
r0
1
= r0 (cos θ + i sin θ) + (cos θ − i sin θ)
r
0
1 1
= r0 + cos θ +i r0 − sin θ,
r0 r0
| {z } | {z }
u v
which gives the following ellipse in the uv-plane (when r0 6= 1):
u2 v2
2 + 2 = 1.
1 1
r0 + r0 r0 − r0
3 3
2 2
1 1
0 0
-1 -1
-2 -2
-3 -3
-3 -2 -1 0 1 2 3 -3 -2 -1 0 1 2 3
Exercise 2.2. Describe and sketch the image of the semi-infinite strip:
Σ = {z ∈ C : 0 ≤ Re(z) ≤ a and Im(z) ≥ 0}
under the map f (z) = z 2 . Here a is a positive real number.
1 Arg(z) 1 Arg(z)+2π 1 Arg(z)+4π
(a) z 7→ Re |z| 3 e 3 i (b) z 7→ Re |z| 3 e 3
i
(c) z 7→ Re |z| 3 e 3
i
Each branch is then a well-defined function (each input gives a unique output). If
we plot all three branches in a single graph, we obtain a beautiful surface below:
Another example of a multi-valued function is the argument map with domain C\{0}.
It is denoted and defined as:
arg(z) = {Arg(z) + 2kπ : k ∈ Z}.
38 2. Holomorphic Functions
1
Figure 2.3. Three branches of z 7→ Re z 3 .
Recall that Arg(z) is defined to be the unique angle θ0 ∈ (−π, π] such that z = |z| eiθ0 .
An element in the set arg(z) is any angle θ such that z = |z| eiθ . For example,
n π π π π o
arg(i) = . . . , − 2π, , + 2π, + 4π, . . .
n 2 2 2 2 o
π π π π
arg(1 + i) = . . . , − 2π, , + 2π, + 4π, . . .
4 4 4 4
For each fixed k ∈ Z, we regard z 7→ Arg(z) + 2kπi as a branch of the multi-valued
function arg(z). Below is the graph of five of its branches:
Solution
(z + w)n − z n
f 0 (z) = lim
w→0 w
(z+w)n
z }| {
z n + C1n z n−1 w + C2n z n−2 w2 + . . . + Cn−1
n
zwn−1 + wn −z n
= lim
w→0 w
C1n z n−1 w + C2n z n−2 w2 + . . . + Cn−1
n
zwn−1 + wn
= lim
w→0 w
n n−1 n n−2 n
= lim C1 z + C2 z w + . . . + Cn−1 zwn−2 + wn−1
w→0
= C1n z n−1 + 0 + . . . + 0
= nz n−1 .
A natural question: How is it the different from single-variable calculus? The key
distinction is that (2.2) is a multivariable limit since w is a complex number! By writing
w = h + ki, z = x + yi and f (x + yi) = u(x, y) + iv(x, y) where u, v, x, y, h, k are real,
the limit in (2.2) can be rewritten as:
f (z+w)=f ((x+yi)+(h+ki)) f (z)
z }| { z }| {
(u(x + h, y + k) + iv(x + h, y + k)) − (u(x, y) + iv(x, y))
f 0 (x + yi) = lim .
(h,k)→(0,0) h + ki
40 2. Holomorphic Functions
Past experience in MATH 2023/3033/3043 tell us that a multivariable limit exists less
likely than a single variable limit, since we require the limit not only exist, but also equal
when (h, k) approaches (0, 0) in all possible directions in the hk-plane. Therefore, there
are many complex-valued functions, which look simple and elementary, are not complex
2
differentiable. For instance, let f (z) = |z| = x2 + y 2 , then:
In general, 2x 6= −2yi (unless x = y = 0) and so f 0 (z) does not exist for almost all
x + yi ∈ C. Nonetheless, f (z) = x2 + y 2 is a polynomial of x and y and so it is real
differentiable everywhere on the xy-plane R2 . From this example, we can see that complex
differentiability is much more restrictive than real differentiability!
Exercise 2.5. Show that f (z) = 1/z is complex differentiable at any z ∈ C\{0}.
and the values of the limit are equal as (h, k) approaches (0, 0) from any direction.
2.2. Complex Differentiability 41
In particular, along the path {k = 0}, the above limit equals to:
f (x0 + y0 i + h) − f (x0 + y0 i)
lim
h→0 h
u(x0 + h, y0 ) + iv(x0 + h, y0 ) − u(x0 , y0 ) − iv(x0 , y0 )
= lim
h→0 h
u(x0 + h, y0 ) − u(x0 , y0 ) v(x0 + h, y0 ) − v(x0 , y0 )
= lim +i
h→0 h h
∂u ∂v
= (x0 , y0 ) + i (x0 , y0 ).
∂x ∂x
Along the path {h = 0}, the limit (2.4) equals to:
f (x0 + y0 i + ki) − f (x0 + y0 i)
lim
k→0 ki
u(x0 , y0 + k) + iv(x0 , y0 + k) − u(x0 , y0 ) − iv(x0 , y0 )
= lim
k→0 ki
u(x0 , y0 + k) − u(x0 , y0 ) v(x0 , y0 + k) − v(x0 , y0 )
= lim +i
k→0 ki ki
1 ∂u ∂v
= (x0 , y0 ) + i (x0 , y0 )
i ∂y ∂y
∂v ∂u
= (x0 , y0 ) − i (x0 , y0 )
∂y ∂y
Therefore, we get
E1 (x, y) = u(x, y) − u(x0 , y0 ) − ux (x0 , y0 ) (x − x0 ) − uy (x0 , y0 ) (y − y0 )
E2 (x, y) = v(x, y) − v(x0 , y0 ) − vx (x0 , y0 ) (x − x0 ) − vy (x0 , y0 ) (y − y0 ).
The fact that |Ei (x, y)| ≤ o(|z − z0 |) is equivalent to u and v being (real) differentiable
at (x0 , y0 ) (simply by the definition of real differentiability).
It is notable that the converse of Proposition 2.1 is not true. If f (z) = u(x, y)+iv(x, y)
satisfies the Cauchy-Riemann equations (2.3) at (x, y), it may not be true that f is complex
differentiable at z = x + yi. Here is one counter-example:
p p
f (z) = |Re(z)Im(z)| = |xy|.
p
Then u(x, y) = |xy| and v(x, y) = 0. We claim that ∂u ∂u
∂x (0, 0) = ∂y (0, 0) = 0:
∂u u(0 + h, 0) − u(0, 0)
(0, 0) = lim
∂x h→0 h
0−0
= lim = 0.
h→0 h
Similarly, one can also show ∂u
∂y (0, 0) = 0. It is obvious that
∂v
∂x = ∂v
∂y = 0. Therefore, the
Cauchy-Riemann equation (2.3) holds at (x, y) = (0, 0).
However, the function f (z) is not complex differentiable at z = 0 + 0i. It is because
when computing the limit (2.4) at (x, y) = (0, 0), we get:
p
f ((0 + 0i) + (h + ki)) − f (0 + 0i) |hk|
lim = lim .
(h,k)→(0,0) h + ki (h,k)→(0,0) h + ki
Along the path {h = 0}, the limit equals to 0. However, along the path {h = k}, the limit
becomes:
√
h2 |h|
lim = lim
h→0 h + hi h→0 h(1 + i)
1 1
which does not exist as it approaches 1+i as h → 0+ , while it approaches − 1+i as
h → 0− .
This f (z) is one example that the Cauchy-Riemann equation holds at a point, but
the function is not complex differentiable at that point. Fortunately, if we assume further
that u(x, y) and v(x, y) are (real) differentiable functions, then the Cauchy-Riemann
equations imply complex differentiability.
Proof. We define:
∂u ∂u
E1 (x, y) = u(x, y) − u(x0 , y0 ) − (x − x0 ) − (y − y0 )
∂x (x0 ,y0 ) ∂y (x0 ,y0 )
∂v ∂v
E2 (x, y) = v(x, y) − v(x0 , y0 ) − (x − x0 ) − (y − y0 ).
∂x (x0 ,y0 ) ∂y (x0 ,y0 )
2.2. Complex Differentiability 43
Since both u and v are differentiable at (x0 , y0 ), we have (by definition of differentiabil-
ity):
E1 (x, y) E2 (x, y)
lim = lim = 0.
(x,y)→(x0 ,y0 ) |z − z0 | (x,y)→(x0 ,y0 ) |z − z0 |
In short, combining Propositions 2.1 and 2.2, for any functions f = u + iv such
that u and v are (real) differentiable on an open domain Ω, complex differentiability
is equivalent to the Cauchy-Riemann equations. Many functions we will encounter are
(real) differentiable. Therefore, to show such a function is complex differentiable it
suffices to verify the Cauchy-Riemann equations. Let’s look at some examples.
44 2. Holomorphic Functions
Example 2.6. Determine all z at which the following functions are complex differ-
entiable.
(a) f (z) = z
2
(b) f (z) = |z|
(c) f (z) = ez
Solution
(a) f (z) = z = x − yi. Hence u = x and v = −y, which are clearly C 1 .
∂u ∂u
=1 =0
∂x ∂y
∂v ∂v
=0 = −1
∂x ∂y
Cauchy-Riemann equations do not hold at every point, hence f is not complex
differentiable at any point z ∈ C.
2
(b) f (z) = |z| = x2 + y 2 . Hence u = x2 + y 2 and v = 0, which are clearly C 1 .
∂u ∂u
= 2x = 2y
∂x ∂y
∂v ∂v
=0 =0
∂x ∂y
Cauchy-Riemann equations hold if and only if (x, y) = (0, 0). Therefore, f is
complex differentiable at z = 0 only.
(c) f (z) = ez = ex cos y + iex sin y. Hence u = ex cos y and v = ex sin y, which are
C 1 functions
∂u ∂u
= ex cos y = −ex sin y
∂x ∂y
∂v ∂v
= ex sin y = ex cos y
∂x ∂y
which are all continuous. Clearly the Cauchy-Riemann equations hold at every
(x, y), hence by Proposition 2.2, the function f (z) = ez is complex differentiable
at every z ∈ C.
From now on, we will use a more professional term to describe complex differentiable
functions:
Solution
(a) Given that f is holomorphic on C, the Cauchy-Riemann equations show:
∂f ∂f
=i .
∂y ∂x
Similarly, f is also holomorphic on C, we have:
∂f ∂f ∂f ∂f ∂f ∂f
=i =⇒ =i =⇒ = −i .
∂y ∂x ∂y ∂x ∂y ∂x
Combining both results, we have:
∂f ∂f ∂f
i = −i =⇒ = 0.
∂x ∂x ∂x
Since f 0 (z) = ∂f 0
∂x for holomorphic functions, we have f (z) = 0 for any z ∈ C.
Therefore, f (z) is a constant function.
(b) Denote h(z) = U (x, y) + iV (x, y). Since h(z) = g(z) = u(x, −y) − iv(x, −y), we
get:
U (x, y) = u(x, −y) and V (x, y) = −v(x, −y).
By Proposition 2.1, g is holomorphic on C implies u and v are real differentiable on
R2 , and hence U and V are real differentiable on R2 as well. Hence, according to
Proposition 2.2, it suffices to show h satisfies the Cauchy-Riemann equations:
∂h ∂ ∂ ∂ ∂
(z) = U (x, y) + i V (x, y) = u(x, −y) + i (−v(x, −y))
∂x ∂x ∂x ∂x ∂x
∂u ∂v ∂u ∂u
= (x, −y) − i (x, −y) = (x, −y) + i (x, −y)
∂x ∂x ∂x ∂y
∂h ∂ ∂ ∂ ∂
(z) = U (x, y) + i V (x, y) = u(x, −y) − i v(x, −y)
∂y ∂y ∂y ∂y ∂y
∂u ∂v ∂u ∂u
= (x, −y) + i (x, −y) = − (x, −y) + i (x, −y)
∂y ∂y ∂y ∂x
∂u ∂u ∂h
=i (x, −y) + i (x, −y) = i (z).
∂x ∂y ∂x
Therefore, h satisfies the Cauchy-Riemann equations and that U and V are real
differentiable on R2 , by Proposition 2.2, h is holomorphic on C.
Exercise 2.7. Determine all z’s in the complex plane at which the following functions
are complex differentiable:
(a) f (z) = 1/z
2
(b) f (z) = z |z|
(c) f (z) = z 2
46 2. Holomorphic Functions
Exercise 2.9. Suppose f (z) is complex differentiable at every z ∈ C. Prove that any
one of the following conditions imply that f is constant:
(i) Re(f ) is constant.
(ii) Im(f ) is constant.
(iii) |f | is constant.
Exercise 2.11. Show that the Cauchy-Riemann equations (2.3) for the function
f = u + iv is equivalent to:
∂f ∂f
=i .
∂y ∂x
Give a geometric interpretation of this result.
Exercise 2.13. Fix a complex number w such that |w| < 1. Consider the map
f : B1 (0) → C defined by:
w−z
f (z) = .
1 − wz
Show that:
(a) f is well-defined on B1 (0);
(b) f holomorphic on B1 (0);
(c) The image of f is B1 (0);
(d) f is bijective as a map f : B1 (0) → B1 (0).
Example 2.8. Find the largest open subset Ω ⊂ C such that the function:
f (x + yi) = x2 − y 2 + 2xyi
is holomorphic on Ω.
Solution
First we divide the complex plane C into open regions:
U = {x + yi : x2 − y 2 > 0} and V = {x + yi : x2 − y 2 < 0}.
Im
Re
Exercise 2.15. Find the largest open subset Ω ⊂ C such that the function:
f (x + yi) = |x| + |y| i
is holomorphic on Ω.
Therefore, ∇u and ∇v are perpendicular to each other provided that they are both non-
zero. Geometrically speaking, it means that the level sets of u and v are perpendicular!
There are more we can say about ∇u and ∇v. Define the matrix:
0 −1
J= .
1 0
48 2. Holomorphic Functions
Given any vector x ∈ R2 , the product J(x) is the vector obtained by rotating x counter-
clockwisely by π2 . The Cauchy-Riemann equations can be rewritten as:
vx −uy 0 −1 ux
= = .
vy ux 1 0 uy
As a result, we have ∇v = J(∇u). It means that for a holomorphic function f , the vector
∇v can be obtained by rotating ∇u counter-clockwisely by π2 .
Proposition 2.6. Let γk (t) = xk (t) + iyk (t) where k = 1, 2 be two C 1 curves in C which
intersect at t = 0 at the point z0 ∈ C. Suppose f : Bε (z0 ) → C is complex differentiable at
z0 and f 0 (z0 ) 6= 0 then we have:
(f ◦ γ1 )0 (0) γ10 (0)
= .
(f ◦ γ2 )0 (0) γ20 (0)
As a result, the angle between γ1 and γ2 at z0 is preserved under f . See Figure 2.4.
Proof. Let f = u + iv (as usual). Recall from Exercise 2.11 that the Cauchy-Riemann
equation is equivalent to
∂f ∂f
=i .
∂y ∂x
Using this and the chain rule, we have:
d ∂f dx1 ∂f dy1
(f ◦ γ1 ) = +
dt ∂x dt ∂y dt
∂f ∂f
(f ◦ γ1 )0 (0) = x0 (0) + i y10 (0)
∂x z0 1 ∂x z0
∂f
= γ10 (0).
∂x z0
Likewise, we have:
∂f
(f ◦ γ2 )0 (0) = γ20 (0)
∂x z0
2.2. Complex Differentiability 49
∂f ∂f
By Exercise 2.12, |f 0 (z0 )| = (z0 ) . Given f 0 (z0 ) 6= 0, we also have (z0 ) 6= 0 and
∂x ∂x
hence:
(f ◦ γ1 )0 (0) fx (z0 ) γ10 (0) γ 0 (0)
0
= 0 = 10
(f ◦ γ2 ) (0) fx (z0 ) γ2 (0) γ2 (0)
as desired.
As such, the complex logarithm is multi-valued and so for each z 6= 0, log z is a set
rather than a single number.
Let’s derive a general formula for log z. Given any z 6= 0, we first express it in polar
form z = |z| eiArg(z) . Then, w = u + iv satisfies ew = z if and only if eu+iv = |z| eiArg(z) :
As a result, we have w = ln |z| + i(Arg(z) + 2kπ) for any k ∈ Z, and using set notations:
Remark 2.8. To avoid confusion, from now on we will denote ln as the real, single-valued
logarithm, and use log for the complex multi-valued logarithm.
Recall that for the argument function arg(z), we define the principal argument Arg(z)
as the unique angle θ ∈ (−π, π] such that z = |z| eiθ . Similarly, we define the principal
logarithm to be:
Definition 2.10 (Principal Logarithm). For any z 6= 0, we define its principal logarithm
to be:
Log(z) := ln |z| + iArg(z).
Exercise 2.17. Find log(z) and Log(z) for each of the following z’s:
(a) z = −2
√
(b) z = 1 + 3i
(c) z = −i
Exercise 2.19. For any two subsets S and T of C, we define the sum of the two sets
to be:
S + T := {s + t : s ∈ S and t ∈ T }.
Show that for any z1 , z2 6= 0, we have:
arg(z1 z2 ) = arg(z1 ) + arg(z2 )
log(z1 z2 ) = log(z1 ) + log(z2 ).
52 2. Holomorphic Functions
d
We next compute log(z). Note that log(z) is multi-valued. Differentiating log(z)
dz
basically means differentiating each branch of log(z) individually. We first verify that
every branch of log(z) is holomorphic on C\{x + 0i : x ≤ 0}:
log(z) = {ln |z| + i(Arg(z) + 2kπ) : k ∈ Z}
1
Let u(x, y) = ln |z| = ln(x2 + y 2 ) and v(x, y) = Arg(z) + 2kπ. We leave it as an exercise
2
for readers to verify that on C\(−∞, 0]:
∂u x ∂u y
= 2 = 2
∂x x + y2 ∂y x + y2
∂v y ∂v x
=− 2 = 2 .
∂x x + y2 ∂y x + y2
Clearly, they are all continuous and satisfy the Cauchy-Riemann equation. By Proposition
2.2, each branch ln |z| + i(Arg(z) + 2kπ) is holomorphic on C\{x + 0i : x ≤ 0}, and so:
d ∂u ∂v x − yi z 1
(ln |z| + i(Arg(z) + 2kπ)) = +i = 2 = = .
dz ∂x ∂x x + y2 zz z
Therefore, every branch of log(z) has the same complex derivative, so we may simply
write:
d 1
(2.6) log(z) = .
dz z
Using (2.7), one can recover the n-th root formula stated in Definition 1.9. For any
complex number a 6= 0, according to (2.7), we have:
1 1
a n = e n log a
n 1 o
= e n (ln|a|+(Arg(a)+2kπ)i) : k ∈ Z
n 1 Arg(a)+2kπ
o
= e n ln|a| e n i
:k∈Z
p
n Arg(a) + 2kπ Arg(a) + 2kπ
= |a| cos + i sin :k∈Z
n n
p
n Arg(a) + 2kπ Arg(a) + 2kπ
= |a| cos + i sin : k = 0, 1, 2, . . . , n − 1
n n
Using the chain rule, one can derive the differentiation formula for complex powers.
Regard w as a fixed complex numbers, we can derive:
d w d w log z d
(2.8) z = e = ew log z w log z
dz dz dz
w
= z w · = wz w−1
z
d z d z log w d
(2.9) w = e = ez log w z log w
dz dz dz
= wz log w
Exercise 2.23. Fix a non-zero complex number w. On what domains are the
following functions holomorphic?
(a) f (z) = wz
(b) g(z) = z w
For complex logarithms, we use log z and Log z to distinguish the multi-valued
logarithm, or the principal branch of logarithm. Similarly, for a complex power z w , which
is defined via logarithms, we can also define its principal branch by ewLog(z) . For instance,
π
the principal branch of ii is e− 2 . Unlike logarithms and arguments, we do not introduce
a new symbol to denote the principal branch of z w . With abuse of notations, we may
π π
sometimes simply write, for instance ii = e− 2 , to mean e− 2 is the principal value of ii .
It may cause ambiguity, but such an ambiguity will cause less nuisance if we state clearly
in the context whether z w means a multi-valued power or the principal branch.
for any x ∈ R. Hence by writing sin x and cos x in terms of the exponentials, we get:
When z is a complex number, it does not make sense to define sin z and cos z by
regarding z as an “angle”. Thanks for the above identities, we define the complex sin and
cos functions as:
d iz d −iz
e = ieiz e = −ie−iz ,
dz dz
and so it is easy to verify that both sin z and cos z are entire functions and their derivatives
are:
d d
(2.10) sin z = cos z cos z = − sin z.
dz dz
The other trigonometric functions are defined similarly as in the real case:
sin z 1
tan z = sec z =
cos z cos z
cos z 1
cot z = csc z =
sin z sin z
Using the product and quotient rules, one can easily derive that:
d d
(2.11) tan z = sec2 z sec z = sec z tan z
dz dz
d d
(2.12) cot z = − csc2 z csc z = − csc z cot z
dz dz
As a result, we have:
1/2
sin−1 z = −i log iz + 1 − z 2 .
For example,
2.3.6. Mapping Properties. Recall that given any z = x+yi, we have ez = ex eiy =
e cos y + iex sin y. Therefore, along the straight-path z = x0 + yi (where x0 is fixed and
x
y varies), the image under the map z 7→ ez is given by (u, v) = ex0 (cos y + i sin y), which
is a circle with radius ex0 centered at the origin. In other words, the complex exponential
maps vertical lines in the xy-plane to concentric circles in the uv-plane. On the other hand,
along the horizontal path z = x + y0 i, the image is given by (u, v) = ex (cos y0 + i sin y0 ),
which is a half-line from the origin (but not passing through it).
1.0
2.0
0.8
1.5
0.6
0.4 1.0
0.2
0.5
0.0
0.0 0.2 0.4 0.6 0.8 1.0 1.0 1.5 2.0 2.5
Exercise 2.28. Show that the complex sine function f (z) = sin z maps horizontal
lines in the xy-plane to ellipses in the uv-plane, and maps vertical lines in the
xy-plane to hyperbolas in the uv-plane.
Chapter 3
Contour Integrals
We will learn how they are defined and how they can be computed soon. In the first
glance, it appears quite similar to line integrals in Multivariable Calculus. However,
when combining with properties of holomorphic functions, there are many beautiful
and amazing results concerning complex contour integrals which did not appear in line
integrals. One notable result is Cauchy’s integral formula, an elegant theorem which
leads to many important results in Complex Analysis and beyond.
57
58 3. Contour Integrals
I
Example 3.1. Compute the line integral f (z) dz for each of the functions below.
γ
Here γ is the circle with radius 2 centered at the origin.
(a) f (z) = z 2
1
(b) f (z) =
z
(c) f (z) = z
Solution
γ can be parametrized by:
z(t) = 2eit , t ∈ [0, 2π].
Therefore dz = 2ieit dt.
(a)
I Z 2π Z 2π
2
it 2 it
z dz = 2e · 2ie dt = 8ie3it dt
γ 0 | {z } | {z } 0
dz
z2
t=2π
1 3it
= 8i e
3i t=0
8 6πi 8
= e − e0 = (1 − 1) = 0.
3 3
(b)
I Z 2π Z t=2π
1 1
dz = · 2ieit dt = i dt = 2πi.
γ z 0 2eit t=0
(c)
I Z 2π Z 2π
−it it
z dz = 2e · 2e dt = 4 dt = 8πi.
γ 0 0
d 1 3it
Remark 3.3. In part (a) of the above example, we have used the fact that dt 3i e =
e3it , and also Fundamental Theorem of Calculus. In general, just like in the real case,
if F (t) is a differentiable function of t on [a, b] such that F 0 (t) = ϕ(t) on [a, b], then we
have
Z t=b
ϕ(t) dt = F (b) − F (a).
t=a
3.1. Complex Integrations 59
However, we sometimes need to be more careful when applying this. Try to find out
what’s wrong with the calculation below:
I Z 2π 2π
1 1 it 1 it
dz = ie dt = − Log(1 − 2e )
|z|=1 1 − 2z 0 1 − 2eit 2 0
1
= − (Log(−1) − Log(−1)) = 0???
2
Solution
First we parametrize L:
z(t) = (1 − t)z1 + tz2 , t ∈ [0, 1].
Then, we have dz = (z2 − z1 ) dt, and so:
Z Z 1
z
e dz = ez1 +t(z2 −z1 ) · (z2 − z1 ) dt
L 0
1
1
= ez1 +t(z2 −z1 ) · (z2 − z1 )
z2 − z1 0
= ez2 − ez1 .
3
Γ2
L1
1 Γ1
x
1 L2 3
Proof. First assume that γ is C 1 . Suppose the path γ can be parametrized by:
z(t) = x(t) + iy(t), t ∈ [a, b].
0
Then, we have dz = z (t) dt, and hence:
Z Z b
f (z) dz = f (z(t)) · z 0 (t) dt
γ a
Z b
= F 0 (z(t)) ·z 0 (t) dt
a | {z }
f (z(t))
Z b
d
= F (z(t)) dt (chain rule)
a dt
= F (z(b)) − F (z(a))
= F (z2 ) − F (z1 ).
If γ is only piecewise C 1 , we can decompose γ = γ1 + · · · + γk so that each γi is C 1 .
Then, one can argue as above for each γi , and finally obtain the desired result by adding
a telescope sum.
Remark 3.5. If such an F (z) in Theorem 3.4 exists, then we call F (z) a primitive function
of f (z).
The above theorem is particularly useful when the anti-derivative of f is easy to find.
For example, if γ is any continuous piecewise C 1 path from z1 to z2 , we can find easily
that:
Z 3 z2
z z 3 − z13
2
z dz = = 2
γ 3 z1 3
Z
z
ez dz = [ez ]z21 = ez2 − ez1 .
γ
Exercise 3.2. Let γ1 be the path which starts from (0, 0), first to (1, 1), then to (0, 2).
Let γ2 be the path which starts from (0, 0), then straight to (0, 2). Verify the following
by direct computations:
Z Z
πz πz
cos dz = cos dz.
γ1 2 γ2 2
Then, verify that Theorem 3.4 gives the same result.
However, it is important to note that Theorem 3.4 requires the curve γ to be inside
Ω (on which F 0 (z) = f (z) holds). Let’s consider the function f (z) = z1 . Although we
3.1. Complex Integrations 61
d
usually simply write dz Log(z) = z1 , it is only true for z ∈ C\{x + 0i : x ≤ 0} since Log(z)
is not continuous on the negative x-axis.
Therefore, we can only apply Theorem 3.4 when the curve γ lies inside Ω :=
C\{x + 0i : x ≤ 0}. For instance, we still have
I
1
=0
γ1 z
where γ1 is the unit circle centered at 2 + 0i with radius 1. This closed curve γ1 is
contained inside Ω.
H
However, it is incorrect to claim γ2 z1 = 0 where γ2 is the unit circle centered at the
origin. The reason is that this closed curve passes through the negative x-axis (hence not
contained inside Ω). In fact we can directly verify that:
I
1
= 2πi.
γ2 z
y
γ2
γ1
Fortunately, we can still apply Theorem 3.4 on f (z) = z12 when the integration curve
γ does not pass through the origin. The reason is that F (z) = − z1 is a primitive function
for f such that F 0 (z) = f (z) holds on C\{0}. Therefore, we have:
I
1
2
=0
γ z
for any closed curve γ not passing through the origin. Also, for a path L in C\{0}
connecting z1 to z2 , we have:
Z z
1 1 2 1 1
2
dz = − = − .
L z z z1 z1 z2
Z
Exercise 3.4. Evaluate the integral |z| dz where γ is each of the following:
γ
Proof. There is a nice trick in the proof that readers are recommended to learn. Let
Z
I= f (z) dz.
γ
Express I in polar form: I = |I| eiθ , then we have e−iθ I = |I| which is real! Suppose γ is
parametrized by z(t) = x(t) + iy(t) where a ≤ t ≤ b, then:
Z Z
−iθ −iθ
e I=e f (z) dz = e−iθ f (z) dz
γ γ
Z b
= Re e−iθ f (z) + iIm e−iθ f (z) (x0 (t) + iy 0 (t)) dt
a
Z b
= Re e−iθ f (z) x0 (t) − Im e−iθ f (z) y 0 (t) dt.
a
The last equality above follows from the fact that e−iθ I is real.
Then, we use Cauchy-Schwarz’s inequality to bound the integrand:
Re e−iθ f (z) x0 (t) − Im e−iθ f (z) y 0 (t)
q q
2 2 2 2
≤ (Re (e−iθ f (z))) + (Im (e−iθ f (z))) (x0 (t)) + (y 0 (t))
= e−iθ f (z) |z 0 (t)| = |f (z)| |z 0 (t)| ≤ M |z 0 (t)| .
Finally, we get:
Z b
e−iθ I ≤ M |z 0 (t)| = M L,
a
and hence |I| ≤ M L, completing the proof.
3.1. Complex Integrations 63
Z
Remark 3.7. If we estimate the integral f (z) dz in a more direct way by writing
γ
f = u + iv and then consider the following:
Z Z Z b
f (z) dz = (u + iv)(dx + idy) = (ux0 − vy 0 ) + i(vx0 + uy 0 ) dt
γ γ a
v !2 !2
u Z Z
u b b
= t 0 0
(ux − vy ) dt + (vx0 + uy 0 ) dt .
a a
then after applying Cauchy-Schwarz’s inequality to each integral, the best we can achieve
is Z √
f (z) dz ≤ 2M L,
γ
which is weaker than the result in Lemma 3.6.
Solution
For any z ∈ C such that |z| = 1, we have:
x y
1 −i x2 +y
e z = e x2 +y2 2
= ex−iy = ex e−iy ,
1
e z = ex ≤ e1 = e.
Here we have used the fact that −1 ≤ x ≤ 1 along the curve |z| = 1.
Therefore, by Lemma 3.6, we have:
I
1
e z dz ≤ |{z}
2π |{z}
e .
|z|=1
L M
Solution
We are interested in the limit when R → +∞, so we can assume R > 1 so that
the contour circle |z| = R does not pass through 1 (at which the integrand is
undefined).
On the contour |z| = R, we have |z − 1| ≥ R − 1 (draw a diagram to convince
yourself on that), so we have:
1 1 1
(z − 1)2
= 2 ≤ (R − 1)2 on |z| = R.
|z − 1|
| {z }
M
64 3. Contour Integrals
The length of the contour |z| = R is 2πR. Hence, by Lemma 3.6, we get
I
1 1
2
dz ≤ 2πR · .
|z|=R (z − 1) (R − 1)2
2πR
From elementary calculus, we have lim = 0, and the desired result
R→+∞ (R − 1)2
follows from the squeeze theorem.
Solution
For |z| = R 1, we have
Log(z) ln(R) + iArg(z) 1 1
= ≤ 2 (|ln(R)| + |Arg(z)|) ≤ 2 (ln(R) + π).
z2 R2 R R
Using ML inequality
I
Log(z) Log(z) 2π
2
dz ≤ 2
· 2πR ≤ (ln(R) + π) → 0 as R → +∞.
|z|=R z z R
Here we have used the l’Hospital’s rule. We conclude:
I
Log(z)
lim dz = 0.
R→+∞ |z|=R z2
Using Cauchy-Goursat’s Theorem, we can immediately conclude that all the integrals
below over any closed curve γ ∈ C are zero, without performing any calculation:
I I I
ez dz, sin z dz, z 2 dz, etc.
γ γ γ
1. Goursat’s theorem 35
(1)
T3
T (1)
T2
(1) (1)
T1 T4
(0)
Figure Figure
3.2. Divide1.
theBisection of 4Ttriangles
contour T into
For simplicity, denote the numerator by E(z) := f (z) − f (z0 ) − f 0 (z0 )(z − z0 ), then we
have:
E(z)
(3.2) lim = 0.
z→z0 z − z0
0
Since the function f (z0 ) + f 0 (z0 )(z − z0 ) has a primitive function zf (z0 ) + f (z0 )
2 (z − z0 )2
(note that z0 is a fixed point), we have
I I I
E(z) dz = [f (z) − f (z0 ) − f 0 (z0 )(z − z0 )] dz = f (z) dz.
T (n) T (n) T (n)
Therefore, to bound the RHS of (3.1), we can consider the integral of E(z) instead,
which is very small according to (3.2).
Now, given any ε > 0, by (3.2), there exists δ > 0 such that whenever z ∈ Bδ (z0 ),
E(z)
we have < ε. Recall that {∆(n) }∞ n=0 is a strictly decreasing sequence of triangles
z − z0
“converging” to the point z0 . Hence, for sufficiently large n, ∆(n) must lie inside the ball
Bδ (z0 ), and so |E(z)| < ε |z − z0 | for any z ∈ ∆(n) ⊂ Bδ (z0 ).
Recall that |z − z0 | is the distance between z and z0 , both of which are in ∆(n) . By
elementary geometry, the distance between any two points in a triangle must be bounded
by the perimeter of the triangle. Hence, we have for any z ∈ ∆(n) ,
εL0
(3.3) |E(z)| < ε |z − z0 | ≤ εLn =
2n
where Ln denotes the perimeter of the triangle T (n) .
Using (3.3), we can apply Lemma 3.6 to show:
I
εL0 εL2
E(z) dz ≤ n · Ln = n0 .
T (n) 2 4
Finally, by considering (3.1), we have proved:
I I I
n n εL20
f (z) dz ≤ 4 f (z) dz = 4 E(z) dz ≤ 4n · = εL20 .
T T (n) T (n) 4n
Since ε > 0 is arbitrary, by letting ε → 0+ , we get:
I
f (z) dz = 0,
T
completing Step 1.
Exercise 3.9. Using the result proved so far, show that Cauchy-Goursat’s Theorem
holds for any closed polygon γ.
Exercise 3.10. Show that if 4ABC is contained inside a simply-connected open set
Ω on which f is holomorphic, then we have:
Z Z Z
f (z) dz = f (z) dz + f (z) dz.
L(A,C) L(A,B) L(B,C)
Exercise 3.11. Which part in the proof of Step 1 will break down if f is not
holomorphic? Also, why will the proof break down if Ω is not simply-connected?
68 3. Contour Integrals
F (z + w) − F (z)
We claim that F 0 (z) = f (z) by showing that the quotient tends to f (z)
w
as w → 0.
z+w
z
z0
Z
Applying Lemma 3.6 on the integral (f (ξ) − f (z)) dξ, we have:
L(z,z+w)
Z
(f (ξ) − f (z)) dξ ≤ ε · |w| ,
L(z,z+w) |{z}
length of contour
Z
1
which implies (f (ξ) − f (z)) dξ ≤ ε (whenever 0 < |w| < δ), or equivalently,
w L(z,z+w)
Z
1
lim (f (ξ) − f (z)) dξ = 0.
w→0 w L(z,z+w)
where γ(z0 , z) is any polygonal path from z0 to z. Can we still claim that F 0 (z) = f (z)
with the same proof? If not, where does the proof break down?
3.2.3. Step 3: Completion of the Proof. We have by far proved that Cauchy-
Goursat’s Theorem holds when at least one of the conditions holds:
(i) γ is a closed polygon; or
(ii) Ω is convex.
Now we deduce the general case based on these special cases.
Given any simply-connected domain Ω and any closed piecewise
I C 1 curve γ ⊂ Ω,
and a holomorphic function f : Ω → C, the key idea to show f (z) dz = 0 is to break
γ
the region enclosed by γ into small rectangles {Rj }N M
j=1 and “partial rectangles” {γk }k=1
(see Figure 3.3). By breaking the region into small enough of these rectangles and partial
rectangles, we may assume that these partial rectangles are contained inside an convex
subset of Ω. This is intuitively true, but the proof involves some deep knowledge on
analysis and topology beyond the scope of this course. Interested readers may read
70 3. Contour Integrals
Stein’s book Appendix B. For most of our applications, we will deal with closed curves
whose enclosed regions are as simple as rectangles and balls, and the union of them.
Figure 3.3
For each rectangle Rj and partial rectangle γk , results from Steps 1 and 2 show
I I
f (z) dz = f (z) dz = 0.
Rj γk
It could be to many people’s surprise that complex analysis in fact has real applications
such as evaluating a complicated real integral which cannot be easily done by finding
anti-derivatives.
Example 3.6. Show that Z ∞
1 − cos x
dx = π
−∞ x2
1 − eiz
by considering the contour below and the function f (z) = .
z2
y
CR
−Cε
x
−R −ε ε R
Solution
iz
The function 1−ez2 is holomorphic on (for instance) {z ∈ C : |z| > 0, Arg(z) 6= 3π 2 }
which is a simply connected open domain containing γR,ε . Using Cauchy-Goursat’s
Theorem: I
1 − eiz
dz = 0 ∀0 < ε < R
γR,ε z2
Therefore, we have:
Z Z Z Z
1 − eiz 1 − eiz 1 − eiz 1 − eiz
2
dz − 2
dz + 2
dz + dz = 0.
CR z Cε z [−R,−ε] z [ε,R] z2
We will then take R → +∞ and ε → 0+ .
We first estimate the integral over CR . For all z ∈ CR , we have Im(z) ≥ 0 =⇒
Re(iz) = −Im(z) ≤ 0 =⇒ eiz = eRe(iz) ≤ 1. Therefore, the numerator satisfies
1 − eiz ≤ 1 + eiz ≤ 2. Using ML inequality,
Z
1 − eiz 2 4π
2
dz ≤ 2 · 2πR = → 0 as R → +∞.
CR z R R
This shows Z
1 − eiz
lim dz = 0.
R→+∞ CR z2
For the integral over Cε , we study the behavior near z = 0. By Taylor’s
expansion of eiz , we have
∞
1 − eiz 1
1 − (1 + iz + 2! (iz)2 + · · · ) i X (iz)n
= = − − .
z2 z2 z n=2 z 2 n!
Since
X∞ X∞ n−2 X∞
(iz)n |z| 1
2 n!
≤ ≤ ≤ e,
n=2
z n=2
n! n=2
n!
we have
1 − eiz i
+ < e.
z2 z
Using ML inequality,
Z Z
1 − eiz i + 1 − eiz i
2
+ dz ≤ e · π → 0 as ε → 0 =⇒ lim 2
+ dz = 0
Cε z z ε→0 +
Cε z z
Using the parametrization zε (t) = εeit where t ∈ [0, π], we have:
Z Z π Z π
i i
lim+ dz = lim+ zε0 (t) dt = lim+ ie−it ieit dt = −π
ε→0 Cε z ε→0 0 zε (t) ε→0 0
Combine the above results
Z Z Z
1 − eiz 1 − eiz i i
lim dz = lim+ + dz − dz = 0 − (−π) = π.
ε→0+ Cε z2 ε→0 Cε z2 z Cε z
γ2
γ1
since z1 is holomorphic on the green region. However, we cannot do the same for γ2 . Any
simply-connected region containing γ2 must contain 0 at which z1 is undefined. In this
section,
I we will introduce Cauchy’s integral formula to deal with contour integrals of the
f (z)
form dz.
γ z−α
For instance, given an entire function f : C → C, a point α, and two closed curves γ1
and γ2 below. Cauchy’s Integral Formula asserts that:
I I
f (z) f (z)
dz = 0 whereas dz = 2πif (α).
γ1 z − α γ2 z −α
γ2
γ1
3.3.1. Elementary Examples. We first illustrate the use of Cauchy’s integral for-
mula by a toy example:
I I (
1 1 2πi · 1 = 2πi if γ encloses 0
dz = dz =
γ z γ z−0 0 if γ does not enclose 0
Here we take f (z) = 1 which is an entire function on C.
Example 3.7. Evaluate the following contour integrals:
I
z
(a) dz
|z|=2 (z + 3i)(z − i)
I
z
(b) dz
|z|=4 (z + 3i)(z − i)
I
ez
(c) 2
dz
|z|=2 z + 1
Solution
(a) The integrand has two singularities: z = −3i and z = i. First observe that
the curve |z| = 2 enclose i only, and hence near the simply-connected region
z
|z| ≤ 2, the function f (z) := z+3i is holomorphic. Apply Cauchy’s integral
formula with this f , we get:
I I z
z z+3i z
dz = dz = 2πi ·
|z|=2 (z + 3i)(z − i) |z|=2 z − i z + 3i z=i
i πi
= 2πi · = .
i + 3i 2
(b) Note that the curve |z| = 4 enclose both singularities −3i and i of the integrand.
We cannot apply Cauchy’s integral formula by writing the integrand as either:
z z
z+3i z−i
or .
z−i z + 3i
The way out is to do partial fractions for the denominator. Let A and B be
complex numbers such that:
1 A B
= + .
(z + 3i)(z − i) z + 3i z − i
We need to solve for A and B:
1 A(z − i) + B(z + 3i)
=
(z + 3i)(z − i) (z + 3i)(z − i)
1 = (A + B)z + (−Ai + 3Bi)
Equating coefficients, we need A + B = 0 and (−Ai + 3Bi) = 1. Solving these
equations, we get A = 41 i and B = − 41 i, and hence:
1 1
1 i i
= 4 − 4 .
(z + 3i)(z − i) z + 3i z − i
74 3. Contour Integrals
Exercise 3.13. Use Cauchy’s integral formula to evaluate the following contour
integrals:
I
1
(a) 2+i
dz
|z|=2 z
I
1
(b) 2
dz
|z−eπi/4 |=1 z + i
I
1
(c) 3−1
dz
|z|=2 z
Try to do the problems in a rather tedious way using partial fractions. We will provide
another approach soon.
The result in Exercise 3.14 implies an important property about holomorphic func-
tions:
Lemma 3.12. If f is holomorphic on an open ball BR (α) such that |f (z)| ≤ |f (α)| for
any z ∈ BR (α), then f is constant on BR (α).
Proof of Proposition 3.11. Consider the set S := {z ∈ Ω : |f (z)| = |f (z0 )|}. By conti-
nuity of |f (z)|, the set S is obviously closed. We will argue that S is also open. Then, by
connectedness of Ω, we must have S = Ω (note that z0 ∈ S so S is non-emtpy).
Take any α ∈ S, so that |f (α)| = |f (z0 )|. Since Ω is open, there exists δ > 0 such
that Bδ (α) ⊂ Ω. Now it is given that |f (z0 )| is the maximum value of |f (z)| among all
z ∈ Ω. In particular, we have |f (α)| = |f (z0 )| ≥ |f (z)| for any z ∈ Bδ (α). Using Lemma
3.12, we conclude that f is constant on Bδ (α), which implies |f (z)| = |f (α)| = |f (z0 )|
for any z ∈ Bδ (α). In other words, Bδ (α) ⊂ S. This shows S is open.
3.3.2. Proof of Cauchy’s Integral Formula. The proof of Cauchy’s integral for-
mula is a reminiscence of the proof of generalized (i.e. with holes) Green’s Theorem in
Multivariable Calculus. Fix α ∈ C and consider a simple closed curve γ enclosing α. We
want to find out the value of the integral:
I
f (z)
dz.
γ z −α
We drill a circular hole near α in the region enclosed by γ, so that the following “key-hole”
contour Γε is produced.
−∂Bε (α)
ε
α
L −L
76 3. Contour Integrals
3.3.3. Cauchy’s Integral Formula with Multiple Holes. We have seen how to
1
apply Cauchy’s integral formula on fractions such as 2 which is not defined on z = i
z +1
and z = −i. If a simple closed contour γ enclosesboth singularities,
then we performed
1 1 1
partial fractions so that the fraction becomes − .
2i z − i z + i
Sometimes, partial fractions can be time-consuming especially when there are many
singularities. However, using the hole-drilling technique demonstrated in the proof of
Cauchy’s integral formula, we can break down the contour integral into a sum of several
contour integrals, each of which is over a contour that encloses only one singularity. Let’s
look at some examples.
Example 3.8. Evaluate the contour integral:
I
z
dz
|z|=4 (z + 3i)(z − i)
without using partial fractions.
Solution
The two singularities are z = −3i and z = i, both are contained inside the contour
|z| = 4. Draw two little circles with small radii ε around each singularity and
consider the key-hole contour:
Γ = γ1 + L1 − ∂Bε (−3i) − L1 + γ2 + L2 − ∂Bε (i) − L2
−L2 L2
ε
i
γ1 γ2
ε
−3i
L1 −L1
Then, the key hole contour Γ encloses a simply-connected region not containing
any singularity of the integrand. Therefore, Cauchy-Goursat’s Theorem asserts that
I
z
dz = 0.
Γ (z + 3i)(z − i)
On the other hand, by cancellation of the common sides, we have:
I Z Z I I I I I
= + − − = − − .
Γ γ1 γ2 |z+3i|=ε |z−i|=ε |z|=4 |z+3i|=ε |z−i|=ε
78 3. Contour Integrals
Therefore,
I
z
0= dz
(z + 3i)(z − i)
IΓ I
z z
= dz − dz
|z|=4 (z + 3i)(z − i) |z+3i|=ε (z + 3i)(z − i)
I
z
− dz.
|z−i|=ε (z + 3i)(z − i)
Therefore, we can break the required integral into the sum of two integrals:
I I I
z z z
dz = dz + dz
|z|=4 (z + 3i)(z − i) |z+3i|=ε (z + 3i)(z − i) |z−i|=ε (z + 3i)(z − i)
z
Since ε is very small, the function is holomorphic on |z + 2i| < ε, and so
z−i
Cauchy’s integral formula asserts that:
I I z
z z−i −3i 3πi
dz = dz = 2πi · = .
|z+3i|=ε (z + 3i)(z − i) |z+3i|=ε z − (−3i) −3i − i 2
For the second integral, we have:
I I z
z z+3i i πi
dz = dz = 2πi · =
|z−i|=ε (z + 3i)(z − i) |z−i|=ε z − i i + 3i 2
Adding up the results, we get:
I
z 3πi πi
dz = + = 2πi.
|z|=4 (z + 3i)(z − i) 2 2
Solution
First factorize the integrand:
1 1
3
=
z −1 (z − 1)(z − ω)(z − ω 2 )
2πi
where ω := e 3 is the cubic root of unity. There are three singularities, namely 1, ω
and ω 2 , all are enclosed by the given contour |z| = 2. By mimicking the hole-drilling
3.3. Cauchy’s Integral Formula I 79
Exercise 3.16. Let n be a positive integer, and ω := e2πi/n denote the n-th root of
unity. Express the contour integral:
I
1
n−1
dz
|z|=2 z
in terms of ω.
Exercise 3.17. Given any real constant a ∈ R, by considering the contour integral
I
eaz
dz, prove the following integration formula:
|z|=1 z
Z π
ea cos θ cos(a sin θ) dθ = π.
0
80 3. Contour Integrals
The corollary is a very remarkable and surprising result. In Real Analysis, there are
many functions which are differentiable for one time but not the second time or further.
However, this theorem and the corollary assert that once f is complex differentiable on a
simply-connected domain (say an open ball), then it is infinitely differentiable on that
domain!
3.4.1. Elementary Examples. Again, we will first see some examples of using the
higher-order Cauchy’s integral formula, then we will give a proof for it. As a quick
example: I
1
2
dz.
|z|=1 z
One way of evaluating it is to argue that its primitive function is − z1 , which is well defined
and holomorphic near the contour |z| = 1. Then by Proposition 3.4, the contour integral
is 0.
Let’s see how to obtain the same result using Theorem 3.13 (with n = 1, and
f (z) ≡ 1): I
1 1 d
1+1
dz = 1 = 0.
2πi |z|=1 z dz z=0
82 3. Contour Integrals
Example 3.10. Evaluate the contour integral using higher-order Cauchy’s integral
formula: I
e2z
3
dz.
|z|=1 z
Solution
In practice, it may be helpful to write the higher-order Cauchy’s integral formula as:
I
f (z) 2πi (n)
n+1
dz = f (α).
γ (z − α) n!
Let f (z) = e2z which is entire, then f 0 (z) = 2e2z and f 00 (z) = 4e2z . By Theorem
3.13 (with n = 2), we get:
I 2z I
e e2z
3
dz = 2+1
dz
γ z γ (z − 0)
2πi 00 2πi
= f (0) = · 4 = 4πi.
2! 2
Solution
The contour |z| = 3 encloses two singularities of the integrand, namely −i and 2i.
By the hole-drilling technique, we can pick a small ε > 0 such that:
I I I !
1 1
2 (z − 2i)3
dz = + 2 (z − 2i)3
dz.
|z|=3 (z + i) |z+i|=ε |z−2i|=ε (z + i)
Then we calculate each integral on the RHS individually:
I I 1
1 (z−2i)3
2 3
dz = 1+1
dz
|z+i|=ε (z + i) (z − 2i) |z+i|=ε (z + i)
2πi d 1 2πi
= =− 3
1! dz z=−i (z − 2i)3 3
I I 1
1 (z+i)2
dz = dz
|z−2i|=ε (z + i)2 (z − 2i)3 |z−2i|=ε (z − 2i)2+1
2
2πi d 1
=
2! dz z=2i (z + i)2
2
6
= πi ·
(z + i)4 z=2i
2πi
= 3
3
Therefore, I
1 2πi 2πi
dz = − 3 + 3 = 0.
|z|=3 (z + i)2 (z − 2i)3 3 3
3.4. Cauchy’s Integral Formula II 83
The higher-order Cauchy’s Integral Formula has some interesting applications on the
proof of some combinatorial identities. It is because for any non-negative integers n and
k (where n ≥ k), we have
I
1 zn 1 dk n n · (n − 1) · · · (n − k + 1)
k+1
dz = z = = Ckn
2πi |z−1|=r (z − 1) k! dz k z=1 k!
for any r > 0. Using this, one can derive some combinatorial identities in the example
below:
Example 3.12. Prove that
(a) C0n − C1n + C2n − · · · + (−1)n Cnn = 0 (where n ≥ 1).
m+n+1
(b) Cnn + Cnn+1 + Cnn+2 + · · · + Cnn+m = Cn+1
Solution
3.4.2. Proof of Higher Order Cauchy’s Integral Formula. Now we discuss the
proof of Theorem 3.13. From the (zeroth order) Cauchy’s integral formula, we know:
I
1 f (z)
f (α) = dz,
2πi γ z − α
where α is a point on the domain Ω, and γ is a simple closed curve in Ω enclosing α.
Note that if w ∈ C is very small, α + w will still be enclosed by γ, and so we have:
I
1 f (z)
f (α + w) = dz.
2πi γ (z − α − w)
Our first goal is to show Theorem 3.13 holds for f 0 (α), i.e. n = 1. Recall that:
f (α + w) − f (α)
f 0 (α) = lim .
w→0 w
We will use the zeroth order Cauchy’s integral formula to evaluate such a limit:
I I
1 1 f (z) f (z)
(3.6) f 0 (α) = lim dz − dz
2πi w→0 w γ z−α−w γ z−α
I
1 1 1 1
= lim f (z) · − dz.
2πi w→0 γ w z−α−w z−α
By straight-forward computation, we get:
1 1 1 1
− = .
w z−α−w z−α (z − α − w)(z − α)
3.4. Cauchy’s Integral Formula II 85
f (z)
The integrand of (3.6) becomes , which is bounded as z ∈ γ is away
(z − α − w)(z − α)
from α and α + w when w is small, and that the holomorphic function f is bounded on γ
by Extreme-Value Theorem. The length of γ is also bounded. Using Lebesgue Dominated
Covergence Theorem (commonly called LDCT in short), we can switch the limit and the
integral sign of (3.6), and get:
I
1 1 1 1
f 0 (α) = lim f (z) · − dz
2πi γ w→0 w z−α−w z−α
I
1 f (z)
= lim dz
2πi γ w→0 (z − α − w)(z − α)
I
1 f (z)
= dz,
2πi γ (z − α)2
proving Theorem 3.13 when n = 1.
The second and higher order cases of Theorem 3.13 can be proved by induction.
Assume the theorem holds for some integer n:
I
n! f (z)
f (n) (α) = dz
2πi γ (z − α)n+1
for any α enclosed by γ. When w is very small, α + w is also enclosed by γ, hence it is
also true that: I
n! f (z)
f (n) (α + w) = dz.
2πi γ (z − α − w)n+1
Our next goal is to determine f (n+1) (α) from the definition:
f (n) (α + w) − f (n) (α)
f (n+1) (α) = lim .
w→0 w
We leave it as an exercise:
Exercise 3.23. Follow the outline listed below, and complete the inductive proof of
Theorem 3.13:
(a) Show that:
1 1 1
−
w (z − α − w)n+1 (z − α)n+1
Xn
1 1
=
(z − α − w)(z − α) j=0
(z − α − w)j (z − α)n−j
(b) Using the induction assumption and LDCT, show that
f (n+1) (α)
I Xn
n! f (z) 1
= lim dz.
2πi γ w→0 (z − α − w)(z − α) j=0 (z − α − w)j (z − α)n−j
(c) Finally, complete the proof.
Theorem 3.15 (Liouville’s Theorem). Any bounded entire function must be constant.
86 3. Contour Integrals
Proof #1. The proof is a consequence of 1st-order Cauchy’s integral formula. Suppose
f : C → C is a entire function and that there exists M > 0 such that |f (z)| ≤ M for any
z ∈ C.
Take an arbitrary α ∈ C, and consider the contour |z − α| = R. By Theorem 3.13
with n = 1, we know: I
1 f (z)
f 0 (α) = dz.
2πi |z−α|=R (z − α)2
Then on the contour, we have:
f (z) M
≤ 2,
(z − α)2 R
and by Lemma 3.6, we can estimate that:
I
f (z) M 2πM
2
dz ≤ 2πR · 2 = .
|z−α|=R (z − α) R R
Therefore, we have for any α ∈ C and R > 0:
I
0 1 f (z) 1 2πM
|f (α)| = 2
dz ≤ · → 0 as R → +∞.
2πi |z−α|=R (z − α) 2π R
This shows f 0 ≡ 0, and hence f is a constant function.
Exercise 3.24. Prove the following general version of Liouville’s Theorem: Suppose
f : C → C is an entire function, and there exists M > 0 and a nonnegative integer k
such that:
k
|f (z)| ≤ M |z| for any z ∈ C.
Show that f is a polynomial of degree at most k.
Liouville’s Theorem is a “luxury" for holomorphic functions. There are many non-
constant bounded functions f : R → R that are (real) differentiable everywhere, while
Liouville’s Theorem says there is no non-constant bounded functions f : C → C which
are complex differentiable everywhere.
1
Proof #1. We prove by contradiction. If p(z) has no root, then p(z) is an entire function.
1 1
Note that |p(z)| → ∞ as |z| → ∞, we have: p(z) → 0 as |z| → ∞. In particular, p(z) is
1
bounded. By Liouville’s Theorem, p(z) is constant, which is a contradiction.
1 1
Proof #2. Suppose p(z) has no root, so that p(z) is an entire function. Since p(z) →0
1 1
as |z| → ∞, there exists R > 0 such that p(z) < p(0) whenever |z| ≥ R. Therefore,
1
supz∈C p(z) is achieved at some point in BR (0). By Maximum Principle (Proposition
1
3.11), the function p(z) must be a constant, which is a contradiction.
0
Proof #3. Suppose p(z) has no root, then pp(z)
(z)
is an entire function. For any R > 0, we
have I 0
p (z) n
− dz = 0 − 2πni = −2πni.
|z|=R p(z) z
p0 (z) n
On the other hand, one can easily compute that p(z) − z is a rational function of the
deg≤n−1
form deg=n+1 . Hence, by Exercise 3.18, we have
I 0
p (z) n
− dz = 0 for large R.
|z|=R p(z) z
It leads to a contradiction.
Remark 3.17. There are many proofs of Fundamental Theorem of Algebra, at least one
in almost all important fields in mathematics. We will give one more using Rouche’s
Theorem later. There is one in Topology using the concept of homotopy. There is even
one geometric proof using Gauss-Bonnet’s Theorem in Differential Geometry! Ironically,
despite the name of the theorem, a purely algebraic proof has not yet been found. The
most purest algebraic proof uses Galois Theory, but that proof is based on the fact that
88 3. Contour Integrals
every real number has a real cubic root (which has to be justified using Intermediate-Value
Theorem in Real Analysis).
Exercise 3.26. Using Liouville’s Theorem, show that if the image of an entire
function f : C → C is disjoint from an open ball Bδ (z0 ), then f is a constant
function.
The Liouville’s Theorem implies that if f : C → C is entire and the image omits a
set with non-empty interior, then it must be a constant. In fact, there is a much stronger
result called the Picard’s Little Theorem, which says that a non-constant entire function
f : C → C must have its image omitting at most 1 point in C.
3.5. Morera’s Theorem 89
Remark 3.19. Although convexity of the domain is needed in Step 2 of the proof of
Cauchy-Goursat’s Theorem, we do not need to assume Ω is convex when using Morera’s
Theorem. It is because complex differentiability is a local property. One can first restrict f
on an open ball Bε (z0 ) which is convex, then prove f is holomorphic on Bε (z0 ). Simply
repeat the same argument on all other open balls in the domain. It will show f is
holomorphic on the whole Ω.
I
In practice, it seems more difficult to verify f (z) dz = 0 for any triangle T than to
T
show f is holomorphic directly. Nonetheless, Morera’s Theorem can come in handy if we
want to show holomorphicity of a function which is not quite explicit. In the last chapter,
we may encounter functions defined in an integral form, such as the Gamma’s function:
Z ∞
Γ(z) = tz−1 e−t dt.
0
It is almost impossible to find an explicit, integral-free expression. Nonetheless, it is
possible toIshow it is a holomorphic function using Morera’s Theorem. The key idea is to
show that Γ(z) dz = 0 for any triangle T in the domain under consideration.
T
90 3. Contour Integrals
Solution
First we show that f is defined on Ω: for any z ∈ Ω and t ∈ [0, ∞), we have:
ezt
≤ ezt ≤ ext
t+1
(as usual, we denote z = x + yi). Note that:
Z ∞ ∞
1 xt 1
ext dt = e = − < ∞.
0 x 0 x
Z ∞ zt
e
Hence, dt is integrable.
0 t+1
It is quite difficult to find an explicit formula for f (z), let alone its derivative.
To show it is holomorphic, we are going to use Morera’s Theorem.
One requirement of applying Morera’s Theorem is that f needs to be a con-
tinuous function. To verify this, we fix z0 ∈ Ω and consider a sequence zj ∈ Ω
converging to z0 . We need to show f (zj ) → f (z0 ) as j → ∞. Since f (z) is defined
using an improper integral, we cannot use uniform convergence to justify the swap-
ping of limit and integral signs. Instead, we use LDCT. As −µ := Re(z0 ) < 0, we
can assume that for sufficiently large j, we have Re(zj ) < − µ2 < 0. Then, for any
t ∈ [0, ∞), and sufficiently large j, we have
µ
ezj t etRe(zj ) e− 2 t
≤= ≤ .
1+t 1+t 1+t
µ
− t µ
As e1+t2
≤ o(e− 4 t ) as t → +∞, and it is a continuous function on [0, ∞), there
exists a constant C > 0 such that
µ
e− 2 t µ
≤ Ce− 4 t for t ∈ [0, ∞).
1+t
µ
e− 2 t
Hence, is integrable over [0, ∞). By LDCT, we have
1+t
Z ∞ zj t Z ∞ Z ∞ z0 t
e ezj t e
lim f (zj ) = lim dt = lim dt = dt = f (z0 ).
j→+∞ j→+∞ 0 1+t 0 j→+∞ 1 + t 0 1 +t
Since this holds for any sequence zj converging to z0 , f is continuous at z0 . Also,
since z0 is arbitrary in Ω, we complete the verification that f is continuous on Ω.
Next we check the major
Z condition of Morera’s Theorem: given any triangle T
in Ω, we want to show f (z) dz = 0.
T
3.5. Morera’s Theorem 91
Z Z Z ∞
ezt
f (z) dz = dt dz
T T 0 t+1
Z ∞Z
ezt
= dz dt (Fubini’s Theorem)
0 T t+1
Z ∞
ezt
= 0 dt (since is holomorphic)
0 t+1
=0
To justify the legitimacy of using Fubini’s Theorem, the integral
Z Z ∞
ezt
dt |dz|
T 0 t+1
Z ∞
ezt 1
needs to be finite. To verify this, we consider dt ≤ − , so that
0 t + 1 x
Z Z ∞ Z
ezt 1
dt |dz| ≤ − |dz|,
T 0 t+1 T x
which is finite since x is away from 0 when z is on any triangle T ⊂ Ω.
Hence, by Morera’s Theorem, f is holomorphic on Ω.
Solution
First we show that Γ is well-defined (i.e. the integral is finite) on Ω by showing that
both Z 1 Z ∞
tz−1 e−t dt and tz−1 e−t dt
0 1
are integrable if Re(z) > 0. For the first integral, we observe that tz−1 e−t =
tx−1 e−t ≤ tx−1 for any t ∈ [0, 1]. Whenever x > 0, we have
Z 1
1
tx−1 dt = < ∞,
0 x
and hence the first integral is integrable if x := Re(z) > 0. For the second integral,
we observe that:
tz−1 e−t
lim = lim tx−1 e−t/2 = 0.
t→+∞ e−t/2 t→+∞
−t/2
Since e is integrable on [1, ∞), we conclude that the second integral is also
integrable on [1, ∞). This shows Γ is well-defined on Re(z) > 0. Also, since
Z ∞
fn (z) := χ[ n1 ,∞) tz−1 e−t dt
0
where χ[ n1 ,∞) is the indicator function of the set [ n1 , ∞). By the fact that
Let ε > 0 such that Re(α) > ε. For sufficiently large j, we have |zj − α| < 2ε , and
hence Re(zj ) > 2ε and Re(zj ) < Re(α) + 2ε . In order to apply LDCT on fn (zj ), we
need to bound the integrand. For any t ∈ [ n1 , ∞), we have for any large j, and
t ∈ [0, ∞):
ε ε
tzj −1 e−t = χ[ n1 ,1) txj −1 e−t + χ[1,∞) txj −1 e−t ≤ χ(0,1) t 2 −1 + χ[1,∞) tRe(α)+ 2 −1 e−t .
| {z }
=:g(t)
ε
2 −1 Re(α)+ 2ε −1
As χ(0,1) t is integrable on (0, 1) and χ[1,∞) t e−t is integrable on [1, ∞),
the function g(t) is integrable on (0, +∞). By LDCT, we conclude:
Z ∞ Z ∞
lim fn (zj ) = lim tzj −1 e−t dt = lim tzj −1 e−t dt = fn (α).
j→∞ j→∞ 1 1 j→∞
n n
By compactness of T , there exists δ, M > 0 such that x := Re(z) ∈ (δ, M ) for any
z ∈ T . Hence,
Z 1 Z 1
1 − n1δ
tx−1 e−t dt ≤ tδ−1 dt = < ∞,
1/n 1/n δ
and Z ∞ Z ∞
x−1 −t
t e dt ≤ tM −1 e−t dt < ∞.
1 1
These combine to show:
I Z ∞
tx−1 e−t dt |dz| < ∞,
T 1/n
3.5. Morera’s Theorem 93
H R∞
hence the switching of T
and 1/n
is justified. By Morera’s Theorem, fn is holo-
morphic on Ω.
Finally, by showing that fn converges uniformly to Γ on Ωε := {z : Re(z) > ε}
for any ε > 0, it will complete the proof that Γ is holomorphic on Ω. We show
uniform convergence from the definition. First we estimate |fn (z) − Γ(z)|. For any
z ∈ Ωε , we have:
Z 1/n
|fn (z) − Γ(z)| = tz−1 e−t dt
0
Z 1/n
≤ tz−1 e−t dt
0
Z 1/n Z 1/n
x−1 −t
= t e dt ≤ tx−1 dt
0 0
1 1
=< ε.
xnx εn
The last step follows from x := Re(z) > ε for any z ∈ Ωε . Taking supremum over
all z ∈ Ωε , we get
1
sup{|fn (z) − Γ(z)| : z ∈ Ωε } ≤ ε
εn
for any n ≥ 1. By squeeze theorem, we conclude:
1
lim sup{|fn (z) − Γ(z)| : z ∈ Ωε } ≤ lim = 0,
n→∞ n→∞ εnε
concluding that fn converges uniformly to Γ on Ωε for any fixed ε > 0.
By Morera’s Theorem (and its corollary), fn being holomorphic on every Ωε
S Γ is holomorphic on every Ωε (where ε > 0). This shows Γ is holomorphic
implies
on ε>0 Ωε (which is Ω).
Proof. Given any z ∈ BR (z0 ), we let ε > 0 be small enough so that |z − z0 | < R − ε. For
simplicity, denote R0 = R − ε.
95
96 4. Taylor and Laurent Series
Then, the contour |z − z0 | = R0 lies inside the open ball BR (z0 ). The key trick to prove
1
the Taylor Theorem is rewriting as a geometric series. Recall that:
ξ−z
1
= 1 + w + w2 + . . . whenever |w| < 1.
1−w
1
We first rewrite into this form:
ξ−z
1 1 1 1
= = ·
ξ−z (ξ − z0 ) − (z − z0 ) ξ − z0 1 − z−z0
ξ−z0
∞ n
1 X z − z0
=
ξ − z0 n=0 ξ − z0
z−z0
Here we have used the fact that ξ−z0 < 1. See the diagram below. The yellow ball is
BR (z0 ), and the red circle is |ξ − z0 | = R0 .
|ξ − z0 | = R0
z0
|z − z0 | < R0
z
Then, whenever z ∈ BR0 (z0 ), the function f (z) can be expressed as:
I
1 1
(4.1) f (z) = f (ξ) · dξ
2πi |ξ−z0 |=R0 ξ−z
I ∞ n
1 f (ξ) X z − z0
= dξ
2πi |ξ−z0 |=R0 ξ − z0 n=0 ξ − z0
I ∞
X
1 f (ξ) (z − z0 )n
= dξ
2πi |ξ−z0 |=R0 n=0 (ξ − z0 )n+1
I
Next we want to see whether we can switch the integral sign and the infinite
|ξ−z0 |=R0
∞
X
summation . For this we need to show uniform convergence of the series below.
n=0
X∞
f (ξ) (z − z0 )n
.
n=0
(ξ − z0 )n+1
4.1. Taylor Series 97
2
Example 4.3. The Taylor series for some composite functions, such as ez , can
be derived by substitution instead of deducing the general n-th derivative of the
function. For example:
z2 z3 z4
ez = 1 + z + + + + ...
2! 3! 4!
2 (z 2 )2 (z 2 )3 (z 2 )4
ez = 1 + z2 + + + + ...
2! 3! 4!
z4 z6 z8
= 1 + z2 + + + + ...
2! 3! 4!
2
Since the series for ez converges for any z ∈ C, the series for ez converges for any
z ∈ C as well.
Similarly, by replacing z by 1 − z in the Taylor series for Log(z), we get:
1 1 1
Log(1 − z) = −z − z 2 − z 3 − z 4 − . . .
2 3 4
The series for Log(z) about 1 converges when |z − 1| < 1, and so the above series
for Log(1 − z) converges when |(1 − z) − 1| < 1, i.e. |z| < 1.
Apart from using Theorem 4.1 to find the Taylor series of a given holomorphic
function, we can also make use of the geometric series formula directly:
1
= 1 + w + w2 + . . . where |w| < 1.
1−w
This method is particularly useful for functions whose n-th derivatives are tedious to
compute.
4.1. Taylor Series 99
Exercise 4.1. Derive the Taylor series of each function below about the given center
z0 :
2π
(a) f (z) = sin 2z; z0 = 3
(b) f (z) = cos 3z; z0 = π
−z 3
(c) f (z) = e ; z0 = 0
(d) f (z) = Log(3 − 2z); z0 = 1
Exercise 4.2. Find the Taylor series about 0 of the functions below up to the z 4
term:
(a) f (z) = e−z cos z
(b) f (z) = Log(1 − ez )
Exercise 4.3. Find the Taylor series about z0 of the function below without using
Theorem 4.1. State its radius of convergence.
1
(a) f (z) = , z0 = 0
(z − 1)(z − 2)
1
(b) f (z) = , z0 = i
(z − 1)(z − 2)
Exercise 4.4. Let α, β and z0 be three distinct complex numbers. Consider the
function
1
f (z) = .
(z − α)(z − β)
Find the Taylor series about z0 of the above function, and state its radius of conver-
gence.
Exercise 4.5. Let α be a fixed non-zero complex number. Consider the principal
branch of (1 + z)α :
(1 + z)α := eαLog(1+z) .
Show that its Taylor series about 0 is given by:
X∞
α(α − 1) · (α − n + 1) n
(1 + z)α = 1 + z .
n=1
n!
State its radius of convergence.
4.1.2. Taylor Series with Remainder Term. In Real Analysis, the Taylor Theorem
with a remainder term asserts that for any smooth (C ∞ ) function f : R → R, we have:
N
X −1 Z x
f (n) (a) 1
f (x) = (x − a)n + (x − t)N −1 f (N ) (t) dt
n! (N − 1)! a
n=0 | {z }
=:RN (x)
The last integral term, commonly denoted as RN (x), measures how fast the Taylor
series converges to f (x) as N → ∞. If lim RN (x) → 0 for any x in an interval I, then
N →∞
X∞
f (n) (a)
the Taylor series (x − a)n converges (pointwise) to f (x) for any x ∈ I. If
n=0
n!
4.1. Taylor Series 101
furthermore, we have:
lim sup |RN (x)| → 0,
N →∞ x∈I
then the Taylor series converges uniformly to f on I. However, it is often not easy to
show RN → 0 as the N -th derivative f (N ) may not be easy to find.
Back to Complex Analysis, we will soon derive the remainder term for the Taylor
series for holomorphic functions. One good thing about the complex version is that the
remainder involves only f , but not its derivatives, making it much easier to handle the
convergence issue of complex Taylor series. It again thanks to Cauchy’s integral formula.
Proposition 4.2. Let f be a holomorphic function defined on BR (z0 ), then for any
z ∈ BR (z0 ), and any simple closed curve γ in BR (z0 ) enclosing both z and z0 , we have:
N
X −1 (n) I N
f (z0 ) n 1 f (ξ) z − z0
f (z) = (z − z0 ) + dξ
n=0
n! 2πi γ ξ − z ξ − z0
| {z }
=:RN (z)
Proof. We only outline the proof since it is modified from the proof of Theorem 4.1.
Using Cauchy’s integral formula, we first have:
I
1 f (ξ)
f (z) = dξ.
2πi γ ξ − z
The key step in the proof of Theorem 4.1 is to write:
1 1 1 1
= = · ,
ξ−z (ξ − z0 ) − (z − z0 ) ξ − z0 1 − z−z
ξ−z0
0
z − z0
so that when < 1, we have:
ξ − z0
∞
X n
1 z − z0
z−z0 = .
1− ξ−z0 n=0
ξ − z0
Now, to prove this proposition, we modify the above key step a bit, by considering:
N
1 − z−z 0 N
X −1 n
ξ−z0 z − z0
= .
1 − z−z
ξ−z0
0
n=0
ξ − z0
We leave the rest of the proof for readers (which is a good exercise to test your under-
standing of the proof of Theorem 4.1).
Exercise 4.8. Let f be a holomorphic function on BR (z0 ). Using this estimate ob-
X∞
f (n) (z0 )
tained in Exercise 4.7, deduce that the Taylor series (z − z0 )n converges
n=0
n!
uniformly to f (z) on any smaller ball Br (z0 ) where 0 < r < R.
X∞
f (n) (z0 )
Remark 4.3. The uniform convergence of (z − z0 )n has many remarkable
n=0
n!
consequences as discussed in MATH 3033/2043. For instance, one can integrate a Taylor
series term-by-term.
4.2.1. Examples of Laurent Series. While a Taylor series gives an analytic expres-
sion for a holomorphic function on a ball, a Laurent series gives an analytic expression
for a function that has a singularity at the center of a ball. Before we discuss a general
theorem about Laurent series, let’s first look at some examples of writing a function as a
Laurent series:
Example 4.5. Consider the function f : C\{1, 2} → C defined by:
1
f (z) = .
(z − 1)(z − 2)
It is holomorphic on its domain C\{1, 2}. Let’s express the above function as a
Laurent series about 1:
1 1 1
= =−
z−2 (z − 1) − 1 1 − (z − 1)
X∞
=− (z − 1)n where |z − 1| < 1.
n=0
Hence, on ∈ A1,0 (1), i.e. the green annulus in the figure below, we have:
∞ ∞
1 1 1 X X
f (z) = · =− (z − 1)n = − (z − 1)n−1
z−1 z−2 z − 1 n=0 n=0
1
=− − 1 − (z − 1) − (z − 1)2 + · · · .
z−1
104 4. Taylor and Laurent Series
x
1 2
However, the green annulus A1,0 (1) is not the only annulus centered at 1 on
which f is holomorphic. There is another one A∞,1 (1) = {1 < |z − 1|} centered at
1, i.e. the yellow annulus in the above figure, on which f is also holomorphic. It is
also possible to express f as a Laurent series on this yellow annulus:
1 1 1 1
= = · 1
z−2 (z − 1) − 1 z − 1 1 − z−1
∞ n
1 X 1
= (where |z − 1| > 1)
z − 1 n=0 z − 1
∞
X 1
=
n=0
(z − 1)n+1
Hence, on the yellow annulus A∞,1 (1), the function f can be expressed as the
following Laurent series:
∞
X
1 1 1 1
f (z) = = · = .
(z − 1)(z − 2) z − 1 z − 2 n=0 (z − 1)n+2
Solution
First recall that the Taylor series for ew is:
X∞
wn
ew = for any w ∈ C.
n=0
n!
Substitute w = z1 , where z 6= 0, we get:
∞
X
1 1
ez = ,
n=0
n!z n
4.2. Laurent Series 105
and hence:
1
f (z) = z 2 e z
∞
X 1
= z2
n=0
n!z n
∞
X 1
= n−2
n=0
n!z
1 1 1
= z2 + z + + + + ···
2 3!z 4!z 2
Exercise 4.12. Find all possible Laurent (or Taylor) series about 1 for the function:
1
f (z) = 2 .
z − 2z
For each series, state the annulus or ball on which it converges.
Exercise 4.13. Find all possible Laurent (or Taylor) series about each z0 below for
1
the function f (z) = .
z
(a) z0 = 0
(b) z0 = 1
(c) z0 = i
For each series, state the annulus or ball on which it converges.
Exercise 4.14. Show that for any w such that |w| < 1, we have:
X∞
1 n(n − 1) n−2
3
= (−1)n w .
(1 + w) n=2
2
[Hint: use Exercise 4.5]
Hence, find all possible Laurent or Taylor series about i for the function:
1
f (z) = 3 .
z
For each series, state the annulus or ball on which it converges.
106 4. Taylor and Laurent Series
Exercise 4.15. Find the Laurent series about 1 on the annulus A∞,0 (1) for the
functions:
1 1
f (z) = sin and g(z) = cos .
z−1 z−1
Hence, find the Laurent series about 1 on A∞,0 (1) for:
z
h(z) = sin .
z−1
Proof. The proof is similar to that of Taylor’s series, but is a bit trickier since an annulus
is not simply-connected and so Cauchy’s integral formula cannot be applied directly.
Fix z ∈ AR,r (z0 ), we first consider a simple closed curve Γ in AR,r (z0 ) which encloses
both z and z0 (just like in the proof of Taylor’s Theorem). However, we cannot apply
Cauchy’s integral formula on the integral:
I
1 f (ξ)
dξ
2πi Γ ξ − z0
C =Γ+L−γ−L
where −γ is the clockwise circle, and L is a straight-path as shown in the figure below.
We can pick Γ to be the circle with radius slightly smaller than R, and γ with radius
slightly bigger than r so that C encloses z.
4.2. Laurent Series 107
Γ
R z
r
z0
−γ
L −L
The key idea of the proof is to express the integral over Γ as a series of non-negative
powers, and the integral over γ as a series of negative powers.
When ξ ∈ Γ, we have |z − z0 | < |ξ − z0 |, so:
∞ n
1 1 1 1 1 X z − z0
= = · =
ξ−z (ξ − z0 ) − (z − z0 ) ξ − z0 1 − z−z
ξ−z0
0
ξ − z0 n=0 ξ − z0
In order to switch the infinite summation and the integral sign, we justify that the series
converges uniformly on ξ ∈ Γ. Suppose Γ has radius R0 , then for any ξ ∈ Γ:
n n
f (ξ) z − z0 1 |z − z0 |
≤ 0 sup |f | .
ξ − z0 ξ − z0 R R0 Γ
Note that supΓ |f | is finite by compactness of Γ. Since |z − z0 | < R0 , the geometric series
X∞ n
1 |z − z0 |
sup |f |
n=0
R0 R0 Γ
converges uniformly on ξ ∈ Γ, so one can switch the summation and integral signs and
get:
I X∞ I
1 f (ξ) 1 f (ξ)
(4.3) dξ = n+1
dξ (z − z0 )n .
2πi Γ ξ−z n=0
2πi Γ (ξ − z0 )
The second integral can be handled similarly. The difference is that when ξ ∈ γ, we
have |ξ − z0 | < |z − z0 | instead. We instead write:
∞ n
1 1 1 1 1 X ξ − z0
= =− · ξ−z0
= −
ξ−z (ξ − z0 ) − (z − z0 ) z − z0 1 − z−z z − z0 n=0 z − z0
0
Hence,
I I X ∞ n
1 f (ξ) 1 f (ξ) ξ − z0
dξ = − dξ.
2πi γ ξ−z 2πi γ n=0 z − z0 z − z0
We leave it as an exercise for readers to argue that the series converges uniformly on
ξ ∈ γ so that we can switch the integral and summations signs:
I X∞ I
1 f (ξ) 1 n 1
(4.4) dξ = − f (ξ)(ξ − z0 ) dξ n+1
.
2πi γ ξ−z n=0
2πi γ (z − z0)
X∞ I
1 1
f (z) = f (ξ)(ξ − z0 )n−1 dξ
n=1
2πi γ (z − z0 )n
| {z }
(4.4)
X∞ I
1 f (ξ)
+ n+1
dξ (z − z0 )n
n=0
2πi Γ (ξ − z0 )
| {z }
(4.3)
Exercise 4.17. Justify the claim in the above proof that the series:
∞
X n
f (ξ) ξ − z0
n=0
z − z0 z − z0
converges uniformly on ξ ∈ γ (and z, z0 are considered to be fixed).
Remark 4.5. Although from the proof of Theorem 4.4 one can express the coefficient cn ’s
of a Laurent series in terms of contour integrals, we do not usually find the coefficients
this way since these contour integrals may not be easy to compute.
4.2. Laurent Series 109
4.2.3. Laurent Series with Remainders. Similar to Taylor series, one can refine
Theorem 4.4 a bit by deriving the remainder terms. Using the remainder terms, one can
argue that for a holomorphic function f defined on an annulus AR,r (z0 ), the Laurent
series converges uniformly to f on every smaller annulus AR0 ,r0 (z0 ) (where r < r0 < R0 <
R). This result is remarkable as it allows us to integrate a Laurent’s series term-by-term.
Proposition 4.6. Let f be a holomorphic function on the annulus AR,r (z0 ), where 0 ≤
r < R ≤ ∞. Then, for each positive integer N and z ∈ AR,r (z0 ), we have:
XN I I N
1 1 1 f (ξ) ξ − z0
f (z) = f (ξ)(ξ − z0 )n−1 dξ + dξ
n=1
2πi γ (z − z0 )n 2πi γ z − ξ z − z0
| {z }
=:rN (z)
N
X −1 I I N
1 f (ξ) 1 f (ξ) z − z0
+ dξ (z − z0 )n + dξ
2πi Γ (ξ − z0 )n+1 2πi Γ ξ−z ξ − z0
n=0 | {z }
=:RN (z)
where Γ and γ are any pair of circles in AR,r (z0 ) centered at z0 such that z is bounded
between Γ and γ.
Proof. We leave the proof of Proposition 4.6 as an exercise. It is very similar to the
proof of Proposition 4.2 for Taylor series. Readers should first digest the whole proof of
Proposition 4.2, then write up a coherent proof for this proposition.
Exercise 4.18. Prove Proposition 4.6. Using this, show that the Laurent series
about z0 for f converges uniformly on every smaller annulus AR0 ,r0 (z0 ) where
r < r0 < R0 < R. Give the proof in two different ways:
(a) Using Weierstrass’ M-test
(b) By estimating the reminder terms RN and rN directly.
1
For the function g(z) = , all singularities are isolated as depicted in the diagram
sin z
below:
y
1
0 is a pole of order 1 for the function . Hence, one can express this function as a
sin z
Laurent series on a small annulus Aε,0 (0):
1 c−1
= + c0 + c1 z + c2 z 2 + · · ·
sin z z
Multiplying z on both sides, we get:
z
= c−1 + c0 z + c1 z 2 + c2 z 3 + · · ·
sin z
and by letting z → 0, we can also conclude that c−1 = 1. Therefore, if γ is a simple close
curve enclosing 0 in this small annulus Aε,0 (0), then we have:
I I I
1 1 1
dz = + c0 + c1 z + c2 z 2 + · · · dz = dz + 0 + 0 + · · · = 2πi.
γ sin z γ z γ z
Exercise 4.19. Find all isolated singularities of each function below, and classify
the nature of these singularities. For poles, state also their orders.
ez − 1
(a) f (z) =
z
Log(z)
(b) g(z) =
(z − 3)5
(c) h(z) = z 4023 cos z1
112 4. Taylor and Laurent Series
Definition 4.8 (Residues). Let z0 be an isolated singularity of f (z) such that the Laurent
series about z0 for f on some annulus Aε,0 (z0 ) is given by:
∞
X
f (z) = cn (z − z0 )n ,
n=−∞
Solution
The denominator has roots −1, 2i and −2i, hence they are isolated singularities of
f . In this solution, we will decompose f (z) into partial fractions. It may not be a
pleasant way finding residues, but we will later provide an easier way.
Note that f is a rational function, we can break it into partial fractions:
A B C D
f (z) = + + + .
(z + 1)2 z + 1 z − 2i z + 2i
We leave it as an exercise for readers to determine the value of A, B, C and D. One
should be able to get:
3
5 − 14
15
7+i 7−i
f (z) = + + 25 + 25 .
(z + 1)2 z + 1 z − 2i z + 2i
| {z }
holomorphic near −1
7+i 7−i
25 25
On a small annulus Aε,0 (−1) about −1, the last two terms are +
z − 2i z + 2i
holomorphic. Therefore, if one express them as a Laurent series about −1, only
1
non-negative powers of z + 1 will appear, and the coefficient of z+1 will not be
affected. Hence, we have:
14
Res(f, −1) = − .
15
By a similar reason, we have:
7+i 7−i
Res(f, 2i) = and Res(f, −2i) = .
25 25
4.3. Residue Calculus 113
It is no doubt that partial fraction decompositions are time-consuming and not fun
(it may remind you the computational nightmare you might have encountered in MATH
1014). Fortunately, there is a better way for finding residues for poles (does not work for
essential singularity).
If we know already that z0 is a pole of order 1 (i.e. simple pole) of a function f (z),
then
c−1
f (z) = + c0 + c1 (z − z0 ) + c2 (z − z0 )2 + · · ·
z − z0
It is then easy to see that
c−1 = lim (z − z0 )f (z).
z→z0
Therefore, in order to find Res(f, z0 ) for a simple pole z0 , we simply need to compute the
above limit.
Now consider the case if z0 is a pole of order k for f , then its Laurent series about z0
is given by:
c−k c−(k−1) c−1
f (z) = k
+ k−1
+ ··· + + c0 + c1 (z − z0 ) + c2 (z − z0 )2 + · · ·
(z − z0 ) (z − z0 ) z − z0
Our goal is to find c−1 . By multiplying both sides by (z − z0 )k , we can get:
(z − z0 )k f (z) = c−k + c−(k−1) (z − z0 ) + · · · + c−1 (z − z0 )k−1 + c0 (z − z0 )k + · · ·
By differentiating both sides for k − 1 times, all terms involving (z − z0 )n with n < k − 1
will disappear:
dk−1
(z − z0 )k f (z) = c−1 (k − 1)! + e c1 (z − z0 )2 + · · ·
c0 (z − z0 ) + e
dz k−1
dk−1
We have used the fact that k−1 (z − z0 )k−1 = (k − 1)!, and e c0 , e
c1 , . . . are some complex
dz
numbers (which we do not need to know their values).
By letting z → z0 , we get:
dk−1
lim (z − z0 )k f (z) = c−1 (k − 1)!
z→z0 dz k−1
which provides a good way to find c−1 without expanding a Laurent series:
Proposition 4.9. Suppose z0 is a pole of order k < ∞ for a function f , then we have:
1 dk−1
Res(f, z0 ) = lim (z − z0 )k f (z).
(k − 1)! z→z0 dz k−1
In particular, for a simple pole z0 , we have:
Res(f, z0 ) = lim (z − z0 )f (z).
z→z0
Solution
As discussed before, the isolated singularities are −1, 2i and −2i. Observe that:
z 2 − 2z 3
lim (z + 1)2 f (z) = lim 2
= 6= 0.
z→−1 z→−1 z + 4 5
Hence −1 is a pole of order 2. From Proposition 4.9, we have:
1 d2−1
Res(f, −1) = lim (z + 1)2 f (z)
(2 − 1)! z→−1 dz 2−1
d z 2 − 2z
= lim
z→−1 dz z 2 + 4
2z 2 + 8z − 8 14
= lim =− .
z→−1 (z 2 + 4)2 25
Both 2i and −2i are simple poles, so we have:
z 2 − 2z 7+i
Res(f, 2i) = lim (z − 2i)f (z) = lim =
z→2i z→2i (z + 1)2 (z + 2i) 25
2
z − 2z 7−i
Res(f, −2i) = lim (z + 2i)f (z) = lim =
z→−2i z→−2i (z + 1)2 (z − 2i) 25
Solution
For each function, we first determine whether 0 is a pole, and find out its order.
For f (z), we consider:
z
lim zf (z) = lim · ez = 1 · e0 = 1 6= 0.
z→0 sin z
z→0
Hence 0 is a simple pole for f , and Res(f, 0) = 1.
For g(z), note that:
z2 z3 z z2
z+ 2! + 3! + ··· 1+ 2! + 3! + ···
lim g(z) = lim z3 z5
= lim z2 z4
= 1 < ∞.
z→0 z→0 z− 3! + 5! − ··· z→0 1− 3! + 5! − ···
Hence 0 is a removable singularity of g(z), and there is no z1 -term in the Laurent
series, and so Res(g, 0) = 0.
For h(z):
z 2
lim z 2 h(z) = lim ez = 1 6= 0.
z→0 z→0 sin z
4.3. Residue Calculus 115
Exercise 4.21. For each function below, find its residue at each isolated singularity
using any method:
z2 − 1 1 z 1997 − 1
z 3 (z 2 + 1) 2
6z + 8z + 9 z 2047 − 1
1 e2zi e2zi − 1
z
e −1 sin z sin z
1 z2 sin z
z sin z e1/z z 2 (z − π)3
where each γj is a small simple-closed curve so that zj is the only singularity it encloses.
Then, each γj -integral is given by 2πi Res(f, zj ), and hence we have the following
theorem:
116 4. Taylor and Laurent Series
Proof. Let ε > 0 be sufficiently such that each circle {|z − zj | = ε}, denoted by γj ,
encloses zj as the only singularity of f (see figure below).
z1
γ1
γ2 γ3
z2 z3
Each γj encloses zj as the only singularity of f . Express f as a Laurent series on Aε,0 (zj ):
∞
X
f (z) = cn (z − zj )n .
n=−∞
I
Recall that (z − zj )n dz 6= 0 only when n = −1, and by uniform convergence of
γj
Laurent series, we get:
I
f (z) = 2πic−1 = 2πi Res(f, zj ).
γj
Therefore, we have:
I N
X
f (z) dz = 2πi Res(f, zj ),
γ j=1
completing the proof.
4.3. Residue Calculus 117
Solution
z 2 − 2z
Denote f (z) = . The singularities of f are −1, 2i and −2i. We have
(z + 1)2 (z 2 + 4)
calculated in Example 4.8 that
14 7+i 7−i
Res(f, −1) = − Res(f, 2i) = Res(f, −2i) =
15 25 25
y
2i
−2i
(a) When 0 < R < 1, the circle |z| = R does not enclose any singularities, hence
I
f (z) dz = 0.
|z|=R
(b) When 1 < R < 1, the circle |z| = R encloses the singularity −1 only, hence
I
28πi
f (z) dz = 2πi Res(f, −1) = −
|z|=R 15
(c) When R > 2, the circle |z| = R encloses all three singularities, hence
I
f (z) dz = 2πi (Res(f, −1) + Res(f, 2i) + Res(f, −2i))
|z|=R
14 7 + i 7 − i 56πi
= 2πi − + + =− .
15 25 25 75
118 4. Taylor and Laurent Series
Example 4.11. Let N be a positive integer and γN be the square contour with
vertices ±(N + 21 ) ± (N + 12 )i. Use Residue Theorem to show:
XN I
1 π2 1 π
2 2
= + cot πz dz.
n=1
n 3 2πi γN z 2
Hence, deduce that:
1 1 π2
1+ + + · · · = .
22 32 6
1 1
The problem about the exact value of 1 + 22 + 32 + · · · was famously known as
“Basel’s Problem”. It was first solved by Euler in 1734.
Solution
π π cos πz
Denote f (z) := cot πz = 2 . Its singularities are the set of all integers n.
z2 z sin πz
First we observe that
πz
lim z 3 f (z) = lim · cos πz = 1,
z→0 z→0 sin πz
so 0 is a pole of order 3 for f . By Proposition 4.9, its residue is given by:
1 d2 1 d2
Res(f, 0) = lim 2 z 3 f (z) = lim 2 πz cot πz
2! z→0 dz 2 z→0 dz
π 2 (πz cos πz − sin πz)
= lim
z→0
sin3 πz
π2 z2 π3 z3
π 2 πz(1 − 2! + · · · ) − (πz − 3! + ···)
= lim
z→0 sin3 πz
3 3
− π 3z + higher-order terms π2
= lim π 2 · 3 =− .
z→0 sin πz 3
For any non-zero integer n, observe that:
π(z − n) cos πz (−1)n 1
lim (z − n)f (z) = lim 2 n = 2 = 2.
z→n z→n z · (−1) sin(π(z − n)) n · (−1)n n
| {z }
=sin πz
1
Hence, n is a simple pole of f (for any n 6= 0), and Res(f, n) = .
n2
Now consider the contour γN . The singularities it encloses are:
0, ±1, ±2 · · · , ±N.
4.3. Residue Calculus 119
y
(N + 21 )i
x
−(N + 12 ) N+ 1
2
−(N + 21 )i
Exercise 4.24. Show, using residues, that for any positive integer k, the Riemann’s
zeta function takes the value:
(−1)k+1 (2π)2k
ζ(2k) = · B2k
2 · (2k)!
where B2k is the 2k-th Bernoulli’s number defined as follows:
The Bernoulli’s numbers (where n ∈ N) are defined to be coefficient of the series:
X∞
z Bn n
= z .
ez − 1 n=0 n!
z 2 3
Note that has a removable singularity at 0 (since ez − 1 = z + z2! + z3! + · · · ).
ez − 1
Therefore, its Laurent series about 0 is indeed a Taylor’s series. The general formula
of these Bernoulli’s numbers can hardly be found, but by direct computations, one
can find:
1 1 1 1
B0 = 1 B1 = − B2 = B3 = 0 B4 = − B5 = 0 B6 = ···
2 6 30 42
π
Hint: consider the function 2k cot πz, and mimic the previous example (which
z
is the case of ζ(2)).
Exercise 4.25. Use Residue Theorem to evaluate the following contour integrals:
I
1
(a) 2+1
dz
|z|=3 z
I
z 3 + 3z + 1
(b) 4 2
dz
|z|=2 z − 5z
I
ez + z
(c) 4
dz
|z−i|=2 (z − 1)
I
sin z
(d) 4023
dz
|z−i|=2 (z − i)
I
(e) tan πz dz where γ is the rectangle contour with vertices:
γ
Exercise 4.26. Let γN be the square contour with vertices ±(N + 21 )π ± (N + 12 )πi
where N is a positive integer. Show that:
I N
1 1 1 2 X (−1)n
csc z dz = + .
2πi γN z 2 6 π 2 n=1 n2
Hence, show that:
1 1 1 π2
1− + − + · · · = .
22 32 42 12
4.3. Residue Calculus 121
Exercise 4.27. Determine the residues of all isolated singularities of the function:
1
f (z) = .
(2z − 1) sin πz
By considering a suitable contour integral of f , show that:
1 1 1 π
1 − + − + ··· = .
3 5 7 4
4.3.4. Evaluation of Real Integrals. Residues are often used to evaluate some
difficult real integrals that physicists and engineers may encounter.
Solution
The key trick is to express the real integral as a complex integral of the circle contour
|z| = 1, which is parametrized by z = eiθ where 0 ≤ θ ≤ 2π.
When z = eiθ is on the contour {|z| = 1}, we have
eiθ + e−iθ 1 1
cos θ = = z+
2 2 z
1 1
dz = ieiθ dθ =⇒ dθ = iθ dz = dz
ie iz
Therefore, the real integral can be written as a complex integral as:
Z 2π I I
1 1 1 1
dθ = b 1
· dz = 2i 2 − 2az + b
dz.
0 a − b cos θ |z|=1 a − 2 z + z
iz |z|=1 bz
We can then use residue theory to evaluate the complex integral. The singularities
of the integrand are roots of the quadratic equation bz 2 − 2az + b = 0, which are:
√ √
a − a 2 − b2 a + a2 − b2
ω1 = and ω2 = .
b b
Note that a > b, so both roots are real. We further observe that:
a+0
ω2 > > 1 and ω1 ω2 = 1,
b
and so |ω1 | < 1. Therefore, ω1 is the only singularity enclosed by the contour |z| = 1.
As ω1 and ω2 are distinct, they are simple poles, and so the contour integral is given
by:
I I
1 1
2
dz = dz
|z|=1 bz − 2az + b |z|=1 b(z − ω1 )(z − ω2 )
I 1
b(z−ω2 ) 1
= dz = 2πi
|z|=1 z − ω1 b(z − ω2 ) z=ω1
2πi πi
= = −√ .
b(ω1 − ω2 ) a2 − b2
Hence, the real integral is given by:
Z 2π
1 πi 2π
dθ = −2i · √ =√ .
0 a − b cos θ a2 − b2 a2 − b2
122 4. Taylor and Laurent Series
Before we proceed to the next example, let’s first prove a useful observation which
will come in handy later on.
Exercise 4.28. Show that the function eiz is bounded on the upper-half plane, i.e.
there exists M > 0 such that eiz ≤ M whenever Im(z) ≥ 0. On the other hand,
show that the function = cos z is unbounded on the upper-half plane.
Solution
Let’s consider the following semi-circle contour:
y
CR
i
x
−R LR R
CR
i
Cε
x
−R −ε ε R
Example 4.14. Let α be a real constant in (0, 1). Evaluate the real integral:
Z ∞
1
I := dx.
0 x (1 + x2 )
α
Solution
First observe that for any x > 0, we have:
xα = eα ln x = eαLog−π/2 (x)
where Log−π/2 is the special branch of logarithm defined in (4.5). It prompts us to
consider a contour integral of the function:
1
f (z) = αLog (z)
.
e −π/2 (1 + z 2 )
We pick a contour as shown in Figure 4.1, where CR and Cε are semi-circles with
radii R and ε respectively. Since the closed contour γR,ε := [−R, −ε] + Cε + [ε, R] +
CR lies completely inside the domain of Log−π/2 (z), by Residue Theorem, we have:
I
1 1
αLog−π/2 (z)
dz = 2πi Res(f, i) = 2πi · αLog
2 −π/2 (z) (z + i)
γR,ε e (1 + z ) e z=i
2πi π
= π = απ i .
2ieα(ln|i|+ 2 i) e2
On the other hand, the γR,ε -integral can break down into:
(4.6)
I Z −ε Z Z R Z !
1 1
αLog−π/2 (z)
dz = + + + αLog−π/2 (z)
dz
2 (1 + z 2 )
γR,ε e (1 + z ) −R Cε ε CR e
When z = x + 0i ∈ [ε, R], the integrand is simply:
1 1
= α .
e αLog−π/2 (x) 2
(1 + x ) x (1 + x2 )
Hence,
Z R Z ∞
1 1
lim lim dz = dx =: I.
ε→0 R→∞ ε eαLog−π/2 (z) (1 + z 2 ) 0 xα (1 + x2 )
4.3. Residue Calculus 125
We are left to analyze the two semi-circular integrals. We will show that they
tend to 0 as ε → 0 and R → ∞.
When z ∈ Cε , we have:
1 1 e−α ln ε ε−α
≤ e−α(ln|z|+iθ(z)) · = = 2
= .
e αLog−π/2 (z)
(1 + z 2 ) 1 − |z|
2 1−ε 1 − ε2
1 1 R−α
≤ e−α(ln|z|+iθ(z)) · = .
eαLog−π/2 (z) (1 + z 2 ) 1 − |z|
2 R2 − 1
Exercise 4.32. Evaluate the following real integrals using residue methods:
Z 2π
1
(a) dθ where a > b > 0.
0 (a + b cos θ)2
126 4. Taylor and Laurent Series
Z 2π
1
(b) dθ where a ∈ R and a 6= ±1.
0 1 − 2a cos θ + a2
Z ∞
x2
(c) dx where a > 0.
(x + a2 )2
2
Z0 ∞
1
(d) dx where n ∈ N.
(x2 + 1)n
Z0 ∞
cos ax
(e) dx where a and b are positive real numbers
x 2 + b2
Z0 ∞
sin ax
(f) dx where a is a positive real number.
x(x2 + 1)
Z0 ∞
ln x
(g) dx where a > 0.
x2 + a2
Z0 ∞
1
(h) dx where α ∈ (0, 1).
0 xα (1 + x4 )
4.3.5. Argument Principle and Rouche’s Theorem. The Argument Principle deals
with functions of the following type:
as one can parametrize f ◦ γ by w = f (γ(eit )), t ∈ [0, 2π]. This integral (up to a factor of
2πi) measures the winding number of the closed curve f ◦ γ around 0.
The Argument Principle is a nice formula relating this integral and the orders of
poles and zeros of f enclosed by γ. The order of a zero z0 of f is defined to be the order
of the pole z0 of f1 .
4.3. Residue Calculus 127
One can check that F is holomorphic on the region enclosed by γ, by showing that all
αj ’s and βk ’s are removable singularities of F :
QN
(z − βk )ord(βk )
F (z) =f (z) · Qk=1
M
j=1 (z − αj )ord(αj )
QN
(αm − βk )ord(βk ) f (z)
lim F (z) = QM k=1 · lim ord(αm )
z→αm
j=1,j6=m (αm − αj )ord(αj ) z→αm (z − αm )
QN ord(βk )
k=1,k6=n (βn − βk ) ord(βn )
lim F (z) = QM · lim f (z) · (z − βn ) .
ord(αj )
j=1 (βn − αj )
z→βn z→βn
M
X ord(αm )
=
j=1
z − αm
128 4. Taylor and Laurent Series
Similarly, we have:
N
h0 (z) X ord(βk )
= .
h(z) z − βk
k=1
Combining both results, we conclude that:
M N
f 0 (z) X ord(αj ) X ord(βk ) F 0 (z)
= − +
f (z) j=1
z − αj z − βk F (z)
k=1
Finally, the desired result (4.7) follows from the Cauchy’s Integral Formula and the fact
0
that FF is holomorphic on the region enclosed by γ.
Remark 4.13. The Argument Principle formula (4.7) can be generalized to the famous
Riemann-Roch’s Theorem for Riemann surfaces.
Example 4.15. Show that all zeros of z 5 + 3z + 1 are inside the circle |z| = 2.
Solution
The dominant term of z 5 + 3z + 1 is z 5 , so we let f (z) = z 5 and h(z) = 3z + 1. On
the circle |z| = 2, we have
5
|f (z)| = |z| = 25
whereas
|h(z)| = 3 |z| + 1 = 7 < 25 = |f (z)| .
By Rouche’s Theorem, we have
I I
f 0 (z) f 0 (z) + h0 (z)
dz = dz
|z|=2 f (z) |z|=2 f (z) + h(z)
Exercise 4.34. Show that z + ez + 3 has exactly one root with negative real part.
Exercise 4.35. Show that all roots of the equation z 6 + (1 + i)z + 1 = 0 lie in the
annulus 12 ≤ |z| < 54 .
The technique above gives another (yes, another one) proof of the Fundamental
Theorem of Algebra.
X∞ X∞
1 1
ζ(z) := z
= z ln n
.
n=1
n n=1
e
131
132 5. What is the Riemann Hypothesis?
Furthermore, the same functional equation holds for fˆ. Let’s verify this. For any z
such that Re(z) > 0, we have:
1 1
fˆ(z + 1) − 2fˆ(z) = f (z + 2) − 2 · f (z + 1)
2 2
1
= (f (z + 2) − 2f (z + 1)) = 0.
2
Now that fˆ is holomorphic on {z : Re(z) > 0} and satisfies the functional equation
fˆ(z + 1) = 2fˆ(z).
One can then repeat the same procedure as before to extend fˆ to a holomorphic function
ˆ
fˆ defined on {z : Re(z) > −1}, which is given by:
ˆ 1
fˆ(z) = fˆ(z + 1), for any z ∈ {Re(z) > −1}.
2
Inductively, we can repeat the same procedure over and over again, and extend f to
a function F : C → C that is holomorphic on the whole complex plane C.
y y
z z+1
1 f (z
fˆ(z) 2
+ 1)
x x
1
fˆ(z) = f (z)
f is originally defined here fˆ(z) = 21 f (z + 1)
x
−1 0 1
z z+1
Γ1 (z) 1 Γ(z + 1)
z
x x
−1
Γ1 (z) = Γ(z)
Γ is originally defined here Γ1 (z) = z1 Γ(z + 1)
136 5. What is the Riemann Hypothesis?
The functional equation for Γ then induces a new functional equation for Γ1 . When-
ever Re(z) > −1, we have:
1
Γ1 (z + 1) = Γ(z + 2) (Definition of Γ1 )
z+1
1
= · (z + 1)Γ(z + 1) (Functional equation for Γ)
z+1
= Γ(z + 1) = zΓ1 (z) (Definition of Γ1 ).
Therefore, one can define:
1
Γ2 (z) := Γ1 (z + 1)
z
for any z ∈ C such that z + 1 is in the domain of Γ1 , i.e. Re(z) > −2 and z 6= −1. As
such, Γ2 is an analytic continuation of Γ1 (and hence of Γ) on {Re(z) > −2}\{0, −1}.
Repeat the above process indefinitely, one can define analytic continuations Γm on
{Re(z) > −m}\{0, −1, −2, · · · , −(m − 1)}, and eventually an analytic continuation Γ̂ of
Γ on the domain C\{0, −1, −2, −3, · · · }.
y
x
−2 −1
Exercise 5.7. Show that for any integer m ≥ 1 and z in the domain of Γm , we have:
Γ(z + m)
Γm (z) =
z(z + 1) · · · (z + m − 1)
Exercise 5.8. Show that each non-positive integer −n is a simple pole of Γ̂, and
that:
(−1)n
Res(Γ̂, −n) = .
n!
Here is a summary of facts about the Gamma function:
• Γ is a priori defined on {z : Re(z) > 0}.
• By the relation Γ(z + 1) = zΓ(z), one can define an analytic continuation Γ̂ of Γ on
C\{0, −1, −2, −3, · · · }.
(−1)n
• Each non-positive integer −n is a simple pole of Γ̂, with residue equal to n! .
Recall that since the analytic continuation must be unique, some textbooks denote the
analytic continuation by simply Γ.
5.2. Riemann ζ Functions 137
Exercise 5.9. Show that when Re(z) > 0, the Gamma function can be decomposed
into: Z ∞
X∞
(−1)n
Γ(z) = + e−t tz−1 dt.
n=0
n!(z + n) 1
X∞
(−1)n
Show also that the infinite sum converges for any z 6= 0, −1, −2, −3, · · ·
n=0
n!(z + n)
Z ∞
and the integral e−t tz−1 dt is an entire function of z.
1
5.2.2. Relation between Γ and ζ. Recall that the Riemann zeta function ζ : {z :
Re(z) > 1} → C is defined by the infinite series:
X∞
1
ζ(z) =
n=1
nz
which converges when Re(z) > 1. The following lemma shows a relation between Γ and
ζ.
Proof. The key step of the proof is the change of variables t = nτ in the integral that
defines Γ:
Z ∞ Z ∞
Γ(z) = tz−1 e−t dt = (nτ )z−1 e−nτ d(nτ )
0 0
Z ∞
= nz τ z−1 e−nτ dτ
0
Z ∞
1
Γ(z) = tz−1 e−nt dt.
nz 0
Here we have used the fact that τ is a dummy variable. Summing up over n, we get:
X∞ X∞ Z ∞
1
(5.2) z
Γ(z) = tz−1 e−nt dt.
n=1
n n=1 0
Next we want to switch the integral and summation signs. It has to be justified using
LDCT. Consider:
tz−1 e−nt ≤ tx−1 e−nt
for any t ∈ [0, ∞). Note that:
∞ Z
X ∞ X∞
1
tx−1 e−nt dt = Γ(x)
n=1 0 n=1
nx
which converges since x > 1. Hence, LDCT shows we can switch the summation and
integral signs of (5.2), and it yields:
X∞ Z ∞ X∞
1
z
Γ(z) = tz−1
e−nt dt.
n=1
n 0 n=1
138 5. What is the Riemann Hypothesis?
∞
X ∞
X
−nt
n
Observing that e = e−t is a geometric series, we get:
n=1 n=1
∞
X e−t 1
e−nt = = t .
n=1
1 − e−t e −1
X ∞
f (n) (0) n
Recall from Exercise 4.8 that the Taylor’s series w converges uniformly
n=0
n!
on every ball B2π−ε (0) slightly smaller than B2π (0), say B2 (0). In particular, since
X∞
f (n) (0) n
[0, 1] ⊂ B2 (0), the convergence of the series t is also uniform on [0, 1]. When
n=0
n!
z is a fixed complex number such that Re(z) > 1, we have tz−1 ≤ tx−1 ≤ 1. Therefore,
X∞
f (n) (0) n+z−1
the series t also converges uniformly on t ∈ [0, 1] regarding z as fixed.
n=0
n!
Using the fact, one can write the first integral as:
Z 1 z−1 Z 1 ∞
!
t 1 X f (n)
(0)
t
dt = tz−1 + tn dt
0 e −1 0 t n=0 n!
Z 1 ∞
!
X f (n) (0) z+n−1
z−2
= t + t dt
0 n=0
n!
z−1 t=1 X ∞ t=1
t f (n) (0) tz+n
= +
z − 1 t=0 n=0 n! z + n t=0
∞
X f (n) (0)
1 1
= + · .
z − 1 n=0 n! z+n
Since Γ has an analytic continuation Γ̂ on C\{0, −1, −2, −3, · · · }. From the above
relation, we can then define an analytic continuation of ζ on C\{1, 0, −1, −2, −3, · · · } as:
∞ Z ∞ z−1 !
1 1 X f (n) (0) 1 t
(5.3) ζ̂(z) = + · + dt
Γ̂(z) z − 1 n=0 n! z+n 1 et − 1
X∞ Z ∞ z−1
1 f (n) (0) 1 1 t
= + + dt.
(z − 1)Γ̂(z) n=0 n! (z + n)Γ̂(z) Γ̂(z) 1 et − 1
It appears (5.3) has singularities at every 1, 0, −1, −2, −3, · · · , yet we can show
0, −1, −2, −3, · · · are all removable. It is because Γ̂ has a simple pole at every of
1/Γ̂
{0, −1, −2, −3, · · · }, so they are zeros of 1/Γ̂. Therefore, z+n has a removable singularity
at −n. Precisely, for any integers m, n ∈ {0, 1, 2, 3, · · · } we have:
(
1
1 if m = n
lim = Res(Γ̂,−n)
z→−m (z + n)Γ̂(z) 0 if m 6= n
This shows {0, 1, 2, 3, · · · } are all removable singularities of ζ̂ since the following limit is
finite for any m = 0, 1, 2, 3, · · ·
f (m) (0) 1
lim ζ̂(z) =
z→−m m! Res(Γ̂, −m).
Therefore, ζ(z) can be holomorphically defined on C\{1} by declaring that
f (m) (0) 1
ζ̂(−m) :=
m! Res(Γ̂, −m)
5.2.4. Special Values of ζ̂. We will determine the value of ζ̂ at some special z ∈ C.
When z = −m where −m is a non-positive integer, then we have already discussed that
f (m) (0) 1
ζ̂(−m) = .
m! Res(Γ̂, −m)
Here f is the function:
1 1
f (w) = − .
ew − 1 w
(−1)m
We have already figured out that Res(Γ̂, −m) = , so ζ̂(−m) = (−1)m f (m) (0).
m!
However, it is not straight-forward to find a general expression for f (m) (0), but by direct
computations one can verify that the first few terms of f (m) (0) are given as follows:
1 1 1
f (0) = − f 0 (0) = f 00 (0) = 0 f (3) (0) = − .
2 12 120
Therefore, the extended Riemann zeta function ζ̂ takes the following values:
1 1
ζ̂(0) = − ζ̂(−1) = −
2 12
1
ζ̂(−2) = 0 ζ̂(−3) =
120
1
To many people’s surprise, the fact that ζ̂(−1) = − 12 is used in String Theory!
However, many “muggles” misunderstand the meaning of it. They misinterpret it as
5.2. Riemann ζ Functions 141
P∞ 1 1
P∞
n=1 n−1 = − 12 , which is mathematically wrong as ζ̂(z) = n=1 n1z only when Re(z) >
1. It would lead to the following awkward and non-sense expression:
1
1 + 2 + 3 + 4 + ··· = − .
12
P∞
Similarly, those “muggles” also mix up ζ̂(0) = − 21 with n=1 n10 = − 12 , and ζ̂(−2) = 0
P∞ 1
with n=1 n−2 = 0, both would lead to awkward expressions:
1
1 + 1 + 1 + 1 + ··· = − .
2
12 + 22 + 32 + 42 + · · · = 0.
Exercise 5.10. Prove that for all positive integers n, one has
Bn+1
ζ̂(−n) = (−1)n
n+1
where Bn ’s are Bernoulli’s numbers defined as in Exercise 4.24.
1 1
where f (w) = − . It is not difficult to show that f (m) (0) = 0 for any even
ew − 1 w
integer m:
x
1
2
The zeros of ζ̂ have deep connections with the distribution of prime numbers. The
renowned Prime Number Theorem asserts that:
π(x)
lim =1
x→∞ x/ ln x
where π(x) is the number of positive prime numbers less than or equal to x. A corollary
of the theorem is that the n-th prime number pn is approximately equal to n ln n. The
proof of Prime Number Theorem relies surprisingly on the fact that there is no zero of
ζ̂ with real part equal to 1. If the Riemann Hypothesis is proven to be true, then the
Prime Number Theorem can be substantially improved, and many mysteries about the
distribution of primes will be revealed.
As of today (February 29, 2020), the mathematics community widely believes that the
conjecture is unsolved. It is one of the most important open problem in Pure Mathematics
nowadays. In 2000, the Clay Mathematics Institute compiled a list of 7 problems, called
Millennium Prize Problems. For each problem in the list, the institute promises to award
US$1,000,000 to the first person who solves or disproves it. Riemann Hypothesis is
one of the problems in the list. The other 6 problems are: P versus NP Problem, Hodge
Conjecture, Poincaré Conjecture, Yang-Mills Existence and Mass Gap, Navier-Stokes
Existence and Smoothness, and Birch Swinnerton-Dyer Conjecture. The only Millennium
Prize Problem that was solved is the Poincaré Conjecture, by Grigori Perelman in 2002-03
using the idea of Ricci flow developed by Richard Hamilton in 1982.
* End of MATH 4023 *
** I hope you have learned a lot and/or enjoyed the course. **
Appendix A
In this appendix we list some important concepts and theorems from MATH 2043/3033
that we will use frequently in this course. Proofs are all omitted since they are essentially
the same as in the real case. This appendix is intended to be brief (no worked example
here). For detail, please consult Chapter 10 of MATH 3033, or Chapter 4 in MATH 2043.
143
144 A. Results from MATH 2043/3033
Proposition A.3. Suppose fn (z) converges uniformly on Ω to the limit function f (z),
then:
• If fn are continuous on Ω for all n, then f is also continuous on Ω.
• For any α ∈ Ω, we have
lim lim fn (z) = lim lim fn (z).
z→α n→∞ n→∞ z→α
• Let [a, b] be a bounded interval in R, and fn (t)’s be integrable functions on [a, b] If
fn (t) converges uniformly to f (t) on [a, b], then
Z b Z b
lim fn (t) dt = lim fn (t) dt.
n→∞ a a n→∞
• Let γ be a curve in C of finite length, and fn (z)’s be integrable functions on γ. If fn (z)
converges uniformly to f (z) on Ω, then
Z Z
lim fn (z) dz = lim fn (z) dz.
n→∞ γ γ n→∞
∞
X
Analogous results hold for uniform convergence series. For instance, if fn (z) converges
n=1
uniformly on Ω, then for any curve γ in Ω of finite length, we have:
Z X∞ X∞ Z
fn (z) dz = fn (z) dz.
γ n=1 n=1 γ
In the above proposition, the conditions that [a, b] is a finite interval, and γ is
a curve of finite length are necessary. For unbounded intervals or curves, uniform
convergence is not sufficient to guarantee the switching of the integral and summation
signs! Fortunately, there is another tool to deal with improper integrals, namely Lebesgue
Dominated Convergence Theorem (LDCT), which stems from measure theory:
then we have: Z ∞
bX ∞ Z
X b
gn (t) dt = gn (t) dt.
a n=1 n=1 a
A. Results from MATH 2043/3033 145
Recall from MATH 3033/2043 that even if fn (x) converges uniformly on (a, b) to
f (x), the derivatives fn0 (x) may not converge to f 0 . Sometimes,
P∞ the limit of fn0 may not
even be differentiable. Likewise, even when the sum n=1 fn (x) converges uniformly
on (a, b), term-by-term differentiation
∞ ∞
d X X
fn (x) = fn0 (x)
dx n=1 n=1
may not hold. Instead, what is needed is pointwise convergence of fn (z), and uniform
convergence of fn0 (z):
Then, we have
∞ ∞
d X X d
fn (z) = fn (z).
dz n=1 n=1
dz
When using Morera’s Theorem, we often consider a double integral of the form:
I Z b
f (z, t) dtdz.
T a
It we can switch the two integral signs, and it happens that f (z, t) is a holomorphic
function for each fixed t ∈ [a, b], then we have:
I Z b Z bI Z b
f (z, t) dtdz = f (z, t) dzdt = 0 dt = 0.
T a a T a
The question is whether we can switch the two integral signs. It thanks for the following
(special case) of Fubini’s Theorem