# Cauchy's integral formula

In mathematics, Cauchy's integral formula, named after Augustin-Louis Cauchy, is a central statement in complex analysis. It expresses the fact that a holomorphic function defined on a disk is completely determined by its values on the boundary of the disk, and it provides integral formulas for all derivatives of a holomorphic function. Cauchy's formula shows that, in complex analysis, "differentiation is equivalent to integration": complex differentiation, like integration, behaves well under uniform limits – a result that does not hold in real analysis.

## Theorem

Let U be an open subset of the complex plane C, and suppose the closed disk D defined as $D={\bigl \z:|z-z_{0}|\leq r{\bigr \}}$ is completely contained in U. Let f : U → C be a holomorphic function, and let γ be the circle, oriented counterclockwise, forming the boundary of D. Then for every a in the interior of D, $f(a)={\frac {1}{2\pi i}}\oint _{\gamma }{\frac {f(z)}{z-a}}\,dz.\,$
The proof of this statement uses the Cauchy integral theorem and like that theorem, it only requires f to be complex differentiable. Since $1/(z-a)$ can be expanded as a power series in the variable $a$ ${\frac {1}{z-a}}={\frac {1+{\frac {a}{z}}+\left({\frac {a}{z}}\right)^{2}+\cdots }{z}}$ it follows that holomorphic functions are analytic, i.e. they can be expanded as convergent power series. In particular f is actually infinitely differentiable, with $f^{(n)}(a)={\frac {n!}{2\pi i}}\oint _{\gamma }{\frac {f(z)}{\left(z-a\right)^{n+1}}}\,dz.$
This formula is sometimes referred to as Cauchy's differentiation formula.
The theorem stated above can be generalized. The circle γ can be replaced by any closed rectifiable curve in U which has winding number one about a. Moreover, as for the Cauchy integral theorem, it is sufficient to require that f be holomorphic in the open region enclosed by the path and continuous on its closure.
Note that not every continuous function on the boundary can be used to produce a function inside the boundary that fits the given boundary function. For instance, if we put the function f(z) = ⁠1/z⁠, defined for |z| = 1, into the Cauchy integral formula, we get zero for all points inside the circle. In fact, giving just the real part on the boundary of a holomorphic function is enough to determine the function up to an imaginary constant — there is only one imaginary part on the boundary that corresponds to the given real part, up to addition of a constant. We can use a combination of a Möbius transformation and the Stieltjes inversion formula to construct the holomorphic function from the real part on the boundary. For example, the function f(z) = i − iz has real part Re f(z) = Im z. On the unit circle this can be written ⁠⁠i/z⁠ − iz/2⁠. Using the Möbius transformation and the Stieltjes formula we construct the function inside the circle. The ⁠i/z⁠ term makes no contribution, and we find the function −iz. This has the correct real part on the boundary, and also gives us the corresponding imaginary part, but off by a constant, namely i.

## Proof sketch

By using the Cauchy integral theorem, one can show that the integral over C (or the closed rectifiable curve) is equal to the same integral taken over an arbitrarily small circle around a. Since f(z) is continuous, we can choose a circle small enough on which f(z) is arbitrarily close to f(a). On the other hand, the integral $\oint _{C}{\frac {1}{z-a}}\,dz=2\pi i,$ over any circle C centered at a. This can be calculated directly via a parametrization (integration by substitution) z(t) = a + εe$^{it}$ where 0 ≤ t ≤ 2π and ε is the radius of the circle.
Letting ε → 0 gives the desired estimate ${\begin{aligned}\left|{\frac {1}{2\pi i}}\oint _{C}{\frac {f(z)}{z-a}}\,dz-f(a)\right|&=\left|{\frac {1}{2\pi i}}\oint _{C}{\frac {f(z)-f(a)}{z-a}}\,dz\right|\\[1ex]&=\left|{\frac {1}{2\pi i}}\int _{0}^{2\pi }\left({\frac {f{\biglz(t){\bigr )}-f(a)}{\varepsilon e^{it}}}\cdot \varepsilon e^{it}i\right)\,dt\right|\\[1ex]&\leq {\frac {1}{2\pi }}\int _{0}^{2\pi }{\frac {\left|f{\biglz(t){\bigr )}-f(a)\right|}{\varepsilon }}\,\varepsilon \,dt\\[1ex]&\leq \max _{|z-a|=\varepsilon }\left|f(z)-f(a)\right|~~{\xrightarrow[{\varepsilon \to 0}]}~~0.\end{aligned}}$

## Example

Let $g(z)={\frac {z^{2}}{z^{2}+2z+2}},$ and let C be the contour described by |z| = 2 (the circle of radius 2).
To find the integral of g(z) around the contour C, we need to know the singularities of g(z). Observe that we can rewrite g as follows: $g(z)={\frac {z^{2}}{(z-z_{1})(z-z_{2})}}$ where z$_1$ = − 1 + i and z$_2$ = − 1 − i.
Thus, g has poles at z$_1$ and z$_2$. The moduli of these points are less than 2 and thus lie inside the contour. This integral can be split into two smaller integrals by Cauchy–Goursat theorem; that is, we can express the integral around the contour as the sum of the integral around z$_1$ and z$_2$ where the contour is a small circle around each pole. Call these contours C$_1$ around z$_1$ and C$_2$ around z$_2$.
Now, each of these smaller integrals can be evaluated by the Cauchy integral formula, but they first must be rewritten to apply the theorem. For the integral around C$_1$, define f$_1$ as f$_1$(z) = (z − z$_1$)g(z). This is analytic (since the contour does not contain the other singularity). We can simplify f$_1$ to be: $f_{1}(z)={\frac {z^{2}}{z-z_{2}}}$ and now $g(z)={\frac {f_{1}(z)}{z-z_{1}}}.$
Since the Cauchy integral formula says that: $\oint _{C}{\frac {f_{1}(z)}{z-a}}\,dz=2\pi i\cdot f_{1}(a),$ we can evaluate the integral as follows: $\oint _{C_{1}}g(z)\,dz=\oint _{C_{1}}{\frac {f_{1}(z)}{z-z_{1}}}\,dz=2\pi i{\frac {z_{1}^{2}}{z_{1}-z_{2}}}.$
Doing likewise for the other contour: $f_{2}(z)={\frac {z^{2}}{z-z_{1}}},$ we evaluate $\oint _{C_{2}}g(z)\,dz=\oint _{C_{2}}{\frac {f_{2}(z)}{z-z_{2}}}\,dz=2\pi i{\frac {z_{2}^{2}}{z_{2}-z_{1}}}.$
The integral around the original contour C then is the sum of these two integrals: ${\begin{aligned}\oint _{C}g(z)\,dz&=\oint _{C_{1}}g(z)\,dz+\oint _{C_{2}}g(z)\,dz\\[.5em]&=2\pi i\left({\frac {z_{1}^{2}}{z_{1}-z_{2}}}+{\frac {z_{2}^{2}}{z_{2}-z_{1}}}\right)\\[.5em]&=2\pi i(-2)\\[.3em]&=-4\pi i.\end{aligned}}$
An elementary trick using partial fraction decomposition: $\oint _{C}g(z)\,dz=\oint _{C}\left(1-{\frac {1}{z-z_{1}}}-{\frac {1}{z-z_{2}}}\right)\,dz=0-2\pi i-2\pi i=-4\pi i$

## Consequences

The integral formula has broad applications. First, it implies that a function which is holomorphic in an open set is in fact infinitely differentiable there. Furthermore, it is an analytic function, meaning that it can be represented as a power series. The proof of this uses the dominated convergence theorem and the geometric series applied to
$f(\zeta )={\frac {1}{2\pi i}}\int _{C}{\frac {f(z)}{z-\zeta }}\,dz.$
The formula is also used to prove the residue theorem, which is a result for meromorphic functions, and a related result, the argument principle. It is known from Morera's theorem that the uniform limit of holomorphic functions is holomorphic. This can also be deduced from Cauchy's integral formula: indeed the formula also holds in the limit and the integrand, and hence the integral, can be expanded as a power series. In addition the Cauchy formulas for the higher order derivatives show that all these derivatives also converge uniformly.
The analog of the Cauchy integral formula in real analysis is the Poisson integral formula for harmonic functions; many of the results for holomorphic functions carry over to this setting. No such results, however, are valid for more general classes of differentiable or real analytic functions. For instance, the existence of the first derivative of a real function need not imply the existence of higher order derivatives, nor in particular the analyticity of the function. Likewise, the uniform limit of a sequence of (real) differentiable functions may fail to be differentiable, or may be differentiable but with a derivative which is not the limit of the derivatives of the members of the sequence.
Another consequence is that if f(z) = Σ a$_n$ z$^n$ is holomorphic in |z| < R and 0 < r < R then the coefficients a$_n$ satisfy Cauchy's estimate $|a_{n}|\leq r^{-n}\sup _{|z|=r}|f(z)|.$
From Cauchy's estimate, one can easily deduce that every bounded entire function must be constant (which is Liouville's theorem).
The formula can also be used to derive Gauss's Mean-Value Theorem, which states $f(z)={\frac {1}{2\pi }}\int _{0}^{2\pi }f(z+re^{i\theta })\,d\theta .$
In other words, the average value of f over the circle centered at z with radius r is f(z). This can be calculated directly via a parametrization of the circle.

## Generalizations
### Smooth functions

A version of Cauchy's integral formula is the Cauchy–Pompeiu formula, and holds for smooth functions as well, as it is based on Stokes' theorem. Let D be a disc in C and suppose that f is a complex-valued C$^1$ function on the closure of D. Then $f(\zeta )={\frac {1}{2\pi i}}\int _{\partial D}{\frac {f(z)\,dz}{z-\zeta }}-{\frac {1}{\pi }}\iint _{D}{\frac {\partial f}{\partial {\bar {z}}}}(z){\frac {dx\wedge dy}{z-\zeta }}.$
One may use this representation formula to solve the inhomogeneous Cauchy–Riemann equations in D. Indeed, if φ is a function in D, then a particular solution f of the equation is a holomorphic function outside the support of μ. Moreover, if in an open set D, $d\mu ={\frac {1}{2\pi i}}\varphi \,dz\wedge d{\bar {z}}$ for some φ ∈ C$^k$(D) (where k ≥ 1), then f(ζ, ζ) is also in C$^k$(D) and satisfies the equation ${\frac {\partial f}{\partial {\bar {z}}}}=\varphi (z,{\bar {z}}).$
The first conclusion is, succinctly, that the convolution μ ∗ k(z) of a compactly supported measure with the Cauchy kernel $k(z)=\operatorname {p.v.} {\frac {1}{z}}$ is a holomorphic function off the support of μ. Here p.v. denotes the principal value. The second conclusion asserts that the Cauchy kernel is a fundamental solution of the Cauchy–Riemann equations. Note that for smooth complex-valued functions f of compact support on C the generalized Cauchy integral formula simplifies to $f(\zeta )={\frac {1}{2\pi i}}\iint {\frac {\partial f}{\partial {\bar {z}}}}{\frac {dz\wedge d{\bar {z}}}{z-\zeta }},$ and is a restatement of the fact that, considered as a distribution, (πz)$^{−1}$ is a fundamental solution of the Cauchy–Riemann operator ⁠∂/∂z̄⁠.
The generalized Cauchy integral formula can be deduced for any bounded open region X with C$^1$ boundary ∂X from this result and the formula for the distributional derivative of the characteristic function χ$_X$ of X: ${\frac {\partial \chi _{X}}{\partial {\bar {z}}}}={\frac {i}{2}}\oint _{\partial X}\,dz,$ where the distribution on the right hand side denotes contour integration along ∂X.
For $\varphi \in {\mathcal {D}}(X)$ calculate:
    ${\begin{aligned}\left\langle {\frac {\partial }{\partial {\bar {z}}}}\left(\chi _{X}\right),\varphi \right\rangle &=-\int _{X}{\frac {\partial \varphi }{\partial {\bar {z}}}}\mathrm {~d} (x,y)\\&=-{\frac {1}{2}}\int _{X}\left(\partial _{x}\varphi +\mathrm {i} \partial _{y}\varphi \right)\mathrm {d} (x,y).\end{aligned}}$
then traverse $\partial X$ in the anti-clockwise direction. Fix a point $p\in \partial X$ and let $s$ denote arc length on $\partial X$ measured from $p$ anti-clockwise. Then, if $\ell$ is the length of $\partial X,[0,\ell ]\ni s\mapsto (x(s),y(s))$ is a parametrization of $\partial X$. The derivative $\tau =\left(x'(s),y'(s)\right)$ is a unit tangent to $\partial X$ and $\nu :=\left(-y'(s),x'(s)\right)$ is the unit outward normal on $\partial X$. We are lined up for use of the divergence theorem: put $V=(\varphi ,\mathrm {i} \varphi )\in {\mathcal {D}}(X)^{2}$ so that $\operatorname {div} V=\partial _{x}\varphi +\mathrm {i} \partial _{y}\varphi$ and we get
    ${\begin{aligned}-{\frac {1}{2}}\int _{X}\left(\partial _{x}\varphi +\mathrm {i} \partial _{y}\varphi \right)\mathrm {d} (x,y)&=-{\frac {1}{2}}\int _{\partial X}V\cdot \nu \mathrm {d} S\\&=-{\frac {1}{2}}\int _{0}^{\ell }\left(\varphi \nu _{1}+\mathrm {i} \varphi \nu _{2}\right)\mathrm {d} s\\&=-{\frac {1}{2}}\int _{0}^{\ell }\varphi (x(s),y(s))\left(y'(s)-\mathrm {i} x'(s)\right)\mathrm {d} s\\&={\frac {1}{2}}\int _{0}^{\ell }\mathrm {i} \varphi (x(s),y(s))\left(x'(s)+\mathrm {i} y'(s)\right)\mathrm {d} s\\&={\frac {\mathrm {i} }{2}}\int _{\partial X}\varphi \mathrm {d} z\end{aligned}}$
Hence we proved ${\frac {\partial \chi _{X}}{\partial {\bar {z}}}}={\frac {i}{2}}\oint _{\partial X}\,dz$.
Now we can deduce the generalized Cauchy integral formula:
Since $u={\frac {\chi _{X}}{\pi \left(z-z_{0}\right)}}\in \mathrm {L} _{\text{loc}}^{1}(X)$ and since $z_{0}\in X$ this distribution is locally in $X$ of the form "distribution times C$^∞$ function", so we may apply the Leibniz rule to calculate its derivatives:
    ${\frac {\partial u}{\partial {\bar {z}}}}={\frac {\partial }{\partial {\bar {z}}}}\left({\frac {1}{\pi \left(z-z_{0}\right)}}\right)\chi _{X}+{\frac {1}{\pi \left(z-z_{0}\right)}}{\frac {\partial }{\partial {\bar {z}}}}\left(\chi _{X}\right)$
Using that (πz)$^{−1}$ is a fundamental solution of the Cauchy–Riemann operator ⁠∂/∂z̄⁠, we get ${\frac {\partial }{\partial {\bar {z}}}}\left({\frac {1}{\pi \left(z-z_{0}\right)}}\right)=\delta _{z_{0}}$:
    ${\frac {\partial u}{\partial {\bar {z}}}}=\delta _{z_{0}}+{\frac {1}{\pi \left(z-z_{0}\right)}}{\frac {\partial }{\partial {\bar {z}}}}\left(\chi _{X}\right)$
Applying ${\frac {\partial u}{\partial {\bar {z}}}}$ to $\phi \in {\mathcal {D}}(X)$:
    ${\begin{aligned}\left\langle {\frac {\partial }{\partial {\bar {z}}}}\left({\frac {\chi _{X}}{\pi \left(z-z_{0}\right)}}\right),\phi \right\rangle &=\phi \left(z_{0}\right)+\left\langle {\frac {1}{\pi \left(z-z_{0}\right)}}{\frac {\partial }{\partial {\bar {z}}}}\left(\chi _{X}\right),\phi \right\rangle \\&=\phi \left(z_{0}\right)+\left\langle {\frac {\partial }{\partial {\bar {z}}}}\left(\chi _{X}\right),{\frac {\phi }{\pi \left(z-z_{0}\right)}}\right\rangle \\&=\phi \left(z_{0}\right)+{\frac {\mathrm {i} }{2}}\int _{\partial X}{\frac {\phi (z)}{\pi \left(z-z_{0}\right)}}\mathrm {d} z\end{aligned}}$
where ${\frac {\partial \chi _{X}}{\partial {\bar {z}}}}={\frac {i}{2}}\oint _{\partial X}\,dz$ is used in the last line.
Rearranging, we get
    $\phi (z_{0})={\frac {1}{2\pi i}}\int _{\partial X}{\frac {\phi (z)\,dz}{z-z_{0}}}-{\frac {1}{\pi }}\iint _{X}{\frac {\partial \phi }{\partial {\bar {z}}}}(z){\frac {dx\wedge dy}{z-z_{0}}}.$
as desired.

### Several variables

In several complex variables, the Cauchy integral formula can be generalized to polydiscs. Let D be the polydisc given as the Cartesian product of n open discs D$_1$, ..., D$_n$: $D=\prod _{i=1}^{n}D_{i}.$
Suppose that f is a holomorphic function in D continuous on the closure of D. Then $f(\zeta )={\frac {1}{\left(2\pi i\right)^{n}}}\int \cdots \iint _{\partial D_{1}\times \cdots \times \partial D_{n}}{\frac {f(z_{1},\ldots ,z_{n})}{(z_{1}-\zeta _{1})\cdots (z_{n}-\zeta _{n})}}\,dz_{1}\cdots dz_{n}$ where ζ = (ζ$_1$,...,ζ$_n$) ∈ D.

### In real algebras

The Cauchy integral formula is generalizable to real vector spaces of two or more dimensions. The insight into this property comes from geometric algebra, where objects beyond scalars and vectors (such as planar bivectors and volumetric trivectors) are considered, and a proper generalization of Stokes' theorem.
Geometric calculus defines a derivative operator ∇ = ê$_i$ ∂$_i$ under its geometric product — that is, for a k-vector field ψ(r), the derivative ∇ψ generally contains terms of grade k + 1 and k − 1. For example, a vector field (k = 1) generally has in its derivative a scalar part, the divergence (k = 0), and a bivector part, the curl (k = 2). This particular derivative operator has a Green's function: $G\left(\mathbf {r} ,\mathbf {r} '\right)={\frac {1}{S_{n}}}{\frac {\mathbf {r} -\mathbf {r} '}{\left|\mathbf {r} -\mathbf {r} '\right|^{n}}}$ where S$_n$ is the surface area of a unit n-ball in the space (that is, S$_2$ = 2π, the circumference of a circle with radius 1, and S$_3$ = 4π, the surface area of a sphere with radius 1). By definition of a Green's function, $\nabla G\left(\mathbf {r} ,\mathbf {r} '\right)=\delta \left(\mathbf {r} -\mathbf {r} '\right).$
It is this useful property that can be used, in conjunction with the generalized Stokes theorem: $\oint _{\partial V}d\mathbf {S} \;f(\mathbf {r} )=\int _{V}d\mathbf {V} \;\nabla f(\mathbf {r} )$ where, for an n-dimensional vector space, dS is an (n − 1)-vector and dV is an n-vector. The function f(r) can, in principle, be composed of any combination of multivectors. The proof of Cauchy's integral theorem for higher dimensional spaces relies on the using the generalized Stokes theorem on the quantity G(r, r′) f(r′) and use of the product rule: $\oint _{\partial V'}G\left(\mathbf {r} ,\mathbf {r} '\right)\;d\mathbf {S} '\;f\left(\mathbf {r} '\right)=\int _{V}\left(\left[\nabla 'G\left(\mathbf {r} ,\mathbf {r} '\right)\right]f\left(\mathbf {r} '\right)+G\left(\mathbf {r} ,\mathbf {r} '\right)\nabla 'f\left(\mathbf {r} '\right)\right)\;d\mathbf {V}$
When ∇f = 0, f(r) is called a monogenic function, the generalization of holomorphic functions to higher-dimensional spaces — indeed, it can be shown that the Cauchy–Riemann condition is just the two-dimensional expression of the monogenic condition. When that condition is met, the second term in the right-hand integral vanishes, leaving only $\oint _{\partial V'}G\left(\mathbf {r} ,\mathbf {r} '\right)\;d\mathbf {S} '\;f\left(\mathbf {r} '\right)=\int _{V}\left[\nabla 'G\left(\mathbf {r} ,\mathbf {r} '\right)\right]f\left(\mathbf {r} '\right)=-\int _{V}\delta \left(\mathbf {r} -\mathbf {r} '\right)f\left(\mathbf {r} '\right)\;d\mathbf {V} =-i_{n}f(\mathbf {r} )$ where i$_n$ is that algebra's unit n-vector, the pseudoscalar. The result is $f(\mathbf {r} )=-{\frac {1}{i_{n}}}\oint _{\partial V}G\left(\mathbf {r} ,\mathbf {r} '\right)\;d\mathbf {S} \;f\left(\mathbf {r} '\right)=-{\frac {1}{i_{n}}}\oint _{\partial V}{\frac {\mathbf {r} -\mathbf {r} '}{S_{n}\left|\mathbf {r} -\mathbf {r} '\right|^{n}}}\;d\mathbf {S} \;f\left(\mathbf {r} '\right)$
Thus, as in the two-dimensional (complex analysis) case, the value of an analytic (monogenic) function at a point can be found by an integral over the surface surrounding the point, and this is valid not only for scalar functions but vector and general multivector functions as well.