13.07.2015 Views

Mathematical Methods for Physics and Engineering - Matematica.NET

Mathematical Methods for Physics and Engineering - Matematica.NET

Mathematical Methods for Physics and Engineering - Matematica.NET

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

24.13 DEFINITE INTEGRALS USING CONTOUR INTEGRATIONWe note that result (24.60) is a special case of (24.63) in which θ 2 is equal toθ 1 +2π.24.13 Definite integrals using contour integrationThe remainder of this chapter is devoted to methods of applying contour integration<strong>and</strong> the residue theorem to various types of definite integral. However, threeapplications of contour integration, in which obtaining a value <strong>for</strong> the integral isnot the prime purpose of the exercise, have been postponed until chapter 25. Theyare the location of the zeros of a complex polynomial, the evaluation of the sumsof certain infinite series <strong>and</strong> the determination of inverse Laplace trans<strong>for</strong>ms.For the integral evalations considered here, not much preamble is given since,<strong>for</strong> this material, the simplest explanation is felt to be via a series of workedexamples that can be used as models.Suppose that an integral of the <strong>for</strong>m24.13.1 Integrals of sinusoidal functions∫ 2π0F(cos θ, sin θ) dθ (24.64)is to be evaluated. It can be made into a contour integral around the unit circleC by writing z =expiθ, <strong>and</strong> hencecos θ = 1 2 (z + z−1 ), sin θ = − 1 2 i(z − z−1 ), dθ = −iz −1 dz. (24.65)This contour integral can then be evaluated using the residue theorem, providedthe trans<strong>for</strong>med integr<strong>and</strong> has only a finite number of poles inside the unit circle<strong>and</strong>noneonit.◮Evaluate∫ 2πcos 2θI =dθ, b > a > 0. (24.66)0 a 2 + b 2 − 2ab cos θBy de Moivre’s theorem (section 3.4),cos nθ = 1 2 (zn + z −n ). (24.67)Using n = 2 in (24.67) <strong>and</strong> straight<strong>for</strong>ward substitution <strong>for</strong> the other functions of θ in(24.66) givesI =i ∮z 4 +12ab C z 2 (z − a/b)(z − b/a) dz.Thus there are two poles inside C, a double pole at z = 0 <strong>and</strong> a simple pole at z = a/b(recall that b>a).We could find the residue of the integr<strong>and</strong> at z = 0 by exp<strong>and</strong>ing the integr<strong>and</strong> as aLaurent series in z <strong>and</strong> identifying the coefficient of z −1 . Alternatively, we may use the861

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!