9.1 Evaluation of Real Integrals Contour integration is a powerful tool for evaluating definite integrals.
9.2 Integrals of Rational Functions over the Real Line Theorem 9.1. If f ( x ) = P ( x ) / Q ( x ) f(x) = P(x)/Q(x) f ( x ) = P ( x ) / Q ( x ) where deg ( Q ) ≥ deg ( P ) + 2 \deg(Q) \geq \deg(P) + 2 deg ( Q ) ≥ deg ( P ) + 2 and Q Q Q has no real roots, Then
∫ − ∞ ∞ f ( x ) d x = 2 π i ∑ I m ( z k ) > 0 R e s ( f , z k ) \int_{-\infty}^{\infty} f(x)\, dx = 2\pi i \sum_{\mathrm{Im}(z_k) > 0} \mathrm{Res}(f, z_k) ∫ − ∞ ∞ f ( x ) d x = 2 π i ∑ Im ( z k ) > 0 Res ( f , z k )
Where the sum is over poles in the upper half-plane.
Proof. Integrate f ( z ) f(z) f ( z ) over the semicircular contour γ R \gamma_R γ R consisting of [ − R , R ] [-R, R] [ − R , R ] on the Real axis and the semicircle ∣ z ∣ = R |z| = R ∣ z ∣ = R in the upper half-plane. As R → ∞ R \to \infty R → ∞ The integral over The semicircle vanishes (since ∣ f ( z ) ∣ ≤ M / R 2 |f(z)| \leq M/R^2 ∣ f ( z ) ∣ ≤ M / R 2 and the length is π R \pi R π R ). ■ \blacksquare ■
9.3 Worked Example Problem. Evaluate ∫ − ∞ ∞ d x x 2 + 1 \int_{-\infty}^{\infty} \frac{dx}{x^2 + 1} ∫ − ∞ ∞ x 2 + 1 d x .
Solution. f ( z ) = 1 z 2 + 1 f(z) = \frac{1}{z^2 + 1} f ( z ) = z 2 + 1 1 has simple poles at z = ± i z = \pm i z = ± i .
Only z = i z = i z = i is in the upper half-plane.
R e s ( 1 z 2 + 1 , i ) = 1 2 z ∣ z = i = 1 2 i \mathrm{Res}\left(\frac{1}{z^2 + 1}, i\right) = \frac{1}{2z}\Big|_{z = i} = \frac{1}{2i} Res ( z 2 + 1 1 , i ) = 2 z 1 z = i = 2 i 1 .
∫ − ∞ ∞ d x x 2 + 1 = 2 π i ⋅ 1 2 i = π \int_{-\infty}^{\infty} \frac{dx}{x^2 + 1} = 2\pi i \cdot \frac{1}{2i} = \pi ∫ − ∞ ∞ x 2 + 1 d x = 2 π i ⋅ 2 i 1 = π . ■ \blacksquare ■
9.4 Integrals Involving Trigonometric Functions For integrals of the form ∫ 0 2 π R ( cos θ , sin θ ) d θ \int_0^{2\pi} R(\cos\theta, \sin\theta)\, d\theta ∫ 0 2 π R ( cos θ , sin θ ) d θ Substitute z = e i θ z = e^{i\theta} z = e i θ So d z = i z d θ dz = iz\, d\theta d z = i z d θ , cos θ = z + z − 1 2 \cos\theta = \frac{z + z^{-1}}{2} cos θ = 2 z + z − 1 sin θ = z − z − 1 2 i \sin\theta = \frac{z - z^{-1}}{2i} sin θ = 2 i z − z − 1 .
The integral becomes ∫ ∣ z ∣ = 1 f ( z ) d z \int_{|z|=1} f(z)\, dz ∫ ∣ z ∣ = 1 f ( z ) d z where f ( z ) f(z) f ( z ) is a rational function.
9.5 Worked Example Problem. Evaluate ∫ 0 2 π d θ 2 + cos θ \int_0^{2\pi} \frac{d\theta}{2 + \cos\theta} ∫ 0 2 π 2 + c o s θ d θ .
Solution. Substitute z = e i θ z = e^{i\theta} z = e i θ : d θ = d z i z d\theta = \frac{dz}{iz} d θ = i z d z cos θ = z + 1 / z 2 \cos\theta = \frac{z + 1/z}{2} cos θ = 2 z + 1/ z .
∫ ∣ z ∣ = 1 d z i z ( 2 + z + 1 / z 2 ) = ∫ ∣ z ∣ = 1 2 d z i ( z 2 + 4 z + 1 ) \int_{|z|=1} \frac{dz}{iz\left(2 + \frac{z + 1/z}{2}\right)} = \int_{|z|=1} \frac{2\, dz}{i(z^2 + 4z + 1)} ∫ ∣ z ∣ = 1 i z ( 2 + 2 z + 1/ z ) d z = ∫ ∣ z ∣ = 1 i ( z 2 + 4 z + 1 ) 2 d z
Poles: z 2 + 4 z + 1 = 0 ⇒ z = − 2 ± 3 z^2 + 4z + 1 = 0 \Rightarrow z = -2 \pm \sqrt{3} z 2 + 4 z + 1 = 0 ⇒ z = − 2 ± 3 .
∣ z 1 ∣ = ∣ − 2 + 3 ∣ = 2 − 3 < 1 |z_1| = |-2 + \sqrt{3}| = 2 - \sqrt{3} \lt 1 ∣ z 1 ∣ = ∣ − 2 + 3 ∣ = 2 − 3 < 1 (inside). ∣ z 2 ∣ = ∣ − 2 − 3 ∣ = 2 + 3 > 1 |z_2| = |-2 - \sqrt{3}| = 2 + \sqrt{3} \gt 1 ∣ z 2 ∣ = ∣ − 2 − 3 ∣ = 2 + 3 > 1 (outside).
R e s ( 1 z 2 + 4 z + 1 , z 1 ) = 1 2 3 \mathrm{Res}\left(\frac{1}{z^2 + 4z + 1}, z_1\right) = \frac{1}{2\sqrt{3}} Res ( z 2 + 4 z + 1 1 , z 1 ) = 2 3 1 .
∫ 0 2 π d θ 2 + cos θ = 2 i ⋅ 2 π i ⋅ 1 2 3 = 2 π 3 \int_0^{2\pi} \frac{d\theta}{2 + \cos\theta} = \frac{2}{i} \cdot 2\pi i \cdot \frac{1}{2\sqrt{3}} = \frac{2\pi}{\sqrt{3}} ∫ 0 2 π 2 + c o s θ d θ = i 2 ⋅ 2 π i ⋅ 2 3 1 = 3 2 π . ■ \blacksquare ■
9.6 Jordan”s Lemma Theorem 9.2 (Jordan’s Lemma). If f ( z ) → 0 f(z) \to 0 f ( z ) → 0 uniformly as ∣ z ∣ → ∞ |z| \to \infty ∣ z ∣ → ∞ in the upper Half-plane and a > 0 a \gt 0 a > 0 Then
lim R → ∞ ∫ C R e i a z f ( z ) d z = 0 \lim_{R \to \infty} \int_{C_R} e^{iaz}f(z)\, dz = 0 lim R → ∞ ∫ C R e ia z f ( z ) d z = 0
Where C R C_R C R is the upper semicircle ∣ z ∣ = R |z| = R ∣ z ∣ = R , I m ( z ) ≥ 0 \mathrm{Im}(z) \geq 0 Im ( z ) ≥ 0 .
This allows evaluation of integrals of the form ∫ − ∞ ∞ f ( x ) cos ( a x ) d x \int_{-\infty}^{\infty} f(x)\cos(ax)\, dx ∫ − ∞ ∞ f ( x ) cos ( a x ) d x and ∫ − ∞ ∞ f ( x ) sin ( a x ) d x \int_{-\infty}^{\infty} f(x)\sin(ax)\, dx ∫ − ∞ ∞ f ( x ) sin ( a x ) d x .
9.7 Fourier-Type Integrals Solution Problem. Evaluate ∫ − ∞ ∞ cos x x 2 + 1 d x \int_{-\infty}^{\infty} \frac{\cos x}{x^2 + 1}\, dx ∫ − ∞ ∞ x 2 + 1 c o s x d x .
Consider ∫ − ∞ ∞ e i x x 2 + 1 d x = 2 π i ⋅ R e s ( e i z z 2 + 1 , i ) \int_{-\infty}^{\infty} \frac{e^{ix}}{x^2 + 1}\, dx = 2\pi i \cdot \mathrm{Res}\!\left(\frac{e^{iz}}{z^2+1}, i\right) ∫ − ∞ ∞ x 2 + 1 e i x d x = 2 π i ⋅ Res ( z 2 + 1 e i z , i ) .
R e s ( e i z z 2 + 1 , i ) = e i ⋅ i 2 i = e − 1 2 i \mathrm{Res}\!\left(\frac{e^{iz}}{z^2+1}, i\right) = \frac{e^{i \cdot i}}{2i} = \frac{e^{-1}}{2i} Res ( z 2 + 1 e i z , i ) = 2 i e i ⋅ i = 2 i e − 1 .
∫ − ∞ ∞ e i x x 2 + 1 d x = 2 π i ⋅ e − 1 2 i = π e \int_{-\infty}^{\infty} \frac{e^{ix}}{x^2 + 1}\, dx = 2\pi i \cdot \frac{e^{-1}}{2i} = \frac{\pi}{e} ∫ − ∞ ∞ x 2 + 1 e i x d x = 2 π i ⋅ 2 i e − 1 = e π .
Taking real parts: ∫ − ∞ ∞ cos x x 2 + 1 d x = π e \int_{-\infty}^{\infty} \frac{\cos x}{x^2 + 1}\, dx = \frac{\pi}{e} ∫ − ∞ ∞ x 2 + 1 c o s x d x = e π .
Problem. Evaluate ∫ − ∞ ∞ x sin x x 2 + a 2 d x \int_{-\infty}^{\infty} \frac{x \sin x}{x^2 + a^2}\, dx ∫ − ∞ ∞ x 2 + a 2 x s i n x d x for a > 0 a \gt 0 a > 0 .
Consider ∫ − ∞ ∞ z e i z z 2 + a 2 d z \int_{-\infty}^{\infty} \frac{z\, e^{iz}}{z^2 + a^2}\, dz ∫ − ∞ ∞ z 2 + a 2 z e i z d z . Only z = i a z = ia z = ia is in the upper half-plane.
R e s ( z e i z z 2 + a 2 , i a ) = i a ⋅ e i ⋅ i a 2 i a = e − a 2 \mathrm{Res}\!\left(\frac{ze^{iz}}{z^2 + a^2}, ia\right) = \frac{ia \cdot e^{i \cdot ia}}{2ia} = \frac{e^{-a}}{2} Res ( z 2 + a 2 z e i z , ia ) = 2 ia ia ⋅ e i ⋅ ia = 2 e − a .
∫ − ∞ ∞ x e i x x 2 + a 2 d x = 2 π i ⋅ e − a 2 = π i e − a \int_{-\infty}^{\infty} \frac{x\, e^{ix}}{x^2 + a^2}\, dx = 2\pi i \cdot \frac{e^{-a}}{2} = \pi i\, e^{-a} ∫ − ∞ ∞ x 2 + a 2 x e i x d x = 2 π i ⋅ 2 e − a = π i e − a .
Taking imaginary parts: ∫ − ∞ ∞ x sin x x 2 + a 2 d x = π e − a \int_{-\infty}^{\infty} \frac{x \sin x}{x^2 + a^2}\, dx = \pi\, e^{-a} ∫ − ∞ ∞ x 2 + a 2 x s i n x d x = π e − a .
Problem. Evaluate ∫ 0 2 π cos 2 θ 5 + 4 cos θ d θ \int_0^{2\pi} \frac{\cos 2\theta}{5 + 4\cos\theta}\, d\theta ∫ 0 2 π 5 + 4 c o s θ c o s 2 θ d θ .
Substitute z = e i θ z = e^{i\theta} z = e i θ : cos θ = ( z + z − 1 ) / 2 \cos\theta = (z + z^{-1})/2 cos θ = ( z + z − 1 ) /2 , cos 2 θ = ( z 2 + z − 2 ) / 2 \cos 2\theta = (z^2 + z^{-2})/2 cos 2 θ = ( z 2 + z − 2 ) /2 .
I = 1 2 i ∫ ∣ z ∣ = 1 z 4 + 1 z 2 ( 2 z + 1 ) ( z + 2 ) d z I = \frac{1}{2i}\int_{|z|=1} \frac{z^4 + 1}{z^2(2z + 1)(z + 2)}\, dz I = 2 i 1 ∫ ∣ z ∣ = 1 z 2 ( 2 z + 1 ) ( z + 2 ) z 4 + 1 d z .
Poles inside ∣ z ∣ = 1 |z| = 1 ∣ z ∣ = 1 : z = 0 z = 0 z = 0 (order 2 2 2 ) and z = − 1 / 2 z = -1/2 z = − 1/2 (simple).
At z = 0 z = 0 z = 0 : R e s = d d z [ z 4 + 1 ( 2 z + 1 ) ( z + 2 ) ] z = 0 = − 5 4 \mathrm{Res} = \frac{d}{dz}\left[\frac{z^4 + 1}{(2z+1)(z+2)}\right]_{z=0} = -\frac{5}{4} Res = d z d [ ( 2 z + 1 ) ( z + 2 ) z 4 + 1 ] z = 0 = − 4 5 .
At z = − 1 / 2 z = -1/2 z = − 1/2 : R e s = 17 / 16 3 / 4 = 17 12 \mathrm{Res} = \frac{17/16}{3/4} = \frac{17}{12} Res = 3/4 17/16 = 12 17 .
I = 1 2 i ⋅ 2 π i ( − 5 4 + 17 12 ) = π 6 I = \frac{1}{2i} \cdot 2\pi i \left(-\frac{5}{4} + \frac{17}{12}\right) = \frac{\pi}{6} I = 2 i 1 ⋅ 2 π i ( − 4 5 + 12 17 ) = 6 π .
9.8 Improper Integrals and Principal Value For integrals where the integrand has poles on the real axis, we use the Cauchy principal value :
P V ∫ − ∞ ∞ f ( x ) d x = lim ε → 0 + ( ∫ − ∞ a − ε f ( x ) d x + ∫ a + ε ∞ f ( x ) d x ) \mathrm{PV}\!\int_{-\infty}^{\infty} f(x)\, dx = \lim_{\varepsilon \to 0^+} \left(\int_{-\infty}^{a-\varepsilon} f(x)\, dx + \int_{a+\varepsilon}^{\infty} f(x)\, dx\right) PV ∫ − ∞ ∞ f ( x ) d x = lim ε → 0 + ( ∫ − ∞ a − ε f ( x ) d x + ∫ a + ε ∞ f ( x ) d x )
Solution Problem. Evaluate P V ∫ − ∞ ∞ sin x x d x \mathrm{PV}\!\int_{-\infty}^{\infty} \frac{\sin x}{x}\, dx PV ∫ − ∞ ∞ x s i n x d x .
Consider ∮ γ e i z z d z \oint_\gamma \frac{e^{iz}}{z}\, dz ∮ γ z e i z d z where γ \gamma γ consists of [ − R , − ε ] [-R, -\varepsilon] [ − R , − ε ] [ ε , R ] [\varepsilon, R] [ ε , R ] on the real axis, small upper semicircle C ε C_\varepsilon C ε around 0 0 0 And large Upper semicircle C R C_R C R .
No poles inside the contour, so the integral is 0 0 0 .
On C R C_R C R : vanishes as R → ∞ R \to \infty R → ∞ by Jordan’s lemma. On C ε C_\varepsilon C ε (indenting above): ∫ C ε e i z z d z → − i π \int_{C_\varepsilon} \frac{e^{iz}}{z}\, dz \to -i\pi ∫ C ε z e i z d z → − iπ as ε → 0 \varepsilon \to 0 ε → 0 (half residue contribution).
0 = P V ∫ − ∞ ∞ e i x x d x + ( − i π ) 0 = \mathrm{PV}\!\int_{-\infty}^{\infty} \frac{e^{ix}}{x}\, dx + (-i\pi) 0 = PV ∫ − ∞ ∞ x e i x d x + ( − iπ ) .
P V ∫ − ∞ ∞ e i x x d x = i π \mathrm{PV}\!\int_{-\infty}^{\infty} \frac{e^{ix}}{x}\, dx = i\pi PV ∫ − ∞ ∞ x e i x d x = iπ .
Taking imaginary parts: P V ∫ − ∞ ∞ sin x x d x = π \mathrm{PV}\!\int_{-\infty}^{\infty} \frac{\sin x}{x}\, dx = \pi PV ∫ − ∞ ∞ x s i n x d x = π .