Chapter 26Cauchy's Integral Formula
In the previous chapter, we proved Cauchy's integral theorem—the statement that if a holomorphic function (a differentiable complex function) is integrated once around a simple closed curve, the result is always $0$. As a corollary, we obtained the “contour integral formula” $\displaystyle\oint_C\frac{1}{z-\alpha}\,dz=2\pi i$ (when $\alpha$ lies inside $C$). In this chapter we combine these two tools and arrive at one of the most striking theorems of complex analysis—Cauchy's integral formula.
In a single sentence, the theorem claims the following: the value of a holomorphic function at a point in the interior of a domain is completely determined by the values of the function on the boundary of that domain (on its perimeter) alone. In the world of real functions, nothing so convenient happens. For example, no matter how well we know the values $f(0),f(1)$ at the two endpoints of the interval $[0,1]$, we can learn nothing about the value $f(0.5)$ at an interior point (if $f$ is merely differentiable). For complex functions, however, from the single condition of being holomorphic—that is, of satisfying the Cauchy–Riemann equations—all the values in the interior can be recovered from boundary information alone. This property, “the boundary determines everything,” is the starting point of nearly all the deep results of complex analysis, such as Taylor expansions, the residue theorem, and the theory of harmonic functions, which we will study later. In Yasuhide Mochizuki's “Mathematics Notes: Complex Analysis,” the source of this chapter, this theorem appears on pp. 34, 39 (the heading, and the figure plus the text of the Theorem), but the proof breaks off with only the English word “Proof” and the ellipsis mark “…”. In this chapter, bringing to bear all the tools assembled in the preceding chapters (Green's theorem, the principle of path deformation, and the contour integral formula), we complete this proof to the end with our own hands.
Before we reach the proof, this chapter makes two detours. The first is to get a feel for Green's theorem, whose proof alone was left as homework in the previous chapter, by verifying it for a concrete function (Section 26.1). The second is to explore, with a keyhole contour, what happens when we integrate $\log z$—a function for which “the answer is not uniquely determined for a single $z$” (a multivalued function)—once around a domain with a hole punched out so as to avoid the origin (a multiply connected domain) (Section 26.2). In fact, this computation is a treacherous investigation: unless we pay the utmost attention to how the branch is chosen, the result appears at first glance to contradict Cauchy's integral theorem. In this chapter, taking care to keep the branch consistent from beginning to end, we compute this integral correctly all the way and confirm that it does not contradict Cauchy's integral theorem. This experience—that “handling a multivalued function requires care in how the path is chosen”—also serves as foreshadowing for the proof of Cauchy's integral formula in Section 26.3.
- Verifying Green's theorem $\displaystyle\int_C f\,dx=-\iint_D\pdiff{f}{y}\,dx\,dy$ for the concrete function $f(x,y)=x^2+y^2$ and the upper half-disk, by computing both the line integral over the boundary and the double integral over the interior
- The idea of a branch and a branch point of the multivalued function $\log z$, and the computation of the integral once around a multiply connected domain with the origin removed, using a keyhole contour
- Integrating $\log z$ in a multiply connected domain, computing correctly while paying attention to the choice of branch and to the sign of the antiderivative, and confirming that the integrals add up to exactly $0$
- The statement of Cauchy's integral formula $\displaystyle f(z)=\frac{1}{2\pi i}\oint_{\partial D}\frac{f(\zeta)}{\zeta-z}\,d\zeta$, and its complete proof in three steps: “shrink to a small circle,” “split $f$ into a constant part and an error part,” and “bound the error with the $ML$ inequality”
- The technique of using Cauchy's integral formula to evaluate contour integrals of the form $\displaystyle\oint\frac{f(z)}{z-\alpha}\,dz$ in an instant, without clearing the denominator of the integrand
- As corollaries, the fact that the value of a holomorphic function can be written as the mean value over a circle surrounding the point (the mean value property), and a gateway to the astonishing fact that a holomorphic function can in fact be differentiated any number of times
Source notes: Yasuhide Mochizuki, Mathematics Notes: Complex Analysis, pp. 34–40.
26.1 A Concrete Example of Green's Theorem — Verification on the Upper Half-Disk
Let us first recall Green's theorem, which we stated without derivation in Section 25.3. For a region $D$ in the plane and its boundary $C=\partial D$ (oriented counterclockwise, that is, in the direction in which $D$ is always on our left as we move along the boundary), if $f(x,y)$ has continuous partial derivatives on a slightly larger set containing $D$ and its boundary, then
$$ \int_C f(x,y)\,dx=-\iint_D \pdiff{f}{y}\,dx\,dy $$holds; this is (part of) Green's theorem. The claim that a “one-dimensional integral over the boundary” equals a “two-dimensional integral over the interior” seems mysterious at first. In this section we choose one specific $f$ and $D$, compute both sides head-on, and see with our own eyes that they really agree. This is exactly the computation on p. 035 of the source notes.
26.1.1 Setting up the region and its boundary
Let $r\gt0$ be a constant, and set
$$ D=\{(x,y)\mid x^2+y^2\le r^2,\ y\ge0\},\qquad C=\{(x,y)\mid x^2+y^2=r^2,\ y\ge0\}\cup\{(x,y)\mid -r\le x\le r,\ y=0\} $$That is, $D$ is the upper half of the disk of radius $r$ (the upper half-disk), and $C$ is its boundary—the arc part and the diameter part on the $x$-axis—joined together into one curve. The function we use for the check is
$$ f(x,y)=x^2+y^2 $$This is the square of the distance from the origin, and it is simple and easy to compute with on $D$. Incidentally, the source notes draw this region as a perspective view with three-dimensional depth, in which the $xy$-plane is seen obliquely from above. However, the problem itself concerns the plane (two dimensions), so in this book we redraw it, as in Figure 26.1, as a two-dimensional picture seen from directly above.
Intuition: Getting the direction of the boundary wrong changes the sign
The left-hand side $\int_C f\,dx$ of Green's theorem changes sign depending on which way we traverse the boundary $C$ (for the same reason as property (iii) in Section 25.1: $\int_{-C}f\,dz=-\int_C f\,dz$). “Counterclockwise” means the direction in which the interior of $D$ is on our left as we face the direction of travel. As in Figure 26.1, the positive direction for this upper half-disk is leftward along the arc (the direction in which $\theta$ increases from $0$ to $\pi$) and from left to right along the diameter. In Section 26.1.3 we check the same function on the lower half-disk and see how the sign changes when the positive direction of the boundary is reversed relative to that for the upper half-disk.
26.1.2 Computing both sides
We split the boundary $C$ into the arc part ($C_{\text{arc}}$) and the diameter part ($C_{\text{diam}}$), and compute $\int_C f\,dx=\int_{C_{\text{arc}}}f\,dx+\int_{C_{\text{diam}}}f\,dx$. First we parametrize the arc part by polar coordinates:
$$ x=r\cos\theta,\qquad y=r\sin\theta\qquad(0\le\theta\le\pi) $$As $\theta$ increases from $0$ to $\pi$, the point $(x,y)$ traces exactly the arc in Figure 26.1 counterclockwise. Then
$$ dx=\diff{x}{\theta}\,d\theta=-r\sin\theta\,d\theta $$(Differentiating $x=r\cos\theta$ with respect to $\theta$ gives $-r\sin\theta$; since $r$ is a constant, it can be taken outside.) Also, on the arc we have $x^2+y^2=r^2$ (because the points lie on the circle), so $f(x,y)=r^2$ is constant on the arc. Hence
$$ \int_{C_{\text{arc}}}f\,dx=\int_0^\pi r^2\cdot(-r\sin\theta)\,d\theta=-r^3\int_0^\pi\sin\theta\,d\theta=-r^3\Big[-\cos\theta\Big]_0^\pi=-r^3\big(-(-1)-(-1)\big)=-r^3\cdot2=-2r^3 $$(We used $\int_0^\pi\sin\theta\,d\theta=[-\cos\theta]_0^\pi=-\cos\pi-(-\cos0)=-(-1)+1=2$.)
Next, the diameter part. Here we stay on the line $y=0$, and $x$ moves from $-r$ to $r$. Since $f(x,0)=x^2+0^2=x^2$,
$$ \int_{C_{\text{diam}}}f\,dx=\int_{-r}^{r}x^2\,dx=\left[\frac13x^3\right]_{-r}^{r}=\frac13r^3-\left(-\frac13r^3\right)=\frac23r^3 $$Adding these together, the line integral over the boundary (the left-hand side) is
$$ \int_C f\,dx=-2r^3+\frac23r^3=-\frac63r^3+\frac23r^3=-\frac43r^3 $$Next we compute the right-hand side, that is, the double integral over the interior. We have $\pdiff{f}{y}=\pdiff{}{y}(x^2+y^2)=2y$ (we regarded $x$ as a constant and differentiated with respect to $y$ alone). Expressing the upper half-disk in polar coordinates $x=r'\cos\theta,\,y=r'\sin\theta$ (with radial coordinate $r'$, $0\le r'\le r$, $0\le\theta\le\pi$), the area element is $dx\,dy=r'\,dr'\,d\theta$ (the Jacobian, from Chapter 7, Section 7.4), so
$$ \iint_D\pdiff{f}{y}\,dx\,dy=\int_0^\pi\int_0^r 2y\cdot r'\,dr'\,d\theta=\int_0^\pi\int_0^r 2(r'\sin\theta)\cdot r'\,dr'\,d\theta=\int_0^\pi\int_0^r 2r'^2\sin\theta\,dr'\,d\theta $$(We substituted $y=r'\sin\theta$.) The integrals in $\theta$ and $r'$ are independent of each other and can be separated (one factor for each variable), so
$$ =\left(\int_0^r 2r'^2\,dr'\right)\left(\int_0^\pi\sin\theta\,d\theta\right)=\left[\frac23r'^3\right]_0^r\cdot\underbrace{\Big[-\cos\theta\Big]_0^\pi}_{=2}=\frac23r^3\cdot2=\frac43r^3 $$Therefore the right-hand side is $-\dfrac43r^3$, which agrees exactly with the value $-\dfrac43r^3$ of the left-hand side. Thus, for a specific function and region, we have seen with our own eyes that Green's theorem
$$ \int_C f\,dx=-\iint_D\pdiff{f}{y}\,dx\,dy $$really does hold.
Illustration 26.1: Verification of Green's theorem ($f=x^2+y^2$, upper half-disk)
Let $f(x,y)=x^2+y^2$ and let $D$ be the upper half-disk of radius $r$. Verify that $\displaystyle\int_C f\,dx=-\iint_D\pdiff{f}{y}\,dx\,dy$ holds.
Solution. As computed in the text, the left-hand side $=-\dfrac43r^3$ and the right-hand side $=-\dfrac43r^3$, so they agree (verified with Python (SymPy): verify/univmath/ch26_green_theorem.py).
26.1.3 Example: What happens on the lower half-disk?
Example 26.1: Green's theorem on the lower half-disk
For the same function $f(x,y)=x^2+y^2$, consider $D'=\{(x,y)\mid x^2+y^2\le r^2,\ y\le0\}$ (the lower half-disk) and its boundary $C'$ (traversed in the positive direction). Compute $\displaystyle\int_{C'}f\,dx$ and $\displaystyle-\iint_{D'}\pdiff{f}{y}\,dx\,dy$ separately, and state how the sign changes compared with the result for the upper half-disk in Section 26.1.2 ($-\dfrac43r^3$).
Solution. First we confirm the positive direction of the boundary $C'$ (the direction in which $D'$ is on our left). The lower half-disk lies on the side $y\le0$, so to traverse the boundary counterclockwise, we go along the arc part from the point $(-r,0)$ through the point $(0,-r)$ to the point $(r,0)$ (the direction in which the angle $\theta$ increases from $-\pi$ to $0$), and then come back along the diameter part from the point $(r,0)$ to the point $(-r,0)$ (the direction in which $x$ decreases from $r$ to $-r$), which closes the loop. Note that both the arc and the diameter are traversed in the opposite direction from the upper half-disk case.
The arc part: $x=r\cos\theta,\,y=r\sin\theta$ ($\theta:-\pi\to0$), and since $f=x^2+y^2=r^2$ on the arc,
$$ \int_{C'_{\text{arc}}}f\,dx=\int_{-\pi}^{0}r^2\cdot(-r\sin\theta)\,d\theta=-r^3\Big[-\cos\theta\Big]_{-\pi}^{0}=-r^3\big(-1-1\big)=2r^3 $$(The value is not $-\cos0-(-\cos(-\pi))=-1-1=-2$ but rather $\big[-\cos\theta\big]_{-\pi}^0=(-\cos0)-(-\cos(-\pi))=(-1)-(1)=-2$. Hence $-r^3\times(-2)=2r^3$.) The diameter part: since $x$ moves from $r$ to $-r$,
$$ \int_{C'_{\text{diam}}}f\,dx=\int_{r}^{-r}x^2\,dx=\left[\frac13x^3\right]_r^{-r}=-\frac13r^3-\frac13r^3=-\frac23r^3 $$Summing up, the left-hand side is
$$ \int_{C'}f\,dx=2r^3-\frac23r^3=\frac43r^3 $$On the other hand, for the double integral on the right-hand side, the only change is that the range of $\theta$ becomes $-\pi\le\theta\le0$, and we get
$$ \iint_{D'}\pdiff{f}{y}\,dx\,dy=\int_{-\pi}^0\int_0^r 2r'^2\sin\theta\,dr'\,d\theta=\frac23r^3\cdot\Big[-\cos\theta\Big]_{-\pi}^0=\frac23r^3\cdot(-2)=-\frac43r^3 $$Hence $-\iint_{D'}\pdiff{f}{y}\,dx\,dy=\dfrac43r^3$, which agrees with the left-hand side $\dfrac43r^3$. Whereas the result for the upper half-disk was $-\dfrac43r^3$, for the lower half-disk the sign has reversed and the value is $+\dfrac43r^3$. This is because, while $f=x^2+y^2$ is symmetric with respect to $x,y$ (its value does not change under $y\to-y$), the positive direction of the boundary runs in opposite senses for the upper and lower half-disks (that is, on the arc, the direction in which $\theta$ moves is exactly reversed between the upper half-disk and the lower half-disk); this illustrates well that the left-hand side of Green's theorem is a quantity that depends on orientation (verified with Python: verify/univmath/ch26_green_theorem.py).
Why is this section placed before Sections 26.2 and 26.3?
Green's theorem is the tool that served as the foundation when we proved Cauchy's integral theorem (if $f$ is holomorphic, then $\oint_C f(z)\,dz=0$) in Section 25.3 of the previous chapter. By checking this “foundation” for ourselves in this section, the idea that “the integral over the boundary is tied to information about the interior” should already be familiar when we prove Cauchy's integral formula in Section 26.3. Moreover, the path integral of $\log z$ treated in Section 26.2 also uses exactly the same technique as this section—computing an integral along a boundary made by joining several arcs and line segments (divide the boundary into pieces, parametrize each piece, and integrate).
26.2 Path Integrals of $\log z$ in a Multiply Connected Domain and Branch Points
In Section 25.2, we saw that computing $\displaystyle\oint_{|z|=r}\log z\,dz$ along the unit circle (the circle centered at the origin) does not give $0$. The reason was that $\log z$ is not holomorphic at $z=0$ (or rather, cannot even be defined at $z=0$). So what happens if we avoid the origin itself and integrate $\log z$ once around a closed path in a region from which only the neighborhood around the origin has been removed? In this section we follow this investigation on p. 036 of the source notes—the computation using a keyhole-shaped contour—through to the end. In fact, this computation hides a small “trap” concerning the choice of branch, and the source notes end with the remark “It doesn't come out to 0…!?” still written at the end. The goal of this section is to avoid this “trap,” keep the branch consistent from beginning to end, and arrive at the correct answer (the result is $0$).
26.2.1 The multivalued function $\log z$ and branches
Write a complex number $z\ne0$ in polar form $z=|z|(\cos\theta+i\sin\theta)=|z|e^{i\theta}$ (Chapter 21, Section 21.2). Then $\theta=\arg z$ (the argument of $z$) is ambiguous up to integer multiples of $2\pi$: both $\theta$ and $\theta+2\pi$ represent the same point $z$. The complex logarithm $\log z$ is defined as follows, so that the property $e^{\log z}=z$ of the real logarithm continues to hold, unchanged, for complex numbers as well:
$$ \log z=\log|z|+i\theta\qquad(\theta\text{ is any one of }z\text{'s arguments}) $$This is the definition given in Chapter 23. Here $\log|z|$ is the ordinary natural logarithm of a real number. The problem is that there are infinitely many candidates $\theta_0,\theta_0+2\pi,\theta_0-2\pi,\dots$ for $\theta$, so $\log z$ has no single determined value for a single $z$. Such a function is called a multivalued function.
Intuition: Following the argument continuously causes a shift
Suppose a point $z=e^{i\theta}$ starts from $\theta=0$ (that is, $z=1$) on the unit circle centered at the origin, goes once around counterclockwise, and returns to $z=1$. The argument at the start was $\theta=0$, but when the point has gone once around and come back to the same point $z=1$, the result of “continuing to increase the argument continuously” is $\theta=2\pi$. Even though it is the same point $z=1$, the value of the argument (and hence the value of $\log z=\log1+i\theta=i\theta$) is shifted by $2\pi i$ between the start and the end. A point such that going once around it (here, $z=0$) does not bring the value of the function back to its original value is called a branch point, and $z=0$ is a branch point of $\log z$. What we obtain by choosing a region that does not enclose a branch point and, in addition, fixing the range of $\theta$ to some interval (for example $-\pi\lt\theta\le\pi$ or $0\le\theta\lt2\pi$) so that a single value is determined is called one branch of $\log z$. We may fix the branch wherever we like, but within a single computation we must use the same branch consistently from beginning to end—this is the most important caution in this section.
26.2.2 Setting up the keyhole contour
Fix positive real numbers $\varepsilon,R$ ($0\lt\varepsilon\lt R$) and consider the closed curve $C=C_1+C_2+C_3+C_4$ consisting of the following four curves (exactly the setup of the source notes).
$$ \begin{aligned} C_1&=\{(x,y)\mid x^2+y^2=R^2,\ y\ge0\} &&\text{(upper semicircular arc of radius }R\text{, counterclockwise)}\\ C_2&=\{(x,y)\mid -R\le x\le-\varepsilon,\ y=0\} &&\text{(on the negative real axis, oriented from the side where }x\text{ decreases toward the origin)}\\ C_3&=\{(x,y)\mid x^2+y^2=\varepsilon^2,\ y\ge0\} &&\text{(upper semicircular arc of radius }\varepsilon\text{, clockwise)}\\ C_4&=\{(x,y)\mid \varepsilon\le x\le R,\ y=0\} &&\text{(on the positive real axis, oriented away from the origin)} \end{aligned} $$Connecting these four curves in the order $C_1\to C_2\to C_3\to C_4$ gives a closed curve shaped just like a keyhole (Figure 26.2). The domain $D_{\varepsilon,R}=\{z\mid \varepsilon\le|z|\le R,\ \operatorname{Im}z\ge0\}\setminus(\text{neighborhood of the origin})$ enclosed by this closed curve does not contain the origin $z=0$ (the branch point of $\log z$). Therefore, if we use the branch that fixes the range of $\theta$ in this domain to $0\le\theta\le\pi$—that is, the branch with $\theta=0$ on the positive real axis ($C_4$) and $\theta=\pi$ on the negative real axis ($C_2$)—then $\log z=\log|z|+i\theta$ is continuous and holomorphic on the whole domain. Particular care is needed with the point that “on $C_2$ we use $\theta=\pi$” (we confirm this in detail in Section 26.2.4).
26.2.3 The integral along each of the four paths
At last, we compute $\displaystyle\int_C\log z\,dz=\int_{C_1}+\int_{C_2}+\int_{C_3}+\int_{C_4}$ carefully, one path at a time. On every path we use the branch fixed earlier (it starts from $\theta=0$ on $C_4$, passes through the upper half-plane, and reaches $\theta=\pi$ on $C_2$).
Computation: $C_1$ (the large arc of radius $R$)
Let $z=Re^{i\theta}$ ($\theta:0\to\pi$). Then $\log z=\log R+i\theta$ and $dz=iRe^{i\theta}\,d\theta$ (differentiating $z=Re^{i\theta}$ with respect to $\theta$ gives $iRe^{i\theta}$). Hence
$$ \int_{C_1}\log z\,dz=\int_0^\pi(\log R+i\theta)\,iRe^{i\theta}\,d\theta =iR\log R\int_0^\pi e^{i\theta}\,d\theta-R\int_0^\pi\theta e^{i\theta}\,d\theta $$(We expanded $(\log R+i\theta)\cdot iRe^{i\theta}=iR\log R\,e^{i\theta}+i^2R\theta e^{i\theta}=iR\log R\,e^{i\theta}-R\theta e^{i\theta}$ and split it into terms that can be integrated with respect to $\theta$.) We start with the first integral:
$$ \int_0^\pi e^{i\theta}\,d\theta=\left[\frac{e^{i\theta}}{i}\right]_0^\pi $$Here we first confirm that $\dfrac1i=-i$ (multiplying both sides of $1/i=i^{-1}$ by $i$, we get $i\cdot i^{-1}=1$ and $i\cdot(-i)=-i^2=1$, so $i^{-1}=-i$). Therefore $\dfrac{e^{i\theta}}{i}=-ie^{i\theta}$, and
$$ \int_0^\pi e^{i\theta}\,d\theta=\Big[-ie^{i\theta}\Big]_0^\pi=-ie^{i\pi}-(-ie^0)=-i(-1)-(-i)(1)=i+i=2i $$(We used $e^{i\pi}=-1,\ e^0=1$.) The second integral is computed by integration by parts ($u=\theta,\ dv=e^{i\theta}d\theta$, $du=d\theta,\ v=\dfrac{e^{i\theta}}{i}=-ie^{i\theta}$).
$$ \int_0^\pi\theta e^{i\theta}\,d\theta=\Big[\theta\cdot(-ie^{i\theta})\Big]_0^\pi-\int_0^\pi(-ie^{i\theta})\,d\theta =\big(\pi\cdot(-i)(-1)-0\big)+i\int_0^\pi e^{i\theta}\,d\theta $$(The $[uv]$ term is $\pi\cdot(-i)e^{i\pi}=\pi\cdot(-i)(-1)=i\pi$ at $\theta=\pi$ and $0$ at $\theta=0$. Into the remaining term $-\int(-ie^{i\theta})d\theta=+i\int e^{i\theta}d\theta$ we substitute $\int_0^\pi e^{i\theta}d\theta=2i$, which we found just above.)
$$ \int_0^\pi\theta e^{i\theta}\,d\theta=i\pi+i\cdot2i=i\pi+2i^2=i\pi-2 $$(We used $i\cdot2i=2i^2=-2$.) This result has been verified exactly with SymPy, using the verification script ch26_logz_keyhole.py that accompanies this book. Using this result, we obtain
(We used $iR\log R\cdot2i=2i^2R\log R=-2R\log R$.)
Caution: Always write out and check the sign when taking the reciprocal of $i$
In the computation above we used the fact $1/i=-i$ many times. Even though it looks like a simple formula, it is easy to get a sign wrong in the middle of an expansion, so in complex analysis computations the fastest way, in the end, is to write out $i^2=-1,\ 1/i=-i,\ i\cdot(-i)=1$ each time as you go along. In the computations that follow in this book we will not repeat the same note, for lack of space, but we would like the reader to reproduce the same computations by hand and to check the signs one by one.
Computation: $C_3$ (the small arc of radius $\varepsilon$, traversed in the opposite direction)
$C_3$ traverses the arc of radius $\varepsilon$ in the direction $\theta:\pi\to0$, opposite to $C_1$ (Figure 26.2). Since $z=\varepsilon e^{i\theta}$, $\log z=\log\varepsilon+i\theta$, and $dz=i\varepsilon e^{i\theta}d\theta$, it suffices to carry out the same computation as for $C_1$ with $R\to\varepsilon$ and with the interval of integration $\pi\to0$.
$$ \int_{C_3}\log z\,dz=\int_\pi^0(\log\varepsilon+i\theta)\,i\varepsilon e^{i\theta}\,d\theta =i\varepsilon\log\varepsilon\int_\pi^0 e^{i\theta}\,d\theta-\varepsilon\int_\pi^0\theta e^{i\theta}\,d\theta $$Since the interval of integration is reversed, $\displaystyle\int_\pi^0 e^{i\theta}d\theta=-\int_0^\pi e^{i\theta}d\theta=-2i$ and $\displaystyle\int_\pi^0\theta e^{i\theta}d\theta=-\int_0^\pi\theta e^{i\theta}d\theta=-(i\pi-2)=-i\pi+2$. Hence
$$ \int_{C_3}\log z\,dz=i\varepsilon\log\varepsilon\cdot(-2i)-\varepsilon(-i\pi+2)=2\varepsilon\log\varepsilon+i\pi\varepsilon-2\varepsilon $$(We used $i\cdot(-2i)=-2i^2=2$.)
Computation: $C_4$ (the positive real axis, the branch with $\theta=0$)
On $C_4$ we have $z=x$ (where $x$ is a positive real number moving from $\varepsilon$ to $R$) and the argument is $\theta=0$, so $\log z=\log x$ (which agrees with the real logarithm function). Since $dz=dx$, this is an ordinary real integral. By integration by parts ($u=\log x,\,dv=dx,\,du=dx/x,\,v=x$),
$$ \int_{C_4}\log z\,dz=\int_\varepsilon^R\log x\,dx=\Big[x\log x\Big]_\varepsilon^R-\int_\varepsilon^R x\cdot\frac1x\,dx=\Big[x\log x\Big]_\varepsilon^R-\Big[x\Big]_\varepsilon^R $$(Since $x\cdot\frac1x=1$, the remaining integral is $\int_\varepsilon^R1\,dx=[x]_\varepsilon^R$.) Hence
$$ \int_{C_4}\log z\,dz=(R\log R-R)-(\varepsilon\log\varepsilon-\varepsilon) $$This result also agrees with the computation on p. 036 of the source notes.
26.2.4 Computing $C_2$
In computing the last remaining piece, $C_2$ (the part on the negative real axis), we must take particular care in choosing the argument. What we must watch is the argument at a point $z=x$ on $C_2$ (where $x$ is a negative real number, $-R\le x\le-\varepsilon$). The branch fixed in Section 26.2.2 was this: “set $\theta=0$ on the positive real axis, and increase $\theta$ continuously as we move counterclockwise through the upper half-plane along $C_1$.” If we stay on this branch and arrive at $C_2$ (the negative real axis), then $\theta$ has increased to exactly $\pi$. That is, the expression
$$ \log z=\log|x|+i\pi\qquad(z=x\lt0,\ \text{branch reached from the upper half-plane}) $$must be used. If we carelessly reuse the same expression $\log z=\log|x|$ (the branch with $\theta=0$) as on $C_4$, we drop this contribution $+i\pi$, so be careful.
Caution: For a multivalued function, do not reuse an expression of “the same form” in a different place
Since $C_2$ and $C_4$ both lie on the $x$-axis, at first glance one is tempted to use the same expression $\log z=\log|x|$ for both. But as we saw in Section 26.2.1, $\log z$ is a multivalued function, and which branch is in use is determined along the path—the argument that we have followed along $C_1\to C_2$ has reached $\theta=\pi$, unlike $\theta=0$ on $C_4$. Rather than reusing an expression on the mere appearance that “it is on the $x$-axis,” we should make a habit of tracking, each time, “which branch we are on now” along the path.
With this in mind, let us compute the integral over $C_2$. We have $z=x$ (with $x:-R\to-\varepsilon$), $\log z=\log|x|+i\pi=\log(-x)+i\pi$ (since $x\lt0$, we have $|x|=-x$), and $dz=dx$.
$$ \int_{C_2}\log z\,dz=\int_{-R}^{-\varepsilon}\big(\log(-x)+i\pi\big)\,dx=\int_{-R}^{-\varepsilon}\log(-x)\,dx+i\pi\int_{-R}^{-\varepsilon}dx $$The second term is easy: $\displaystyle i\pi\int_{-R}^{-\varepsilon}dx=i\pi\big((-\varepsilon)-(-R)\big)=i\pi(R-\varepsilon)$. For the first term, $\displaystyle\int_{-R}^{-\varepsilon}\log(-x)\,dx$, we must be careful in constructing the antiderivative. Integrating by parts with $u=\log(-x),\ dv=dx$, we have $v=x$ (an antiderivative of $dv=dx$ is $x$, not $-x$), and since $\dfrac{d}{dx}\log(-x)=\dfrac{1}{-x}\cdot(-1)=\dfrac1x$ (by the chain rule: differentiating the outer function $\log(\cdot)$ gives $\frac{1}{-x}$, differentiating the inner function $-x$ gives $-1$, and multiplying these gives $\frac{1}{-x}\cdot(-1)=\frac1x$), we get $du=\dfrac1x\,dx$. Therefore, we obtain
$$ \int\log(-x)\,dx=x\log(-x)-\int x\cdot\frac1x\,dx=x\log(-x)-\int1\,dx=x\log(-x)-x $$(note that an antiderivative of $dv=dx$ is $v=x$, not $v=-x$—if we carelessly took $v=-x$, the signs in what follows would be off). Using this antiderivative $x\log(-x)-x$, we get
$$ \int_{-R}^{-\varepsilon}\log(-x)\,dx=\Big[x\log(-x)-x\Big]_{-R}^{-\varepsilon} =\big((-\varepsilon)\log\varepsilon-(-\varepsilon)\big)-\big((-R)\log R-(-R)\big) $$ $$ =(-\varepsilon\log\varepsilon+\varepsilon)-(-R\log R+R)=R\log R-R-\varepsilon\log\varepsilon+\varepsilon $$(We expanded term by term, keeping in mind that $-(-x)=x$. When $x=-\varepsilon$, we have $-x=\varepsilon$, so $\log(-x)=\log\varepsilon$; when $x=-R$, we have $-x=R$, so $\log(-x)=\log R$.)
From the above, the integral over $C_2$ is
$$ \int_{C_2}\log z\,dz=\big(R\log R-R-\varepsilon\log\varepsilon+\varepsilon\big)+i\pi(R-\varepsilon) $$This is the result.
26.2.5 The sum of the four integrals — it comes out to a clean $0$
We list all four integrals found in Sections 26.2.3 and 26.2.4.
$$ \begin{aligned} \int_{C_1}\log z\,dz&=-2R\log R-i\pi R+2R\\ \int_{C_2}\log z\,dz&=R\log R-R-\varepsilon\log\varepsilon+\varepsilon+i\pi R-i\pi\varepsilon\\ \int_{C_3}\log z\,dz&=2\varepsilon\log\varepsilon+i\pi\varepsilon-2\varepsilon\\ \int_{C_4}\log z\,dz&=R\log R-R-\varepsilon\log\varepsilon+\varepsilon \end{aligned} $$Collecting the $R\log R$ terms gives $-2R\log R+R\log R+R\log R=0$. Collecting the $\varepsilon\log\varepsilon$ terms gives $-\varepsilon\log\varepsilon+2\varepsilon\log\varepsilon-\varepsilon\log\varepsilon=0$. Collecting the $R$ terms (those containing no logarithm) gives $2R-R+R=2R$. “Hmm, this does not vanish…” you might think, but since $C_2$ has a $-R$ and $C_4$ has a $-R$, the total is $2R-R-R=0$ (compare the four expressions carefully and be sure to pick up every like term without missing any). Similarly, the $\varepsilon$ terms give $\varepsilon-2\varepsilon+\varepsilon=0$. Finally, the terms containing $i\pi$ give $-i\pi R+i\pi R-i\pi\varepsilon+i\pi\varepsilon=0$. All the terms cancel out neatly, giving
This is an equality that holds exactly for any $0\lt\varepsilon\lt R$ (even without taking the limits $\varepsilon\to0$ or $R\to\infty$). The function $\log z$ is holomorphic in the keyhole domain $D_{\varepsilon,R}$ (once the branch is fixed), and $D_{\varepsilon,R}$ can be regarded as a simply connected domain (Section 25.4: if we regard $C_2,C_4$ as two walls with opposite orientations along the “slit” of the keyhole, the whole becomes a single simply connected domain). Therefore, by Cauchy's integral theorem (Section 25.3), we should have $\oint_C\log z\,dz=0$, and Eq. \eqref{eq:26-logz-total} agrees perfectly with this expectation. We have now confirmed that, as long as the choice of branch and the sign of the antiderivative are handled correctly, there is no problem at all with $\log z$ itself or with Cauchy's integral theorem.
Math Note: The divergent terms cancel out exactly
As $R\to\infty$ we have $R\log R\to\infty$; as $\varepsilon\to0$ we have $\varepsilon\log\varepsilon\to0$ (this can be shown, for example, by l'Hôpital's rule); and as $R\to\infty$ we have $i\pi R\to\infty$. Thus $C_1$ and $C_2$, viewed individually, contain divergent terms. Nevertheless, when the four integrals are summed, all the divergent terms cancel and the total is exactly $0$—this is no accident but a consequence of the fact that $\log z$ is holomorphic on the whole keyhole domain. The phenomenon in which the integrals along individual paths appear to diverge while the whole is finite (indeed, $0$) appears repeatedly throughout complex analysis (for example, in residue calculations, where a real integral is evaluated by replacing it with a complex integral).
Example 26.2: Numerical check of the keyhole integral
Let $\varepsilon=0.1,\ R=10$. Compute numerically the value of the integral over each of $C_1,C_2,C_3,C_4$, and verify that the sum is exactly $0$ (to within rounding error).
Solution. We simply substitute the numerical values into the formulas of Sections 26.2.3 and 26.2.4 (alternatively, one may parametrize each path and integrate numerically). The results computed with Python (mpmath) are
$$ \int_{C_1}\approx-26.05-31.42i,\quad \int_{C_2}\approx13.36+31.10i,\quad \int_{C_3}\approx-0.66+0.31i,\quad \int_{C_4}\approx13.36 $$and the four add up to $\approx(-1.6\times10^{-30})+(-6.4\times10^{-31})i$, that is, exactly $0$ to within the rounding error of the computer (verify/univmath/ch26_logz_keyhole.py). Note that if we mistakenly computed $C_2$ with the branch left at $\theta=0$ (the same expression as on $C_4$), the sum would be $i\pi(\varepsilon-R)\approx-31.10i$, far from $0$. This shows clearly how much the result changes depending on the choice of branch alone.
Application: Keyhole contours are used to evaluate real integrals
The technique practiced here—building a keyhole-shaped contour that avoids the origin (or a branch point) and estimating the contributions of the large and small arcs—will later form the backbone of the standard method of using complex analysis to evaluate real improper integrals (for example, integrals containing logarithms or fractional powers, such as $\int_0^\infty\frac{\log x}{1+x^2}\,dx$). This book does not go as far as the residue theorem, but the keyhole idea itself is an important tool that appears repeatedly in physics as well (for example, in phase calculations in scattering theory).
26.3 Cauchy's Integral Formula — Statement and Proof
Up to the previous section, we have become familiar with four tools—Green's theorem, Cauchy's integral theorem, the contour integral formula, and the principle of path deformation—by working through concrete examples. Now we put all of them to work to state and prove, in complete form, the main theorem of this chapter—the theorem that is written as “Theorem” on p. 039 of the source notes, where its proof breaks off at “Proof…”.
26.3.1 Statement of the theorem
Theorem 26.1: Cauchy's integral formula
Let $D$ be a bounded domain whose boundary $\partial D$ is a piecewise smooth simple closed curve. Suppose that the function $w=f(z)$ is holomorphic on an open set that contains $D$ and its boundary $\partial D$ and is slightly larger than $D$. Then, for every point $z$ in the interior of $D$,
\begin{align} f(z)=\frac{1}{2\pi i}\oint_{\partial D}\frac{f(\zeta)}{\zeta-z}\,d\zeta \label{eq:26-cif} \end{align}holds. Here the integral is taken around $\partial D$ counterclockwise (the positive direction), and the integration variable is denoted by the Greek letter $\zeta$ (zeta) to distinguish it from $z$.
Intuition: The boundary information alone determines the interior
What Theorem 26.1 asserts is this: “as long as we know the values of $f$ on the boundary $\partial D$ of $D$, we can compute the value $f(z)$ at any point $z$ in the interior of $D$ from the integral in Eq. \eqref{eq:26-cif}.” A holomorphic function is so “rigid” (has so few degrees of freedom) that the boundary information alone determines everything in the interior. This resembles the situation with an X-ray photograph, in which information visible from the outside alone lets us reproduce the internal structure completely. Such a property never arises from the differentiability of a real function.
26.3.2 Proof strategy
The idea of the proof consists of the following three steps.
- First, using the principle of path deformation (Section 25.4), we replace the integral along the boundary $\partial D$ by an integral along a small circle $C_\rho:\ |\zeta-z|=\rho$ centered at the point $z$.
- We split the numerator $f(\zeta)$ of the integrand into two parts: the constant $f(z)$ and the difference $f(\zeta)-f(z)$. By the contour integral formula (Section 25.5), the constant part gives exactly $2\pi i\,f(z)$.
- For the remaining difference part, we use the continuity of $f$ and the $ML$ inequality (Section 25.1) to show that it can be made as small as we like by making $\rho$ small enough. Moreover, the value of this difference part itself must be a constant independent of $\rho$ (by the principle of path deformation), so by the reasoning that “the only constant that can be made arbitrarily small is $0$,” we conclude that the difference part is exactly $0$.
26.3.3 Proof
Proof
Step 1: Deformation to a small circle. Since $z$ is an interior point of $D$, we can choose $\rho_0\gt0$ so that a sufficiently small disk $|\zeta-z|\le\rho_0$ centered at $z$ is contained entirely in the interior of $D$. Then, for any $\rho$ satisfying $0\lt\rho\le\rho_0$, if we consider the domain $D\setminus\{|\zeta-z|\lt\rho\}$ obtained by removing the disk $|\zeta-z|\lt\rho$ from $D$, the integrand $\dfrac{f(\zeta)}{\zeta-z}$ is holomorphic on this domain (because the denominator is nowhere $0$). Therefore, by the principle of path deformation from Section 25.4 (the method of joining the outer boundary $\partial D$ and the inner circle $C_\rho:\ |\zeta-z|=\rho$ by a thin corridor, as in Figure 26.4, and regarding the whole as a single simply connected domain),
$$ \oint_{\partial D}\frac{f(\zeta)}{\zeta-z}\,d\zeta=\oint_{C_\rho}\frac{f(\zeta)}{\zeta-z}\,d\zeta $$holds for all $\rho$ satisfying $0\lt\rho\le\rho_0$ (the two integrals along the corridor, traversed back and forth, run in opposite directions and cancel each other, so they make no contribution—the same argument as in Section 25.4). Since the left-hand side is a quantity independent of $\rho$, the right-hand side also has the same value regardless of the choice of $\rho$. We will use this fact later.
Step 2: Splitting the numerator into a constant part and a difference. We rewrite the numerator $f(\zeta)$ of the integrand on the right-hand side as $f(\zeta)=f(z)+\big(f(\zeta)-f(z)\big)$ (note that $f(z)$ is a constant with respect to $\zeta$). Then
$$ \oint_{C_\rho}\frac{f(\zeta)}{\zeta-z}\,d\zeta =f(z)\oint_{C_\rho}\frac{d\zeta}{\zeta-z}+\oint_{C_\rho}\frac{f(\zeta)-f(z)}{\zeta-z}\,d\zeta $$(The integral distributes over sums—property (i) in Section 25.1. Also, since $f(z)$ is a constant, it can be taken outside the integral—property (ii) of the same section.) The integral in the first term on the right is exactly the contour integral formula (Section 25.5) itself (the case $\alpha=z$), and $\displaystyle\oint_{C_\rho}\frac{d\zeta}{\zeta-z}=2\pi i$. Therefore we can write
$$ \oint_{C_\rho}\frac{f(\zeta)}{\zeta-z}\,d\zeta=2\pi i\,f(z)+J(\rho),\qquad J(\rho):=\oint_{C_\rho}\frac{f(\zeta)-f(z)}{\zeta-z}\,d\zeta $$Combining this with the result of Step 1, we obtain $\displaystyle\oint_{\partial D}\frac{f(\zeta)}{\zeta-z}\,d\zeta-2\pi i\,f(z)=J(\rho)$, and since the left-hand side does not depend on $\rho$, we see that the value of $J(\rho)$ itself is a constant independent of $\rho$. Let us simply write $J$ for this constant.
Step 3: Showing that $J=0$. The function $f$ is continuous at $z$ (since it is holomorphic, it is in particular continuous). By the definition of continuity (the $\varepsilon$-$\delta$ argument of Chapter 1), for every $\varepsilon\gt0$ there exists a $\delta\gt0$ such that
$$ |\zeta-z|\lt\delta\ \Longrightarrow\ |f(\zeta)-f(z)|\lt\varepsilon $$holds. Now we choose $\rho$ so that $\rho\lt\min(\delta,\rho_0)$ (since $J(\rho)$ has the same value $J$ regardless of $\rho$, we may choose any $\rho$ we like). Every point $\zeta$ on the circle $C_\rho$ satisfies $|\zeta-z|=\rho\lt\delta$, so $|f(\zeta)-f(z)|\lt\varepsilon$ holds. Also, since $|\zeta-z|=\rho$ on $C_\rho$, the magnitude of the integrand is bounded by
$$ \left|\frac{f(\zeta)-f(z)}{\zeta-z}\right|=\frac{|f(\zeta)-f(z)|}{|\zeta-z|}\lt\frac{\varepsilon}{\rho} $$Since the circumference of the circle $C_\rho$ is $2\pi\rho$, we apply the $ML$ inequality (property (v) in Section 25.1: $\left|\int_C g\,d\zeta\right|\le M\cdot L$, where $M$ is the maximum of $|g|$ on the path and $L$ is the length of the path) with $g(\zeta)=\dfrac{f(\zeta)-f(z)}{\zeta-z}$, $M=\dfrac{\varepsilon}{\rho}$, and $L=2\pi\rho$, and obtain
$$ |J|=|J(\rho)|=\left|\oint_{C_\rho}\frac{f(\zeta)-f(z)}{\zeta-z}\,d\zeta\right|\le\frac{\varepsilon}{\rho}\cdot2\pi\rho=2\pi\varepsilon $$However, $\varepsilon\gt0$ is a number that can be chosen arbitrarily (as small as we like) from the start, and $J$ itself is a fixed constant that does not depend on $\varepsilon$ at all. If the statement “the constant $|J|$, which is greater than or equal to $0$, is smaller than the positive number $2\pi\varepsilon$, however small it is” can be made for every $\varepsilon\gt0$, then we must have $|J|=0$, that is, $J=0$ (if we had $|J|\gt0$, then taking $\varepsilon=|J|/(4\pi)$ would give $2\pi\varepsilon=|J|/2\lt|J|$, contradicting the claim that “$|J|\le2\pi\varepsilon$ holds for every $\varepsilon$”).
Conclusion. Since we have shown that $J=0$, returning to the equation of Step 2 gives
$$ \oint_{\partial D}\frac{f(\zeta)}{\zeta-z}\,d\zeta=2\pi i\,f(z)+0=2\pi i\,f(z) $$Dividing both sides by $2\pi i$, we obtain the desired formula, Eq. \eqref{eq:26-cif}:
$$ f(z)=\frac{1}{2\pi i}\oint_{\partial D}\frac{f(\zeta)}{\zeta-z}\,d\zeta $$∎
Caution: Three checkpoints when using the theorem
When using Cauchy's integral formula, always check the following three points. (1) Is $f$ holomorphic on a slightly larger region containing $\partial D$ (it must be $f$ itself that is holomorphic, not $\dfrac{f(\zeta)}{\zeta-z}$)? (2) Is $z$ a point in the interior of $D$, rather than a point on the boundary or in the exterior (if $z$ lies on the boundary, the integral on the right-hand side of the formula is not even meaningful, and if $z$ lies in the exterior, then by the argument for the contour integral formula in Section 25.5 the integral is instead $0$)? (3) Is $\partial D$ traversed counterclockwise? If the formula is applied mechanically in a situation that does not satisfy these conditions, a completely different (in many cases incorrect) value results.
26.3.4 Examples: Applying the formula directly
Example 26.3: $\displaystyle\oint_{|z|=1}\frac{e^z}{z}\,dz$
Using Cauchy's integral formula, compute the contour integral $\displaystyle\oint_{|z|=1}\frac{e^z}{z}\,dz$ along the circle $|z|=1$ (counterclockwise).
Solution. We view the integrand in the form $\dfrac{f(z)}{z-0}$; that is, we let $f(z)=e^z$ and $\alpha=0$. The function $f(z)=e^z$ is holomorphic on all of $\C$, and also on a region containing $D=\{|z|\le1\}$ together with its boundary. The point $\alpha=0$ is an interior point of $|z|\le1$. Therefore, by Theorem 26.1 (the case $z=0$), we have
$$ \oint_{|z|=1}\frac{e^z}{z}\,dz=\oint_{|z|=1}\frac{f(z)}{z-0}\,dz=2\pi i\,f(0)=2\pi i\,e^0=2\pi i $$(We used $e^0=1$.) There is no need to clear the denominator and compute the integral directly; all we have to do is “evaluate the numerator function at the position of the singularity and multiply by $2\pi i$,” and this is what makes the formula so valuable. A check by numerical integration (verify/univmath/ch26_cauchy_examples.py) also gives the same value $2\pi i\approx6.283i$.
Example 26.4: $\displaystyle\oint_{|z|=2}\frac{z^2}{z-1}\,dz$
Compute the contour integral $\displaystyle\oint_{|z|=2}\frac{z^2}{z-1}\,dz$ along the circle $|z|=2$ (counterclockwise).
Solution. Let $f(z)=z^2$ (holomorphic on all of $\C$) and $\alpha=1$. The point $\alpha=1$ lies inside $|z|\le2$ ($|1|=1\lt2$). By Theorem 26.1, we have
$$ \oint_{|z|=2}\frac{z^2}{z-1}\,dz=2\pi i\,f(1)=2\pi i\cdot1^2=2\pi i $$Had the singularity $\alpha=1$ been outside the path of integration (for example, had the integral been taken along $|z|=0.5$), the integrand $z^2/(z-1)$ would be holomorphic throughout the region, and so the integral would immediately be $0$ by Cauchy's integral theorem (Section 25.3). Note that the form of the answer changes completely depending on “whether or not the singularity lies inside.”
26.4 Putting Cauchy's Integral Formula to Work — Examples and Corollaries
In this section, we put Cauchy's integral formula (Theorem 26.1) to further use and look at (1) how to handle the case of several singularities and (2) two important corollaries of the formula (the mean value property and the infinite differentiability of holomorphic functions). The source notes contain nothing beyond this point (the notes end at p. 040), but these results follow naturally as applications of Theorem 26.1 and are also a gateway to later developments in complex analysis (Taylor expansion and the residue theorem), so this book supplies them on its own.
26.4.1 The case of several singularities
When the denominator of the integrand is an expression of degree two or higher that becomes $0$ at several points, it suffices to use partial fraction decomposition (the same technique used for rational expressions in high school mathematics, extended to the complex numbers) and apply Cauchy's integral formula to each singularity one by one.
Example 26.5: $\displaystyle\oint_{|z-2|=1}\frac{\sin z}{z-2}\,dz$
Compute the integral $\displaystyle\oint_{|z-2|=1}\frac{\sin z}{z-2}\,dz$ along the circle $|z-2|=1$ (center $2$, radius $1$, counterclockwise).
Solution. The function $f(z)=\sin z$ is holomorphic on all of $\C$, and $\alpha=2$ is the center of the circle $|z-2|=1$. The distance from the center to the center itself is $0$, which is less than the radius $1$ ($0\lt1$), so $\alpha=2$ is an interior point. By Theorem 26.1, we have
$$ \oint_{|z-2|=1}\frac{\sin z}{z-2}\,dz=2\pi i\,f(2)=2\pi i\sin2 $$(Here $\sin2$ is the sine of $2$ measured in radians; it is a real number that cannot be simplified any further. Numerically, $\sin2\approx0.9093$.) Numerical integration also agrees, giving $2\pi i\sin2\approx5.713i$ (verify/univmath/ch26_cauchy_examples.py).
Example 26.6: $\displaystyle\oint_{|z|=3}\frac{1}{z^2+1}\,dz$ — the sum of the contributions of two singularities
Compute the integral $\displaystyle\oint_{|z|=3}\frac{1}{z^2+1}\,dz$ along the circle $|z|=3$ (counterclockwise) as the sum of the contributions of each of the singularities $z=i,-i$.
Solution. Factoring the denominator gives $z^2+1=(z-i)(z+i)$ (because $i^2=-1$, we have $(z-i)(z+i)=z^2-i^2=z^2+1$). There are two singularities, $z=i$ and $z=-i$, and both lie inside $|z|=3$ ($|i|=|-i|=1\lt3$). As in Figure 26.5, the circle $|z|=3$ can be deformed into the sum of a small circle $C_i$ surrounding $i$ and a small circle $C_{-i}$ surrounding $-i$ (the same idea as in the application to multiply connected domains in Section 25.4).
On $C_i$, we view $\dfrac{1}{z^2+1}=\dfrac{1}{(z-i)(z+i)}=\dfrac{f_1(z)}{z-i}$ (where $f_1(z):=\dfrac{1}{z+i}$ is holomorphic near $z=i$), and by Theorem 26.1 we have
$$ \oint_{C_i}\frac{1}{z^2+1}\,dz=2\pi i\,f_1(i)=2\pi i\cdot\frac{1}{i+i}=2\pi i\cdot\frac{1}{2i}=\pi $$($\dfrac{2\pi i}{2i}=\pi$: the $i$ cancels and disappears.) Similarly, on $C_{-i}$ we view $\dfrac{1}{z^2+1}=\dfrac{f_2(z)}{z-(-i)}$ (where $f_2(z):=\dfrac{1}{z-i}$), and obtain
$$ \oint_{C_{-i}}\frac{1}{z^2+1}\,dz=2\pi i\,f_2(-i)=2\pi i\cdot\frac{1}{-i-i}=2\pi i\cdot\frac{1}{-2i}=-\pi $$Adding the two, we get
$$ \oint_{|z|=3}\frac{1}{z^2+1}\,dz=\pi+(-\pi)=0 $$The fact that the contribution $+\pi$ at $z=i$ and the contribution $-\pi$ at $z=-i$ cancel each other exactly reflects the fact that the two singularities are placed symmetrically about the origin. Numerical integration also agrees with $0$, to within rounding error (verify/univmath/ch26_cauchy_examples.py).
Example 26.7: $\displaystyle\oint_{|z-1|=0.5}\frac{z^3-1}{z-1}\,dz$
For $f(z)=z^3-1$, compute the integral $\displaystyle\oint_{|z-1|=0.5}\frac{z^3-1}{z-1}\,dz$ along the circle $|z-1|=0.5$ (counterclockwise).
Solution. The function $f(z)=z^3-1$ is holomorphic on all of $\C$, and $\alpha=1$ is the center of the circle $|z-1|=0.5$, so it lies inside. By Theorem 26.1, we have
$$ \oint_{|z-1|=0.5}\frac{z^3-1}{z-1}\,dz=2\pi i\,f(1)=2\pi i\,(1^3-1)=2\pi i\cdot0=0 $$At first glance, the denominator $z-1$ seems to create a singularity, but the numerator $f(z)=z^3-1$ happens to become $0$ exactly at $z=1$ ($1^3-1=0$), so $f(1)=0$ and the whole integral is $0$. The strength of this formula is that the answer comes out just by “evaluating the numerator at the point where the denominator becomes $0$”; note that we can carry out the computation without even noticing that $\dfrac{z^3-1}{z-1}=z^2+z+1$ (the numerator can be factored and the fraction reduced).
26.4.2 Corollary: A holomorphic function can be written as an average over a circle
In Theorem 26.1, let us take the boundary $\partial D$ to be the circle $C_\rho:\ \zeta=z+\rho e^{i\theta}$ ($0\le\theta\le2\pi$) itself, with center $z$ and radius $\rho$ (this is nothing other than Theorem 26.1 for the case where $D$ itself is the disk $|\,\zeta-z\,|\le\rho$). Substituting $d\zeta=i\rho e^{i\theta}d\theta$, we get
$$ f(z)=\frac{1}{2\pi i}\int_0^{2\pi}\frac{f(z+\rho e^{i\theta})}{\rho e^{i\theta}}\cdot i\rho e^{i\theta}\,d\theta =\frac{1}{2\pi i}\int_0^{2\pi}f(z+\rho e^{i\theta})\cdot i\,d\theta $$(The $\rho e^{i\theta}$ in the numerator and denominator cancel, and only $i$ remains.) Since $i/i=1$, we obtain
This means that “the value of a holomorphic function at the center $z$ equals the average of its values on any circle centered at that point,” and it is called the mean value property. For example, for $f(z)=e^z$, $z=1+i$, $\rho=0.7$, averaging the values $f(z+\rho e^{i\theta})$ on the circle from $\theta=0$ to $2\pi$ returns exactly $f(1+i)=e^{1+i}$, which can also be confirmed numerically (verify/univmath/ch26_cauchy_examples.py).
Why “mean value”?
Even in the world of real functions, there are functions with a similar property: the temperature distribution after heat has spread uniformly once enough time has passed (a harmonic function) has the property that the value at every point equals the average of the temperatures on a circle centered at that point. In fact, the real part $u$ and the imaginary part $v$ of a holomorphic function $f=u+iv$ both turn out to be harmonic functions (this follows from the Cauchy–Riemann equations, Chapter 24, but we do not go into it in this book), so the mean value property of Eq. \eqref{eq:26-mean-value} is also connected with the physical intuition that “the temperature inside can be determined from the temperature on the boundary.”
26.4.3 Corollary: A holomorphic function can be differentiated any number of times
We introduce another surprising consequence of Theorem 26.1, without going into the details of the proof. Let us formally differentiate the right-hand side of Cauchy's integral formula
$$ f(z)=\frac{1}{2\pi i}\oint_{\partial D}\frac{f(\zeta)}{\zeta-z}\,d\zeta $$with respect to $z$. Differentiating $\dfrac{1}{\zeta-z}$ inside the integral sign with respect to $z$ gives $\dfrac{1}{(\zeta-z)^2}$ ($\dfrac{d}{dz}(\zeta-z)^{-1}=-(\zeta-z)^{-2}\cdot(-1)=(\zeta-z)^{-2}$). It can be shown that the order of integration and differentiation may be interchanged (we omit the proof in this book, but it follows from the fact that the integrand is continuous in $\zeta$ and smooth in $z$), so we have
$$ f'(z)=\frac{1}{2\pi i}\oint_{\partial D}\frac{f(\zeta)}{(\zeta-z)^2}\,d\zeta $$Repeating the same operation, we find that in general the following formula holds for the $n$th derivative as well:
($n!$ is the factorial of $n$; for $n=0$ we have $0!=1$, and we return to Eq. \eqref{eq:26-cif}.) For example, for $f(z)=z^4$, $z=1$, $\rho=0.5$, $n=2$, computing the right-hand side of Eq. \eqref{eq:26-higher-deriv} by numerical integration gives $12$, which agrees with the value $f''(1)=12$ computed directly from $f''(z)=12z^2$ (verify/univmath/ch26_cauchy_examples.py).
Math Note: This is astonishing
For real functions, “being differentiable once (of class $C^1$)” and “being differentiable any number of times (of class $C^\infty$)” are entirely different properties, and one can construct plenty of real functions that are differentiable once but cannot be differentiated a second time. But what Eq. \eqref{eq:26-higher-deriv} shows is the fact that as long as a complex function is differentiable just once on a domain (that is, holomorphic), it is automatically differentiable any number of times in the interior of that domain. This is one of the ways in which complex analysis differs essentially from real analysis, and it also underlies the fact that a holomorphic function can be expanded in a power series (a Taylor expansion). This book stops here, but in textbooks on complex analysis (see the references at the end of the chapter), the story continues from this point on to Taylor expansions, Laurent expansions, and the residue theorem.
Application: Cauchy's integral formula and Green's functions
The structure of Cauchy's integral formula, in which “the values on the boundary are integrated to reproduce the value in the interior,” can be regarded as the complex-analytic version of the idea behind the Green's function in physics—the method of finding the field inside a region (an electric potential, a wave function, a temperature distribution, and so on) by multiplying the information about the boundary conditions and sources by an integral kernel (here $1/(\zeta-z)$) and integrating. The idea of “reconstructing the values in the interior from an integral of the boundary or initial conditions” appears again and again in many areas of science and engineering, for instance in potential problems in electrostatics and in the calculation of propagators in quantum mechanics.
26.5 Summary and Exercises
26.5.1 Summary
- We checked Green's theorem $\displaystyle\int_C f\,dx=-\iint_D\pdiff{f}{y}\,dx\,dy$ for $f(x,y)=x^2+y^2$ on the upper half-disk and on the lower half-disk, and confirmed that reversing the orientation of the boundary (the positive direction is the one in which the region lies on your left) reverses the sign.
- A multivalued function such as $\log z$ must be considered in a region that avoids the branch point ($z=0$ for $\log z$), and one branch, with the range of $\theta$ (the argument) fixed to a single choice, must be used consistently along the entire path. If the correct branch is used, the integral of $\log z$ along the keyhole contour $C_1+C_2+C_3+C_4$ is exactly $0$ regardless of the values of $\varepsilon,R$ (consistent with Cauchy's integral theorem).
- Cauchy's integral formula $\displaystyle f(z)=\frac{1}{2\pi i}\oint_{\partial D}\frac{f(\zeta)}{\zeta-z}\,d\zeta$ (Eq. \eqref{eq:26-cif})—the value at an interior point $z$ of a function $f$ holomorphic on $D$ is determined by the values on the boundary $\partial D$ alone.
- Outline of the proof: ① deform the path of integration to a small circle $C_\rho$ by the principle of path deformation; ② split the numerator into $f(z)$ (a constant) and $f(\zeta)-f(z)$, where the former gives $2\pi i\,f(z)$ by the contour integral formula; ③ for the latter, continuity and the $ML$ inequality show that, although it is a constant independent of $\rho$, it can be made arbitrarily small, and hence we conclude that it is $0$.
- A tip for computation: $\displaystyle\oint\frac{f(z)}{z-\alpha}\,dz=2\pi i\,f(\alpha)$ (when $\alpha$ lies in the interior and $f$ is holomorphic near it). When there are several singularities, decompose into partial fractions and apply the formula at each singularity.
- Corollaries: the mean value property $\displaystyle f(z)=\frac{1}{2\pi}\int_0^{2\pi}f(z+\rho e^{i\theta})\,d\theta$ (Eq. \eqref{eq:26-mean-value}), and the infinite differentiability of holomorphic functions $\displaystyle f^{(n)}(z)=\frac{n!}{2\pi i}\oint_{\partial D}\frac{f(\zeta)}{(\zeta-z)^{n+1}}\,d\zeta$ (Eq. \eqref{eq:26-higher-deriv}).
26.5.2 Exercises
Exercise 26.1: Checking in two ways
Compute $\displaystyle\oint_{|z|=2}\frac{z+1}{z(z-1)}\,dz$ in two ways, (a) using the partial fraction decomposition $\dfrac{z+1}{z(z-1)}=\dfrac{A}{z}+\dfrac{B}{z-1}$, and (b) using Cauchy's integral formula at each singularity, and verify that the two results agree.
Hint: Find $A,B$ by substituting $z=0,1$ into $z+1=A(z-1)+Bz$. In (b), regard the integrand as $f_1(z)=(z+1)/(z-1)$ near $z=0$ and as $f_2(z)=(z+1)/z$ near $z=1$.
Exercise 26.2: The meaning of Step 1 of the proof
Explain how the operation used in Step 1 of the proof of Theorem 26.1, “deforming the integral along the boundary $\partial D$ into an integral along a small circle $C_\rho$ centered at $z$,” is related to the principle of path deformation in Section 25.4. In particular, describe the role played by the “thin corridor” in Figure 26.4, relating it to the discussion of multiply connected domains in Section 25.4.
Hint: Recall what happens to the integrals along the corridor on the way out and on the way back in the example of the triply connected domain in Section 25.4 (the bridging segments $C_3,C_4$).
Exercise 26.3: Behavior of the keyhole integral in the limit
For the keyhole integral in Section 26.2, consider the limit $\varepsilon\to0$, $R\to\infty$. Investigate how each of the integrals along $C_1,C_2,C_3,C_4$ behaves (does it diverge, approach a finite value, or approach $0$?), and explain why the total is nevertheless always $0$.
Hint: In the four expressions found in Section 26.2.5, write out and compare, path by path, the coefficients of the terms $R\log R$, $\varepsilon\log\varepsilon$, $R$, $\varepsilon$, $i\pi R$, and $i\pi\varepsilon$.
Exercise 26.4: Cauchy's integral formula cannot be used for a function that is not holomorphic
Compute $\displaystyle\oint_{|z|=1}\bar z\,dz$ directly by setting $z=e^{i\theta}$. Also state why Cauchy's integral formula (or Cauchy's integral theorem from the previous chapter) cannot be applied to this function $f(z)=\bar z$.
Hint: Substituting $\bar z=e^{-i\theta}$ and $dz=ie^{i\theta}d\theta$ makes the integrand a constant. Check the condition for holomorphy (the Cauchy–Riemann equations; Chapter 24) for $\bar z=x-iy$.
Exercise 26.5: Verification of Green's theorem ($f=x^2-y^2$)
Redo the computation of Section 26.1 with $f(x,y)=x^2-y^2$ (on the upper half-disk), and verify that $\displaystyle\int_C f\,dx=-\iint_D\pdiff{f}{y}\,dx\,dy$ holds.
Hint: The computation is easier to follow if you notice that on the arc one can write $f=x^2-y^2=r^2\cos2\theta$ (the double-angle formula). Note also that $\pdiff{f}{y}=-2y$.
Exercise 26.6: The meaning of the mean value property
Retrace once more, in your own words, the steps by which the mean value property of Section 26.4.2 (Eq. \eqref{eq:26-mean-value}) is derived from Cauchy's integral formula (Eq. \eqref{eq:26-cif}), and explain, as if you were explaining it to a high-school student, what it means to say that “the value of a holomorphic function at a point is given as the average of its values on a circle surrounding that point.”
Hint: Follow how $d\zeta$ changes when we set $\zeta=z+\rho e^{i\theta}$, and which part of the integrand cancels and disappears.
26.5.3 References
- Yasuhide Mochizuki, Mathematics Notes: Complex Analysis (handwritten lecture notes, in Japanese). The source of this chapter.
- Reiji Takahashi, Fukuso Kaiseki [Complex Analysis], University of Tokyo Press, 1990 (in Japanese).
- Lars V. Ahlfors, Fukuso Kaiseki [Complex Analysis], translated by Kenkichi Kasahara, Gendai Sūgakusha, 1982 (Japanese translation of Complex Analysis; in Japanese).
- Teiji Takagi, Kaiseki Gairon [Outline of Analysis], 3rd revised ed., Iwanami Shoten, 1983 (in Japanese).