Contour integration

Last updated

In the mathematical field of complex analysis, contour integration is a method of evaluating certain integrals along paths in the complex plane. [1] [2] [3]

Contents

Contour integration is closely related to the calculus of residues, [4] a method of complex analysis.

One use for contour integrals is the evaluation of integrals along the real line that are not readily found by using only real variable methods. [5]

Contour integration methods include:

One method can be used, or a combination of these methods, or various limiting processes, for the purpose of finding these integrals or sums.

Curves in the complex plane

In complex analysis a contour is a type of curve in the complex plane. In contour integration, contours provide a precise definition of the curves on which an integral may be suitably defined. A curve in the complex plane is defined as a continuous function from a closed interval of the real line to the complex plane: .

This definition of a curve coincides with the intuitive notion of a curve, but includes a parametrization by a continuous function from a closed interval. This more precise definition allows us to consider what properties a curve must have for it to be useful for integration. In the following subsections we narrow down the set of curves that we can integrate to include only those that can be built up out of a finite number of continuous curves that can be given a direction. Moreover, we will restrict the "pieces" from crossing over themselves, and we require that each piece have a finite (non-vanishing) continuous derivative. These requirements correspond to requiring that we consider only curves that can be traced, such as by a pen, in a sequence of even, steady strokes, which stop only to start a new piece of the curve, all without picking up the pen. [6]

Directed smooth curves

Contours are often defined in terms of directed smooth curves. [6] These provide a precise definition of a "piece" of a smooth curve, of which a contour is made.

A smooth curve is a curve with a non-vanishing, continuous derivative such that each point is traversed only once (z is one-to-one), with the possible exception of a curve such that the endpoints match (). In the case where the endpoints match, the curve is called closed, and the function is required to be one-to-one everywhere else and the derivative must be continuous at the identified point (). A smooth curve that is not closed is often referred to as a smooth arc. [6]

The parametrization of a curve provides a natural ordering of points on the curve: comes before if . This leads to the notion of a directed smooth curve. It is most useful to consider curves independent of the specific parametrization. This can be done by considering equivalence classes of smooth curves with the same direction. A directed smooth curve can then be defined as an ordered set of points in the complex plane that is the image of some smooth curve in their natural order (according to the parametrization). Note that not all orderings of the points are the natural ordering of a smooth curve. In fact, a given smooth curve has only two such orderings. Also, a single closed curve can have any point as its endpoint, while a smooth arc has only two choices for its endpoints.

Contours

Contours are the class of curves on which we define contour integration. A contour is a directed curve which is made up of a finite sequence of directed smooth curves whose endpoints are matched to give a single direction. This requires that the sequence of curves be such that the terminal point of coincides with the initial point of for all such that . This includes all directed smooth curves. Also, a single point in the complex plane is considered a contour. The symbol is often used to denote the piecing of curves together to form a new curve. Thus we could write a contour that is made up of curves as

Contour integrals

The contour integral of a complex function is a generalization of the integral for real-valued functions. For continuous functions in the complex plane, the contour integral can be defined in analogy to the line integral by first defining the integral along a directed smooth curve in terms of an integral over a real valued parameter. A more general definition can be given in terms of partitions of the contour in analogy with the partition of an interval and the Riemann integral. In both cases the integral over a contour is defined as the sum of the integrals over the directed smooth curves that make up the contour.

For continuous functions

To define the contour integral in this way one must first consider the integral, over a real variable, of a complex-valued function. Let be a complex-valued function of a real variable, . The real and imaginary parts of are often denoted as and , respectively, so that Then the integral of the complex-valued function over the interval is given by

Now, to define the contour integral, let be a continuous function on the directed smooth curve . Let be any parametrization of that is consistent with its order (direction). Then the integral along is denoted and is given by [6]

This definition is well defined. That is, the result is independent of the parametrization chosen. [6] In the case where the real integral on the right side does not exist the integral along is said not to exist.

As a generalization of the Riemann integral

The generalization of the Riemann integral to functions of a complex variable is done in complete analogy to its definition for functions from the real numbers. The partition of a directed smooth curve is defined as a finite, ordered set of points on . The integral over the curve is the limit of finite sums of function values, taken at the points on the partition, in the limit that the maximum distance between any two successive points on the partition (in the two-dimensional complex plane), also known as the mesh, goes to zero.

Direct methods

Direct methods involve the calculation of the integral through methods similar to those in calculating line integrals in multivariate calculus. This means that we use the following method:

Example

A fundamental result in complex analysis is that the contour integral of 1/z is i, where the path of the contour is taken to be the unit circle traversed counterclockwise (or any positively oriented Jordan curve about 0). In the case of the unit circle there is a direct method to evaluate the integral

In evaluating this integral, use the unit circle |z| = 1 as a contour, parametrized by z(t) = eit, with t ∈ [0, 2π], then dz/dt = ieit and

which is the value of the integral. This result only applies to the case in which z is raised to power of -1. If the power is not equal to -1, then the result will always be zero.

Applications of integral theorems

Applications of integral theorems are also often used to evaluate the contour integral along a contour, which means that the real-valued integral is calculated simultaneously along with calculating the contour integral.

Integral theorems such as the Cauchy integral formula or residue theorem are generally used in the following method:

Example 1

Consider the integral

To evaluate this integral, we look at the complex-valued function

which has singularities at i and i. We choose a contour that will enclose the real-valued integral, here a semicircle with boundary diameter on the real line (going from, say, a to a) will be convenient. Call this contour C.

There are two ways of proceeding, using the Cauchy integral formula or by the method of residues:

Using the Cauchy integral formula

Note that: thus

Furthermore, observe that

Since the only singularity in the contour is the one at i, then we can write

which puts the function in the form for direct application of the formula. Then, by using Cauchy's integral formula,

We take the first derivative, in the above steps, because the pole is a second-order pole. That is, (zi) is taken to the second power, so we employ the first derivative of f(z). If it were (zi) taken to the third power, we would use the second derivative and divide by 2!, etc. The case of (zi) to the first power corresponds to a zero order derivative—just f(z) itself.

We need to show that the integral over the arc of the semicircle tends to zero as a → ∞, using the estimation lemma

where M is an upper bound on |f(z)| along the arc and L the length of the arc. Now, So

Using the method of residues

Consider the Laurent series of f(z) about i, the only singularity we need to consider. We then have

(See the sample Laurent calculation from Laurent series for the derivation of this series.)

It is clear by inspection that the residue is i/4, so, by the residue theorem, we have

Thus we get the same result as before.

Contour note

As an aside, a question can arise whether we do not take the semicircle to include the other singularity, enclosing i. To have the integral along the real axis moving in the correct direction, the contour must travel clockwise, i.e., in a negative direction, reversing the sign of the integral overall.

This does not affect the use of the method of residues by series.

Example 2 – Cauchy distribution

The integral

the contour ContourDiagram.png
the contour

(which arises in probability theory as a scalar multiple of the characteristic function of the Cauchy distribution) resists the techniques of elementary calculus. We will evaluate it by expressing it as a limit of contour integrals along the contour C that goes along the real line from a to a and then counterclockwise along a semicircle centered at 0 from a to a. Take a to be greater than 1, so that the imaginary unit i is enclosed within the curve. The contour integral is

Since eitz is an entire function (having no singularities at any point in the complex plane), this function has singularities only where the denominator z2 + 1 is zero. Since z2 + 1 = (z + i)(zi), that happens only where z = i or z = −i. Only one of those points is in the region bounded by this contour. The residue of f(z) at z = i is

According to the residue theorem, then, we have

The contour C may be split into a "straight" part and a curved arc, so that and thus

According to Jordan's lemma, if t > 0 then

Therefore, if t > 0 then

A similar argument with an arc that winds around i rather than i shows that if t < 0 then and finally we have this:

(If t = 0 then the integral yields immediately to real-valued calculus methods and its value is π.)

Example 3 – trigonometric integrals

Certain substitutions can be made to integrals involving trigonometric functions, so the integral is transformed into a rational function of a complex variable and then the above methods can be used in order to evaluate the integral.

As an example, consider

We seek to make a substitution of z = eit. Now, recall and

Taking C to be the unit circle, we substitute to get:

The singularities to be considered are at Let C1 be a small circle about and C2 be a small circle about Then we arrive at the following:

Example 3a – trigonometric integrals, the general procedure

The above method may be applied to all integrals of the type

where P and Q are polynomials, i.e. a rational function in trigonometric terms is being integrated. Note that the bounds of integration may as well be π and −π, as in the previous example, or any other pair of endpoints 2π apart.

The trick is to use the substitution z = eit where dz = ieit dt and hence

This substitution maps the interval [0, 2π] to the unit circle. Furthermore, and so that a rational function f(z) in z results from the substitution, and the integral becomes which is in turn computed by summing the residues of f(z)1/iz inside the unit circle.

TrigonometricToComplex.png

The image at right illustrates this for which we now compute. The first step is to recognize that

The substitution yields

The poles of this function are at 1 ± 2 and −1 ± 2. Of these, 1 + 2 and −1 − 2 are outside the unit circle (shown in red, not to scale), whereas 1 − 2 and −1 + 2 are inside the unit circle (shown in blue). The corresponding residues are both equal to i2/16, so that the value of the integral is

Example 4 – branch cuts

Consider the real integral

We can begin by formulating the complex integral

Keyhole contour.svg

We can use the Cauchy integral formula or residue theorem again to obtain the relevant residues. However, the important thing to note is that z1/2 = e(Log z)/2, so z1/2 has a branch cut. This affects our choice of the contour C. Normally the logarithm branch cut is defined as the negative real axis, however, this makes the calculation of the integral slightly more complicated, so we define it to be the positive real axis.

Then, we use the so-called keyhole contour, which consists of a small circle about the origin of radius ε say, extending to a line segment parallel and close to the positive real axis but not touching it, to an almost full circle, returning to a line segment parallel, close, and below the positive real axis in the negative sense, returning to the small circle in the middle.

Note that z = −2 and z = −4 are inside the big circle. These are the two remaining poles, derivable by factoring the denominator of the integrand. The branch point at z = 0 was avoided by detouring around the origin.

Let γ be the small circle of radius ε, Γ the larger, with radius R, then

It can be shown that the integrals over Γ and γ both tend to zero as ε → 0 and R → ∞, by an estimation argument above, that leaves two terms. Now since z1/2 = e(Log z)/2, on the contour outside the branch cut, we have gained 2π in argument along γ. (By Euler's identity, eiπ represents the unit vector, which therefore has π as its log. This π is what is meant by the argument of z. The coefficient of 1/2 forces us to use 2π.) So

Therefore:

By using the residue theorem or the Cauchy integral formula (first employing the partial fractions method to derive a sum of two simple contour integrals) one obtains

Example 5 – the square of the logarithm

KeyholeContourLeftTikz.tif

This section treats a type of integral of which is an example.

To calculate this integral, one uses the function and the branch of the logarithm corresponding to −π < arg z ≤ π.

We will calculate the integral of f(z) along the keyhole contour shown at right. As it turns out this integral is a multiple of the initial integral that we wish to calculate and by the Cauchy residue theorem we have

Let R be the radius of the large circle, and r the radius of the small one. We will denote the upper line by M, and the lower line by N. As before we take the limit when R → ∞ and r → 0. The contributions from the two circles vanish. For example, one has the following upper bound with the ML lemma:

In order to compute the contributions of M and N we set z = −x + on M and z = −x on N, with 0 < x < ∞:

which gives

Example 6 – logarithms and the residue at infinity

ContourLogs.png

We seek to evaluate

This requires a close study of

We will construct f(z) so that it has a branch cut on [0, 3], shown in red in the diagram. To do this, we choose two branches of the logarithm, setting and

The cut of z34 is therefore (−∞, 0] and the cut of (3 − z)1/4 is (−∞, 3]. It is easy to see that the cut of the product of the two, i.e. f(z), is [0, 3], because f(z) is actually continuous across (−∞, 0). This is because when z = −r < 0 and we approach the cut from above, f(z) has the value

When we approach from below, f(z) has the value

But

so that we have continuity across the cut. This is illustrated in the diagram, where the two black oriented circles are labelled with the corresponding value of the argument of the logarithm used in z34 and (3 − z)1/4.

We will use the contour shown in green in the diagram. To do this we must compute the value of f(z) along the line segments just above and just below the cut.

Let z = r (in the limit, i.e. as the two green circles shrink to radius zero), where 0 ≤ r ≤ 3. Along the upper segment, we find that f(z) has the value and along the lower segment,

It follows that the integral of f(z)/5 − z along the upper segment is iI in the limit, and along the lower segment, I.

If we can show that the integrals along the two green circles vanish in the limit, then we also have the value of I, by the Cauchy residue theorem. Let the radius of the green circles be ρ, where ρ < 0.001 and ρ → 0, and apply the ML inequality. For the circle CL on the left, we find

Similarly, for the circle CR on the right, we have

Now using the Cauchy residue theorem, we have where the minus sign is due to the clockwise direction around the residues. Using the branch of the logarithm from before, clearly

The pole is shown in blue in the diagram. The value simplifies to

We use the following formula for the residue at infinity:

Substituting, we find and where we have used the fact that −1 = eπi for the second branch of the logarithm. Next we apply the binomial expansion, obtaining

The conclusion is that

Finally, it follows that the value of I is which yields

Evaluation with residue theorem

Using the residue theorem, we can evaluate closed contour integrals. The following are examples on evaluating contour integrals with the residue theorem.

Using the residue theorem, let us evaluate this contour integral.

Recall that the residue theorem states

where is the residue of , and the are the singularities of lying inside the contour (with none of them lying directly on ).

has only one pole, . From that, we determine that the residue of to be

Thus, using the residue theorem, we can determine:

Multivariable contour integrals

To solve multivariable contour integrals (i.e. surface integrals, complex volume integrals, and higher order integrals), we must use the divergence theorem. For now, let be interchangeable with . These will both serve as the divergence of the vector field denoted as . This theorem states:

In addition, we also need to evaluate where is an alternate notation of . The divergence of any dimension can be described as

Example 1

Let the vector field and be bounded by the following

The corresponding double contour integral would be set up as such:

OiintLaTeX.svg

We now evaluate . Meanwhile, set up the corresponding triple integral:

Example 2

Let the vector field , and remark that there are 4 parameters in this case. Let this vector field be bounded by the following:

To evaluate this, we must utilize the divergence theorem as stated before, and we must evaluate . Let

OiiintLaTeX.svg

Thus, we can evaluate a contour integral with . We can use the same method to evaluate contour integrals for any vector field with as well.

Integral representation

An integral representation of a function is an expression of the function involving a contour integral. Various integral representations are known for many special functions. Integral representations can be important for theoretical reasons, e.g. giving analytic continuation or functional equations, or sometimes for numerical evaluations.

Hankel's contour Hankel contour.png
Hankel's contour

For example, the original definition of the Riemann zeta function ζ(s) via a Dirichlet series,

is valid only for Re(s) > 1. But

where the integration is done over the Hankel contour H, is valid for all complex s not equal to 1.

See also

Related Research Articles

<span class="mw-page-title-main">Bessel function</span> Families of solutions to related differential equations

Bessel functions, first defined by the mathematician Daniel Bernoulli and then generalized by Friedrich Bessel, are canonical solutions y(x) of Bessel's differential equation for an arbitrary complex number , which represents the order of the Bessel function. Although and produce the same differential equation, it is conventional to define different Bessel functions for these two values in such a way that the Bessel functions are mostly smooth functions of .

<span class="mw-page-title-main">Cauchy's integral theorem</span> Theorem in complex analysis

In mathematics, the Cauchy integral theorem in complex analysis, named after Augustin-Louis Cauchy, is an important statement about line integrals for holomorphic functions in the complex plane. Essentially, it says that if is holomorphic in a simply connected domain Ω, then for any simply closed contour in Ω, that contour integral is zero.

<span class="mw-page-title-main">Cauchy's integral formula</span> Provides integral formulas for all derivatives of a holomorphic function

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.

<span class="mw-page-title-main">Residue (complex analysis)</span> Attribute of a mathematical function

In mathematics, more specifically complex analysis, the residue is a complex number proportional to the contour integral of a meromorphic function along a path enclosing one of its singularities. Residues can be computed quite easily and, once known, allow the determination of general contour integrals via the residue theorem.

<span class="mw-page-title-main">Residue theorem</span> Concept of complex analysis

In complex analysis, the residue theorem, sometimes called Cauchy's residue theorem, is a powerful tool to evaluate line integrals of analytic functions over closed curves; it can often be used to compute real integrals and infinite series as well. It generalizes the Cauchy integral theorem and Cauchy's integral formula. The residue theorem should not be confused with special cases of the generalized Stokes' theorem; however, the latter can be used as an ingredient of its proof.

<span class="mw-page-title-main">Stirling's approximation</span> Approximation for factorials

In mathematics, Stirling's approximation is an asymptotic approximation for factorials. It is a good approximation, leading to accurate results even for small values of . It is named after James Stirling, though a related but less precise result was first stated by Abraham de Moivre.

<span class="mw-page-title-main">Error function</span> Sigmoid shape special function

In mathematics, the error function, often denoted by erf, is a function defined as:

<span class="mw-page-title-main">Inverse trigonometric functions</span> Inverse functions of sin, cos, tan, etc.

In mathematics, the inverse trigonometric functions are the inverse functions of the trigonometric functions, under suitably restricted domains. Specifically, they are the inverses of the sine, cosine, tangent, cotangent, secant, and cosecant functions, and are used to obtain an angle from any of the angle's trigonometric ratios. Inverse trigonometric functions are widely used in engineering, navigation, physics, and geometry.

<span class="mw-page-title-main">Trigonometric substitution</span> Technique of integral evaluation

In mathematics, a trigonometric substitution replaces a trigonometric function for another expression. In calculus, trigonometric substitutions are a technique for evaluating integrals. In this case, an expression involving a radical function is replaced with a trigonometric one. Trigonometric identities may help simplify the answer. Like other methods of integration by substitution, when evaluating a definite integral, it may be simpler to completely deduce the antiderivative before applying the boundaries of integration.

The Basel problem is a problem in mathematical analysis with relevance to number theory, concerning an infinite sum of inverse squares. It was first posed by Pietro Mengoli in 1650 and solved by Leonhard Euler in 1734, and read on 5 December 1735 in The Saint Petersburg Academy of Sciences. Since the problem had withstood the attacks of the leading mathematicians of the day, Euler's solution brought him immediate fame when he was twenty-eight. Euler generalised the problem considerably, and his ideas were taken up more than a century later by Bernhard Riemann in his seminal 1859 paper "On the Number of Primes Less Than a Given Magnitude", in which he defined his zeta function and proved its basic properties. The problem is named after Basel, hometown of Euler as well as of the Bernoulli family who unsuccessfully attacked the problem.

<span class="mw-page-title-main">Gaussian integral</span> Integral of the Gaussian function, equal to sqrt(π)

The Gaussian integral, also known as the Euler–Poisson integral, is the integral of the Gaussian function over the entire real line. Named after the German mathematician Carl Friedrich Gauss, the integral is

In mathematics, Laplace's method, named after Pierre-Simon Laplace, is a technique used to approximate integrals of the form

<span class="mw-page-title-main">Argument principle</span> Theorem in complex analysis

In complex analysis, the argument principle is a theorem relating the difference between the number of zeros and poles of a meromorphic function to a contour integral of the function's logarithmic derivative.

<span class="mw-page-title-main">Voigt profile</span> Probability distribution

The Voigt profile is a probability distribution given by a convolution of a Cauchy-Lorentz distribution and a Gaussian distribution. It is often used in analyzing data from spectroscopy or diffraction.

<span class="mw-page-title-main">Dirichlet integral</span> Integral of sin(x)/x from 0 to infinity.

In mathematics, there are several integrals known as the Dirichlet integral, after the German mathematician Peter Gustav Lejeune Dirichlet, one of which is the improper integral of the sinc function over the positive real number line.

In probability theory, calculation of the sum of normally distributed random variables is an instance of the arithmetic of random variables.

In mathematics, the Nørlund–Rice integral, sometimes called Rice's method, relates the nth forward difference of a function to a line integral on the complex plane. It commonly appears in the theory of finite differences and has also been applied in computer science and graph theory to estimate binary tree lengths. It is named in honour of Niels Erik Nørlund and Stephen O. Rice. Nørlund's contribution was to define the integral; Rice's contribution was to demonstrate its utility by applying saddle-point techniques to its evaluation.

<span class="mw-page-title-main">Reciprocal gamma function</span> Mathematical function

In mathematics, the reciprocal gamma function is the function

In complex analysis, Jordan's lemma is a result frequently used in conjunction with the residue theorem to evaluate contour integrals and improper integrals. The lemma is named after the French mathematician Camille Jordan.

The Egorychev method is a collection of techniques introduced by Georgy Egorychev for finding identities among sums of binomial coefficients, Stirling numbers, Bernoulli numbers, Harmonic numbers, Catalan numbers and other combinatorial numbers. The method relies on two observations. First, many identities can be proved by extracting coefficients of generating functions. Second, many generating functions are convergent power series, and coefficient extraction can be done using the Cauchy residue theorem.

References

  1. Stalker, John (1998). Complex Analysis: Fundamentals of the Classical Theory of Functions. Springer. p. 77. ISBN   0-8176-4038-X.
  2. Bak, Joseph; Newman, Donald J. (1997). "Chapters 11 & 12". Complex Analysis. Springer. pp. 130–156. ISBN   0-387-94756-6.
  3. Krantz, Steven George (1999). "Chapter 2". Handbook of Complex Variables. Springer. ISBN   0-8176-4011-8.
  4. Mitrinović, Dragoslav S.; Kečkić, Jovan D. (1984). "Chapter 2". The Cauchy Method of Residues: Theory and Applications. Springer. ISBN   90-277-1623-4.
  5. Mitrinović, Dragoslav S.; Kečkić, Jovan D. (1984). "Chapter 5". The Cauchy Method of Residues: Theory and Applications. Springer. ISBN   90-277-1623-4.
  6. 1 2 3 4 5 Saff, Edward B.; Snider, Arthur David (2003). "Chapter 4". Fundamentals of Complex Analysis with Applications to Engineering, Science, and Mathematics (3rd ed.). Prentice Hall. ISBN   0-1390-7874-6.

Further reading