书架/Thomas' Calculus

Chapter 8: Techniques of Integration

教材插图

OVERVIEW The Fundamental Theorem tells us how to evaluate a definite integral once we have an antiderivative for the integrand function. However, finding antiderivatives (or indefinite integrals) is not as straightforward as finding derivatives. In this chapter we study a number of important techniques that apply to finding integrals for specialized classes of functions such as trigonometric functions, products of certain functions, and rational functions. Since we cannot always find an antiderivative, we develop numerical methods for calculating definite integrals. We also study integrals for which the domain or range is infinite, called improper integrals.

8.1 Using Basic Integration Formulas

Table 8.1 summarizes the indefinite integrals of many of the functions we have studied so far, and the substitution method helps us use the table to evaluate more complicated functions involving these basic ones. In this section we combine the Substitution Rules (studied in Chapter 5) with algebraic methods and trigonometric identities to help us use Table 8.1. A more extensive Table of Integrals is given at the back of the chapter, and we discuss its use in Section 8.6.

Sometimes we have to rewrite an integral to match it to a standard form of the type displayed in Table 8.1. We start with an example of this procedure.

EXAMPLE 1 Evaluate the integral

∫532𝑥−3√𝑥2−3𝑥+1𝑑𝑥.

Solution We rewrite the integral and apply the Substitution Rule for Definite Integrals presented in Section 5.6, to find

∫532𝑥−3√𝑥2−3𝑥+1𝑑𝑥=∫111𝑑𝑢√𝑢𝑢=𝑥2−3𝑥+1,𝑑𝑢=(2𝑥−3)𝑑𝑥;=∫111𝑢−1/2𝑑𝑢=2√𝑢]111=2(√11−1)≈4.63.(Table8.1,Formula2)

TABLE 8.1 Basic integration formulas

  1. ∫𝑘𝑑𝑥 =𝑘𝑥 +𝐶 (any number 𝑘 )

  2. ∫𝑥𝑛𝑑𝑥 =𝑥𝑛+1𝑛+1 +𝐶(𝑛 ≠ −1)

  3. ∫𝑑𝑥𝑥 =ln⁡|𝑥| +𝐶

  4. ∫𝑒𝑥𝑑𝑥 =𝑒𝑥 +𝐶

  5. ∫𝑎𝑥𝑑𝑥 =𝑎𝑥ln⁡𝑎 +𝐶(𝑎 >0,𝑎 ≠1)

  6. ∫sin⁡𝑥𝑑𝑥 = −cos⁡𝑥 +𝐶

  7. ∫cos⁡𝑥𝑑𝑥 =sin⁡𝑥 +𝐶

  8. ∫sec2⁡𝑥𝑑𝑥 =tan⁡𝑥 +𝐶

  9. ∫csc2⁡𝑥𝑑𝑥 = −cot⁡𝑥 +𝐶

  10. ∫sec⁡𝑥tan⁡𝑥𝑑𝑥 =sec⁡𝑥 +𝐶

  11. ∫csc⁡𝑥cot⁡𝑥𝑑𝑥 = −csc⁡𝑥 +𝐶

  12. ∫tan⁡𝑥𝑑𝑥 =ln⁡|sec⁡𝑥| +𝐶

  13. ∫cot⁡𝑥𝑑𝑥 =ln⁡|sin⁡𝑥| +𝐶

  14. ∫sec⁡𝑥𝑑𝑥 =ln⁡|sec⁡𝑥 +tan⁡𝑥| +𝐶

  15. ∫csc⁡𝑥𝑑𝑥 = −ln⁡|csc⁡𝑥 +cot⁡𝑥| +𝐶

  16. ∫sinh⁡𝑥𝑑𝑥 =cosh⁡𝑥 +𝐶

  17. ∫cosh⁡𝑥𝑑𝑥 =sinh⁡𝑥 +𝐶

  18. ∫𝑑𝑥√𝑎2−𝑥2 =arcsin⁡(𝑥𝑎) +𝐶

  19. ∫𝑑𝑥𝑎2+𝑥2 =1𝑎arctan⁡(𝑥𝑎) +𝐶

  20. ∫𝑑𝑥𝑥√𝑥2−𝑎2 =1𝑎arcsec∣𝑥𝑎∣ +𝐶

  21. ∫𝑑𝑥√𝑎2+𝑥2 =sinh−1⁡(𝑥𝑎) +𝐶(𝑎 >0)

  22. ∫𝑑𝑥√𝑥2−𝑎2 =cosh−1⁡(𝑥𝑎) +𝐶(𝑥 >𝑎 >0)

EXAMPLE 2 Complete the square to evaluate

∫𝑑𝑥√8𝑥−𝑥2.

Solution We complete the square to simplify the denominator:

8𝑥−𝑥2=−(𝑥2−8𝑥)=−(𝑥2−8𝑥+16−16)=−(𝑥2−8𝑥+16)+16=16−(𝑥−4)2.

Then

∫𝑑𝑥√8𝑥−𝑥2=∫𝑑𝑥√16−(𝑥−4)2=∫𝑑𝑢√𝑎2−𝑢2𝑎=4,𝑢=(𝑥−4),=arcsin⁡(𝑢𝑎)+𝐶𝑑𝑢=𝑑𝑥=arcsin⁡(𝑥−44)+𝐶.(Table8.1,Formula18)

EXAMPLE 3 Evaluate the integral

∫(cos⁡𝑥sin⁡2𝑥+sin⁡𝑥cos⁡2𝑥)𝑑𝑥.

Solution We can replace the integrand with an equivalent trigonometric expression using the Sine Addition Formula to obtain a simple substitution:

∫(cos⁡𝑥sin⁡2𝑥+sin⁡𝑥cos⁡2𝑥)𝑑𝑥=∫sin⁡(𝑥+2𝑥)𝑑𝑥=∫sin⁡3𝑥𝑑𝑥=∫13sin⁡𝑢𝑑𝑢𝑢=3𝑥,𝑑𝑢=3𝑑𝑥=−13cos⁡3𝑥+𝐶. Table 8.1, Formula 6 

In Section 5.5 we found the indefinite integral of the secant function by multiplying it by a fractional form equal to one, and then integrating the equivalent result. We can use that same procedure in other instances as well, as we illustrate next.

EXAMPLE 4 Find ∫𝜋/40𝑑𝑥1−sin⁡𝑥.

Solution We multiply the numerator and denominator of the integrand by 1 +sin⁡𝑥 . This procedure transforms the integral into one we can evaluate:

∫𝜋/40𝑑𝑥1−sin⁡𝑥=∫𝜋/4011−sin⁡𝑥⋅1+sin⁡𝑥1+sin⁡𝑥𝑑𝑥 Multiply and divide by  conjugate. =∫𝜋/401+sin⁡𝑥1−sin2⁡𝑥𝑑𝑥 Simplify. =∫𝜋/401+sin⁡𝑥cos2⁡𝑥𝑑𝑥1−sin2⁡𝑥=cos2⁡𝑥=∫𝜋/40(sec2⁡𝑥+sec⁡𝑥tan⁡𝑥)𝑑𝑥 Use Table 8.1,  Formulas 8 and 10 =[tan⁡𝑥+sec⁡𝑥]𝜋/40=(1+√2−(0+1))=√2.

EXAMPLE 5 Evaluate

∫3𝑥2−7𝑥3𝑥+2𝑑𝑥. 𝑥−33𝑥+2――――――)3𝑥2−7𝑥3𝑥2+2𝑥−9𝑥−9𝑥−6+6

Solution The integrand is an improper fraction since the degree of the numerator is greater than the degree of the denominator. To integrate it, we perform long division to obtain a quotient plus a remainder that is a proper fraction:

3𝑥2−7𝑥3𝑥+2=𝑥−3+63𝑥+2.

Therefore,

∫3𝑥2−7𝑥3𝑥+2𝑑𝑥=∫(𝑥−3+63𝑥+2)𝑑𝑥=𝑥22−3𝑥+2ln⁡|3𝑥+2|+𝐶.

Reducing an improper fraction by long division (Example 5) does not always lead to an expression we can integrate directly. We see what to do about that in Section 8.5.

EXAMPLE 6 Evaluate

∫3𝑥+2√1−𝑥2𝑑𝑥.

Solution We first separate the integrand to get

∫3𝑥+2√1−𝑥2𝑑𝑥=3∫𝑥𝑑𝑥√1−𝑥2+2∫𝑑𝑥√1−𝑥2.

In the first of these new integrals, we substitute

𝑢=1−𝑥2,𝑑𝑢=−2𝑥𝑑𝑥, so 𝑥𝑑𝑥=−12𝑑𝑢.

Then we obtain

3∫𝑥𝑑𝑥√1−𝑥2=3∫(−1/2)𝑑𝑢√𝑢=−32∫𝑢−1/2𝑑𝑢=−32⋅𝑢1/21/2+𝐶1=−3√1−𝑥2+𝐶1.

The second of the new integrals is a standard form,

2∫𝑑𝑥√1−𝑥2=2arcsin⁡𝑥+𝐶2. Table 8.1, Formula 18 

Combining these results and renaming 𝐶1 +𝐶2 as C gives

∫3𝑥+2√1−𝑥2𝑑𝑥=−3√1−𝑥2+2arcsin⁡𝑥+𝐶.

The question of what to substitute for in an integrand is not always quite so clear. Sometimes we simply proceed by trial-and-error, and if nothing works out, we then try another method altogether. The next several sections of the text present some of these new methods, but substitution works in the following example.

EXAMPLE 7 Evaluate

∫𝑑𝑥(1+√𝑥)3.

Solution We might try substituting for the term √𝑥 , but the derivative factor 1/√𝑥 is missing from the integrand, so this substitution will not help. The other possibility is to substitute for (1 +√𝑥) , and it turns out this works:

∫𝑑𝑥(1+√𝑥)3=∫2(𝑢−1)𝑑𝑢𝑢3𝑢=1+√𝑥,𝑑𝑢=12√𝑥𝑑𝑥;𝑑𝑥=2√𝑥𝑑𝑢=2(𝑢−1)𝑑𝑢=∫(2𝑢2−2𝑢3)𝑑𝑢 =−2𝑢+1𝑢2+𝐶 =1−2𝑢𝑢2+𝐶 =1−2(1+√𝑥)(1+√𝑥)2+𝐶 =𝐶−1+2√𝑥(1+√𝑥)2.

When evaluating definite integrals, a property of the integrand may help us in calculating the result.

EXAMPLE 8 Evaluate

∫𝜋/2−𝜋/2𝑥3cos⁡𝑥𝑑𝑥.

Solution No substitution or algebraic manipulation is clearly helpful here. But we observe that the interval of integration is the symmetric interval [ −𝜋/2,𝜋/2] . Moreover, the factor 𝑥3 is an odd function, and cos⁡𝑥 is an even function, so their product is odd. Therefore,

∫𝜋/2−𝜋/2𝑥3cos⁡𝑥𝑑𝑥 =0. Theorem 8, Section 5.6

EXERCISES 8.1

Assorted Integrations

The integrals in Exercises 1–44 are in no particular order. Evaluate each integral using any algebraic method, trigonometric identity, or substitution you think is appropriate.

  1. ∫1016𝑥8𝑥2+2𝑑𝑥

  2. ∫𝑥2𝑥2+1𝑑𝑥

  3. ∫(sec⁡𝑥 −tan⁡𝑥)2𝑑𝑥

  4. ∫𝜋/3𝜋/4𝑑𝑥cos2⁡𝑥tan⁡𝑥

  5. ∫1−𝑥√1−𝑥2𝑑𝑥

  6. ∫𝑑𝑥𝑥−√𝑥

  7. ∫𝑒−cot⁡𝑧sin2⁡𝑧𝑑𝑧

  8. ∫2ln⁡𝑧316𝑧𝑑𝑧

  9. ∫𝑑𝑧𝑒𝑧+𝑒−𝑧

  10. ∫218𝑑𝑥𝑥2−2𝑥+2

  11. ∫0−14𝑑𝑥1+(2𝑥+1)2

  12. ∫3−14𝑥2−72𝑥+3𝑑𝑥

  13. ∫𝑑𝑡1−sec⁡𝑡

  14. ∫csc⁡𝑡sin⁡3𝑡𝑑𝑡

  15. ∫𝜋/401+sin⁡𝜃cos2⁡𝜃𝑑𝜃

  16. ∫𝑑𝜃√2𝜃−𝜃2

  17. ∫ln⁡𝑦𝑦+4𝑦ln2⁡𝑦𝑑𝑦

  18. ∫2√𝑦𝑑𝑦2√𝑦

  19. ∫𝑑𝜃sec⁡𝜃+tan⁡𝜃

  20. ∫𝑑𝑡𝑡√3+𝑡2

  21. ∫4𝑡3−𝑡2+16𝑡𝑡2+4𝑑𝑡

  22. ∫𝑥+2√𝑥−12𝑥√𝑥−1𝑑𝑥

  23. ∫𝜋/20√1−cos⁡𝜃𝑑𝜃

  24. ∫(sec⁡𝑡 +cot⁡𝑡)2𝑑𝑡

  25. ∫𝑑𝑦√𝑒2𝑦−1

  26. ∫6𝑑𝑦√𝑦(1+𝑦)

  27. ∫2𝑑𝑥𝑥√1−4ln2⁡𝑥

  28. ∫𝑑𝑥(𝑥−2)√𝑥2−4𝑥+3

  29. ∫(csc⁡𝑥 −sec⁡𝑥)(sin⁡𝑥 +cos⁡𝑥)𝑑𝑥

  30. ∫3sinh⁡(𝑥2+ln⁡5)𝑑𝑥

  31. ∫3√22𝑥3𝑥2−1𝑑𝑥

  32. ∫1−1√1+𝑥2sin⁡𝑥𝑑𝑥

  33. ∫0−1√1+𝑦1−𝑦𝑑𝑦

  34. ∫𝑒𝑧+𝑒𝑧𝑑𝑧

  35. ∫7𝑑𝑥(𝑥−1)√𝑥2−2𝑥−48

  36. ∫𝑑𝑥(2𝑥+1)√4𝑥+4𝑥2

  37. ∫2𝜃3−7𝜃2+7𝜃2𝜃−5𝑑𝜃

  38. ∫𝑑𝜃cos⁡𝜃−1

  39. ∫𝑑𝑥1+𝑒𝑥

  40. ∫√𝑥1+𝑥3𝑑𝑥

Hint: Use long division.

 Hint: Let 𝑢=𝑥3/2.
  1. ∫𝑒3𝑥𝑒𝑥+1𝑑𝑥

  2. ∫2𝑥−13𝑥𝑑𝑥

  3. ∫1√𝑥(1+𝑥)𝑑𝑥

  4. ∫tan⁡𝜃+3sin⁡𝜃𝑑𝜃

Theory and Examples

  1. Area Find the area of the region bounded above by 𝑦 =2cos⁡𝑥 and below by 𝑦 =sec⁡𝑥, −𝜋/4 ≤𝑥 ≤𝜋/4 .

  2. Volume Find the volume of the solid generated by revolving the region in Exercise 45 about the x-axis.

  3. Arc length Find the length of the curve 𝑦 =ln⁡(cos⁡𝑥) , 0 ≤𝑥 ≤𝜋/3 .

  4. Arc length Find the length of the curve 𝑦 =ln⁡(sec⁡𝑥) , 0 ≤𝑥 ≤𝜋/4 .

  5. Centroid Find the centroid of the region bounded by the 𝑥 -axis, the curve 𝑦 =sec⁡𝑥 , and the lines 𝑥 = −𝜋/4 , 𝑥 =𝜋/4 .

  6. Centroid Find the centroid of the region bounded by the 𝑥 -axis, the curve 𝑦 =csc⁡𝑥 , and the lines 𝑥 =𝜋/6 , 𝑥 =5𝜋/6 .

  7. The functions 𝑦 =𝑒𝑥3 and 𝑦 =𝑥3𝑒𝑥3 do not have elementary antiderivatives, but 𝑦 =(1 +3𝑥3)𝑒𝑥3 does. Evaluate

∫(1+3𝑥3)𝑒𝑥3𝑑𝑥.52.$𝑈𝑠𝑒𝑡ℎ𝑒𝑠𝑢𝑏𝑠𝑡𝑖𝑡𝑢𝑡𝑖𝑜𝑛$𝑢=tan⁡𝑥$𝑡𝑜𝑒𝑣𝑎𝑙𝑢𝑎𝑡𝑒𝑡ℎ𝑒𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙$∫𝑑𝑥1+sin2⁡𝑥.53.$𝑈𝑠𝑒𝑡ℎ𝑒𝑠𝑢𝑏𝑠𝑡𝑖𝑡𝑢𝑡𝑖𝑜𝑛$𝑢=𝑥4+1$𝑡𝑜𝑒𝑣𝑎𝑙𝑢𝑎𝑡𝑒𝑡ℎ𝑒𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙$∫𝑥7√𝑥4+1𝑑𝑥.54.$𝑈𝑠𝑖𝑛𝑔𝑑𝑖𝑓𝑓𝑒𝑟𝑒𝑛𝑡𝑠𝑢𝑏𝑠𝑡𝑖𝑡𝑢𝑡𝑖𝑜𝑛𝑠𝑆ℎ𝑜𝑤𝑡ℎ𝑎𝑡𝑡ℎ𝑒𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙$∫((𝑥2−1)(𝑥+1))−2/3𝑑𝑥

can be evaluated with any of the following substitutions.

a. 𝑢 =1/(𝑥 +1)

b. 𝑢 =((𝑥−1)/(𝑥+1))𝑘 for k = 1, 1/2, 1/3, -1/3, -2/3, and -1

c. 𝑢 =arctan⁡𝑥

d. 𝑢 =tan−1⁡√𝑥

e. 𝑢 =tan−1⁡((𝑥−1)/2)

f. 𝑢 =arccos⁡𝑥

g. 𝑢 =cosh−1⁡𝑥

8.2 Integration by Parts

What is the value of the integral?

Integration by parts is a technique for simplifying integrals of the form

∫𝑢(𝑥)𝑣′(𝑥)𝑑𝑥.

It is useful when 𝑢 can be differentiated repeatedly and 𝑣′ can be integrated repeatedly without difficulty. The integrals

∫𝑥cos⁡𝑥𝑑𝑥 and ∫𝑥2𝑒𝑥𝑑𝑥

are such integrals because 𝑢(𝑥) =𝑥 or 𝑢(𝑥) =𝑥2 can be differentiated repeatedly, and 𝑣′(𝑥) =cos⁡𝑥 or 𝑣′(𝑥) =𝑒𝑥 can be integrated repeatedly without difficulty. Integration by parts also applies to integrals like

∫ln⁡𝑥𝑑𝑥 and ∫𝑒𝑥cos⁡𝑥𝑑𝑥.

In the first case, the integrand ln⁡𝑥 can be rewritten as (ln⁡𝑥)(1) , and 𝑢(𝑥) =ln⁡𝑥 is easy to differentiate while 𝑣′(𝑥) =1 easily integrates to x. In the second case, each part of the integrand appears again after repeated differentiation or integration.

Product Rule in Integral Form

If 𝑢 and 𝑣 are differentiable functions of 𝑥 , the Product Rule says that

𝑑𝑑𝑥[𝑢(𝑥)𝑣(𝑥)]=𝑢′(𝑥)𝑣(𝑥)+𝑢(𝑥)𝑣′(𝑥).

In terms of indefinite integrals, this equation becomes

∫𝑑𝑑𝑥[𝑢(𝑥)𝑣(𝑥)]𝑑𝑥=∫[𝑢′(𝑥)𝑣(𝑥)+𝑢(𝑥)𝑣′(𝑥)]𝑑𝑥

or

∫𝑑𝑑𝑥[𝑢(𝑥)𝑣(𝑥)]𝑑𝑥=∫𝑢′(𝑥)𝑣(𝑥)𝑑𝑥+∫𝑢(𝑥)𝑣′(𝑥)𝑑𝑥.

Rearranging the terms of this last equation, we get

∫𝑢(𝑥)𝑣′(𝑥)𝑑𝑥=∫𝑑𝑑𝑥[𝑢(𝑥)𝑣(𝑥)]𝑑𝑥−∫𝑣(𝑥)𝑢′(𝑥)𝑑𝑥,

leading to the following integration by parts formula.

Integration by Parts Formula

∫𝑢(𝑥)𝑣′(𝑥)𝑑𝑥=𝑢(𝑥)𝑣(𝑥)−∫𝑣(𝑥)𝑢′(𝑥)𝑑𝑥(1)

This formula allows us to exchange the problem of computing the integral ∫𝑢(𝑥)𝑣′(𝑥)𝑑𝑥 for the problem of computing a different integral, ∫𝑣(𝑥)𝑢′(𝑥)𝑑𝑥 . In many cases, we can choose the functions u and v so that the second integral is easier to compute than the first. There can be many choices for u and v, and it is not always clear which choice works best, so sometimes we need to try several.

The formula is often given in differential form. With 𝑣′(𝑥) 𝑑𝑥 =𝑑𝜈 and 𝑢′(𝑥) 𝑑𝑥 =𝑑𝑢 , the integration by parts formula becomes

Integration by Parts Formula—Differential Version

∫𝑢𝑑𝑣=𝑢𝑣−∫𝑣𝑑𝑢(2)

The next examples illustrate the technique.

EXAMPLE 1 Find

∫𝑥cos⁡𝑥𝑑𝑥.

Solution There is no obvious antiderivative of 𝑥cos⁡𝑥 , so we use the integration by parts formula

∫𝑢(𝑥)𝑣′(𝑥)𝑑𝑥=𝑢(𝑥)𝑣(𝑥)−∫𝑣(𝑥)𝑢′(𝑥)𝑑𝑥

to change this expression to one that is easier to integrate. We first decide how to choose the functions 𝑢(𝑥) and 𝑣(𝑥) . There is more than one way to do this, but here we choose to factor the expression 𝑥cos⁡𝑥 into

𝑢(𝑥)=𝑥 and 𝑣′(𝑥)=cos⁡𝑥.

Next we differentiate 𝑢(𝑥) and find an antiderivative of 𝑣′(𝑥) ,

𝑢′(𝑥)=1 and 𝑣(𝑥)=sin⁡𝑥.

When finding an antiderivative for 𝑣′(𝑥) , we have a choice of how to pick a constant of integration C. We choose the constant C = 0, since that makes this antiderivative as simple as possible. We now apply the integration by parts formula:

∫𝑢(𝑥)𝑥cos⁡𝑥𝑑𝑥=𝑥sin⁡𝑥−∫𝑣(𝑥)sin⁡𝑥(1)𝑑𝑥 Integration by parts formula =𝑥sin⁡𝑥+cos⁡𝑥+𝐶 Integrate and simplify. 

and we have found the integral of the original function.

There are at least four apparent choices available for 𝑢(𝑥) and 𝑣′(𝑥) in Example 1:

𝟏 . 𝐋 𝐞 𝐭𝑢(𝑥)=1 and 𝑣′(𝑥)=𝑥cos⁡𝑥.𝟐 . 𝐋 𝐞 𝐭𝑢(𝑥)=𝑥 and 𝑣′(𝑥)=cos⁡𝑥.𝟑 . 𝐋 𝐞 𝐭𝑢(𝑥)=𝑥cos⁡𝑥 and 𝑣′(𝑥)=1.𝟒 . 𝐋 𝐞 𝐭𝑢(𝑥)=cos⁡𝑥 and 𝑣′(𝑥)=𝑥.

We used choice 2 in Example 1. The other three choices lead to integrals that we do not know how to evaluate. For instance, Choice 3, with 𝑢′(𝑥) =cos⁡𝑥 −𝑥sin⁡𝑥 , leads to the integral

∫(𝑥cos⁡𝑥−𝑥2sin⁡𝑥)𝑑𝑥.

The goal of integration by parts is to go from an integral ∫𝑢(𝑥)𝑣′(𝑥)𝑑𝑥 that we don’t see how to evaluate to an integral ∫𝑣(𝑥)𝑢′(𝑥)𝑑𝑥 that we can evaluate. Generally, we choose 𝑣′(𝑥) first to be as much of the integrand as we can readily integrate; then we let 𝑢(𝑥) be the leftover part. When finding 𝑣(𝑥) from 𝑣′(𝑥) , any antiderivative will work, and we usually pick the simplest one. In particular, no arbitrary constant of integration is needed in 𝑣(𝑥) because it would simply cancel out of the right-hand side of Equation (2).

 **EXAMPLE 2**  Find ∫ln⁡𝑥𝑑𝑥.

Solution We have not yet seen how to find an antiderivative for ln⁡𝑥 . If we set 𝑢(𝑥) =ln⁡𝑥 , then 𝑢′(𝑥) is the simpler function 1/x. It may not appear that a second function 𝑣′(𝑥) is multiplying 𝑢(𝑥) =ln⁡𝑥 , but we can choose 𝑣′(𝑥) to be the constant function 𝑣′(𝑥) =1 . We use the integration by parts formula given in Equation (1), with

𝑢(𝑥)=ln⁡𝑥 and 𝑣′(𝑥)=1.

We differentiate 𝑢(𝑥) and find an antiderivative of 𝑣′(𝑥) ,

𝑢′(𝑥)=1𝑥 and 𝑣(𝑥)=𝑥.

Then

∫ln⁡𝑥⋅1𝑑𝑥=(ln⁡𝑥)𝑥−∫𝑥1𝑥𝑑𝑥𝑢(𝑥)𝑣′(𝑥)=𝑢(𝑥)𝑣(𝑥)𝑣(𝑥)𝑢′(𝑥)=𝑥ln⁡𝑥−∫1𝑑𝑥=𝑥ln⁡𝑥−𝑥+𝐶 Integration by parts formula  Simplify and integrate. 

In the following examples we use the differential form to indicate the process of integration by parts. The computations are the same, with du and dv providing shorter expressions for 𝑢′(𝑥) dx and 𝑣′(𝑥) dx.

Sometimes we have to use integration by parts more than once, as in the next example.

EXAMPLE 3 Evaluate

∫𝑥2𝑒𝑥𝑑𝑥.

Solution We use the integration by parts formula given in Equation (1), with

𝑢(𝑥)=𝑥2 and 𝑣′(𝑥)=𝑒𝑥.

We differentiate 𝑢(𝑥) and find an antiderivative of 𝑣′(𝑥) ,

𝑢′(𝑥)=2𝑥 and 𝑣(𝑥)=𝑒𝑥.

We summarize this choice by setting 𝑑𝑢 =𝑢′(𝑥) 𝑑𝑥 and 𝑑𝑣 =𝑣′(𝑥) 𝑑𝑥 , so

𝑑𝑢=2𝑥𝑑𝑥 and 𝑑𝑣=𝑒𝑥𝑑𝑥.

We then have

∫𝑢𝑥2𝑒𝑥𝑑𝑥⏟𝑑𝑣=𝑥2𝑢𝑒𝑥−∫𝑣𝑒𝑥2𝑥𝑑𝑥⏟𝑑𝑢=𝑥2𝑒𝑥−2∫𝑥𝑒𝑥𝑑𝑥 Integration by parts formula 

The new integral is less complicated than the original because the exponent on x is reduced by one. To evaluate the integral on the right, we integrate by parts again with u = x, 𝑑𝑣 =𝑒𝑥𝑑𝑥 . Then du = dx, 𝑣 =𝑒𝑥 , and

∫𝑥𝑒𝑥𝑑𝑥⏟𝑢𝑑𝑣=𝑥𝑒𝑥⏟𝑢𝑣−∫𝑒𝑥𝑑𝑥⏟𝑣𝑑𝑢=𝑥𝑒𝑥−𝑒𝑥+𝐶. Integration by parts Equation(2) 𝑢=𝑥,𝑑𝑣=𝑒𝑥𝑑𝑥𝑣=𝑒𝑥,𝑑𝑢=𝑑𝑥

Using this last evaluation, we then obtain

∫𝑥2𝑒𝑥𝑑𝑥=𝑥2𝑒𝑥−2∫𝑥𝑒𝑥𝑑𝑥=𝑥2𝑒𝑥−2𝑥𝑒𝑥+2𝑒𝑥+𝐶,

where the constant of integration is renamed after substituting for the integral on the right.

The technique of Example 3 works for any integral ∫𝑥𝑛𝑒𝑥𝑑𝑥 in which 𝑛 is a positive integer, because differentiating 𝑥𝑛 will eventually lead to a constant, and repeatedly integrating 𝑒𝑥 is easy.

Integrals like the one in the next example occur in electrical engineering. Their evaluation requires two integrations by parts, followed by solving for the unknown integral.

EXAMPLE 4 Evaluate

∫𝑒𝑥cos⁡𝑥𝑑𝑥.

Solution Let 𝑢 =𝑒𝑥 and dv = cos x dx. Then du = 𝑒𝑥 dx, v = sin x, and

∫𝑒𝑥cos⁡𝑥𝑑𝑥=𝑒𝑥sin⁡𝑥−∫𝑒𝑥sin⁡𝑥𝑑𝑥.

The second integral is like the first except that it has sin⁡𝑥 in place of cos⁡𝑥 . To evaluate it, we use integration by parts with

𝑢=𝑒𝑥,𝑑𝑣=sin⁡𝑥𝑑𝑥,𝑣=−cos⁡𝑥,𝑑𝑢=𝑒𝑥𝑑𝑥.

Then

∫𝑒𝑥cos⁡𝑥𝑑𝑥=𝑒𝑥sin⁡𝑥−(−𝑒𝑥cos⁡𝑥−∫(−cos⁡𝑥)(𝑒𝑥𝑑𝑥))=𝑒𝑥sin⁡𝑥+𝑒𝑥cos⁡𝑥−∫𝑒𝑥cos⁡𝑥𝑑𝑥.

The unknown integral now appears on both sides of the equation, but with opposite signs. Adding the integral to both sides and adding the constant of integration gives

2∫𝑒𝑥cos⁡𝑥𝑑𝑥=𝑒𝑥sin⁡𝑥+𝑒𝑥cos⁡𝑥+𝐶1.

Dividing by 2 and renaming the constant of integration then gives

∫𝑒𝑥cos⁡𝑥𝑑𝑥=𝑒𝑥sin⁡𝑥+𝑒𝑥cos⁡𝑥2+𝐶.

EXAMPLE 5 Obtain a formula that expresses the integral

∫cos𝑛⁡𝑥𝑑𝑥

in terms of an integral of a lower power of cos⁡𝑥 .

Solution We may think of cos𝑛⁡𝑥 as cos𝑛−1⁡𝑥 ⋅cos⁡𝑥 . Then we let

𝑢=cos𝑛−1⁡𝑥 and 𝑑𝑣=cos⁡𝑥𝑑𝑥,

so that

𝑑𝑢=(𝑛−1)(cos𝑛−2⁡𝑥)(−sin⁡𝑥𝑑𝑥) and 𝑣=sin⁡𝑥.

Integration by parts then gives

∫cos𝑛⁡𝑥𝑑𝑥=cos𝑛−1⁡𝑥sin⁡𝑥+(𝑛−1)∫sin2⁡𝑥cos𝑛−2⁡𝑥𝑑𝑥=cos𝑛−1⁡𝑥sin⁡𝑥+(𝑛−1)∫(1−cos2⁡𝑥)cos𝑛−2⁡𝑥𝑑𝑥=cos𝑛−1⁡𝑥sin⁡𝑥+(𝑛−1)∫cos𝑛−2⁡𝑥𝑑𝑥−(𝑛−1)∫cos𝑛⁡𝑥𝑑𝑥.

If we add

(𝑛−1)∫cos𝑛⁡𝑥𝑑𝑥

to both sides of this equation, we obtain

𝑛∫cos𝑛⁡𝑥𝑑𝑥=cos𝑛−1⁡𝑥sin⁡𝑥+(𝑛−1)∫cos𝑛−2⁡𝑥𝑑𝑥.

We then divide through by 𝑛 , and the final result is

∫cos𝑛⁡𝑥𝑑𝑥=cos𝑛−1⁡𝑥sin⁡𝑥𝑛+𝑛−1𝑛∫cos𝑛−2⁡𝑥𝑑𝑥.

The formula found in Example 5 is called a reduction formula because it replaces an integral containing some power of a function with an integral of the same form having the power reduced. When n is a positive integer, we may apply the formula repeatedly until the remaining integral is easy to evaluate. For example, the result in Example 5 tells us that

∫cos3⁡𝑥𝑑𝑥=cos2⁡𝑥sin⁡𝑥3+23∫cos⁡𝑥𝑑𝑥=13cos2⁡𝑥sin⁡𝑥+23sin⁡𝑥+𝐶.

Evaluating Definite Integrals by Parts

The integration by parts formula in Equation (1) can be combined with Part 2 of the Fundamental Theorem in order to evaluate definite integrals by parts. Assuming that both 𝑢′ and 𝑣′ are continuous over the interval [𝑎,𝑏] , Part 2 of the Fundamental Theorem gives

Integration by Parts Formula for Definite Integrals

∫𝑏𝑎𝑢(𝑥)𝑣′(𝑥)𝑑𝑥=𝑢(𝑥)𝑣(𝑥)]𝑏𝑎−∫𝑏𝑎𝑣(𝑥)𝑢′(𝑥)𝑑𝑥(3)

教材插图

FIGURE 8.1 The region in Example 6.

EXAMPLE 6 Find the area of the region bounded by the curve 𝑦 =𝑥𝑒−𝑥 and the 𝑥 -axis from 𝑥 =0 to 𝑥 =4 .

Solution The region is shaded in Figure 8.1. Its area is

∫40𝑥𝑒−𝑥𝑑𝑥.

Let 𝑢 =𝑥 , 𝑑𝑣 =𝑒−𝑥𝑑𝑥 , 𝑣 = −𝑒−𝑥 , and 𝑑𝑢 =𝑑𝑥 . Then

∫40𝑥𝑒−𝑥𝑑𝑥 = −𝑥𝑒−𝑥∣40 −∫40( −𝑒−𝑥)𝑑𝑥 Integration by parts Formula (3) =[ −4𝑒−4 −( −0𝑒−0)] +∫40𝑒−𝑥𝑑𝑥 = −4𝑒−4 −𝑒−𝑥∣40 = −4𝑒−4 −(𝑒−4 −𝑒−0) =1 −5𝑒−4 ≈0.91.

EXERCISES 8.2

Integration by Parts

Evaluate the integrals in Exercises 1–24 using integration by parts.

  1. ∫𝑥sin⁡𝑥2𝑑𝑥

  2. ∫𝜃cos⁡𝜋𝜃𝑑𝜃

  3. ∫𝑡2cos⁡𝑡𝑑𝑡

  4. ∫𝑥2sin⁡𝑥𝑑𝑥

  5. ∫21𝑥ln⁡𝑥𝑑𝑥

  6. ∫𝑒1𝑥3ln⁡𝑥𝑑𝑥

  7. ∫𝑥𝑒𝑥𝑑𝑥

  8. ∫𝑥𝑒3𝑥𝑑𝑥

  9. ∫𝑥2𝑒−𝑥𝑑𝑥

  10. ∫(𝑥2 −2𝑥 +1)𝑒2𝑥𝑑𝑥

  11. ∫tan−1⁡𝑦𝑑𝑦

  12. ∫arcsin⁡𝑦𝑑𝑦

  13. ∫𝑥sec2⁡𝑥𝑑𝑥

  14. ∫4𝑥sec2⁡2𝑥𝑑𝑥

  15. ∫𝑥3𝑒𝑥𝑑𝑥

  16. ∫𝑝4𝑒−𝑝𝑑𝑝

  17. ∫(𝑥2 −5𝑥)𝑒𝑥𝑑𝑥

  18. ∫(𝑟2 +𝑟 +1)𝑒𝑟𝑑𝑟

∫√𝑥ln⁡𝑥𝑑𝑥
  1. ∫𝑥5𝑒𝑥𝑑𝑥

  2. ∫𝑡2𝑒4𝑡𝑑𝑡

  3. ∫𝑒𝜃sin⁡𝜃𝑑𝜃

  4. ∫𝑒−𝑦cos⁡𝑦𝑑𝑦

  5. ∫𝑒2𝑥cos⁡3𝑥𝑑𝑥

  6. ∫𝑒−2𝑥sin⁡2𝑥𝑑𝑥

Using Substitution

Evaluate the integrals in Exercises 25–30 by using a substitution prior to integration by parts.

  1. ∫𝑒√3𝑠+9𝑑𝑠

  2. ∫10𝑥√1−𝑥𝑑𝑥

  3. ∫𝜋/30𝑥tan2⁡𝑥𝑑𝑥

  4. ∫ln⁡(𝑥 +𝑥2)𝑑𝑥

  5. ∫sin⁡(ln⁡𝑥)𝑑𝑥

  6. ∫𝑧(ln⁡𝑧)2𝑑𝑧

Evaluating Integrals

Evaluate the integrals in Exercises 31–56. Some integrals do not require integration by parts.

  1. ∫𝑥sec⁡𝑥2𝑑𝑥

  2. ∫cos⁡√𝑥√𝑥𝑑𝑥

  3. ∫𝑥(ln⁡𝑥)2𝑑𝑥

  4. ∫1𝑥(ln⁡𝑥)2𝑑𝑥

  5. ∫ln⁡𝑥𝑥2𝑑𝑥

  6. ∫(ln⁡𝑥)3𝑥𝑑𝑥

  7. ∫𝑥3𝑒𝑥4𝑑𝑥

  8. ∫𝑥5𝑒𝑥3𝑑𝑥

  9. ∫𝑥3√𝑥2+1𝑑𝑥

  10. ∫𝑥2sin⁡𝑥3𝑑𝑥

  11. ∫sin⁡3𝑥cos⁡2𝑥𝑑𝑥

  12. ∫sin⁡2𝑥cos⁡4𝑥𝑑𝑥

  13. ∫𝑒√𝑥√𝑥𝑑𝑥

  14. ∫cos⁡√𝑥𝑑𝑥

  15. ∫√𝑥𝑒√𝑥𝑑𝑥

  16. ∫𝜋/20𝜃2sin⁡2𝜃𝑑𝜃

  17. ∫𝜋/20𝑥3cos⁡2𝑥𝑑𝑥

  18. ∫22/√3𝑡sec−1⁡𝑡𝑑𝑡

  19. ∫1/√202𝑥arcsin⁡(𝑥2)𝑑𝑥

  20. ∫𝑥arctan⁡𝑥𝑑𝑥

  21. ∫𝑥2tan−1⁡𝑥2𝑑𝑥

  22. ∫(1 +2𝑥2)𝑒𝑥2𝑑𝑥

  23. ∫𝑥𝑒𝑥(𝑥+1)2𝑑𝑥

  24. ∫√𝑥(arcsin⁡√𝑥)𝑑𝑥

  25. ∫(sin−1⁡𝑥)2√1−𝑥2𝑑𝑥

Theory and Examples

  1. Finding area Find the area of the region enclosed by the curve 𝑦 =𝑥sin⁡𝑥 and the x-axis (see the accompanying figure) for

a. 0 ≤𝑥 ≤𝜋.

b. 𝜋 ≤𝑥 ≤2𝜋.

c. 2𝜋 ≤𝑥 ≤3𝜋.

d. What pattern do you see here? What is the area between the curve and the x-axis for 𝑛𝜋 ≤𝑥 ≤(𝑛 +1)𝜋 , n an arbitrary nonnegative integer? Give reasons for your answer.

教材插图

  1. Finding area Find the area of the region enclosed by the curve 𝑦 =𝑥cos⁡𝑥 and the x-axis (see the accompanying figure) for

a. 𝜋/2 ≤𝑥 ≤3𝜋/2.

b. 3𝜋/2 ≤𝑥 ≤5𝜋/2 .

c. 5𝜋/2 ≤𝑥 ≤7𝜋/2.

d. What pattern do you see? What is the area between the curve and the x-axis for

(2𝑛−12)𝜋≤𝑥≤(2𝑛+12)𝜋,

n an arbitrary positive integer? Give reasons for your answer.

教材插图

  1. Finding volume Find the volume of the solid generated by revolving the region in the first quadrant bounded by the coordinate axes, the curve 𝑦 =𝑒𝑥 , and the line 𝑥 =ln⁡2 about the line 𝑥 =ln⁡2 .

  2. Finding volume Find the volume of the solid generated by revolving the region in the first quadrant bounded by the coordinate axes, the curve 𝑦 =𝑒−𝑥 , and the line x = 1

a. about the y-axis.

b. about the line x = 1.

  1. Finding volume Find the volume of the solid generated by revolving the region in the first quadrant bounded by the coordinate axes and the curve 𝑦 =cos⁡𝑥 , 0 ≤𝑥 ≤𝜋/2 , about

a. the y-axis.

b. the line 𝑥 =𝜋/2 .

  1. Finding volume Find the volume of the solid generated by revolving the region bounded by the x-axis and the curve 𝑦 =𝑥sin⁡𝑥,0 ≤𝑥 ≤𝜋 , about

a. the y-axis.

b. the line 𝑥 =𝜋 .

(See Exercise 57 for a graph.)

  1. Consider the region bounded by the graphs of 𝑦 =ln⁡𝑥 , 𝑦 =0 , and 𝑥 =𝑒 .

a. Find the area of the region.

b. Find the volume of the solid formed by revolving this region about the x-axis.

c. Find the volume of the solid formed by revolving this region about the line x = -2.

d. Find the centroid of the region.

  1. Consider the region bounded by the graphs of 𝑦 =arctan⁡𝑥 , 𝑦 =0 , and 𝑥 =1 .

a. Find the area of the region.

b. Find the volume of the solid formed by revolving this region about the y-axis.

  1. Average value A retarding force, symbolized by the dashpot in the accompanying figure, slows the motion of the weighted spring so that the mass’s position at time 𝑡 is
𝑦=2𝑒−𝑡cos⁡𝑡,𝑡≥0.

Find the average value of 𝑦 over the interval 0 ≤𝑡 ≤2𝜋 .

教材插图

  1. Average value In a mass-spring-dashpot system like the one in Exercise 65, the mass’s position at time 𝑡 is
𝑦=4𝑒−𝑡(sin⁡𝑡−cos⁡𝑡),𝑡≥0.

Find the average value of 𝑦 over the interval 0 ≤𝑡 ≤2𝜋 .

Reduction Formulas

In Exercises 67–73, use integration by parts to establish the reduction formula.

  1. ∫𝑥𝑛cos⁡𝑥𝑑𝑥 =𝑥𝑛sin⁡𝑥 −𝑛∫𝑥𝑛−1sin⁡𝑥𝑑𝑥

  2. ∫𝑥𝑛sin⁡𝑥𝑑𝑥 = −𝑥𝑛cos⁡𝑥 +𝑛∫𝑥𝑛−1cos⁡𝑥𝑑𝑥

  3. ∫𝑥𝑛𝑒𝑎𝑥𝑑𝑥 =𝑥𝑛𝑒𝑎𝑥𝑎 −𝑛𝑎∫𝑥𝑛−1𝑒𝑎𝑥𝑑𝑥,𝑎 ≠0

  4. ∫(ln⁡𝑥)𝑛𝑑𝑥 =𝑥(ln⁡𝑥)𝑛 −𝑛∫(ln⁡𝑥)𝑛−1𝑑𝑥

71.∫𝑥𝑚(ln⁡𝑥)𝑛𝑑𝑥=𝑥𝑚+1𝑚+1(ln⁡𝑥)𝑛−𝑛𝑚+1∫𝑥𝑚(ln⁡𝑥)𝑛−1𝑑𝑥,𝑚≠−1 72.∫𝑥𝑛√𝑥+1𝑑𝑥=2𝑥𝑛2𝑛+3(𝑥+1)3/2−2𝑛2𝑛+3∫𝑥𝑛−1√𝑥+1𝑑𝑥 73.∫𝑥𝑛√𝑥+1𝑑𝑥=2𝑥𝑛2𝑛+1√𝑥+1−2𝑛2𝑛+1∫𝑥𝑛−1√𝑥+1𝑑𝑥74.$𝑈𝑠𝑒𝐸𝑥𝑎𝑚𝑝𝑙𝑒5𝑡𝑜𝑠ℎ𝑜𝑤𝑡ℎ𝑎𝑡$∫𝜋/20sin𝑛⁡𝑥𝑑𝑥=∫𝜋/20cos𝑛⁡𝑥𝑑𝑥=⎧{ {⎨{ {⎩(𝜋2)1⋅3⋅5⋯(𝑛−1)2⋅4⋅6⋯𝑛,𝑛even2⋅4⋅6⋯(𝑛−1)1⋅3⋅5⋯𝑛,𝑛odd75.$𝑆ℎ𝑜𝑤𝑡ℎ𝑎𝑡$∫𝑏𝑎(∫𝑏𝑥𝑓(𝑡)𝑑𝑡)𝑑𝑥=∫𝑏𝑎(𝑥−𝑎)𝑓(𝑥)𝑑𝑥.76.$𝑈𝑠𝑒𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑡𝑖𝑜𝑛𝑏𝑦𝑝𝑎𝑟𝑡𝑠𝑡𝑜𝑜𝑏𝑡𝑎𝑖𝑛𝑡ℎ𝑒𝑓𝑜𝑟𝑚𝑢𝑙𝑎$∫√1−𝑥2𝑑𝑥=12𝑥√1−𝑥2+12∫1√1−𝑥2𝑑𝑥.

Integrating Inverses of Functions

Integration by parts leads to a rule for integrating inverses that usually gives good results:

∫𝑓−1(𝑥)𝑑𝑥=∫𝑦𝑓′(𝑦)𝑑𝑦𝑦=𝑓−1(𝑥),𝑥=𝑓(𝑦)𝑑𝑥=𝑓′(𝑦)𝑑𝑦=𝑦𝑓(𝑦)−∫𝑓(𝑦)𝑑𝑦 Integration by parts with 𝑢=𝑦,𝑑𝑣=𝑓′(𝑦)𝑑𝑦=𝑥𝑓−1(𝑥)−∫𝑓(𝑦)𝑑𝑦

The idea is to take the most complicated part of the integral, in this case 𝑓−1(𝑥) , and simplify it first. For the integral of ln⁡𝑥 , we get

∫ln⁡𝑥𝑑𝑥=∫𝑦𝑒𝑦𝑑𝑦𝑦=ln⁡𝑥,𝑥=𝑒𝑦𝑑𝑥=𝑒𝑦𝑑𝑦=𝑦𝑒𝑦−𝑒𝑦+𝐶=𝑥ln⁡𝑥−𝑥+𝐶.

For the integral of arccos⁡𝑥 , we get

∫arccos⁡𝑥𝑑𝑥=𝑥arccos⁡𝑥−∫cos⁡𝑦𝑑𝑦=𝑥arccos⁡𝑥−sin⁡𝑦+𝐶=𝑥arccos⁡𝑥−sin⁡(arccos⁡𝑥)+𝐶.𝑦=arccos⁡𝑥

Use the formula

∫𝑓−1(𝑥)𝑑𝑥=𝑥𝑓−1(𝑥)−∫𝑓(𝑦)𝑑𝑦𝑦=𝑓−1(𝑥)(4)

to evaluate the integrals in Exercises 77–80. Express your answers in terms of x.

  1. ∫arcsec⁡𝑥 𝑑𝑥

  2. ∫arctan⁡𝑥 𝑑𝑥

  3. ∫sec−1⁡𝑥 𝑑𝑥

  4. ∫log2⁡𝑥 𝑑𝑥

Another way to integrate 𝑓−1(𝑥) (when 𝑓−1 is integrable) is to use integration by parts with 𝑢 =𝑓−1(𝑥) and dv = dx to rewrite the integral of 𝑓−1 as

∫𝑓−1(𝑥)𝑑𝑥=𝑥𝑓−1(𝑥)−∫𝑥(𝑑𝑑𝑥𝑓−1(𝑥))𝑑𝑥.(5)

Exercises 81 and 82 compare the results of using Equations (4) and (5).

  1. Equations (4) and (5) give different formulas for the integral of arccos⁡𝑥 :
𝐚.∫arccos⁡𝑥𝑑𝑥=𝑥arccos⁡𝑥−sin⁡(arccos⁡𝑥)+𝐶(Eq. (4) ∫arccos⁡𝑥𝑑𝑥=𝑥arccos⁡𝑥−√1−𝑥2+𝐶(Eq. (5)

Can both integrations be correct? Explain.

  1. Equations (4) and (5) lead to different formulas for the integral of arctan⁡𝑥 :
𝐚.∫arctan⁡𝑥𝑑𝑥=𝑥arctan⁡𝑥−ln⁡sec⁡(arctan⁡𝑥)+𝐶(Eq. (4) ∫arctan⁡𝑥𝑑𝑥=𝑥arctan⁡𝑥−ln⁡√1+𝑥2+𝐶(Eq. (5)

Can both integrations be correct? Explain.

Evaluate the integrals in Exercises 83 and 84 with (a) Eq. (4) and (b) Eq. (5). In each case, check your work by differentiating your answer with respect to x.

83.∫sinh−1⁡𝑥𝑑𝑥84.∫tanh−1⁡𝑥𝑑𝑥

8.3 Trigonometric Integrals

Trigonometric integrals involve algebraic combinations of the six basic trigonometric functions. In principle, we can always express such integrals in terms of sines and cosines, but it is often simpler to work with other functions, as in the integral

∫sec2⁡𝑥𝑑𝑥=tan⁡𝑥+𝐶.

The general idea is to use identities to transform the integrals we must find into integrals that are easier to work with.

Products of Powers of Sines and Cosines

We begin with integrals of the form

∫sin𝑚⁡𝑥cos𝑛⁡𝑥𝑑𝑥,

where m and n are nonnegative integers (positive or zero). We can divide the appropriate substitution into three cases according to m and n being odd or even.

Case 1 If 𝑚 is odd in ∫sin𝑚⁡𝑥cos𝑛⁡𝑥𝑑𝑥 , we write 𝑚 as 2𝑘 +1 and use the identity sin2⁡𝑥 =1 −cos2⁡𝑥 to obtain

sin𝑚⁡𝑥=sin2𝑘+1⁡𝑥=(sin2⁡𝑥)𝑘sin⁡𝑥=(1−cos2⁡𝑥)𝑘sin⁡𝑥.(1)

Then we substitute 𝑢 =cos⁡𝑥 and 𝑑𝑢 = −sin⁡𝑥𝑑𝑥 .

Case 2 If 𝑛 is odd in ∫sin𝑚⁡𝑥cos𝑛⁡𝑥𝑑𝑥 , we write 𝑛 as 2𝑘 +1 and use the identity cos2⁡𝑥 =1 −sin2⁡𝑥 to obtain

cos𝑛⁡𝑥=cos2𝑘+1⁡𝑥=(cos2⁡𝑥)𝑘cos⁡𝑥=(1−sin2⁡𝑥)𝑘cos⁡𝑥.

We then substitute 𝑢 =sin⁡𝑥 and 𝑑𝑢 =cos⁡𝑥𝑑𝑥 .

Case 3 If both m and n are even in ∫sin𝑚⁡𝑥cos𝑛⁡𝑥𝑑𝑥 , we substitute

sin2⁡𝑥=1−cos⁡2𝑥2,cos2⁡𝑥=1+cos⁡2𝑥2(2)

to reduce the integrand to one in lower powers of cos⁡2𝑥 .

Here are some examples illustrating each case.

EXAMPLE 1 Evaluate ∫sin3⁡𝑥cos2⁡𝑥𝑑𝑥 .

Solution This is an example of Case 1.

∫sin3⁡𝑥cos2⁡𝑥𝑑𝑥=∫sin2⁡𝑥cos2⁡𝑥sin⁡𝑥𝑑𝑥𝑚 is odd. =∫(1−cos2⁡𝑥)(cos2⁡𝑥)sin⁡𝑥𝑑𝑥sin2⁡𝑥=1−cos2⁡𝑥=∫(1−𝑢2)(𝑢2)(−𝑑𝑢)𝑢=cos⁡𝑥,𝑑𝑢=−sin⁡𝑥𝑑𝑥=∫(𝑢4−𝑢2)𝑑𝑢 Distribute. =𝑢55−𝑢33+𝐶=cos5⁡𝑥5−cos3⁡𝑥3+𝐶◼

EXAMPLE 2 Evaluate

∫cos5⁡𝑥𝑑𝑥.

Solution This is an example of Case 2, where m = 0 is even and n = 5 is odd.

∫cos5⁡𝑥𝑑𝑥=∫cos4⁡𝑥cos⁡𝑥𝑑𝑥=∫(1−sin2⁡𝑥)2cos⁡𝑥𝑑𝑥=∫(1−𝑢2)2𝑑𝑢𝑢=sin⁡𝑥,𝑑𝑢=cos⁡𝑥𝑑𝑥=∫(1−2𝑢2+𝑢4)𝑑𝑢 Square 1−𝑢2.=𝑢−23𝑢3+15𝑢5+𝐶=sin⁡𝑥−23sin3⁡𝑥+15sin5⁡𝑥+𝐶

EXAMPLE 3 Evaluate

∫sin2⁡𝑥cos4⁡𝑥𝑑𝑥.

Solution This is an example of Case 3.

∫sin2⁡𝑥cos4⁡𝑥𝑑𝑥=∫(1−cos⁡2𝑥2)(1+cos⁡2𝑥2)2𝑑𝑥𝑚 and 𝑛 both even =18∫(1−cos⁡2𝑥)(1+2cos⁡2𝑥+cos2⁡2𝑥)𝑑𝑥=18∫(1+cos⁡2𝑥−cos2⁡2𝑥−cos3⁡2𝑥)𝑑𝑥=18[𝑥+12sin⁡2𝑥−∫(cos2⁡2𝑥+cos3⁡2𝑥)𝑑𝑥]

For the term involving cos2⁡2𝑥 , we use

∫cos2⁡2𝑥𝑑𝑥=12∫(1+cos⁡4𝑥)𝑑𝑥 Use the identity cos2⁡𝜃=(1+cos⁡2𝜃)/2, with 𝜃=2𝑥.=12(𝑥+14sin⁡4𝑥)+𝐶1.

For the cos3⁡2𝑥 term, we have

∫cos3⁡2𝑥𝑑𝑥=∫(1−sin2⁡2𝑥)cos⁡2𝑥𝑑𝑥=12∫(1−𝑢2)𝑑𝑢=12(sin⁡2𝑥−13sin3⁡2𝑥)+𝐶2.

Combining everything and simplifying, we get

∫sin2⁡𝑥cos4⁡𝑥𝑑𝑥=116(𝑥−14sin⁡4𝑥+13sin3⁡2𝑥)+𝐶.

Eliminating Square Roots

In the next example, we use the identity cos2⁡𝜃 =(1 +cos⁡2𝜃)/2 to eliminate a square root.

EXAMPLE 4 Evaluate

∫𝜋/40√1+cos⁡4𝑥𝑑𝑥.

Solution To eliminate the square root, we use the identity

cos2⁡𝜃=1+cos⁡2𝜃2 or 1+cos⁡2𝜃=2cos2⁡𝜃.

With 𝜃 =2𝑥 , this becomes

1+cos⁡4𝑥=2cos2⁡2𝑥.

Therefore,

∫𝜋/40√1+cos⁡4𝑥𝑑𝑥=∫𝜋/40√2cos2⁡2𝑥𝑑𝑥=∫𝜋/40√2√cos2⁡2𝑥𝑑𝑥=√2∫𝜋/40|cos⁡2𝑥|𝑑𝑥=√2∫𝜋/40cos⁡2𝑥𝑑𝑥 on [0,𝜋/4]=√2[sin⁡2𝑥2]𝜋/40=√22[1−0]=√22.

Integrals of Powers of tan x and sec x

We know how to integrate the tangent and secant functions and their squares. To integrate higher powers, we use the identities tan2⁡𝑥 =sec2⁡𝑥 −1 and sec2⁡𝑥 =tan2⁡𝑥 +1 , and integrate by parts when necessary to reduce the higher powers to lower powers.

EXAMPLE 5 Evaluate

∫tan4⁡𝑥𝑑𝑥.

Solution

∫tan4⁡𝑥𝑑𝑥=∫tan2⁡𝑥⋅tan2⁡𝑥𝑑𝑥=∫tan2⁡𝑥⋅(sec2⁡𝑥−1)𝑑𝑥tan2⁡𝑥=sec2⁡𝑥−1=∫tan2⁡𝑥sec2⁡𝑥𝑑𝑥−∫tan2⁡𝑥𝑑𝑥=∫tan2⁡𝑥sec2⁡𝑥𝑑𝑥−∫(sec2⁡𝑥−1)𝑑𝑥tan2⁡𝑥=sec2⁡𝑥−1=∫tan2⁡𝑥sec2⁡𝑥𝑑𝑥−∫sec2⁡𝑥𝑑𝑥+∫𝑑𝑥

In the first integral, we let

𝑢=tan⁡𝑥,𝑑𝑢=sec2⁡𝑥𝑑𝑥

and have

∫𝑢2𝑑𝑢=13𝑢3+𝐶1.

The remaining integrals are standard forms, so

∫tan4⁡𝑥𝑑𝑥=13tan3⁡𝑥−tan⁡𝑥+𝑥+𝐶.

EXAMPLE 6 Evaluate

∫sec3⁡𝑥𝑑𝑥.

Solution We integrate by parts using

𝑢=sec⁡𝑥,𝑑𝑣=sec2⁡𝑥𝑑𝑥,𝑣=tan⁡𝑥,𝑑𝑢=sec⁡𝑥tan⁡𝑥𝑑𝑥.

Then

∫sec3⁡𝑥𝑑𝑥=sec⁡𝑥tan⁡𝑥−∫(tan⁡𝑥)(sec⁡𝑥tan⁡𝑥)𝑑𝑥=sec⁡𝑥tan⁡𝑥−∫(sec2⁡𝑥−1)sec⁡𝑥𝑑𝑥=sec⁡𝑥tan⁡𝑥−∫sec3⁡𝑥𝑑𝑥+∫sec⁡𝑥𝑑𝑥 Integrate by parts. tan2⁡𝑥=sec2⁡𝑥−1

Combining the two secant-cubed integrals gives

2∫sec3⁡𝑥𝑑𝑥=sec⁡𝑥tan⁡𝑥+∫sec⁡𝑥𝑑𝑥

and therefore

∫sec3⁡𝑥𝑑𝑥=12sec⁡𝑥tan⁡𝑥+12ln⁡|sec⁡𝑥+tan⁡𝑥|+𝐶.

EXAMPLE 7 Evaluate

∫tan4⁡𝑥sec4⁡𝑥𝑑𝑥.

Solution

∫(tan4⁡𝑥)(sec4⁡𝑥)𝑑𝑥=∫(tan4⁡𝑥)(1+tan2⁡𝑥)(sec2⁡𝑥)𝑑𝑥sec2⁡𝑥=1+tan2⁡𝑥=∫(tan4⁡𝑥+tan6⁡𝑥)(sec2⁡𝑥)𝑑𝑥 Distribute. =∫(𝑢4+𝑢6)𝑑𝑢=𝑢55+𝑢77+𝐶𝑢=tan⁡𝑥,=tan5⁡𝑥5+tan7⁡𝑥7+𝐶𝑑𝑢=sec2⁡𝑥𝑑𝑥

Products of Sines and Cosines

The integrals

∫sin⁡𝑚𝑥sin⁡𝑛𝑥𝑑𝑥,∫sin⁡𝑚𝑥cos⁡𝑛𝑥𝑑𝑥, and ∫cos⁡𝑚𝑥cos⁡𝑛𝑥𝑑𝑥

arise in many applications involving periodic functions. We can evaluate these integrals through integration by parts, but two such integrations are required in each case. It is simpler to use the following identities.

sin⁡𝑚𝑥sin⁡𝑛𝑥=12[cos⁡(𝑚−𝑛)𝑥−cos⁡(𝑚+𝑛)𝑥](3) sin⁡𝑚𝑥cos⁡𝑛𝑥=12[sin⁡(𝑚−𝑛)𝑥+sin⁡(𝑚+𝑛)𝑥](4) cos⁡𝑚𝑥cos⁡𝑛𝑥=12[cos⁡(𝑚−𝑛)𝑥+cos⁡(𝑚+𝑛)𝑥](5)

These identities come from the angle sum formulas for the sine and cosine functions (Section 1.3). They give functions whose antiderivatives are easily found.

EXAMPLE 8 Evaluate

∫sin⁡3𝑥cos⁡5𝑥𝑑𝑥.

Solution From Equation (4) with m = 3 and n = 5, we get

∫sin⁡3𝑥cos⁡5𝑥𝑑𝑥=12∫[sin⁡(−2𝑥)+sin⁡8𝑥]𝑑𝑥=12∫(sin⁡8𝑥−sin⁡2𝑥)𝑑𝑥=−cos⁡8𝑥16+cos⁡2𝑥4+𝐶.

EXERCISES 8.3

Powers of Sines and Cosines

Evaluate the integrals in Exercises 1–22.

  1. ∫cos⁡2𝑥𝑑𝑥

  2. ∫𝜋03sin⁡𝑥3𝑑𝑥

  3. ∫cos3⁡𝑥sin⁡𝑥𝑑𝑥

  4. ∫sin4⁡2𝑥cos⁡2𝑥𝑑𝑥

  5. ∫sin3⁡𝑥𝑑𝑥

  6. ∫cos3⁡4𝑥𝑑𝑥

  7. ∫sin5⁡𝑥𝑑𝑥

  8. ∫𝜋0sin5⁡𝑥2𝑑𝑥

  9. ∫cos3⁡𝑥𝑑𝑥

  10. ∫𝜋/603cos5⁡3𝑥𝑑𝑥

  11. ∫sin3⁡𝑥cos3⁡𝑥𝑑𝑥

  12. ∫cos3⁡2𝑥sin5⁡2𝑥𝑑𝑥

  13. ∫cos2⁡𝑥𝑑𝑥

  14. ∫𝜋/20sin2⁡𝑥𝑑𝑥

  15. ∫𝜋/20sin7⁡𝑦𝑑𝑦

  16. ∫7cos7⁡𝑡𝑑𝑡

  17. ∫𝜋08sin4⁡𝑥𝑑𝑥

  18. ∫8cos4⁡2𝜋𝑥𝑑𝑥

  19. ∫16sin2⁡𝑥cos2⁡𝑥𝑑𝑥

  20. ∫𝜋08sin4⁡𝑦cos2⁡𝑦𝑑𝑦

  21. ∫8cos3⁡2𝜃sin⁡2𝜃𝑑𝜃

  22. ∫𝜋/20sin2⁡2𝜃cos3⁡2𝜃𝑑𝜃

Integrating Square Roots

Evaluate the integrals in Exercises 23–32.

  1. ∫2𝜋0√1−cos⁡𝑥2𝑑𝑥

  2. ∫𝜋0√1−cos⁡2𝑥𝑑𝑥

  3. ∫𝜋0√1−sin2⁡𝑡𝑑𝑡

  4. ∫𝜋0√1−cos2⁡𝜃𝑑𝜃

Exercises 59–64 require the use of various trigonometric identities before you evaluate the integrals.

  1. ∫𝜋/2𝜋/3sin2⁡𝑥√1−cos⁡𝑥𝑑𝑥

  2. ∫𝜋/60√1+sin⁡𝑥𝑑𝑥

  3. ∫𝜋5𝜋/6cos4⁡𝑥√1−sin⁡𝑥𝑑𝑥

  4. ∫3𝜋/4𝜋/2√1−sin⁡2𝑥𝑑𝑥

Assorted Integrations

  1. ∫𝜋/20𝜃√1−cos⁡2𝜃𝑑𝜃

  2. ∫𝜋−𝜋(1 −cos2⁡𝑡)3/2𝑑𝑡

Powers of Tangents and Secants

Evaluate the integrals in Exercises 33–52.

  1. ∫sec2⁡𝑥tan⁡𝑥𝑑𝑥

  2. ∫sec⁡𝑥tan2⁡𝑥𝑑𝑥

  3. ∫sec3⁡𝑥tan⁡𝑥𝑑𝑥

  4. ∫sec3⁡𝑥tan3⁡𝑥𝑑𝑥

Applications

  1. ∫sec2⁡𝑥tan2⁡𝑥𝑑𝑥

  2. ∫sec4⁡𝑥tan2⁡𝑥𝑑𝑥

  3. ∫0−𝜋/32sec3⁡𝑥𝑑𝑥

  4. ∫𝑒𝑥sec3⁡𝑒𝑥𝑑𝑥

  5. ∫sec4⁡𝜃𝑑𝜃

  6. ∫tan4⁡𝑥sec3⁡𝑥𝑑𝑥

  7. ∫𝜋/2𝜋/4csc4⁡𝜃𝑑𝜃

  8. ∫sec6⁡𝑥𝑑𝑥

  9. ∫4tan3⁡𝑥𝑑𝑥

  10. ∫𝜋/4−𝜋/46tan4⁡𝑥𝑑𝑥

  11. ∫tan5⁡𝑥𝑑𝑥

  12. ∫cot6⁡2𝑥𝑑𝑥

  13. ∫𝜋/3𝜋/6cot3⁡𝑥𝑑𝑥

  14. ∫8cot4⁡𝑡𝑑𝑡

  15. ∫𝜋/3𝜋/4tan5⁡𝜃sec4⁡𝜃𝑑𝜃

  16. ∫cot3⁡𝑡csc4⁡𝑡𝑑𝑡

Products of Sines and Cosines

Evaluate the integrals in Exercises 53–58.

  1. ∫sin⁡3𝑥cos⁡2𝑥𝑑𝑥

  2. ∫sin⁡2𝑥cos⁡3𝑥𝑑𝑥

  3. ∫𝜋−𝜋sin⁡3𝑥sin⁡3𝑥𝑑𝑥

  4. ∫𝜋/20sin⁡𝑥cos⁡𝑥𝑑𝑥

  5. ∫cos⁡3𝑥cos⁡4𝑥𝑑𝑥

  6. ∫𝜋/2−𝜋/2cos⁡𝑥cos⁡7𝑥𝑑𝑥

  7. ∫sin2⁡𝜃cos⁡3𝜃𝑑𝜃

  8. ∫cos2⁡2𝜃sin⁡𝜃𝑑𝜃

Hint: Multiply by √1−sin⁡𝑥1−sin⁡𝑥 .

  1. ∫cos3⁡𝜃sin⁡2𝜃𝑑𝜃

  2. ∫sin3⁡𝜃cos⁡2𝜃𝑑𝜃

  3. ∫sin⁡𝜃cos⁡𝜃cos⁡3𝜃𝑑𝜃

  4. ∫sin⁡𝜃sin⁡2𝜃sin⁡3𝜃𝑑𝜃

Use any method to evaluate the integrals in Exercises 65–70.

  1. ∫sec3⁡𝑥tan⁡𝑥𝑑𝑥

  2. ∫sin3⁡𝑥cos4⁡𝑥𝑑𝑥

  3. ∫tan2⁡𝑥csc⁡𝑥𝑑𝑥

  4. ∫cot⁡𝑥cos2⁡𝑥𝑑𝑥

  5. ∫𝑥sin2⁡𝑥𝑑𝑥

  6. ∫𝑥cos3⁡𝑥𝑑𝑥

  7. Arc length Find the length of the curve 𝑦 =ln⁡(sin⁡𝑥), 𝜋6 ≤𝑥 ≤𝜋2.

  8. Center of gravity Find the center of gravity of the region bounded by the 𝑥 -axis, the curve 𝑦 =sec⁡𝑥 , and the lines 𝑥 = −𝜋/4 , 𝑥 =𝜋/4 .

  9. Volume Find the volume generated by revolving one arch of the curve 𝑦 =sin⁡𝑥 about the x-axis.

  10. Area Find the area between the 𝑥 -axis and the curve 𝑦 =√1+cos⁡4𝑥 , 0 ≤𝑥 ≤𝜋 .

  11. Centroid Find the centroid of the region bounded by the graphs of 𝑦 =𝑥 +cos⁡𝑥 and 𝑦 =0 for 0 ≤𝑥 ≤2𝜋 .

  12. Volume Find the volume of the solid formed by revolving the region bounded by the graphs of 𝑦 =sin⁡𝑥 +sec⁡𝑥 , 𝑦 =0 , 𝑥 =0 , and 𝑥 =𝜋/3 about the 𝑥 -axis.

  13. Volume Find the volume of the solid formed by revolving the region bounded by the graphs of 𝑦 =arctan⁡𝑥 , 𝑥 =0 , and 𝑦 =𝜋/4 about the 𝑦 -axis.

  14. Average Value Find the average value of the function 𝑓(𝑥) =11−sin⁡𝜃 on [0,𝜋/6] .

8.4 Trigonometric Substitutions

教材插图

教材插图

教材插图

FIGURE 8.3 The arctangent, arcsine, and arcsecant of x/a, graphed as functions of x/a.

Trigonometric substitutions occur when we replace the variable of integration by a trigonometric function. The most common substitutions are 𝑥 =𝑎tan⁡𝜃 , 𝑥 =𝑎sin⁡𝜃 , and 𝑥 =𝑎sec⁡𝜃 . These substitutions are effective in transforming integrals involving √𝑎2+𝑥2 , √𝑎2−𝑥2 , and √𝑥2−𝑎2 into integrals with respect to 𝜃 , since they come from the reference right triangles in Figure 8.2.

教材插图

FIGURE 8.2 Reference triangles for the three basic substitutions, identifying the sides labeled x and a for each substitution.

With 𝑥 =𝑎tan⁡𝜃 ,

𝑎2+𝑥2=𝑎2+𝑎2tan2⁡𝜃=𝑎2(1+tan2⁡𝜃)=𝑎2sec2⁡𝜃.

With 𝑥 =𝑎sin⁡𝜃

𝑎2−𝑥2=𝑎2−𝑎2sin2⁡𝜃=𝑎2(1−sin2⁡𝜃)=𝑎2cos2⁡𝜃.

With 𝑥 =𝑎sec⁡𝜃

𝑥2−𝑎2=𝑎2sec2⁡𝜃−𝑎2=𝑎2(sec2⁡𝜃−1)=𝑎2tan2⁡𝜃.

We want any substitution we use in an integration to be reversible so that we can change back to the original variable afterward. For example, if 𝑥 =𝑎tan⁡𝜃 , we want to be able to set 𝜃 =arctan⁡(𝑥/𝑎) after the integration takes place. If 𝑥 =𝑎sin⁡𝜃 , we want to be able to set 𝜃 =arcsin⁡(𝑥/𝑎) when we’re done, and similarly for 𝑥 =𝑎sec⁡𝜃 .

As we know from Section 1.5, the functions in these substitutions have inverses only for selected values of 𝜃 (Figure 8.3). For reversibility,

𝑥=𝑎tan⁡𝜃 requires 𝜃=arctan⁡(𝑥𝑎) with −𝜋2<𝜃<𝜋2, 𝑥=𝑎sin⁡𝜃 requires 𝜃=arcsin⁡(𝑥𝑎) with −𝜋2≤𝜃≤𝜋2,

with

{0≤𝜃<𝜋2 if 𝑥𝑎≥1,𝜋2<𝜃≤𝜋 if 𝑥𝑎≤−1.

To simplify calculations with the substitution 𝑥 =𝑎sec⁡𝜃 , we will restrict its use to integrals in which 𝑥/𝑎 ≥1 . This will place 𝜃 in [0,𝜋/2) and make tan⁡𝜃 ≥0 . We will then have √𝑥2−𝑎2 =√𝑎2tan2⁡𝜃 =|𝑎tan⁡𝜃| =𝑎tan⁡𝜃 , free of absolute values, provided 𝑎 >0 .

Procedure for a Trigonometric Substitution

  1. Write down the substitution for x, calculate the differential dx, and specify the selected values of 𝜃 for the substitution.

  2. Substitute the trigonometric expression and the calculated differential into the integrand, and then simplify the results algebraically.

  3. Evaluate the trigonometric integral, keeping in mind the restrictions on the angle 𝜃 for reversibility.

  4. Draw an appropriate reference triangle to reverse the substitution in the integration result and convert it back to the original variable x.

EXAMPLE 1 Evaluate

∫𝑑𝑥√4+𝑥2.

教材插图

FIGURE 8.4 Reference triangle for 𝑥 =2tan⁡𝜃 (Example 1):

tan⁡𝜃=𝑥2

Solution We set

and

sec⁡𝜃=√4+𝑥22.

Then

𝑥=2tan⁡𝜃,𝑑𝑥=2sec2⁡𝜃𝑑𝜃,−𝜋2<𝜃<𝜋2,4+𝑥2=4+4tan2⁡𝜃=4(1+tan2⁡𝜃)=4sec2⁡𝜃. ∫𝑑𝑥√4+𝑥2=∫2sec2⁡𝜃𝑑𝜃√4sec2⁡𝜃=∫sec2⁡𝜃𝑑𝜃|sec⁡𝜃|√sec2⁡𝜃=|sec⁡𝜃|=∫sec⁡𝜃𝑑𝜃sec⁡𝜃>0 for −𝜋2<𝜃<𝜋2=ln⁡|sec⁡𝜃+tan⁡𝜃|+𝐶=ln⁡∣√4+𝑥22+𝑥2∣+𝐶. From Fig.8.4 

Notice how we expressed ln⁡|sec⁡𝜃+tan⁡𝜃| in terms of x: We drew a reference triangle for the original substitution 𝑥 =2tan⁡𝜃 (Figure 8.4) and read the ratios from the triangle.

EXAMPLE 2 Here we find an expression for the inverse hyperbolic sine function in terms of the natural logarithm. Following the same procedure as in Example 1, we find that

∫𝑑𝑥√𝑎2+𝑥2=∫sec⁡𝜃𝑑𝜃=ln⁡|sec⁡𝜃+tan⁡𝜃|+𝐶=ln⁡∣√𝑎2+𝑥2𝑎+𝑥𝑎∣+𝐶𝑥=𝑎tan⁡𝜃,𝑑𝑥=𝑎sec2⁡𝜃𝑑𝜃(Fig.8.2)

From Table 7.9, sinh−1⁡(𝑥/𝑎) is also an antiderivative of 1/√𝑎2+𝑥2 , so the two anti-derivatives differ by a constant, giving

sinh−1⁡𝑥𝑎=ln⁡∣√𝑎2+𝑥2𝑎+𝑥𝑎∣+𝐶.

Setting 𝑥 =0 in this last equation, we find that 0 =ln⁡|1| +𝐶 , so 𝐶 =0 . Since √𝑎2+𝑥2 >|𝑥| , we conclude that √𝑎2+𝑥2𝑎 +𝑥𝑎 >0 , and therefore

sinh−1⁡𝑥𝑎=ln⁡(√𝑎2+𝑥2𝑎+𝑥𝑎)

(See also Exercise 76 in Section 7.3.)

EXAMPLE 3 Evaluate

FIGURE 8.5 Reference triangle for 𝑥 =3sin⁡𝜃 (Example 3):

教材插图

sin⁡𝜃=𝑥3

and

∫𝑥2𝑑𝑥√9−𝑥2.

Solution We set

cos⁡𝜃=√9−𝑥23. 𝑥=3sin⁡𝜃,𝑑𝑥=3cos⁡𝜃𝑑𝜃,−𝜋2<𝜃<𝜋2

Then

9−𝑥2=9−9sin2⁡𝜃=9(1−sin2⁡𝜃)=9cos2⁡𝜃. ∫𝑥2𝑑𝑥√9−𝑥2=∫9sin2⁡𝜃⋅3cos⁡𝜃𝑑𝜃|3cos⁡𝜃|=9∫sin2⁡𝜃𝑑𝜃cos⁡𝜃>0 for −𝜋2<𝜃<𝜋2=9∫1−cos⁡2𝜃2𝑑𝜃sin2⁡𝜃=1−cos⁡2𝜃2=92(𝜃−sin⁡2𝜃2)+𝐶=92(𝜃−sin⁡𝜃cos⁡𝜃)+𝐶sin⁡2𝜃=2sin⁡𝜃cos⁡𝜃=92(arcsin⁡𝑥3−𝑥3⋅√9−𝑥23)+𝐶 From Fig.8.5 =92arcsin⁡𝑥3−𝑥2√9−𝑥2+𝐶.

EXAMPLE 4 Evaluate

∫𝑑𝑥√25𝑥2−4,𝑥>25.

Solution We first rewrite the radical as

√25𝑥2−4=√25(𝑥2−425)=5√𝑥2−(25)2√𝑥2−𝑎2 with 𝑎=25

to put the radicand in the form 𝑥2 −𝑎2 . We then substitute

𝑥=25sec⁡𝜃,𝑑𝑥=25sec⁡𝜃tan⁡𝜃𝑑𝜃,0<𝜃<𝜋2.

We then get

𝑥2−(25)2=425sec2⁡𝜃−425=425(sec2⁡𝜃−1)=425tan2⁡𝜃

教材插图

and

√𝑥2−(25)2=25|tan⁡𝜃|=25tan⁡𝜃.tan⁡𝜃>0 for 0<𝜃<𝜋/2

With these substitutions, we have

FIGURE 8.6 If 𝑥 =(2/5)sec⁡𝜃 ,

0 <𝜃 <𝜋/2 , then 𝜃 =arcsec⁡(5𝑥/2) ,

and we can read the values of the other trigonometric functions of 𝜃 from this right triangle (Example 4).

∫𝑑𝑥√25𝑥2−4 =∫𝑑𝑥5√𝑥2−(4/25) =∫(2/5)sec⁡𝜃tan⁡𝜃𝑑𝜃5⋅(2/5)tan⁡𝜃 =15∫sec⁡𝜃𝑑𝜃 =15ln⁡|sec⁡𝜃 +tan⁡𝜃| +𝐶 =15ln⁡∣5𝑥2+√25𝑥2−42∣ +𝐶. From Fig.8.6

EXERCISES 8.4

Using Trigonometric Substitutions

Evaluate the integrals in Exercises 1–14.

  1. ∫𝑑𝑥√9+𝑥2

  2. ∫3𝑑𝑥√1+9𝑥2

  3. ∫2−2𝑑𝑥4+𝑥2

  4. ∫20𝑑𝑥8+2𝑥2

  5. ∫3/20𝑑𝑥√9−𝑥2

  6. ∫1/2√202𝑑𝑥√1−4𝑥2

  7. ∫√25−𝑡2𝑑𝑡

  8. ∫√1−9𝑡2𝑑𝑡

  9. ∫𝑑𝑥√4𝑥2−49,𝑥 >72

  10. ∫5𝑑𝑥√25𝑥2−9,𝑥 >35

  11. ∫√𝑦2−49𝑦𝑑𝑦,𝑦 >7

  12. ∫√𝑦2−25𝑦3𝑑𝑦,𝑦 >5

  13. ∫𝑑𝑥𝑥2√𝑥2−1,𝑥 >1

  14. ∫2𝑑𝑥𝑥3√𝑥2−1,𝑥 >1

Assorted Integrations

Use any method to evaluate the integrals in Exercises 15–38. Most will require trigonometric substitutions, but some can be evaluated by other methods.

  1. ∫𝑑𝑥𝑥√𝑥2−1

  2. ∫𝑑𝑥1+𝑥2

  3. ∫𝑥𝑑𝑥√𝑥2−1

  4. ∫𝑑𝑥√1−𝑥2

  5. ∫𝑥√9−𝑥2𝑑𝑥

  6. ∫𝑥24+𝑥2𝑑𝑥

  7. ∫𝑥3𝑑𝑥√𝑥2+4

  8. ∫𝑑𝑥𝑥2√𝑥2+1

  9. ∫8𝑑𝑤𝑤2√4−𝑤2

  10. ∫√9−𝑤2𝑤2𝑑𝑤

  11. ∫√𝑥+11−𝑥𝑑𝑥

  12. ∫𝑥√𝑥2−4𝑑𝑥

  13. ∫√3/204𝑥2𝑑𝑥(1−𝑥2)3/2

  14. ∫10𝑑𝑥(4−𝑥2)3/2

  15. ∫𝑑𝑥(𝑥2−1)3/2,𝑥 >1

  16. ∫𝑥2𝑑𝑥(𝑥2−1)5/2,𝑥 >1

  17. ∫(1−𝑥2)3/2𝑥6𝑑𝑥

  18. ∫(1−𝑥2)1/2𝑥4𝑑𝑥

  19. ∫8𝑑𝑥(4𝑥2+1)2

  20. ∫6𝑑𝑡(9𝑡2+1)2

  21. ∫𝑥3𝑑𝑥𝑥2−1

  22. ∫𝑥𝑑𝑥25+4𝑥2

  23. ∫𝑣2𝑑𝑣(1−𝑣2)5/2

  24. ∫(1−𝑟2)5/2𝑟8𝑑𝑟

In Exercises 39–48, use an appropriate substitution and then a trigonometric substitution to evaluate the integrals.

  1. ∫ln⁡40𝑒𝑡𝑑𝑡√𝑒2𝑡+9

  2. ∫ln⁡(4/3)ln⁡(3/4)𝑒𝑡𝑑𝑡(1+𝑒2𝑡)3/2

  3. ∫1/41/122𝑑𝑡√𝑡+4𝑡√𝑡

  4. ∫𝑒1𝑑𝑦𝑦√1+(ln⁡𝑦)2

  5. ∫𝑥𝑑𝑥√1+𝑥4

  6. ∫√1−(ln⁡𝑥)2𝑥ln⁡𝑥𝑑𝑥

  7. ∫√4−𝑥𝑥𝑑𝑥

  8. ∫√𝑥1−𝑥3𝑑𝑥 (Hint: Let 𝑥 =𝑢2 .)

(Hint: Let 𝑢 =𝑥3/2 .)

  1. ∫√𝑥√1−𝑥𝑑𝑥

  2. ∫√𝑥−2√𝑥−1𝑑𝑥

Complete the Square Before Using Trigonometric Substitutions For Exercises 49–52, complete the square before using an appropriate trigonometric substitution. 49. ∫√8−2𝑥−𝑥2𝑑𝑥

  1. ∫1√𝑥2−2𝑥+5𝑑𝑥

  2. ∫√𝑥2+4𝑥+3𝑥+2𝑑𝑥

  3. ∫√𝑥2+2𝑥+2𝑥2+2𝑥+1𝑑𝑥

Initial Value Problems

Initial Value Problems Solve the initial value problems in Exercises 53–56 for y as a function of x. 53. 𝑥𝑑𝑦𝑑𝑥 =√𝑥2−4, 𝑥 ≥2, 𝑦(2) =0

  1. √𝑥2−9𝑑𝑦𝑑𝑥 =1, 𝑥 >3, 𝑦(5) =ln⁡3

  2. (𝑥2 +4)𝑑𝑦𝑑𝑥 =3, 𝑦(2) =0

  3. (𝑥2 +1)2𝑑𝑦𝑑𝑥 =√𝑥2+1, 𝑦(0) =1

Applications and Examples

  1. Area Find the area of the region in the first quadrant that is enclosed by the coordinate axes and the curve 𝑦 =√9−𝑥2/3 .

  2. Area Find the area enclosed by the ellipse 𝑥2𝑎2 +𝑦2𝑏2 =1.

  3. Consider the region bounded by the graphs of 𝑦 =sin−1⁡𝑥 , 𝑦 =0 , and 𝑥 =1/2 .
    a. Find the area of the region.
    b. Find the centroid of the region.

  4. Consider the region bounded by the graphs of 𝑦 =√𝑥arctan⁡𝑥 and y = 0 for 0 ≤𝑥 ≤1 . Find the volume of the solid formed by revolving this region about the x-axis (see accompanying figure).

教材插图

  1. Evaluate ∫𝑥3√1−𝑥2𝑑𝑥 using a. integration by parts. b. a 𝑢 -substitution. c. a trigonometric substitution.

  2. Path of a water skier Suppose that a boat is positioned at the origin with a water skier tethered to the boat at the point (10, 0) on a rope 10m long. As the boat travels along the positive 𝑦 -axis, the skier is pulled behind the boat along an unknown path 𝑦 =𝑓(𝑥) , as shown in the accompanying figure.
    a. Show that 𝑓′(𝑥) =−√100−𝑥2𝑥 .
    (Hint: Assume that the skier is always pointed directly at the boat and the rope is on a line tangent to the path 𝑦 =𝑓(𝑥) .)
    b. Solve the equation in part (a) for 𝑓(𝑥) , using 𝑓(10) =0 .

教材插图

  1. Find the average value of 𝑓(𝑥) =√𝑥+1√𝑥 on the interval [1, 3].
  2. Find the length of the curve 𝑦 =1 −𝑒−𝑥 , 0 ≤𝑥 ≤1 .

8.5 Integration of Rational Functions by Partial Fractions

This section shows how to express a rational function (a quotient of polynomials) as a sum of simpler fractions, called partial fractions, which are more easily integrated. For instance, the rational function (5𝑥 −3)/(𝑥2 −2𝑥 −3) can be rewritten as

5𝑥−3𝑥2−2𝑥−3=2𝑥+1+3𝑥−3.

You can verify this equation algebraically by placing the fractions on the right side over a common denominator (𝑥 +1)(𝑥 −3) . The skill acquired in writing rational functions as such a sum is useful in other settings as well (for instance, when using certain transform methods to solve differential equations). To integrate the rational function (5𝑥 −3)/(𝑥2 −2𝑥 −3) on the left side of the expression we are considering, we simply sum the integrals of the fractions on the right side:

∫5𝑥−3(𝑥+1)(𝑥−3)𝑑𝑥=∫2𝑥+1𝑑𝑥+∫3𝑥−3𝑑𝑥=2ln⁡|𝑥+1|+3ln⁡|𝑥−3|+𝐶.

The method for rewriting rational functions as a sum of simpler fractions is called the method of partial fractions. In the case of our example, it consists of finding constants A and B such that

5𝑥−3𝑥2−2𝑥−3=𝐴𝑥+1+𝐵𝑥−3.(1)

(Pretend for a moment that we do not know that A = 2 and B = 3 will work.) We call the fractions 𝐴/(𝑥 +1) and 𝐵/(𝑥 −3) partial fractions because their denominators are only part of the original denominator 𝑥2 −2𝑥 −3 . We call A and B undetermined coefficients until suitable values for them have been found.

To find 𝐴 and 𝐵 , we first clear Equation (1) of fractions and regroup in powers of 𝑥 , obtaining

5𝑥−3=𝐴(𝑥−3)+𝐵(𝑥+1)=(𝐴+𝐵)𝑥−3𝐴+𝐵.

This will be an identity in 𝑥 if and only if the coefficients of like powers of 𝑥 on the two sides are equal:

𝐴+𝐵=5,−3𝐴+𝐵=−3.

Solving these equations simultaneously gives A = 2 and B = 3.

General Description of the Method

Success in writing a rational function 𝑓(𝑥)/𝑔(𝑥) as a sum of partial fractions depends on three things:

  • The degree of 𝑓(𝑥) must be less than the degree of 𝑔(𝑥) . That is, the fraction must be proper. If it isn’t, divide 𝑓(𝑥) by 𝑔(𝑥) and work with the remainder term. Example 3 of this section illustrates such a case.

  • We must know the factors of 𝑔(𝑥) . In theory, any polynomial with real coefficients can be written as a product of real linear factors and real quadratic factors. In practice, the factors may be hard to find.

  • The values of the undetermined coefficients form a system of 𝑛 linear equations in 𝑛 unknowns. For large 𝑛 , solving such systems may require linear algebra methods (such as Gaussian Elimination).

Here is how we find the partial fractions of a proper fraction 𝑓(𝑥)/𝑔(𝑥) when the factors of g are known. A quadratic polynomial (or factor) is irreducible if it cannot be written as the product of two linear factors with real coefficients. That is, the polynomial has no real roots.

Method of Partial Fractions When 𝑓(𝑥)/𝑔(𝑥) Is Proper

  1. Let 𝑥 −𝑟 be a linear factor of 𝑔(𝑥) . Suppose that (𝑥 −𝑟)𝑚 is the highest power of 𝑥 −𝑟 that divides 𝑔(𝑥) . Then, to this factor, assign the sum of the 𝑚 partial fractions:
𝐴1(𝑥−𝑟)+𝐴2(𝑥−𝑟)2+⋯+𝐴𝑚(𝑥−𝑟)𝑚.

Do this for each distinct linear factor of 𝑔(𝑥) .

  1. Let 𝑥2 +𝑝𝑥 +𝑞 be an irreducible quadratic factor of 𝑔(𝑥) . In this case, 𝑥2 +𝑝𝑥 +𝑞 has no real roots. Suppose that (𝑥2 +𝑝𝑥 +𝑞)𝑛 is the highest power of this factor that divides 𝑔(𝑥) . Then, to this factor, assign the sum of the n partial fractions:
𝐵1𝑥+𝐶1(𝑥2+𝑝𝑥+𝑞)+𝐵2𝑥+𝐶2(𝑥2+𝑝𝑥+𝑞)2+⋯+𝐵𝑛𝑥+𝐶𝑛(𝑥2+𝑝𝑥+𝑞)𝑛.

Do this for each distinct quadratic factor of 𝑔(𝑥) .

  1. Set the original fraction 𝑓(𝑥)/𝑔(𝑥) equal to the sum of all these partial fractions. Clear the resulting equation of fractions.

  2. Find the values of the undetermined coefficients.

There are often multiple ways to find the values of the undetermined coefficients in Step 4. To find the values of the coefficients that satisfy Equation (1), we equated coefficients of like powers of x. In the next example, we instead will assign convenient values of x, leading to simple equations that we can solve for the undetermined coefficients.

EXAMPLE 1 Use partial fractions to evaluate

∫𝑥2+4𝑥+1(𝑥−1)(𝑥+1)(𝑥+3)𝑑𝑥.

Solution Note that each of the factors (𝑥 −1) , (𝑥 +1) , and (𝑥 +3) is raised only to the first power. Therefore, the partial fraction decomposition has the form

𝑥2+4𝑥+1(𝑥−1)(𝑥+1)(𝑥+3)=𝐴𝑥−1+𝐵𝑥+1+𝐶𝑥+3.

To find the values of the undetermined coefficients A, B, and C, we clear fractions and get

𝑥2+4𝑥+1=𝐴(𝑥+1)(𝑥+3)+𝐵(𝑥−1)(𝑥+3)+𝐶(𝑥−1)(𝑥+1).

On the right side, we notice that a factor (𝑥 −1) is present in all terms except for the one containing A. Therefore, letting x = 1 allows us to solve for A.

𝑥=1:12+4(1)+1=𝐴(2)(4)+𝐵(0)+𝐶(0) 6=8𝐴 𝐴=34

In a similar manner, we can let x equal -1 to find B or -3 to find C.

𝑥=−1:(−1)2+4(−1)+1=𝐴(0)+𝐵(−2)(2)+𝐶(0)−2=−4𝐵 𝐵=12 𝑥=−3:(−3)2+4(−3)+1=𝐴(0)+𝐵(0)+𝐶(−4)(−2)−2=8𝐶 𝐶=−14

Hence we have

∫𝑥2+4𝑥+1(𝑥−1)(𝑥+1)(𝑥+3)𝑑𝑥=∫[341𝑥−1+121𝑥+1−141𝑥+3]𝑑𝑥=34ln⁡|𝑥−1|+12ln⁡|𝑥+1|−14ln⁡|𝑥+3|+𝐾,

where K is the arbitrary constant of integration (we call it K here to avoid confusion with the undetermined coefficient we labeled as C).

You can solve for the undetermined coefficients (𝐴,𝐵,etc.) by equating coefficients of like powers of x or by assigning convenient values to x. You should choose the method that is most convenient for the problem at hand.

EXAMPLE 2 Use partial fractions to evaluate

∫6𝑥+7(𝑥+2)2𝑑𝑥.

Solution First we express the integrand as a sum of partial fractions with undetermined coefficients.

6𝑥+7(𝑥+2)2=𝐴𝑥+2+𝐵(𝑥+2)2 Two terms because (𝑥+2) is squared 6𝑥+7=𝐴(𝑥+2)+𝐵 Multiply both sides by (𝑥+2)2.=𝐴𝑥+(2𝐴+𝐵)

Equating coefficients of corresponding powers of x gives

𝐴=6 and 2𝐴+𝐵=12+𝐵=7, or 𝐴=6 and 𝐵=−5.

Therefore,

∫6𝑥+7(𝑥+2)2𝑑𝑥=∫(6𝑥+2−5(𝑥+2)2)𝑑𝑥=6∫𝑑𝑥𝑥+2−5∫(𝑥+2)−2𝑑𝑥=6ln⁡|𝑥+2|+5(𝑥+2)−1+𝐶.

The next example shows how to handle the case when 𝑓(𝑥)/𝑔(𝑥) is an improper fraction. It is a case where the degree of f is larger than the degree of g.

EXAMPLE 3 Use partial fractions to evaluate

∫2𝑥3−4𝑥2−𝑥−3𝑥2−2𝑥−3𝑑𝑥.

Solution First we divide the denominator into the numerator to get a polynomial plus a proper fraction.

𝑥2−2𝑥−3――――――――――)2𝑥3−4𝑥2−𝑥−32𝑥3−4𝑥2−6𝑥――――――――5𝑥−3

Then we write the improper fraction as a polynomial plus a proper fraction.

2𝑥3−4𝑥2−𝑥−3𝑥2−2𝑥−3=2𝑥+5𝑥−3𝑥2−2𝑥−3

We found the partial fraction decomposition of the fraction on the right in the opening example, so

∫2𝑥3−4𝑥2−𝑥−3𝑥2−2𝑥−3𝑑𝑥=∫2𝑥𝑑𝑥+∫5𝑥−3𝑥2−2𝑥−3𝑑𝑥=∫2𝑥𝑑𝑥+∫2𝑥+1𝑑𝑥+∫3𝑥−3𝑑𝑥=𝑥2+2ln⁡|𝑥+1|+3ln⁡|𝑥−3|+𝐶.

EXAMPLE 4 Use partial fractions to evaluate

∫−2𝑥+4(𝑥2+1)(𝑥−1)2𝑑𝑥.

Solution The denominator has an irreducible quadratic factor 𝑥2 +1 as well as a repeated linear factor (𝑥 −1)2 , so we write

−2𝑥+4(𝑥2+1)(𝑥−1)2=𝐴𝑥+𝐵𝑥2+1+𝐶𝑥−1+𝐷(𝑥−1)2.(2)

Clearing the equation of fractions gives

−2𝑥+4=(𝐴𝑥+𝐵)(𝑥−1)2+𝐶(𝑥−1)(𝑥2+1)+𝐷(𝑥2+1)=(𝐴+𝐶)𝑥3+(−2𝐴+𝐵−𝐶+𝐷)𝑥2+(𝐴−2𝐵+𝐶)𝑥+(𝐵−𝐶+𝐷).

Equating coefficients of like terms gives

 Coefficients of 𝑥3:0=𝐴+𝐶  Coefficients of 𝑥2:0=−2𝐴+𝐵−𝐶+𝐷  Coefficients of 𝑥1:−2=𝐴−2𝐵+𝐶  Coefficients of 𝑥0:4=𝐵−𝐶+𝐷

We solve these equations simultaneously to find the values of A, B, C, and D.

−4=−2𝐴,𝐴=2 Subtract fourth equation from second.  𝐶=−𝐴=−2 From the first equation  𝐵=(𝐴+𝐶+2)/2=1 From the third equation and 𝐶=−𝐴 𝐷=4−𝐵+𝐶=1. From the fourth equation 

We substitute these values into Equation (2), obtaining

−2𝑥+4(𝑥2+1)(𝑥−1)2=2𝑥+1𝑥2+1−2𝑥−1+1(𝑥−1)2.

Finally, using the expansion above, we can integrate:

∫−2𝑥+4(𝑥2+1)(𝑥−1)2𝑑𝑥=∫(2𝑥+1𝑥2+1−2𝑥−1+1(𝑥−1)2)𝑑𝑥=∫(2𝑥𝑥2+1+1𝑥2+1−2𝑥−1+1(𝑥−1)2)𝑑𝑥=ln⁡(𝑥2+1)+tan−1⁡𝑥−2ln⁡|𝑥−1|−1𝑥−1+𝐾.

We use the letter K instead of C to represent an arbitrary constant here because we have already used C to represent a variable in the partial fraction representation.

EXAMPLE 5 Use partial fractions to evaluate

∫𝑑𝑥𝑥(𝑥2+1)2.

Solution The form of the partial fraction decomposition is

1𝑥(𝑥2+1)2=𝐴𝑥+𝐵𝑥+𝐶𝑥2+1+𝐷𝑥+𝐸(𝑥2+1)2.

Multiplying by 𝑥(𝑥2 +1)2 , we have

1=𝐴(𝑥2+1)2+(𝐵𝑥+𝐶)𝑥(𝑥2+1)+(𝐷𝑥+𝐸)𝑥=𝐴(𝑥4+2𝑥2+1)+𝐵(𝑥4+𝑥2)+𝐶(𝑥3+𝑥)+𝐷𝑥2+𝐸𝑥=(𝐴+𝐵)𝑥4+𝐶𝑥3+(2𝐴+𝐵+𝐷)𝑥2+(𝐶+𝐸)𝑥+𝐴.

If we equate coefficients, we get the system

𝐴+𝐵=0,𝐶=0,2𝐴+𝐵+𝐷=0,𝐶+𝐸=0,𝐴=1.

Solving this system gives A = 1, B = -1, C = 0, D = -1, and E = 0. Thus,

HISTORICAL BIOGRAPHY Oliver Heaviside (1850–1925)

Heaviside studied electricity and languages on his own. He was able to simplify Maxwell’s 20 equations into the two we now call Maxwell’s equations. Heaviside’s contributions in mathematics are in the areas of vector algebra and vector calculus.

To know more, visit the companion Website.

∫𝑑𝑥𝑥(𝑥2+1)2=∫[1𝑥+−𝑥𝑥2+1+−𝑥(𝑥2+1)2]𝑑𝑥=∫𝑑𝑥𝑥−∫𝑥𝑑𝑥𝑥2+1−∫𝑥𝑑𝑥(𝑥2+1)2=∫𝑑𝑥𝑥−12∫𝑑𝑢𝑢−12∫𝑑𝑢𝑢2𝑢=𝑥2+1,=ln⁡|𝑥|−12ln⁡|𝑢|+12𝑢+𝐾=ln⁡|𝑥|−12ln⁡(𝑥2+1)+12(𝑥2+1)+𝐾=ln⁡|𝑥|√𝑥2+1+12(𝑥2+1)+𝐾.

Determining Coefficients by Differentiating

Another way to determine the constants that appear in partial fractions is to differentiate, as in the next example.

EXAMPLE 6 Find 𝐴,𝐵 , and 𝐶 in the equation

𝑥−1(𝑥+1)3=𝐴𝑥+1+𝐵(𝑥+1)2+𝐶(𝑥+1)3.

Solution We first clear fractions:

𝑥−1=𝐴(𝑥+1)2+𝐵(𝑥+1)+𝐶.

Substituting 𝑥 = −1 shows 𝐶 = −2 . We then differentiate both sides with respect to 𝑥 , obtaining

1=2𝐴(𝑥+1)+𝐵.

Substituting 𝑥 = −1 shows 𝐵 =1 . We differentiate again to get 0 =2𝐴 , which shows 𝐴 =0 . Hence,

𝑥−1(𝑥+1)3=1(𝑥+1)2−2(𝑥+1)3.

EXERCISES 8.5

Expanding Quotients into Partial Fractions

Expand the quotients in Exercises 1–8 by partial fractions.

  1. 5𝑥−13(𝑥−3)(𝑥−2)

  2. 5𝑥−7𝑥2−3𝑥+2

  3. 𝑥+4(𝑥+1)2

  4. 2𝑥+2𝑥2−2𝑥+1

  5. 𝑧+1𝑧2(𝑧−1)

  6. 𝑧𝑧3−𝑧2−6𝑧

  7. 𝑡2+8𝑡2−5𝑡+6

  8. 𝑡4+9𝑡4+9𝑡2

Nonrepeated Linear Factors In Exercises 9–16, express the integrand as a sum of partial fractions and evaluate the integrals. 9. ∫𝑑𝑥1−𝑥2

  1. ∫𝑑𝑥𝑥2+2𝑥

  2. ∫𝑥+4𝑥2+5𝑥−6𝑑𝑥

  3. ∫2𝑥+1𝑥2−7𝑥+12𝑑𝑥

  4. ∫84𝑦𝑑𝑦𝑦2−2𝑦−3

  5. ∫11/2𝑦+4𝑦2+𝑦𝑑𝑦

  6. ∫𝑑𝑡𝑡3+𝑡2−2𝑡

  7. ∫𝑥+32𝑥3−8𝑥𝑑𝑥

Repeated Linear Factors

In Exercises 17–20, express the integrand as a sum of partial fractions and evaluate the integrals.

∫10𝑥3𝑑𝑥𝑥2+2𝑥+1
  1. ∫0−1𝑥3𝑑𝑥𝑥2−2𝑥+1

  2. ∫𝑑𝑥(𝑥2−1)2

  3. ∫𝑥2𝑑𝑥(𝑥−1)(𝑥2+2𝑥+1)

In Exercises 21–32, express the integrand as a sum of partial fractions and evaluate the integrals.

  1. ∫10𝑑𝑥(𝑥+1)(𝑥2+1)

  2. ∫√313𝑡2+𝑡+4𝑡3+𝑡𝑑𝑡

  3. ∫𝑦2+2𝑦+1(𝑦2+1)2𝑑𝑦

  4. ∫8𝑥2+8𝑥+2(4𝑥2+1)2𝑑𝑥

  5. ∫2𝑠+2(𝑠2+1)(𝑠−1)3𝑑𝑠

  6. ∫𝑠4+81𝑠(𝑠2+9)2𝑑𝑠

  7. ∫𝑥2−𝑥+2𝑥3−1𝑑𝑥

  8. ∫1𝑥4+𝑥𝑑𝑥

  9. ∫𝑥2𝑥4−1𝑑𝑥

  10. ∫𝑥2+𝑥𝑥4−3𝑥2−4𝑑𝑥

  11. ∫2𝜃3+5𝜃2+8𝜃+4(𝜃2+2𝜃+2)2𝑑𝜃

  12. ∫𝜃4−4𝜃3+2𝜃2−3𝜃+1(𝜃2+1)3𝑑𝜃

Improper Fractions

In Exercises 33–38, perform long division on the integrand, write the proper fraction as a sum of partial fractions, and then evaluate the integral.

  1. ∫2𝑥3−2𝑥2+1𝑥2−𝑥𝑑𝑥

  2. ∫𝑥4𝑥2−1𝑑𝑥

  3. ∫9𝑥3−3𝑥+1𝑥3−𝑥2𝑑𝑥

  4. ∫16𝑥34𝑥2−4𝑥+1𝑑𝑥

  5. ∫𝑦4+𝑦2−1𝑦3+𝑦𝑑𝑦

  6. ∫2𝑦4𝑦3−𝑦2+𝑦−1𝑑𝑦

Evaluating Integrals

Evaluate the integrals in Exercises 39–54.

  1. ∫𝑒𝑡𝑑𝑡𝑒2𝑡+3𝑒𝑡+2

  2. ∫𝑒4𝑡+2𝑒2𝑡−𝑒𝑡𝑒2𝑡+1𝑑𝑡

  3. ∫cos⁡𝑦𝑑𝑦sin2⁡𝑦+sin⁡𝑦−6

  4. ∫sin⁡𝜃𝑑𝜃cos2⁡𝜃+cos⁡𝜃−2

  5. ∫(𝑥−2)2tan−1⁡(2𝑥)−12𝑥3−3𝑥(4𝑥2+1)(𝑥−2)2𝑑𝑥

  6. ∫(𝑥+1)2tan−1⁡(3𝑥)+9𝑥3+𝑥(9𝑥2+1)(𝑥+1)2𝑑𝑥

  7. ∫1𝑥3/2−√𝑥𝑑𝑥

  8. ∫1(𝑥1/3−1)√𝑥𝑑𝑥 (Hint: Let 𝑥 =𝑢6 .)

  9. ∫√𝑥+1𝑥𝑑𝑥

  10. ∫1𝑥√𝑥+9𝑑𝑥 (Hint: Let 𝑥 +1 =𝑢2 .)

  11. ∫1𝑥(𝑥4+1)𝑑𝑥 (Hint: Multiply by 𝑥3𝑥3 .)

  12. ∫1𝑥6(𝑥5+4)𝑑𝑥

  13. ∫1cos⁡2𝜃sin⁡𝜃𝑑𝜃

  14. ∫1cos⁡𝜃+sin⁡2𝜃𝑑𝜃

  15. ∫√1+√𝑥𝑥𝑑𝑥

  16. ∫√𝑥√2−√𝑥+√𝑥𝑑𝑥

Use any method to evaluate the integrals in Exercises 55–66. 55. ∫𝑥3−2𝑥2−3𝑥𝑥+2𝑑𝑥

  1. ∫𝑥+2𝑥3−2𝑥2−3𝑥𝑑𝑥

  2. ∫2𝑥−2−𝑥2𝑥+2−𝑥𝑑𝑥

  3. ∫2𝑥22𝑥+2𝑥−2𝑑𝑥

  4. ∫1𝑥4−1𝑑𝑥

  5. ∫𝑥4−1𝑥5−5𝑥+1𝑑𝑥

  6. ∫ln⁡𝑥+2𝑥(ln⁡𝑥+1)(ln⁡𝑥+3)𝑑𝑥

  7. ∫2𝑥(ln⁡𝑥−2)3𝑑𝑥

  8. ∫1√𝑥2−1𝑑𝑥

  9. ∫𝑥𝑥+√𝑥2+2𝑑𝑥

  10. ∫𝑥5√𝑥3+1𝑑𝑥

  11. ∫𝑥2√1−𝑥2𝑑𝑥

Initial Value Problems

Solve the initial value problems in Exercises 67–70 for x as a function of t.

  1. (𝑡2 −3𝑡 +2)𝑑𝑥𝑑𝑡 =1 (t>2), 𝑥(3) =0

  2. (3𝑡4 +4𝑡2 +1)𝑑𝑥𝑑𝑡 =2√3 , 𝑥(1) = −𝜋√3/4

  3. (𝑡2 +2𝑡)𝑑𝑥𝑑𝑡 =2𝑥 +2 (t,x>0), 𝑥(1) =1

  4. (𝑡 +1)𝑑𝑥𝑑𝑡 =𝑥2 +1 (t>-1), 𝑥(0) =0

Applications and Examples

In Exercises 71 and 72, find the volume of the solid generated by revolving the shaded region about the indicated axis.

  1. The 𝑥 -axis

教材插图

  1. The 𝑦 -axis

教材插图

  1. Find the length of the curve 𝑦 =ln⁡(1 −𝑥2) , 0 ≤𝑥 ≤12 .

  2. Evaluate ∫sec⁡𝜃𝑑𝜃 by a. multiplying by sec⁡𝜃+tan⁡𝜃sec⁡𝜃+tan⁡𝜃 and then using a 𝑢 -substitution b. writing the integral as ∫1cos⁡𝜃𝑑𝜃 . Then multiply by cos⁡𝜃cos⁡𝜃 , use a trigonometric identity and a 𝑢 -substitution, and finally integrate using partial fractions.

教材插图

T 75. Find, to two decimal places, the 𝑥 -coordinate of the centroid of the region in the first quadrant bounded by the 𝑥 -axis, the curve 𝑦 =arctan⁡𝑥 , and the line 𝑥 =√3 .

T 76. Find the 𝑥 -coordinate of the centroid of this region to two decimal places.

教材插图

  1. Social diffusion Sociologists sometimes use the phrase “social diffusion” to describe the way information spreads through a population. The information might be a rumor, a cultural fad, or news about a technical innovation. In a sufficiently large population, the number of people x who have the information is treated as a differentiable function of time t, and the rate of diffusion, dx/dt, is assumed to be proportional to the number of people who have the information times the number of people who do not. This leads to the equation
𝑑𝑥𝑑𝑡=𝑘𝑥(𝑁−𝑥),

where N is the number of people in the population.

Suppose t is in days, k = 1/250, and two people start a rumor at time t = 0 in a population of N = 1000 people.

a. Find x as a function of t.

b. When will half the population have heard the rumor? (This is when the rumor will be spreading the fastest.)

  1. Second-order chemical reactions Many chemical reactions are the result of the interaction of two molecules that undergo a change to produce a new product. The rate of the reaction typically depends on the concentrations of the two kinds of molecules. If 𝑎 is the amount of substance 𝐴 and 𝑏 is the amount of substance 𝐵 at time 𝑡 =0 , and if 𝑥 is the amount of product at time 𝑡 , then the rate of formation of 𝑥 may be given by the differential equation
𝑑𝑥𝑑𝑡=𝑘(𝑎−𝑥)(𝑏−𝑥),

or

1(𝑎−𝑥)(𝑏−𝑥)𝑑𝑥𝑑𝑡=𝑘,

where k is a constant for the reaction. Integrate both sides of this equation to obtain a relation between x and t (a) if a = b, and (b) if 𝑎 ≠𝑏 . Assume in each case that x = 0 when t = 0.

8.6 Integral Tables and Computer Algebra Systems

In this section we discuss how to use tables and computer algebra systems (CAS) to evaluate integrals.

Integral Tables

A Brief Table of Integrals is provided at the back of the text, after the index. (More extensive tables appear in compilations such as CRC Mathematical Tables, which contain thousands of integrals.) The integration formulas are stated in terms of constants a, b, c, m, n, and so on. These constants can usually assume any real value and need not be integers. Occasional limitations on their values are stated with the formulas. Formula 21 requires 𝑛 ≠ −1 , for example, and Formula 27 requires 𝑛 ≠ −2 .

The formulas also assume that the constants do not take on values that require dividing by zero or taking even roots of negative numbers. For example, Formula 24 assumes that 𝑎 ≠0 , and Formulas 29a and 29b cannot be used unless b is positive.

EXAMPLE 1 Find

∫𝑥(2𝑥+5)−1𝑑𝑥.

Solution We use Formula 24 at the back of the text (not 22, which requires 𝑛 ≠ −1 ):

∫𝑥(𝑎𝑥+𝑏)−1𝑑𝑥=𝑥𝑎−𝑏𝑎2ln⁡|𝑎𝑥+𝑏|+𝐶.

With 𝑎 =2 and 𝑏 =5 , we have

∫𝑥(2𝑥+5)−1𝑑𝑥=𝑥2−54ln⁡|2𝑥+5|+𝐶.

EXAMPLE 2 Find

∫𝑑𝑥𝑥√2𝑥−4.

Solution We use Formula 29b:

∫𝑑𝑥𝑥√𝑎𝑥−𝑏=2√𝑏arctan⁡√𝑎𝑥−𝑏𝑏+𝐶.

With 𝑎 =2 and 𝑏 =4 , we have

∫𝑑𝑥𝑥√2𝑥−4=2√4arctan⁡√2𝑥−44+𝐶=arctan⁡√𝑥−22+𝐶.

EXAMPLE 3 Find

∫𝑥arcsin⁡𝑥𝑑𝑥.

Solution We begin by using Formula 106:

∫𝑥𝑛arcsin⁡𝑎𝑥𝑑𝑥=𝑥𝑛+1𝑛+1arcsin⁡𝑎𝑥−𝑎𝑛+1∫𝑥𝑛+1𝑑𝑥√1−𝑎2𝑥2,𝑛≠−1.

With 𝑛 =1 and 𝑎 =1 , we have

∫𝑥arcsin⁡𝑥𝑑𝑥=𝑥22arcsin⁡𝑥−12∫𝑥2𝑑𝑥√1−𝑥2.

Next we use Formula 49 to find the integral on the right:

∫𝑥2√𝑎2−𝑥2𝑑𝑥=𝑎22arcsin⁡(𝑥𝑎)−12𝑥√𝑎2−𝑥2+𝐶.

With a = 1,

∫𝑥2𝑑𝑥√1−𝑥2=12arcsin⁡𝑥−12𝑥√1−𝑥2+𝐶.

The combined result is

∫𝑥arcsin⁡𝑥𝑑𝑥=𝑥22arcsin⁡𝑥−12(12arcsin⁡𝑥−12𝑥√1−𝑥2+𝐶)=(𝑥22−14)arcsin⁡𝑥+14𝑥√1−𝑥2+𝐶′.

Reduction Formulas

The time required for repeated integrations by parts can sometimes be shortened by applying reduction formulas like the following.

∫tan𝑛⁡𝑥𝑑𝑥=1𝑛−1tan𝑛−1⁡𝑥−∫tan𝑛−2⁡𝑥𝑑𝑥(1) ∫(ln⁡𝑥)𝑛𝑑𝑥=𝑥(ln⁡𝑥)𝑛−𝑛∫(ln⁡𝑥)𝑛−1𝑑𝑥(2) ∫sin𝑛⁡𝑥cos𝑚⁡𝑥𝑑𝑥=−sin𝑛−1⁡𝑥cos𝑚+1⁡𝑥𝑚+𝑛+𝑛−1𝑚+𝑛∫sin𝑛−2⁡𝑥cos𝑚⁡𝑥𝑑𝑥(𝑛≠−𝑚).(3)

By applying such a formula repeatedly, we can eventually express the original integral in terms of a power low enough to be evaluated directly. The next example illustrates this procedure.

EXAMPLE 4 Find

∫tan5⁡𝑥𝑑𝑥.

Solution We apply Equation (1) with n = 5 to get

∫tan5⁡𝑥𝑑𝑥=14tan4⁡𝑥−∫tan3⁡𝑥𝑑𝑥.

We then apply Equation (1) again, with 𝑛 =3 , to evaluate the remaining integral:

∫tan3⁡𝑥𝑑𝑥=12tan2⁡𝑥−∫tan⁡𝑥𝑑𝑥=12tan2⁡𝑥+ln⁡|cos⁡𝑥|+𝐶1.

The combined result is

∫tan5⁡𝑥𝑑𝑥=14tan4⁡𝑥−12tan2⁡𝑥−ln⁡|cos⁡𝑥|+𝐶.

As their form suggests, reduction formulas are derived using integration by parts. (See Example 5 in Section 8.3.)

Integration with a CAS

A powerful capability of computer algebra systems is their ability to integrate symbolically. This is performed with the integrate command specified by the particular system (for example, int in Maple, Integrate in Mathematica).

EXAMPLE 5 Suppose that you want to evaluate the indefinite integral of the function

𝑓(𝑥)=𝑥2√𝑎2+𝑥2.

Using Maple, you first define or name the function:

𝑓:=𝑥∧2∗sqrt⁡(𝑎∧2+𝑥∧2);

Then you use the integrate command on f, identifying the variable of integration:

int⁡(𝑓,𝑥);

Maple returns the answer

𝑥(𝑎2+𝑥2)3/24−𝑎2𝑥√𝑎2+𝑥28−𝑎4ln⁡(𝑥+√𝑎2+𝑥2)8

If you want to see whether the answer can be simplified, enter

 simplify (%);

Maple returns

−𝑎4ln⁡(𝑥+√𝑎2+𝑥2)8+𝑥√𝑎2+𝑥2(𝑎2+2𝑥2)8

If you want the definite integral for 0 ≤𝑥 ≤𝜋/2 , you can use the format

int⁡(𝑓,𝑥=0..Pi/2);

Maple will return the expression

𝑎4ln⁡(𝑎2)16+𝜋(𝜋2+4𝑎2)3/264−𝑎2𝜋√𝜋2+4𝑎232+𝑎4ln⁡(2)8−𝑎4ln⁡(𝜋+√𝜋2+4𝑎2)8

You can also find the definite integral for a particular value of the constant 𝑎 :

>𝑎:=1;>int⁡(𝑓,𝑥=0..1);

Maple returns the numerical answer

38√2+18ln⁡(√2−1).

EXAMPLE 6 Use a CAS to find

∫sin2⁡𝑥cos3⁡𝑥𝑑𝑥.

Solution With Maple, we have the entry

int⁡((sin∧⁡2)(𝑥)∗(cos∧⁡3)(𝑥),𝑥);

with the immediate return

−15sin⁡(𝑥)cos⁡(𝑥)4+115cos⁡(𝑥)2sin⁡(𝑥)+215sin⁡(𝑥).

Computer algebra systems vary in how they process integrations. We used Maple in Examples 5 and 6. Mathematica would have returned somewhat different results:

  1. In Example 5, given
𝐼𝑛[1]:= Integrate [𝑥∧2∗ Sqrt [𝑎∧2+𝑥∧2],𝑥]

Mathematica returns

𝑂𝑢𝑡[1]=18√𝑎2+𝑥2⎛⎜ ⎜ ⎜ ⎜⎝𝑎2𝑥+2𝑥3−𝑎3sinh−1⁡(𝑥𝑎)√1+𝑥2𝑎2⎞⎟ ⎟ ⎟ ⎟⎠

without having to simplify an intermediate result. The answer is different from, but equivalent to, Formula 36 in the integral tables.

  1. The Mathematica answer to the integral
𝐼𝑛[2]:= Integrate [ Sin [𝑥]∧2∗ Cos [𝑥]∧3,𝑥]

in Example 6 is

𝑂𝑢𝑡[2]=sin⁡[𝑥]8−148sin⁡[3𝑥]−180sin⁡[5𝑥]

differing from the Maple answer. Both answers are correct.

Although a CAS is very powerful and can aid us in solving difficult problems, each CAS has its own limitations. There are even situations where a CAS may further complicate a problem (in the sense of producing an answer that is extremely difficult to use or interpret). Note, too, that neither Maple nor Mathematica returns an arbitrary constant +C. On the other hand, a little mathematical thinking on your part may reduce the problem to one that is quite easy to handle. We provide an example in Exercise 67.

Nonelementary Integrals

Many functions have antiderivatives that cannot be expressed using the standard functions that we have encountered, such as polynomials, trigonometric functions, and exponential functions. Integrals of functions that do not have elementary antiderivatives are called nonelementary integrals. These integrals can sometimes be expressed with infinite series (Chapter 9) or approximated using numerical methods (Section 8.7). Examples of nonelementary integrals include the error function (which measures the probability of random errors)

erf⁡(𝑥)=2√𝜋∫𝑥0𝑒−𝑡2𝑑𝑡

and integrals such as

∫sin⁡𝑥2𝑑𝑥 and ∫√1+𝑥4𝑑𝑥

that arise in engineering and physics. These and a number of others, such as

∫𝑒𝑥𝑥𝑑𝑥,∫𝑒(𝑒𝑥)𝑑𝑥,∫1ln⁡𝑥𝑑𝑥,∫√1−𝑘2sin2⁡𝑥𝑑𝑥,0<𝑘<1,∫ln⁡(ln⁡𝑥)𝑑𝑥,∫sin⁡𝑥𝑥𝑑𝑥,0<𝑘<1,

look so easy they tempt us to try them just to see how they turn out. It can be proved, however, that there is no way to express any of these integrals as finite combinations of elementary functions. The same applies to integrals that can be changed into these by substitution. The functions in these integrals all have antiderivatives, as a consequence of the Fundamental Theorem of Calculus, Part 1, because they are continuous. However, none of the antiderivatives are elementary. The integrals you are asked to evaluate in this chapter have elementary antiderivatives.

EXERCISES 8.6

Using Integral Tables

Use the table of integrals at the back of the text to evaluate the integrals in Exercises 1–26.

  1. ∫𝑑𝑥𝑥√𝑥−3

  2. ∫𝑑𝑥𝑥√𝑥+4

  3. ∫𝑥𝑑𝑥√𝑥−2

  4. ∫𝑥𝑑𝑥(2𝑥+3)3/2

  5. ∫𝑥√2𝑥−3𝑑𝑥

  6. ∫𝑥(7𝑥 +5)3/2𝑑𝑥

  7. ∫√9−4𝑥𝑥2𝑑𝑥

  8. ∫𝑑𝑥𝑥2√4𝑥−9

  9. ∫𝑥√4𝑥−𝑥2𝑑𝑥

  10. ∫√𝑥−𝑥2𝑥𝑑𝑥

  11. ∫𝑑𝑥𝑥√7+𝑥2

  12. ∫𝑑𝑥𝑥√7−𝑥2

  13. ∫√4−𝑥2𝑥𝑑𝑥

  14. ∫√𝑥2−4𝑥𝑑𝑥

  15. ∫𝑒2𝑡cos⁡3𝑡𝑑𝑡

  16. ∫𝑒−3𝑡sin⁡4𝑡𝑑𝑡

  17. ∫𝑥arccos⁡𝑥𝑑𝑥

  18. ∫𝑥arctan⁡𝑥𝑑𝑥

  19. ∫𝑥2arctan⁡𝑥𝑑𝑥

  20. ∫tan−1⁡𝑥𝑥2𝑑𝑥

  21. ∫sin⁡3𝑥cos⁡2𝑥𝑑𝑥

  22. ∫sin⁡2𝑥cos⁡3𝑥𝑑𝑥

  23. ∫8sin⁡4𝑡sin⁡𝑡2𝑑𝑡

  24. ∫sin⁡𝑡3sin⁡𝑡6𝑑𝑡

  25. ∫cos⁡𝜃3cos⁡𝜃4𝑑𝜃

  26. ∫cos⁡𝜃2cos⁡7𝜃𝑑𝜃

Substitution and Integral Tables

In Exercises 27–40, use a substitution to change the integral into one you can find in the table. Then evaluate the integral. 27. ∫𝑥3+𝑥+1(𝑥2+1)2𝑑𝑥

  1. ∫𝑥2+6𝑥(𝑥2+3)2𝑑𝑥

  2. ∫arcsin⁡√𝑥𝑑𝑥

  3. ∫cos−1⁡√𝑥√𝑥𝑑𝑥

  4. ∫√𝑥√1−𝑥𝑑𝑥

  5. ∫√2−𝑥√𝑥𝑑𝑥

  6. ∫cot⁡𝑡√1−sin2⁡𝑡𝑑𝑡,0 <𝑡 <𝜋/2

  7. ∫𝑑𝑡tan⁡𝑡√4−sin2⁡𝑡

  8. ∫𝑑𝑦𝑦√3+(ln⁡𝑦)2

  9. ∫tan−1⁡√𝑦𝑑𝑦

  10. ∫1√𝑥2+2𝑥+5𝑑𝑥 (Hint: Complete the square.)

  11. ∫𝑥2√𝑥2−4𝑥+5𝑑𝑥

  12. ∫√5−4𝑥−𝑥2𝑑𝑥

  13. ∫𝑥2√2𝑥−𝑥2𝑑𝑥

Using Reduction Formulas Use reduction formulas to evaluate the integrals in Exercises 41–50. 41. ∫sin5⁡2𝑥𝑑𝑥

  1. ∫8cos4⁡2𝜋𝑡𝑑𝑡

  2. ∫sin2⁡2𝜃cos3⁡2𝜃𝑑𝜃

  3. ∫2sin2⁡𝑡sec4⁡𝑡𝑑𝑡

  4. ∫4tan3⁡2𝑥𝑑𝑥

  5. ∫8cot4⁡𝑡𝑑𝑡

  6. ∫2sec3⁡𝜋𝑥𝑑𝑥

  7. ∫3sec4⁡3𝑥𝑑𝑥

  8. ∫csc5⁡𝑥𝑑𝑥

  9. ∫16𝑥3(ln⁡𝑥)2𝑑𝑥

Evaluate the integrals in Exercises 51–56 by making a substitution (possibly trigonometric) and then applying a reduction formula.

  1. ∫𝑒𝑡sec3⁡(𝑒𝑡 −1)𝑑𝑡

  2. ∫csc3⁡√𝜃√𝜃𝑑𝜃

  3. ∫102√𝑥2+1𝑑𝑥

  4. ∫√3/20𝑑𝑦(1−𝑦2)5/2

  5. ∫21(𝑟2−1)3/2𝑟𝑑𝑟

  6. ∫1/√30𝑑𝑡(𝑡2+1)7/2

Applications

  1. Surface area Find the area of the surface generated by revolving the curve 𝑦 =√𝑥2+2 , 0 ≤𝑥 ≤√2 , about the x-axis.

  2. Arc length Find the length of the curve 𝑦 =𝑥2 , 0 ≤𝑥 ≤√3/2 .

  3. Centroid Find the centroid of the region cut from the first quadrant by the curve 𝑦 =1/√𝑥+1 and the line 𝑥 =3 .

  4. Moment about y-axis A thin plate of constant density 𝛿 =1 occupies the region enclosed by the curve 𝑦 =36/(2𝑥 +3) and the line x = 3 in the first quadrant. Find the moment of the plate about the y-axis.

  5. Use the integral table and a calculator to find, to two decimal places, the area of the surface generated by revolving the curve 𝑦 =𝑥2 , −1 ≤𝑥 ≤1 , about the 𝑥 -axis.

  6. Volume The head of your firm’s accounting department has asked you to find a formula she can use in a computer program to calculate the year-end inventory of gasoline in the company’s tanks. A typical tank is shaped like a right circular cylinder of radius 𝑟 and length 𝐿 , mounted horizontally, as shown in the accompanying figure. The data come to the accounting office as depth measurements taken with a vertical measuring stick marked in centimeters.

a. Show, in the notation of the figure, that the volume of gasoline that fills the tank to a depth 𝑑 is

𝑉=2𝐿∫−𝑟+𝑑−𝑟√𝑟2−𝑦2𝑑𝑦.

b. Evaluate the integral.

教材插图

  1. What is the largest value that
∫𝑏𝑎√𝑥−𝑥2𝑑𝑥

can have for any 𝑎 and 𝑏 ? Give reasons for your answer.

  1. What is the largest value that
∫𝑏𝑎𝑥√2𝑥−𝑥2𝑑𝑥

can have for any a and b? Give reasons for your answer.

COMPUTER EXPLORATIONS

In Exercises 65 and 66, use a CAS to perform the integrations.

  1. Evaluate the integrals

a. ∫𝑥ln⁡𝑥𝑑𝑥 b. ∫𝑥2ln⁡𝑥𝑑𝑥 c. ∫𝑥3ln⁡𝑥𝑑𝑥.

d. What pattern do you see? Predict the formula for ∫𝑥4ln⁡𝑥𝑑𝑥 and then see if you are correct by evaluating it with a CAS.

e. What is the formula for ∫𝑥𝑛ln⁡𝑥𝑑𝑥,𝑛 ≥1 ? Check your answer using a CAS.

  1. Evaluate the integrals

a. ∫ln⁡𝑥𝑥2𝑑𝑥

b. ∫ln⁡𝑥𝑥3𝑑𝑥

∫ln⁡𝑥𝑥4𝑑𝑥.

d. What pattern do you see? Predict the formula for

∫ln⁡𝑥𝑥5𝑑𝑥

and then see if you are correct by evaluating it with a CAS.

e. What is the formula for

∫ln⁡𝑥𝑥𝑛𝑑𝑥,𝑛≥2?

Check your answer using a CAS.

  1. a. Use a CAS to evaluate
∫𝜋/20sin𝑛⁡𝑥sin𝑛⁡𝑥+cos𝑛⁡𝑥𝑑𝑥,

where n is an arbitrary positive integer. Does your CAS find the result?

b. In succession, find the integral when 𝑛 =1,2,3,5 , and 7. Comment on the complexity of the results.

c. Now substitute 𝑥 =(𝜋/2) −𝑢 and add the new and old integrals. What is the value of

∫𝜋/20sin𝑛⁡𝑥sin𝑛⁡𝑥+cos𝑛⁡𝑥𝑑𝑥?

This exercise illustrates how a little mathematical ingenuity can sometimes solve a problem not immediately amenable to solution by a CAS.

8.7 Numerical Integration

The antiderivatives of some functions, like sin⁡(𝑥2) , 1/ln⁡𝑥 , and √1+𝑥4 , have no elementary formulas. When we cannot find a workable antiderivative for a function f that we have to integrate, we can partition the interval of integration, replace f by a closely fitting polynomial on each subinterval, integrate the polynomials, and add the results to approximate the definite integral of f. This procedure is an example of numerical integration. In this section we start by revisiting the Midpoint Rule, which we studied in Section 5.2. We then study two new methods, the Trapezoidal Rule and Simpson’s Rule. A key goal in our analysis is to control the possible error that is introduced when computing an approximation to an integral.

Approximating Integrals with the Midpoint Rule

In Section 5.2 we introduced the Midpoint Rule to approximate a definite integral over an interval [𝑎,𝑏] . The rule is based on subdividing [𝑎,𝑏] into n equal subintervals,

[𝑥0,𝑥1],[𝑥1,𝑥2],…,[𝑥𝑛−1,𝑥𝑛]

each of width

Δ𝑥=𝑏−𝑎𝑛.

We then approximate the integral using n rectangles, where the height of the kth rectangle is the value of f at the midpoint 𝑐𝑘 =(𝑥𝑘−1 +𝑥𝑘)/2 of the kth subinterval [𝑥𝑘−1,𝑥𝑘] .

Midpoint Rule for Approximating a Definite Integral

∫𝑏𝑎𝑓(𝑥)𝑑𝑥≈𝑛∑𝑘=1𝑓(𝑐𝑘)(𝑏−𝑎𝑛)=[𝑓(𝑐1)+𝑓(𝑐2)+⋯+𝑓(𝑐𝑛)](𝑏−𝑎𝑛)

with 𝑐𝑘 =𝑥𝑘−1+𝑥𝑘2 and 𝑥𝑘 =𝑎 +𝑘(𝑏−𝑎𝑛) .

Trapezoidal Approximations

The Trapezoidal Rule for the value of a definite integral is based on approximating the region between a curve and the 𝑥 -axis with trapezoids instead of rectangles, as in Figure 8.7. It is not necessary for the subdivision points 𝑥0,𝑥1,𝑥2,…,𝑥𝑛 in the figure to be evenly spaced, but the resulting formula is simpler if they are. We therefore assume that the length of each subinterval is

Δ𝑥=𝑏−𝑎𝑛.

The length Δ𝑥 =(𝑏 −𝑎)/𝑛 is called the step size or mesh size. The area of the trapezoid that lies above the ith subinterval is

Δ𝑥(𝑦𝑖−1+𝑦𝑖2)=Δ𝑥2(𝑦𝑖−1+𝑦𝑖),

教材插图

FIGURE 8.7 The Trapezoidal Rule approximates short stretches of the curve 𝑦 =𝑓(𝑥) with line segments. To approximate the integral of f from a to b, we add the areas of the trapezoids made by vertically joining the ends of the segments to the x-axis.

where 𝑦𝑖−1 =𝑓(𝑥𝑖−1) and 𝑦𝑖 =𝑓(𝑥𝑖) . (See Figure 8.7.) The area below the curve 𝑦 =𝑓(𝑥) and above the x-axis is then approximated by adding the areas of all the trapezoids:

𝑇=12(𝑦0+𝑦1)Δ𝑥+12(𝑦1+𝑦2)Δ𝑥+…+12(𝑦𝑛−2+𝑦𝑛−1)Δ𝑥+12(𝑦𝑛−1+𝑦𝑛)Δ𝑥=Δ𝑥(12𝑦0+𝑦1+𝑦2+⋯+𝑦𝑛−1+12𝑦𝑛)=Δ𝑥2(𝑦0+2𝑦1+2𝑦2+⋯+2𝑦𝑛−1+𝑦𝑛),

where

𝑦0=𝑓(𝑎),𝑦1=𝑓(𝑥1),…,𝑦𝑛−1=𝑓(𝑥𝑛−1),𝑦𝑛=𝑓(𝑏).

The Trapezoidal Rule says: Use T to estimate the integral of f from a to b.

教材插图

FIGURE 8.8 The trapezoidal approximation of the area under the graph of 𝑦 =𝑥2 from 𝑥 =1 to 𝑥 =2 is a slight overestimate (Example 1).

TABLE 8.2

x𝑦 =𝑥2
11
542516
643616
744916
24

教材插图

FIGURE 8.9 Simpson’s Rule approximates short stretches of the curve with parabolas.

The Trapezoidal Rule To approximate ∫𝑏𝑎𝑓(𝑥)𝑑𝑥 , use [ T = \frac{\Delta x}{2}(y_0 + 2y_1 + 2y_2 + \cdots + 2y_{n-1} + y_n). ] The y’s are the values of f at the partition points 𝑥0 =𝑎,𝑥1 =𝑎 +Δ𝑥,𝑥2 =𝑎 +2Δ𝑥,…,𝑥𝑛−1 =𝑎 +(𝑛 −1)Δ𝑥,𝑥𝑛 =𝑏, where Δ𝑥 =(𝑏 −𝑎)/𝑛.

EXAMPLE 1 Use the Trapezoidal Rule with n = 4 to estimate ∫21𝑥2𝑑𝑥 . Compare the estimate with the exact value.

Solution Partition [1, 2] into four subintervals of equal length (Figure 8.8). Then evaluate 𝑦 =𝑥2 at each partition point (Table 8.2).

Using these 𝑦 -values, 𝑛 =4 , and Δ𝑥 =(2 −1)/4 =1/4 in the Trapezoidal Rule, we have

𝑇=Δ𝑥2(𝑦0+2𝑦1+2𝑦2+2𝑦3+𝑦4)=18(1+2(2516)+2(3616)+2(4916)+4)=7532=2.34375.

Since the parabola is concave up, the approximating segments lie above the curve, giving each trapezoid slightly more area than the corresponding strip under the curve. The exact value of the integral is

∫21𝑥2𝑑𝑥=𝑥33]21=83−13=73.

The 𝑇 approximation overestimates the integral by about half a percent of its true value of 7/3. The percentage error is (2.34375 −7/3)/(7/3) ≈0.00446 , or 0.446% .

Simpson’s Rule: Approximations Using Parabolas

Another rule for approximating the definite integral of a continuous function results from using parabolas instead of the straight-line segments that produced trapezoids. As before, we partition the interval [𝑎,𝑏] into n subintervals of equal length ℎ =Δ𝑥 =(𝑏 −𝑎)/𝑛 , but this time we require that n be an even number. On each consecutive pair of intervals we approximate the curve 𝑦 =𝑓(𝑥) ≥0 by a parabola, as shown in Figure 8.9. A typical parabola passes through three consecutive points (𝑥𝑖−1,𝑦𝑖−1) , (𝑥𝑖,𝑦𝑖) , and (𝑥𝑖+1,𝑦𝑖+1) on the curve.

Let’s calculate the shaded area beneath a parabola passing through three consecutive points. To simplify our calculations, we first take the case where 𝑥0 = −ℎ , 𝑥1 =0 , and 𝑥2 =ℎ (Figure 8.10), where ℎ =Δ𝑥 =(𝑏 −𝑎)/𝑛 . The area under the parabola will be the same if we shift the 𝑦 -axis to the left or right. The parabola has an equation of the form

𝑦=𝐴𝑥2+𝐵𝑥+𝐶,

教材插图

FIGURE 8.10 By integrating from -h to h, we find the shaded area to be

ℎ3(𝑦0+4𝑦1+𝑦2).

so the area under it from x = -h to x = h is

𝐴𝑝=∫ℎ−ℎ(𝐴𝑥2+𝐵𝑥+𝐶)𝑑𝑥=[𝐴𝑥33+𝐵𝑥22+𝐶𝑥]ℎ−ℎ=2𝐴ℎ33+2𝐶ℎ=ℎ3(2𝐴ℎ2+6𝐶).

Since the curve passes through the three points ( −ℎ,𝑦0),(0,𝑦1) , and (ℎ,𝑦2) , we also have

𝑦0=𝐴ℎ2−𝐵ℎ+𝐶,𝑦1=𝐶,𝑦2=𝐴ℎ2+𝐵ℎ+𝐶.

After some algebraic manipulation, we find that

𝐴𝑝=ℎ3(𝑦0+4𝑦1+𝑦2).

Now shifting the parabola horizontally to its shaded position in Figure 8.9 does not change the area under it. Thus the area under the parabola through (𝑥0,𝑦0) , (𝑥1,𝑦1) , and (𝑥2,𝑦2) in Figure 8.9 is still

ℎ3(𝑦0+4𝑦1+𝑦2).

Similarly, the area under the parabola through the points (𝑥2,𝑦2) , (𝑥3,𝑦3) , and (𝑥4,𝑦4) is

ℎ3(𝑦2+4𝑦3+𝑦4).

Computing the areas under all the parabolas and adding the results give the approximation

HISTORICAL BIOGRAPHY

To know more, visit the companion Website.

Thomas Simpson (1720–1761)

Simpson was a successful text writer and did most of his research on probability. Simpson’s rule to approximate definite integrals was developed before he was born. It is another of history’s beautiful quirks that one of the ablest mathematicians of the 18th century is remembered not for his own work but for a rule that was never his, that he never claimed, and that bears his name only because he happened to mention it in one of his books.

∫𝑏𝑎𝑓(𝑥)𝑑𝑥≈ℎ3(𝑦0+4𝑦1+𝑦2)+ℎ3(𝑦2+4𝑦3+𝑦4)+…+ℎ3(𝑦𝑛−2+4𝑦𝑛−1+𝑦𝑛)=ℎ3(𝑦0+4𝑦1+2𝑦2+4𝑦3+2𝑦4+⋯+2𝑦𝑛−2+4𝑦𝑛−1+𝑦𝑛).

The result is known as Simpson’s Rule. The function need not be positive, as in our derivation, but the number 𝑛 of subintervals must be even for us to apply the rule because each parabolic arc uses two subintervals.

Simpson’s Rule

To approximate ∫𝑏𝑎𝑓(𝑥)𝑑𝑥 , use

𝑆=Δ𝑥3(𝑦0+4𝑦1+2𝑦2+4𝑦3+⋯+2𝑦𝑛−2+4𝑦𝑛−1+𝑦𝑛).

The 𝑦 ‘s are the values of 𝑓 at the partition points

𝑥0=𝑎,𝑥1=𝑎+Δ𝑥,𝑥2=𝑎+2Δ𝑥,…,𝑥𝑛−1=𝑎+(𝑛−1)Δ𝑥,𝑥𝑛=𝑏.

The number 𝑛 is even, and Δ𝑥 =(𝑏 −𝑎)/𝑛 .

Note the pattern of the coefficients in the above rule: 1, 4, 2, 4, 2, 4, 2, …, 4, 1.

EXAMPLE 2 Use Simpson’s Rule with 𝑛 =4 to approximate ∫205𝑥4𝑑𝑥 .

TABLE 8.3

x𝑦 =5𝑥4
00
12516
15
3240516
280

Solution Partition [0, 2] into four subintervals and evaluate 𝑦 =5𝑥4 at the partition points (Table 8.3). Then apply Simpson’s Rule with 𝑛 =4 and Δ𝑥 =1/2 :

𝑆=Δ𝑥3(𝑦0+4𝑦1+2𝑦2+4𝑦3+𝑦4)=16(0+4(516)+2(5)+4(40516)+80)=32112.

This estimate differs from the exact value (32) by only 1/12, a percentage error of less than three-tenths of one percent, and this was with just four subintervals.

Error Analysis

Whenever we use an approximation technique, we must consider how accurate the approximation might be. The following theorem gives formulas for estimating the errors when using the Midpoint Rule, the Trapezoidal Rule, and Simpson’s Rule. The error is the difference between the approximation obtained by using the rule and the actual value of the definite integral ∫𝑏𝑎𝑓(𝑥)𝑑𝑥 .

THEOREM 1—Error Estimates in the Midpoint, Trapezoidal, and Simpson’s Rules

If 𝑓″ is continuous and M is any upper bound for the values of |𝑓″| on [a, b], then the error 𝐸𝑀 in the Midpoint Rule approximation of the integral of f from a to b for n steps satisfies the inequality

|𝐸𝑀|≤𝑀(𝑏−𝑎)324𝑛2. Midpoint Rule 

If 𝑓″ is continuous and M is any upper bound for the values of |𝑓″| on [a, b], then the error 𝐸𝑇 in the Trapezoidal Rule approximation of the integral of f from a to b for n steps satisfies the inequality

|𝐸𝑇|≤𝑀(𝑏−𝑎)312𝑛2. Trapezoidal Rule 

If 𝑓(4) is continuous and 𝑀 is any upper bound for the values of |𝑓(4)| on [𝑎,𝑏] , then the error 𝐸𝑆 in the Simpson’s Rule approximation of the integral of 𝑓 from 𝑎 to 𝑏 for 𝑛 steps satisfies the inequality

|𝐸𝑆|≤𝑀(𝑏−𝑎)5180𝑛4. Simpson's Rule 

To give an idea of why Theorem 1 is true in the case of the Trapezoidal Rule, we begin with a result which says that if 𝑓″ is continuous on the interval [𝑎,𝑏] , then

∫𝑏𝑎𝑓(𝑥)𝑑𝑥=𝑇−𝑏−𝑎12⋅𝑓′′(𝑐)(Δ𝑥)2

for some number 𝑐 between 𝑎 and 𝑏 . This result follows from Taylor’s Remainder Theorem, which we discuss in Theorem 24 of Section 9.9. It follows that as Δ𝑥 approaches zero, the error defined by

𝐸𝑇=−𝑏−𝑎12⋅𝑓′′(𝑐)(Δ𝑥)2

approaches zero at the rate of the square of Δ𝑥 .

The inequality

|𝐸𝑇|≤𝑏−𝑎12max|𝑓′′(𝑥)|(Δ𝑥)2,

where “max” refers to the maximum of |𝑓″(𝑥)| over the interval [𝑎,𝑏] , gives an upper bound for the magnitude of the error. In practice, we usually cannot find the exact value of max|𝑓″(𝑥)| and have to estimate an upper bound or “worst case” value for it instead. If M is any upper bound for the values of |𝑓″(𝑥)| on [𝑎,𝑏] , so that |𝑓″(𝑥)| ≤𝑀 for every x in [𝑎,𝑏] , then

|𝐸𝑇|≤𝑏−𝑎12𝑀(Δ𝑥)2.

If we substitute (𝑏 −𝑎)/𝑛 for Δ𝑥 , we get

|𝐸𝑇|≤𝑀(𝑏−𝑎)312𝑛2.

The upper bound for the error in the Midpoint Rule,

|𝐸𝑀|≤𝑀(𝑏−𝑎)324𝑛2,

is based on a similar argument. Taylor’s Remainder Theorem implies that the Midpoint Rule approximation and the definite integral differ by at most

|𝑓′′(𝑐)|(𝑏−𝑎)324𝑛2

for some 𝑐 ∈[𝑎,𝑏] . Taking 𝑀 to be at least as large as the maximum of |𝑓″(𝑐)| on [𝑎,𝑏] , we obtain the error bound in Theorem 1. Note that this is half as large as the error bound for the Trapezoidal Rule. This does not mean that the Midpoint Rule is always more accurate than the Trapezoidal Rule, but rather that the largest possible error that might occur is only half as large as the largest possible error for the Trapezoidal Rule.

To estimate the error in Simpson’s Rule, we start with a result, again following from Taylor’s Remainder Theorem, that says that if the fourth derivative 𝑓(4) is continuous, then

∫𝑏𝑎𝑓(𝑥)𝑑𝑥=𝑆−𝑏−𝑎180⋅𝑓(4)(𝑐)(Δ𝑥)4

for some point c between a and b. Thus, as Δ𝑥 approaches zero, the error,

𝐸𝑆=−𝑏−𝑎180⋅𝑓(4)(𝑐)(Δ𝑥)4,

approaches zero as the fourth power of Δ𝑥 . (This helps to explain why Simpson’s Rule is likely to give better results than the Trapezoidal Rule.)

The inequality

|𝐸𝑆|≤𝑏−𝑎180max∣𝑓(4)(𝑥)∣(Δ𝑥)4,

where “max” refers to the maximum of |𝑓(4)(𝑥)| over the interval [a, b], gives an upper bound for the magnitude of the error. As with max|𝑓″| in the error formula for the Trapezoidal Rule, we usually cannot find the exact value of max|𝑓(4)(𝑥)| and have to replace it with an upper bound. If M is any upper bound for the values of |𝑓(4)(𝑥)| on [a, b], then

|𝐸𝑆|≤𝑏−𝑎180𝑀(Δ𝑥)4.

Substituting (𝑏 −𝑎)/𝑛 for Δ𝑥 in this last expression gives

|𝐸𝑠|≤𝑀(𝑏−𝑎)5180𝑛4.

You might wonder why we don’t just take 𝑀 to be the maximum value of |𝑓″(𝑥)| on [𝑎,𝑏] for the first two rules, or the maximum value of |𝑓(4)(𝑥)| for Simpson’s Rule. The reason is that sometimes this maximum is hard to compute, while a less accurate upper bound is easily found. We can certainly set 𝑀 to the maximum value if we can compute it.

EXAMPLE 3 Find an upper bound for the error in estimating ∫205𝑥4𝑑𝑥 using Simpson’s Rule with 𝑛 =4 (Example 2).

Solution To estimate the error, we first find an upper bound M for the magnitude of the fourth derivative of 𝑓(𝑥) =5𝑥4 on the interval 0 ≤𝑥 ≤2 . Since the fourth derivative has the constant value 𝑓(4)(𝑥) =120 , we take M = 120. With b - a = 2 and n = 4, the error estimate for Simpson’s Rule gives

|𝐸𝑠|≤𝑀(𝑏−𝑎)5180𝑛4=120(2)5180⋅44=112.

This estimate is consistent with the result of Example 2.

Theorem 1 can also be used to estimate the number of subintervals required when using the Trapezoidal or Simpson’s Rule if we specify a certain tolerance for the error.

EXAMPLE 4 Estimate the minimum number of subintervals needed to approximate the integral in Example 3 using Simpson’s Rule with an error of magnitude less than 10−4 .

Solution Using the inequality in Theorem 1, if we choose the number of subintervals 𝑛 to satisfy

𝑀(𝑏−𝑎)5180𝑛4<10−4,

then the error 𝐸𝑆 in Simpson’s Rule satisfies |𝐸𝑆| <10−4 , as required.

From the solution in Example 3, we have 𝑀 =120 and 𝑏 −𝑎 =2 , so we want 𝑛 to satisfy

120(2)5180𝑛4<1104,

or, equivalently,

𝑛4>64⋅1043.

It follows that

𝑛>10(643)1/4≈21.5.

Since 𝑛 must be even in Simpson’s Rule, we estimate the minimum number of subintervals required for the error tolerance to be 𝑛 =22 .

EXAMPLE 5 As we saw in Chapter 7, the value of ln⁡2 can be calculated from the integral

ln⁡2=∫211𝑥𝑑𝑥.

Table 8.4 shows values of 𝑇 and 𝑆 for approximations of ∫21(1/𝑥)𝑑𝑥 using various values of 𝑛 . Notice how Simpson’s Rule dramatically improves over the Trapezoidal Rule.

TABLE 8.4 Trapezoidal Rule approximations (𝑇𝑛) and Simpson’s Rule approximations (𝑆𝑛) of ln⁡2 =∫21(1/𝑥)𝑑𝑥

n𝑇𝑛|Error| less than...𝑆𝑛|Error| less than...
100.69377140320.00062422270.69315023070.0000030502
200.69330338180.00015620130.69314737470.0000001942
300.69321661540.00006943490.69314721900.0000000385
400.69318624000.00003905950.69314719270.0000000122
500.69317217930.00002499880.69314718560.0000000050
1000.69315343050.00000625000.69314718090.0000000004

In particular, notice that when we double the value of n (thereby halving the value of ℎ =Δ𝑥 ), the T error is divided by 2 squared, whereas the S error is divided by 2 to the fourth.

This has a dramatic effect as Δ𝑥 =(2 −1)/𝑛 gets very small. The Simpson approximation for n = 50 rounds accurately to seven places and for n = 100 is accurate to nine decimal places (billionths)!

If 𝑓(𝑥) is a polynomial of degree less than 4, then its fourth derivative is zero, and

𝐸𝑆=−𝑏−𝑎180𝑓(4)(𝑐)(Δ𝑥)4=−𝑏−𝑎180(0)(Δ𝑥)4=0.

Thus, there will be no error in the Simpson approximation of any integral of 𝑓 . In other words, if 𝑓 is a constant, a linear function, or a quadratic or cubic polynomial, Simpson’s Rule will give the value of any integral of 𝑓 exactly, whatever the number of subdivisions. Similarly, if 𝑓 is a constant or a linear function, then its second derivative is zero, and

𝐸𝑇=−𝑏−𝑎12𝑓′′(𝑐)(Δ𝑥)2=−𝑏−𝑎12(0)(Δ𝑥)2=0.

The Trapezoidal Rule will therefore give the exact value of any integral of 𝑓 . This is no surprise, for the trapezoids fit the graph perfectly.

Although decreasing the step size Δ𝑥 reduces the error in the Simpson and Trapezoidal approximations in theory, it may fail to do so in practice. When Δ𝑥 is very small, say Δ𝑥 =10−8 , computer or calculator round-off errors in the arithmetic required to evaluate S and T may accumulate to such an extent that the error formulas no longer describe what is going on. Shrinking Δ𝑥 below a certain size can actually make things worse. You should consult a text on numerical analysis for more sophisticated methods if you are having problems with round-off error using the rules discussed in this section.

教材插图

FIGURE 8.11 The dimensions of the swamp in Example 6.

EXAMPLE 6 A town wants to drain and fill a polluted swamp (Figure 8.11). The swamp averages 1.5 m deep. About how many cubic meters of dirt will it take to fill the area after the swamp is drained?

Solution To calculate the volume of the swamp, we estimate the surface area and multiply by 1.5. To estimate the area, we use Simpson’s Rule with Δ𝑥 =6m , and the 𝑦 s equal to the distances measured across the swamp, as shown in Figure 8.11.

𝑆=Δ𝑥3(𝑦0+4𝑦1+2𝑦2+4𝑦3+2𝑦4+4𝑦5+𝑦6)=63(44+148+46+64+24+36+4)=732

The volume is about (732)(1.5) = 1098 m 4 .

EXERCISES 8.7

For some exercises, a calculator may be helpful for expressing answers in decimal form.

Estimating Definite Integrals

The instructions for the integrals in Exercises 1–10 have three parts, one for the Midpoint Rule, one for the Trapezoidal Rule, and one for Simpson’s Rule.

I. Using the Midpoint Rule

a. Estimate the integral with 𝑛 =4 steps and find an upper bound for |𝐸𝑀| .

b. Evaluate the integral directly and find |𝐸𝑀| .

c. Use the formula (|𝐸𝑀|/(true value)) ×100 to express |𝐸𝑀| as a percentage of the integral’s true value.

II. Using the Trapezoidal Rule

a. Estimate the integral with 𝑛 =4 steps and find an upper bound for |𝐸𝑇| .

b. Evaluate the integral directly and find |𝐸𝑇| .

c. Use the formula (|𝐸𝑇|/(true value)) ×100 to express |𝐸𝑇| as a percentage of the integral’s true value.

III. Using Simpson’s Rule

a. Estimate the integral with 𝑛 =4 steps and find an upper bound for |𝐸𝑆| .

b. Evaluate the integral directly and find |𝐸𝑆| .

c. Use the formula (|𝐸𝑆|/(true value)) ×100 to express |𝐸𝑆| as a percentage of the integral’s true value.

  1. ∫21𝑥𝑑𝑥

  2. ∫31(2𝑥 −1)𝑑𝑥

  3. ∫1−1(𝑥2 +1)𝑑𝑥

  4. ∫0−2(𝑥2 −1)𝑑𝑥

∫211𝑠2𝑑𝑠
  1. ∫20(𝑡3 +𝑡)𝑑𝑡

  2. ∫1−1(𝑡3 +1)𝑑𝑡

  3. ∫421(𝑠−1)2𝑑𝑠

∫𝜋0sin⁡𝑡𝑑𝑡 ∫10sin⁡𝜋𝑡𝑑𝑡

Estimating the Number of Subintervals

In Exercises 11–22, estimate the minimum number of subintervals needed to approximate the integrals with an error of magnitude less than 10−4 by (a) the Trapezoidal Rule and (b) Simpson’s Rule. (The integrals in Exercises 11–18 are the integrals from Exercises 1–8.)

  1. ∫21𝑥𝑑𝑥

  2. ∫31(2𝑥 −1)𝑑𝑥

  3. ∫1−1(𝑥2 +1)𝑑𝑥

  4. ∫0−2(𝑥2 −1)𝑑𝑥

  5. ∫20(𝑡3 +𝑡)𝑑𝑡

  6. ∫1−1(𝑡3 +1)𝑑𝑡

  7. ∫211𝑠2𝑑𝑠

  8. ∫421(𝑠−1)2𝑑𝑠

  9. ∫30√𝑥+1𝑑𝑥

  10. ∫301√𝑥+1𝑑𝑥

  11. ∫20sin⁡(𝑥 +1)𝑑𝑥

  12. ∫1−1cos⁡(𝑥 +𝜋)𝑑𝑥

Estimates with Numerical Data

  1. Volume of water in a swimming pool A rectangular swimming pool is 5 m wide and 10 m long. The accompanying table shows the depth ℎ(𝑥) of the water at 1-m intervals from one end of the pool to the other. Estimate the volume of water in the pool using the Trapezoidal Rule with n = 10 applied to the integral
𝑉=∫1005⋅ℎ(𝑥)𝑑𝑥.
Position (m) xDepth (m) h(x)Position (m) xDepth (m) h(x)
01.2062.30
11.6472.38
21.8282.46
31.9892.54
42.10102.60
52.20
  1. Distance traveled The accompanying table shows time-to-speed data for a car accelerating from rest to 130 km/h. How far had the car traveled by the time it reached this speed? (Use trapezoids to estimate the area under the velocity curve, but be careful: The time intervals vary in length.)
Speed changeTime (s)
Zero to 30 km/h2.2
40 km/h3.2
50 km/h4.5
60 km/h5.9
70 km/h7.8
80 km/h10.2
90 km/h12.7
100 km/h16.0
110 km/h20.6
120 km/h26.2
130 km/h37.1
  1. Wing design The design of a new airplane requires a gasoline tank of constant cross-sectional area in each wing. A scale drawing of a cross-section is shown here. The tank must hold 2000 kg of gasoline, which has a density of 673kg/m3 . Estimate the length of the tank by Simpson’s Rule.

教材插图

𝑦3=0.65r 𝑦4=0.7 𝑦5=𝑦6=0.75
  1. Oil consumption on Pathfinder Island A diesel generator runs continuously, consuming oil at a gradually increasing rate until it must be temporarily shut down to have the filters replaced. Use the Trapezoidal Rule to estimate the amount of oil consumed by the generator during that week.
DayOil consumption rate (liters/hour)
Sun0.019
Mon0.020
Tue0.021
Wed0.023
Thu0.025
Fri0.028
Sat0.031
Sun0.035

Theory and Examples

  1. Usable values of the sine-integral function The sine-integral function,
Si⁡(𝑥)=∫𝑥0sin⁡𝑡𝑡𝑑𝑡, "Sine integral of 𝑥 " 

is one of the many functions in engineering whose formulas cannot be simplified. There is no elementary formula for the anti-derivative of (sin⁡𝑡)/𝑡 . The values of Si(𝑥) , however, are readily estimated by numerical integration.

Although the notation does not show it explicitly, the function being integrated is

𝑓(𝑡)={sin⁡𝑡𝑡,𝑡≠01,𝑡=0,

the continuous extension of (sin⁡𝑡)/𝑡 to the interval [0,𝑥] . The function has derivatives of all orders at every point of its domain. Its graph is smooth, and you can expect good results from Simpson’s Rule.

教材插图

a. Use the fact that |𝑓(4)| ≤1 on [0,𝜋/2] to give an upper bound for the error that will occur if

Si⁡(𝜋2)=∫𝜋/20sin⁡𝑡𝑡𝑑𝑡

is estimated by Simpson’s Rule with 𝑛 =4 .

b. Estimate Si⁡(𝜋/2) by Simpson’s Rule with 𝑛 =4 .

c. Express the error bound you found in part (a) as a percentage of the value you found in part (b).

  1. The error function The error function,
erf⁡(𝑥)=2√𝜋∫𝑥0𝑒−𝑡2𝑑𝑡,

which is important in probability and in the theories of heat flow and signal transmission, must be evaluated numerically because there is no elementary expression for the antiderivative of 𝑒−𝑡2 .

a. Use Simpson’s Rule with 𝑛 =10 to estimate erf (1).

b. In [0, 1],

∣𝑑4𝑑𝑡4(𝑒−𝑡2)∣≤12.

Give an upper bound for the magnitude of the error of the estimate in part (a).

  1. Prove that the sum T in the Trapezoidal Rule for ∫𝑏𝑎𝑓(𝑥)𝑑𝑥 is a Riemann sum for f continuous on [a, b]. (Hint: Use the Intermediate Value Theorem to show the existence of 𝑐𝑘 in the subinterval [𝑥𝑘−1,𝑥𝑘] satisfying 𝑓(𝑐𝑘) =(𝑓(𝑥𝑘−1) +𝑓(𝑥𝑘))/2 .)

  2. Prove that the sum 𝑆 in Simpson’s Rule for ∫𝑏𝑎𝑓(𝑥)𝑑𝑥 is a Riemann sum for 𝑓 continuous on [𝑎,𝑏] . (See Exercise 29.)

  3. Elliptic integrals The length of the ellipse

𝑥2𝑎2+𝑦2𝑏2=1

turns out to be

 Length =4𝑎∫𝜋/20√1−𝑒2cos2⁡𝑡𝑑𝑡,

where 𝑒 =√𝑎2−𝑏2/𝑎 is the ellipse’s eccentricity. The integral in this formula, called an elliptic integral, is nonelementary except when 𝑒 =0 or 1.

a. Use the Trapezoidal Rule with 𝑛 =10 to estimate the length of the ellipse when 𝑎 =1 and 𝑒 =1/2 .

b. Use the fact that the absolute value of the second derivative of 𝑓(𝑡) =√1−𝑒2cos2⁡𝑡 is less than 1 to find an upper bound for the error in the estimate you obtained in part (a).

Applications

  1. The length of one arch of the curve 𝑦 =sin⁡𝑥 is given by
𝐿=∫𝜋0√1+cos2⁡𝑥𝑑𝑥.

Estimate 𝐿 by Simpson’s Rule with 𝑛 =8 .

When solving Exercises 33-40, you may need to use a calculator or a computer.

  1. Your metal fabrication company is bidding for a contract to make sheets of corrugated iron roofing like the one shown here. The cross-sections of the corrugated sheets are to conform to the curve
𝑦=sin⁡3𝜋20𝑥,0≤𝑥≤20cm.

If the roofing is to be stamped from flat sheets by a process that does not stretch the material, how wide should the original material be? To find out, use numerical integration to approximate the length of the sine curve to two decimal places.

教材插图

  1. Your engineering firm is bidding for the contract to construct the tunnel shown here. The tunnel is 90 m long and 15 m wide at the base. The cross-section is shaped like one arch of the curve 𝑦 =7.5cos⁡(𝜋𝑥/15) . Upon completion, the tunnel’s inside surface (excluding the roadway) will be treated with a waterproof sealer that costs $26.11 per square meter to apply. How much will it cost to apply the sealer? (Hint: Use numerical integration to find the length of the cosine curve.)

教材插图

Find, to two decimal places, the areas of the surfaces generated by revolving the curves in Exercises 35 and 36 about the x-axis.

  1. 𝑦 =sin⁡𝑥,0 ≤𝑥 ≤𝜋

  2. 𝑦 =𝑥2/4,0 ≤𝑥 ≤2

  3. Use numerical integration to estimate the value of

arcsin⁡0.6=∫0.60𝑑𝑥√1−𝑥2.

For reference, arcsin 0.6 = 0.64350 to five decimal places.

  1. Use numerical integration to estimate the value of
𝜋=4∫1011+𝑥2𝑑𝑥.39.$𝐷𝑟𝑢𝑔𝑎𝑠𝑠𝑖𝑚𝑖𝑙𝑎𝑡𝑖𝑜𝑛𝐴𝑛𝑎𝑣𝑒𝑟𝑎𝑔𝑒𝑎𝑑𝑢𝑙𝑡𝑢𝑛𝑑𝑒𝑟𝑎𝑔𝑒60𝑦𝑒𝑎𝑟𝑠𝑎𝑠𝑠𝑖𝑚𝑖𝑙𝑎𝑡𝑒𝑠𝑎12−ℎ𝑜𝑢𝑟𝑐𝑜𝑙𝑑𝑚𝑒𝑑𝑖𝑐𝑖𝑛𝑒𝑖𝑛𝑡𝑜ℎ𝑖𝑠𝑜𝑟ℎ𝑒𝑟𝑠𝑦𝑠𝑡𝑒𝑚𝑎𝑡𝑎𝑟𝑎𝑡𝑒𝑚𝑜𝑑𝑒𝑙𝑒𝑑𝑏𝑦$𝑑𝑦𝑑𝑡=6−ln⁡(2𝑡2−3𝑡+3),

where 𝑦 is measured in milligrams and 𝑡 is the time in hours since the medication was taken. What amount of medicine is absorbed into a person’s system over a 12-hour period?

  1. Effects of an antihistamine The concentration of an antihistamine in the bloodstream of a healthy adult is modeled by
𝐶=12.5−4ln⁡(𝑡2−3𝑡+4),

where C is measured in grams per liter and t is the time in hours since the medication was taken. What is the average level of concentration in the bloodstream over a 6-hour period?

8.8 Improper Integrals

Up to now, we have required definite integrals to satisfy two properties. First, the domain of integration [𝑎,𝑏] must be finite. Second, the range of the integrand must be finite on this domain. In practice, we may encounter problems that fail to meet one or both of these conditions. The integral for the area under the curve 𝑦 =(ln⁡𝑥)/𝑥2 from x = 1 to 𝑥 =∞ is an example for which the domain is infinite (Figure 8.12a). The integral for the area under the curve of 𝑦 =1/√𝑥 between x = 0 and x = 1 is an example for which the range of the integrand is infinite (Figure 8.12b). In either case, the integrals are said to be improper and are calculated as limits. We will see in Chapter 9 that improper integrals are useful for investigating the convergence of certain infinite series.

教材插图

教材插图

(b)

教材插图

(a)

教材插图

(b)

FIGURE 8.13 (a) The area in the first quadrant under the curve 𝑦 =𝑒−𝑥/2 . (b) The area is an improper integral of the first type.

(a)

FIGURE 8.12 Are the areas under these infinite curves finite? We will see that the answer is yes for both curves.

Infinite Limits of Integration

Consider the infinite region (unbounded on the right) that lies under the curve 𝑦 =𝑒−𝑥/2 in the first quadrant (Figure 8.13a). You might think this region has infinite area, but we will see that the value is finite. We assign a value to the area in the following way. First find the area 𝐴(𝑏) of the portion of the region that is bounded on the right by x = b (Figure 8.13b).

𝐴(𝑏)=∫𝑏0𝑒−𝑥/2𝑑𝑥=−2𝑒−𝑥/2]𝑏0=−2𝑒−𝑏/2+2=2−2𝑒−𝑏/2,

which is a little less than 2. Then find the limit of 𝐴(𝑏) as 𝑏 →∞ .

lim𝑏→∞𝐴(𝑏)=lim𝑏→∞(2−2𝑒−𝑏/2)=2

Therefore, the value we assign to the area under the curve from 0 to ∞ is

∫∞0𝑒−𝑥/2𝑑𝑥=lim𝑏→∞∫𝑏0𝑒−𝑥/2𝑑𝑥=2.

DEFINITION Integrals with infinite limits of integration are improper integrals of Type I.

  1. If 𝑓(𝑥) is continuous on [𝑎,∞) , then
∫∞𝑎𝑓(𝑥)𝑑𝑥=lim𝑏→∞∫𝑏𝑎𝑓(𝑥)𝑑𝑥.
  1. If 𝑓(𝑥) is continuous on ( −∞,𝑏] , then
∫𝑏−∞𝑓(𝑥)𝑑𝑥=lim𝑎→−∞∫𝑏𝑎𝑓(𝑥)𝑑𝑥.
  1. If 𝑓(𝑥) is continuous on ( −∞,∞) , then
∫∞−∞𝑓(𝑥)𝑑𝑥=∫𝑐−∞𝑓(𝑥)𝑑𝑥+∫∞𝑐𝑓(𝑥)𝑑𝑥,

where 𝑐 is any real number.

In each case, if the limit exists and is finite, we say that the improper integral converges and that the limit is the value of the improper integral. If the limit fails to exist, the improper integral diverges.

教材插图

FIGURE 8.14 The area under this curve is an improper integral (Example 1).

HISTORICAL BIOGRAPHY

Lejeune Dirichlet (1805–1859)

Dirichlet, a German mathematician, investigated the solution and equilibrium of systems of differential equations and discovered many results on the convergence of series. In 1855, Dirichlet succeeded Gauss as the professor of mathematics at Göttingen.

To know more, visit the companion Website.

The choice of 𝑐 in Part 3 of the definition is unimportant. We can evaluate or determine the convergence or divergence of ∫∞−∞𝑓(𝑥)𝑑𝑥 with any convenient choice.

Any of the integrals in the above definition can be interpreted as an area if 𝑓 ≥0 on the interval of integration. For instance, we interpreted the improper integral in Figure 8.13 as an area. In that case, the area has the finite value 2. If 𝑓 ≥0 and the improper integral diverges, we say the area under the curve is infinite.

EXAMPLE 1 Is the area under the curve 𝑦 =(ln⁡𝑥)/𝑥2 from 𝑥 =1 to 𝑥 =∞ finite? If so, what is its value?

Solution We find the area under the curve from x = 1 to x = b and examine the limit as 𝑏 →∞ . If the limit is finite, we take it to be the area under the curve (Figure 8.14). The area from 1 to b is

∫𝑏1ln⁡𝑥𝑥2𝑑𝑥=[(ln⁡𝑥)(−1𝑥)]𝑏1−∫𝑏1(−1𝑥)(1𝑥)𝑑𝑥 Integration by parts with 𝑢=ln⁡𝑥,𝑑𝑣=𝑑𝑥/𝑥2,𝑑𝑢=𝑑𝑥/𝑥,𝑣=−1/𝑥=−ln⁡𝑏𝑏−[1𝑥]𝑏1=−ln⁡𝑏𝑏−1𝑏+1.

The limit of the area as 𝑏 →∞ is

∫∞1ln⁡𝑥𝑥2𝑑𝑥=lim𝑏→∞∫𝑏1ln⁡𝑥𝑥2𝑑𝑥=lim𝑏→∞[−ln⁡𝑏𝑏−1𝑏+1]=−[lim𝑏→∞ln⁡𝑏𝑏]−0+1=−[lim𝑏→∞1/𝑏1]+1=0+1=1.(\text{L'Hôpital's Rule})

Thus, the improper integral converges and the area has finite value 1.

EXAMPLE 2 Evaluate

∫∞−∞𝑑𝑥1+𝑥2.

Solution According to Part 3 of the definition, we can choose c = 0 and write

∫∞−∞𝑑𝑥1+𝑥2=∫0−∞𝑑𝑥1+𝑥2+∫∞0𝑑𝑥1+𝑥2.

Next we evaluate each improper integral on the right side of the equation above.

∫0−∞𝑑𝑥1+𝑥2=lim𝑎→−∞∫0𝑎𝑑𝑥1+𝑥2=lim𝑎→−∞tan−1⁡𝑥]0𝑎=lim𝑎→−∞(tan−1⁡0−tan−1⁡𝑎)=0−(−𝜋2)=𝜋2

教材插图

FIGURE 8.15 The area under this curve is finite (Example 2).

Thus,

∫∞0𝑑𝑥1+𝑥2=lim𝑏→∞∫𝑏0𝑑𝑥1+𝑥2=lim𝑏→∞tan−1⁡𝑥]𝑏0=lim𝑏→∞(tan−1⁡𝑏−tan−1⁡0)=𝜋2−0=𝜋2 ∫∞−∞𝑑𝑥1+𝑥2=𝜋2+𝜋2=𝜋.

Since 1/(1 +𝑥2) >0 , the improper integral can be interpreted as the (finite) area beneath the curve and above the x-axis (Figure 8.15).

The Integral ∫∞1𝑑𝑥𝑥𝑝

The function y = 1/x is the boundary between the convergent and divergent improper integrals with integrands of the form 𝑦 =1/𝑥𝑝 . As the next example shows, the improper integral converges if p > 1 and diverges if 𝑝 ≤1 .

EXAMPLE 3 For what values of p does the integral ∫∞1𝑑𝑥/𝑥𝑝 converge? When the integral does converge, what is its value?

Solution If 𝑝 ≠1 , then

∫𝑏1𝑑𝑥𝑥𝑝=𝑥−𝑝+1−𝑝+1]𝑏1=(11−𝑝)(𝑏−𝑝+1−1)=(11−𝑝)(1𝑏𝑝−1−1).

Thus,

∫∞1𝑑𝑥𝑥𝑝=lim𝑏→∞∫𝑏1𝑑𝑥𝑥𝑝=lim𝑏→∞[(11−𝑝)(1𝑏𝑝−1−1)]={1𝑝−1,𝑝>1∞,𝑝<1

because

lim𝑏→∞1𝑏𝑝−1={0,𝑝>1∞,𝑝<1.

Therefore, the integral converges to the value 1/(𝑝 −1) if 𝑝 >1 , and it diverges if 𝑝 <1 . If 𝑝 =1 , the integral also diverges:

∫∞1𝑑𝑥𝑥𝑝=∫∞1𝑑𝑥𝑥=lim𝑏→∞∫𝑏1𝑑𝑥𝑥=lim𝑏→∞[ln⁡|𝑥|]𝑏1=lim𝑏→∞(ln⁡𝑏−ln⁡1)=∞.

教材插图

FIGURE 8.16 The area under this curve is an example of an improper integral of the second kind.

Integrands with Vertical Asymptotes

Another type of improper integral arises when the integrand has a vertical asymptote—an infinite discontinuity—at a limit of integration or at some point between the limits of integration. If the integrand f is positive over the interval of integration, we can again interpret the improper integral as the area under the graph of f and above the x-axis between the limits of integration.

Consider the region in the first quadrant that lies under the curve 𝑦 =1/√𝑥 from 𝑥 =0 to 𝑥 =1 (Figure 8.12b). First we find the area of the portion from 𝑎 to 1 (Figure 8.16):

∫1𝑎𝑑𝑥√𝑥=2√𝑥∣1𝑎=2−2√𝑎.

Then we find the limit of this area as 𝑎 →0+ :

lim𝑎→0+∫1𝑎𝑑𝑥√𝑥=lim𝑎→0+(2−2√𝑎)=2.

Therefore, the area under the curve from 0 to 1 is finite and is defined to be

∫10𝑑𝑥√𝑥=lim𝑎→0+∫1𝑎𝑑𝑥√𝑥=2.

DEFINITION Integrals of functions that become infinite at a point within the interval of integration are improper integrals of Type II.

  1. If 𝑓(𝑥) is continuous on (𝑎,𝑏] and discontinuous at a, then
∫𝑏𝑎𝑓(𝑥)𝑑𝑥=lim𝑐→𝑎+∫𝑏𝑐𝑓(𝑥)𝑑𝑥.
  1. If 𝑓(𝑥) is continuous on [𝑎,𝑏) and discontinuous at b, then
∫𝑏𝑎𝑓(𝑥)𝑑𝑥=lim𝑐→𝑏−∫𝑐𝑎𝑓(𝑥)𝑑𝑥.
  1. If 𝑓(𝑥) is discontinuous at 𝑐 , where 𝑎 <𝑐 <𝑏 , and continuous on [𝑎,𝑐) ∪(𝑐,𝑏] , then
∫𝑏𝑎𝑓(𝑥)𝑑𝑥=∫𝑐𝑎𝑓(𝑥)𝑑𝑥+∫𝑏𝑐𝑓(𝑥)𝑑𝑥.

In each case, if the limit exists and is finite, we say that the improper integral converges and that the limit is the value of the improper integral. If the limit does not exist, the integral diverges.

In Part 3 of the definition, the integral on the left side of the equation converges if both integrals on the right side converge; otherwise, it diverges.

EXAMPLE 4 Investigate the convergence of

∫1011−𝑥𝑑𝑥.

教材插图

FIGURE 8.17 The area beneath the curve and above the x-axis for [0,1) is not a real number (Example 4).

教材插图

FIGURE 8.18 Example 5 shows that the area under the curve exists (so it is a real number).

Solution The integrand 𝑓(𝑥) =1/(1 −𝑥) is continuous on [0, 1) but is discontinuous at x = 1 and becomes infinite as 𝑥 →1− (Figure 8.17). We evaluate the integral as

lim𝑏→1−∫𝑏011−𝑥𝑑𝑥=lim𝑏→1−[−ln⁡|1−𝑥|]𝑏0=lim𝑏→1−[−ln⁡(1−𝑏)+0]=∞.

The limit is infinite, so the integral diverges.

EXAMPLE 5 Evaluate

∫30𝑑𝑥(𝑥−1)2/3.

Solution The integrand has a vertical asymptote at x = 1 and is continuous on [0,1) and (1,3] (Figure 8.18). Thus, by Part 3 of the definition above,

∫30𝑑𝑥(𝑥−1)2/3=∫10𝑑𝑥(𝑥−1)2/3+∫31𝑑𝑥(𝑥−1)2/3.

Next, we evaluate each improper integral on the right-hand side of this equation.

∫10𝑑𝑥(𝑥−1)2/3=lim𝑏→1−∫𝑏0𝑑𝑥(𝑥−1)2/3=lim𝑏→1−3(𝑥−1)1/3]𝑏0=lim𝑏→1−[3(𝑏−1)1/3+3]=3 ∫31𝑑𝑥(𝑥−1)2/3=lim𝑐→1+∫3𝑐𝑑𝑥(𝑥−1)2/3=lim𝑐→1+3(𝑥−1)1/3]3𝑐=lim𝑐→1+[3(3−1)1/3−3(𝑐−1)1/3]=33√2

We conclude that

∫30𝑑𝑥(𝑥−1)2/3=3+33√2.

Improper Integrals with a CAS

Computer algebra systems can evaluate many convergent improper integrals. To evaluate the integral

∫∞2𝑥+3(𝑥−1)(𝑥2+1)𝑑𝑥

(which converges) using Maple, enter

𝑓:=(𝑥+3)/((𝑥−1)∗(𝑥∧2+1));

Then use the integration command

int⁡(𝑓,𝑥=2.. infinity );

Maple returns the answer

ln⁡(5)+arctan⁡(2)−𝜋2

To obtain a numerical result, use the evaluation command evalf and specify the number of digits as follows:

 evalf (%,6);

The symbol % instructs the computer to evaluate the last expression on the screen, in this case ( −1/2)𝜋 +ln⁡(5) +arctan⁡(2) . Maple returns 1.14579.

If you are using Mathematica, entering

 In [1]:= Integrate [(𝑥+3)/((𝑥−1)(𝑥∧2+1)),{𝑥,2, Infinity }]

returns

𝑂𝑢𝑡[1]=−𝜋2+ArcTan⁡[2]+Log⁡[5].

教材插图

FIGURE 8.19 The graph of 𝑒−𝑥2 lies below the graph of 𝑒−𝑥 for x > 1 (Example 6a).

HISTORICAL BIOGRAPHY Karl Weierstrass (1815–1897)

Weierstrass attended the University of Bonn to learn public administration, but he found that his passion was for mathematics. In his Berlin lectures in the 1860s, he also proved several theorems for continuous and complex functions. The standards of rigor that he set greatly affected the future of mathematics.

To obtain a numerical result with six digits, use the command “N[%, 6]”; it also yields 1.14579.

To know more, visit the companion Website.

Tests for Convergence and Divergence

When we cannot evaluate an improper integral directly, we try to determine whether it converges or diverges. If the integral diverges, that’s the end of the story. If it converges, we can use numerical methods to approximate its value. The principal tests for convergence or divergence are the Direct Comparison Test and the Limit Comparison Test.

THEOREM 2—Direct Comparison Test Let f and g be continuous on [𝑎,∞) with 0 ≤𝑓(𝑥) ≤𝑔(𝑥) for all 𝑥 ≥𝑎 . Then

  1. if ∫∞𝑎𝑔(𝑥)𝑑𝑥 converges, then ∫∞𝑎𝑓(𝑥)𝑑𝑥 also converges.
  2. if ∫∞𝑎𝑓(𝑥)𝑑𝑥 diverges, then ∫∞𝑎𝑔(𝑥)𝑑𝑥 also diverges.

Outline of a Proof Assume that ∫∞𝑎𝑔(𝑥)𝑑𝑥 converges. For each number 𝑏 ≥𝑎 , let 𝐼(𝑏) =∫𝑏𝑎𝑓(𝑥)𝑑𝑥 . Since 𝑓(𝑥) is nonnegative, we know that 𝐼(𝑏) is an increasing function of b. Using the completeness property of the real numbers (Appendix A.7), it can be shown that since this integral increases with b, it either has a limit or it diverges to infinity as 𝑏 →∞ . We will not give a proof of this fact here. Since 𝑓(𝑥) ≤𝑔(𝑥) for every x,

𝐼(𝑏)=∫𝑏𝑎𝑓(𝑥)𝑑𝑥≤∫𝑏𝑎𝑔(𝑥)𝑑𝑥≤∫∞𝑎𝑔(𝑥)𝑑𝑥.

Thus, the finite number 𝑀 =∫∞𝑎𝑔(𝑥) 𝑑𝑥 is an upper bound to the values of 𝐼(𝑏) . Therefore 𝐼(𝑏) cannot diverge to infinity, so it must must converge to a finite value as 𝑏 →∞ . Hence ∫∞𝑎𝑓(𝑥) 𝑑𝑥 converges. This establishes that statement 1 holds. Since statement 2 is the contrapositive form of statement 1, it must hold as well.

Although the theorem is stated for Type I improper integrals, a similar result is true for integrals of Type II as well.

EXAMPLE 6 These examples illustrate how we use Theorem 2.

(a) ∫∞1𝑒−𝑥2𝑑𝑥 converges because 0 <𝑒−𝑥2 <𝑒−𝑥 for every 𝑥 ≥1 (Figure 8.19) and

∫∞1𝑒−𝑥𝑑𝑥=lim𝑏→∞∫𝑏1𝑒−𝑥𝑑𝑥=lim𝑏→∞[−𝑒−𝑥𝑑𝑥]𝑏1=lim𝑏→∞(−𝑒−𝑏+𝑒−1)=1𝑒

converges.

(b) ∫∞1sin2⁡𝑥𝑥2𝑑𝑥 converges because

0≤sin2⁡𝑥𝑥2≤1𝑥2 on [1,∞) and ∫∞11𝑥2𝑑𝑥 converges. (Example3) ∫∞11√𝑥2−0.1𝑑𝑥 diverges because (c) 1√𝑥2−0.1≥1𝑥 on [1,∞) and ∫∞11𝑥𝑑𝑥 diverges. (Example3)

(d) ∫𝜋/20cos⁡𝑥√𝑥𝑑𝑥 converges because

0≤cos⁡𝑥√𝑥≤1√𝑥 on [0,𝜋2],0≤cos⁡𝑥≤1 on [0,𝜋2]

and

∫𝜋/20𝑑𝑥√𝑥=lim𝑎→0+∫𝜋/2𝑎𝑑𝑥√𝑥=lim𝑎→0+√4𝑥∣𝜋/2𝑎2√𝑥=√4𝑥=lim𝑎→0+(√2𝜋−√4𝑎)=√2𝜋 converges .

Although we have shown that the integrals in parts (a), (b), and (d) of Example 6 converge, we do not know the exact values of these integrals. For example, our computations in part (d) imply that

0≤∫𝜋/20cos⁡𝑥√𝑥𝑑𝑥≤∫𝜋/20𝑑𝑥√𝑥=√2𝜋,

but unless we do further calculations the most we can say is that the integral is some real number between 0 and √2𝜋 .

THEOREM 3—Limit Comparison Test

If the positive functions 𝑓 and 𝑔 are continuous on [𝑎,∞) , and if

lim𝑥→∞𝑓(𝑥)𝑔(𝑥)=𝐿,0<𝐿<∞,

then

∫∞𝑎𝑓(𝑥)𝑑𝑥 and ∫∞𝑎𝑔(𝑥)𝑑𝑥

either both converge or both diverge.

We omit the proof of Theorem 3, which is similar to that of Theorem 2.

If two functions 𝑓(𝑥) and 𝑔(𝑥) satisfy the hypotheses of Theorem 3 and their improper integrals both converge, it need not be the case that these two integrals have the same value. This is illustrated in the next example.

EXAMPLE 7 Show that

∫∞1𝑑𝑥1+𝑥2

converges by comparison with ∫∞1(1/𝑥2)𝑑𝑥 . Find and compare the values of the two integrals.

教材插图

Solution The functions 𝑓(𝑥) =1/𝑥2 and 𝑔(𝑥) =1/(1 +𝑥2) are positive and continuous on [1,∞) . Also,

lim𝑥→∞𝑓(𝑥)𝑔(𝑥)=lim𝑥→∞1/𝑥21/(1+𝑥2)=lim𝑥→∞1+𝑥2𝑥2=lim𝑥→∞(1𝑥2+1)=0+1=1,

which is a positive finite limit (Figure 8.20). Therefore, ∫∞1𝑑𝑥1+𝑥2 converges because ∫∞1𝑑𝑥𝑥2 converges.

FIGURE 8.20 The functions in Example 7.

The integrals converge to different values, however:

∫∞1𝑑𝑥𝑥2=12−1=1 Example 3 

and

∫∞1𝑑𝑥1+𝑥2=lim𝑏→∞∫𝑏1𝑑𝑥1+𝑥2=lim𝑏→∞[tan−1⁡𝑏−tan−1⁡1]=𝜋2−𝜋4=𝜋4.

TABLE 8.5

b∫𝑏11−𝑒−𝑥𝑥𝑑𝑥
20.5226637569
51.3912002736
102.0832053156
1004.3857862516
10006.6883713446
100008.9909564376
10000011.2935415306

EXAMPLE 8 Investigate the convergence of ∫∞11−𝑒−𝑥𝑥𝑑𝑥 .

Solution The integrand suggests a comparison of 𝑓(𝑥) =(1 −𝑒−𝑥)/𝑥 with 𝑔(𝑥) =1/𝑥 . However, we cannot use the Direct Comparison Test because 𝑓(𝑥) ≤𝑔(𝑥) and the integral of 𝑔(𝑥) diverges. On the other hand, using the Limit Comparison Test, we find that

lim𝑥→∞𝑓(𝑥)𝑔(𝑥)=lim𝑥→∞(1−𝑒−𝑥𝑥)(𝑥1)=lim𝑥→∞(1−𝑒−𝑥)=1,

which is a positive finite limit. Therefore, ∫∞11−𝑒−𝑥𝑥𝑑𝑥 diverges because ∫∞1𝑑𝑥𝑥 diverges. Approximations to the improper integral are given in Table 8.5. Note that the values of these approximations do not appear to approach a fixed finite limit as 𝑏 →∞ .

EXERCISES 8.8

Evaluating Improper Integrals

The integrals in Exercises 1–34 converge. Evaluate the integrals without using tables.

  1. ∫∞0𝑑𝑥𝑥2+1

  2. ∫∞1𝑑𝑥𝑥1.001

  3. ∫10𝑑𝑥√𝑥

  4. ∫40𝑑𝑥√4−𝑥

  5. ∫1−1𝑑𝑥𝑥2/3

  6. ∫1−8𝑑𝑥𝑥1/3

  7. ∫10𝑑𝑥√1−𝑥2

  8. ∫10𝑑𝑟𝑟0.999

  9. ∫−2−∞2𝑑𝑥𝑥2−1

  10. ∫2−∞2𝑑𝑥𝑥2+4

  11. ∫∞22𝑣2−𝑣𝑑𝑣

  12. ∫∞22𝑑𝑡𝑡2−1

  13. ∫∞−∞2𝑥𝑑𝑥(𝑥2+1)2

  14. ∫∞−∞𝑥𝑑𝑥(𝑥2+4)3/2

  15. ∫10𝜃+1√𝜃2+2𝜃𝑑𝜃

  16. ∫20𝑠+1√4−𝑠2𝑑𝑠

  17. ∫∞0𝑑𝑥(1+𝑥)√𝑥

  18. ∫∞11𝑥√𝑥2−1𝑑𝑥

  19. ∫∞0𝑑𝑣(1+𝑣2)(1+tan−1⁡𝑣)

  20. ∫∞016tan−1⁡𝑥1+𝑥2𝑑𝑥

  21. ∫0−∞𝜃𝑒𝜃𝑑𝜃

  22. ∫∞02𝑒−𝜃sin⁡𝜃𝑑𝜃

  23. ∫0−∞𝑒−|𝑥|𝑑𝑥

  24. ∫∞−∞2𝑥𝑒−𝑥2𝑑𝑥

  25. ∫10𝑥ln⁡𝑥𝑑𝑥

  26. ∫10( −ln⁡𝑥)𝑑𝑥

  27. ∫20𝑑𝑠√4−𝑠2

  28. ∫104𝑟𝑑𝑟√1−𝑟4

  29. ∫21𝑑𝑠𝑠√𝑠2−1

  30. ∫42𝑑𝑡𝑡√𝑡2−4

  31. ∫4−1𝑑𝑥√|𝑥|

  32. ∫20𝑑𝑥√|𝑥−1|

  33. ∫∞−1𝑑𝜃𝜃2+5𝜃+6

  34. ∫∞0𝑑𝑥(𝑥+1)(𝑥2+1)

Testing for Convergence

In Exercises 35–68, use integration, the Direct Comparison Test, or the Limit Comparison Test to test the integrals for convergence. If more than one method applies, use whatever method you prefer.

  1. ∫21/2𝑑𝑥𝑥ln⁡𝑥

  2. ∫1−1𝑑𝜃𝜃2−2𝜃

  3. ∫∞1/2𝑑𝑥𝑥(ln⁡𝑥)3

  4. ∫∞0𝑑𝜃𝜃2−1

  5. ∫𝜋/20tan⁡𝜃𝑑𝜃

  6. ∫𝜋/20cot⁡𝜃𝑑𝜃

  7. ∫10ln⁡𝑥𝑥2𝑑𝑥

  8. ∫21𝑑𝑥𝑥ln⁡𝑥

Theory and Examples

  1. ∫ln⁡20𝑥−2𝑒−1/𝑥𝑑𝑥

  2. ∫10𝑒−√𝑥√𝑥𝑑𝑥

  3. ∫𝜋0𝑑𝑡√𝑡+sin⁡𝑡

  4. ∫10𝑑𝑡𝑡−sin⁡𝑡 (Hint: 𝑡 ≥sin⁡𝑡 for 𝑡 ≥0 )

  5. ∫20𝑑𝑥1−𝑥2

  6. ∫20𝑑𝑥1−𝑥

  7. ∫1−1ln⁡|𝑥|𝑑𝑥

diverges and hence that

  1. ∫1−1 −𝑥ln⁡|𝑥|𝑑𝑥
∫∞−∞2𝑥𝑑𝑥𝑥2+1
  1. ∫∞1𝑑𝑥𝑥3+1

  2. ∫∞4𝑑𝑥√𝑥−1

  3. ∫∞2𝑑𝑣√𝑣−1

  4. ∫∞0𝑑𝜃1+𝑒𝜃

  5. ∫∞0𝑑𝑥√𝑥6+1

  6. ∫∞2𝑑𝑥√𝑥2−1

diverges. Then show that

  1. ∫∞1√𝑥+1𝑥2𝑑𝑥

  2. ∫∞2𝑥𝑑𝑥√𝑥4−1

lim𝑏→∞∫𝑏−𝑏2𝑥𝑑𝑥𝑥2+1 =0.

  1. ∫∞𝜋2+cos⁡𝑥𝑥𝑑𝑥

  2. ∫∞𝜋1+sin⁡𝑥𝑥2𝑑𝑥

Exercises 83–86 are about the infinite region in the first quadrant between the curve 𝑦 =𝑒−𝑥 and the x-axis.

  1. ∫∞42𝑑𝑡𝑡3/2−1

  2. ∫∞21ln⁡𝑥𝑑𝑥

  3. ∫∞1𝑒𝑥𝑥𝑑𝑥

  4. ∫∞𝑒𝑒ln⁡(ln⁡𝑥)𝑑𝑥

  5. ∫∞11√𝑒𝑥−𝑥𝑑𝑥

  6. ∫∞11𝑒𝑥−2𝑥𝑑𝑥

  7. ∫∞−∞𝑑𝑥√𝑥4+1

  8. ∫∞−∞𝑑𝑥𝑒𝑥+𝑒−𝑥

In Exercises 69–80, determine whether the improper integral converges or diverges. If it converges, evaluate the integral.

  1. ∫101𝑥√𝑥𝑑𝑥

  2. ∫∞21𝑥√𝑥𝑑𝑥

  3. ∫32015√𝑥𝑑𝑥

  4. ∫∞115√𝑥𝑑𝑥

  5. ∫∞31𝑥4𝑑𝑥

  6. ∫1−21𝑥4𝑑𝑥

  7. ∫∞0𝑥2𝑒𝑥3𝑑𝑥

  8. ∫0−∞𝑥2𝑒𝑥3𝑑𝑥

  9. ∫0−31𝑥2+3𝑥𝑑𝑥

  10. ∫∞11𝑥2+3𝑥𝑑𝑥

  11. ∫4−∞𝑥(𝑥2+9)5/2𝑑𝑥

  12. ∫4−∞𝑥(𝑥2+9)2/5𝑑𝑥

  13. Find the values of 𝑝 for which each integral converges. a. ∫21𝑑𝑥𝑥(ln⁡𝑥)𝑝 b. ∫∞2𝑑𝑥𝑥(ln⁡𝑥)𝑝

  14. ∫∞−∞𝑓(𝑥)𝑑𝑥 may not equal lim𝑏→∞∫𝑏−𝑏𝑓(𝑥)𝑑𝑥. Show that ∫∞02𝑥𝑑𝑥𝑥2+1

  15. Find the area of the region.

  16. Find the centroid of the region.

  17. Find the volume of the solid generated by revolving the region about the y-axis.

  18. Find the volume of the solid generated by revolving the region about the x-axis.

  19. Find the area of the region that lies between the curves 𝑦 =sec⁡𝑥 and 𝑦 =tan⁡𝑥 from 𝑥 =0 to 𝑥 =𝜋/2 .

  20. The region in Exercise 87 is revolved about the 𝑥 -axis to generate a solid.

a. Find the volume of the solid.

b. Show that the inner and outer surfaces of the solid have infinite area.

  1. Consider the infinite region in the first quadrant bounded by the graphs of 𝑦 =1𝑥2 , 𝑦 =0 , and 𝑥 =1 .

a. Find the area of the region.

b. Find the volume of the solid formed by revolving the region (i) about the x-axis; (ii) about the y-axis.

  1. Consider the infinite region in the first quadrant bounded by the graphs of 𝑦 =1√𝑥 , 𝑦 =0 , 𝑥 =0 , and 𝑥 =1 .

a. Find the area of the region.

b. Find the volume of the solid formed by revolving the region (i) about the x-axis; (ii) about the y-axis.

  1. Evaluate the integrals.

a. ∫10𝑑𝑡√𝑡(1+𝑡)

𝐛.∫∞0𝑑𝑡√𝑡(1+𝑡)
  1. Evaluate ∫∞3𝑑𝑥𝑥√𝑥2−9

  2. Estimating the value of a convergent improper integral whose domain is infinite

a. Show that

∫∞3𝑒−3𝑥𝑑𝑥=13𝑒−9<0.000042,

and hence that ∫∞3𝑒−𝑥2𝑑𝑥 <0.000042 . Explain why this means that ∫∞0𝑒−𝑥2𝑑𝑥 can be replaced by ∫30𝑒−𝑥2𝑑𝑥 without introducing an error of magnitude greater than 0.000042.

T b. Evaluate ∫30𝑒−𝑥2𝑑𝑥 numerically.

  1. The infinite paint can or Gabriel’s horn As Example 3 shows, the integral ∫∞1(𝑑𝑥/𝑥) diverges. This means that the integral
∫∞12𝜋1𝑥√1+1𝑥4𝑑𝑥,

which measures the surface area of the solid of revolution traced out by revolving the curve 𝑦 =1/𝑥,1 ≤𝑥 , about the x-axis, diverges also. By comparing the two integrals, we see that, for every finite value b > 1,

∫𝑏12𝜋1𝑥√1+1𝑥4𝑑𝑥>2𝜋∫𝑏11𝑥𝑑𝑥.

教材插图

However, the integral

∫∞1𝜋(1𝑥)2𝑑𝑥

for the volume of the solid converges.

a. Calculate it.

b. This solid of revolution is sometimes described as a can that does not hold enough paint to cover its own interior. Think about that for a moment. It is common sense that a finite

amount of paint cannot cover an infinite surface. But if we fill the horn with paint (a finite amount), then we will have covered an infinite surface. Explain the apparent contradiction.

  1. Sine-integral function The integral
Si⁡(𝑥)=∫𝑥0sin⁡𝑡𝑡𝑑𝑡,

called the sine-integral function, has important applications in optics.

T a. Plot the integrand (sin⁡𝑡)/𝑡 for t > 0. Is the sine-integral function everywhere increasing or decreasing? Do you think Si(𝑥) =0 for 𝑥 ≥0 ? Check your answers by graphing the function Si(𝑥) for 0 ≤𝑥 ≤25 .

b. Explore the convergence of

∫∞0sin⁡𝑡𝑡𝑑𝑡.

If it converges, what is its value?

  1. Error function The function
erf⁡(𝑥)=∫𝑥02𝑒−𝑡2√𝜋𝑑𝑡,

called the error function, has important applications in probability and statistics.

T a. Plot the error function for 0 ≤𝑥 ≤25 .

b. Explore the convergence of

∫∞02𝑒−𝑡2√𝜋𝑑𝑡.

If it converges, what appears to be its value? You will see how to confirm your estimate in Section 14.4, Exercise 41.

  1. Normal probability distribution The function
𝑓(𝑥)=1𝜎√2𝜋𝑒−12(𝑥−𝜇𝜎)2

is called the normal probability density function with mean 𝜇 and standard deviation 𝜎 . The number 𝜇 tells where the distribution is centered, and 𝜎 measures the “scatter” around the mean.

From the theory of probability, it is known that

∫∞−∞𝑓(𝑥)𝑑𝑥=1.

In what follows, let 𝜇 =0 and 𝜎 =1 .

T a. Draw the graph of f. Find the intervals on which f is increasing, the intervals on which f is decreasing, and any local extreme values and where they occur.

b. Evaluate

∫𝑛−𝑛𝑓(𝑥)𝑑𝑥

for n = 1, 2, and 3.

c. Give a convincing argument that

∫∞−∞𝑓(𝑥)𝑑𝑥=1.

(Hint: Show that 0 <𝑓(𝑥) <𝑒−𝑥/2 for 𝑥 >1 , and for 𝑏 >1 ,

∫∞𝑏𝑒−𝑥/2𝑑𝑥→0 as 𝑏→∞.)
  1. Show that if 𝑓(𝑥) is integrable on every interval of real numbers, and if a and b are real numbers with a < b, then

a. ∫𝑎−∞𝑓(𝑥)𝑑𝑥 and ∫∞𝑎𝑓(𝑥)𝑑𝑥 both converge if and only if

∫𝑏−∞𝑓(𝑥)𝑑𝑥 and ∫∞𝑏𝑓(𝑥)𝑑𝑥 both converge.

𝐛.∫𝑎−∞𝑓(𝑥)𝑑𝑥+∫∞𝑎𝑓(𝑥)𝑑𝑥=∫𝑏−∞𝑓(𝑥)𝑑𝑥+∫∞𝑏𝑓(𝑥)𝑑𝑥

when the integrals involved converge.

COMPUTER EXPLORATIONS

In Exercises 99–102, use a CAS to explore the integrals for various values of p (include noninteger values). For what values of p does the

  1. ∫∞𝑒𝑥𝑝ln⁡𝑥𝑑𝑥

integral converge? What is the value of the integral when it does converge? Plot the integrand for various values of p.

  1. ∫∞0𝑥𝑝ln⁡𝑥𝑑𝑥

  2. ∫∞−∞𝑥𝑝ln⁡|𝑥|𝑑𝑥

Use a CAS to evaluate the integrals.

  1. ∫2/𝜋0sin⁡1𝑥𝑑𝑥

  2. ∫2/𝜋0𝑥sin⁡1𝑥𝑑𝑥

CHAPTER 8 Questions to Guide Your Review

  1. What is the formula for integration by parts? Where does it come from? Why might you want to use it?

  2. When applying the formula for integration by parts, how do you choose the u and dv? How can you apply integration by parts to an integral of the form ∫𝑓(𝑥)𝑑𝑥 ?

  3. If an integrand is a product of the form sin𝑛⁡𝑥cos𝑚⁡𝑥 , where m and n are nonnegative integers, how do you evaluate the integral? Give a specific example of each case.

  4. What substitutions are made to evaluate integrals of sin⁡𝑚𝑥 sin nx, sin⁡𝑚𝑥cos⁡𝑛𝑥 , and cos⁡𝑚𝑥cos⁡𝑛𝑥 ? Give an example of each case.

  5. ∫𝑒0𝑥𝑝ln⁡𝑥𝑑𝑥

  6. What substitutions are sometimes used to transform integrals involving √𝑎2−𝑥2 , √𝑎2+𝑥2 , and √𝑥2−𝑎2 into integrals that can be evaluated directly? Give an example of each case.

  7. What restrictions can you place on the variables involved in the three basic trigonometric substitutions to make sure the substitutions are reversible (have inverses)?

  8. What is the goal of the method of partial fractions?

  9. When the degree of a polynomial 𝑓(𝑥) is less than the degree of a polynomial 𝑔(𝑥) , how do you write 𝑓(𝑥)/𝑔(𝑥) as a sum of partial fractions if 𝑔(𝑥)

a. is a product of distinct linear factors?

b. consists of a repeated linear factor?

c. contains an irreducible quadratic factor?

What do you do if the degree of 𝑓 is not less than the degree of 𝑔 ?

  1. How are integral tables typically used? What do you do if a particular integral you want to evaluate is not listed in the table?

  2. What is a reduction formula? How are reduction formulas used? Give an example.

  3. How would you compare the relative merits of the Midpoint Rule, the Trapezoidal Rule, and Simpson’s Rule?

  4. What is an improper integral of Type I? Type II? How are the values of various types of improper integrals defined? Give examples.

  5. What tests are available for determining the convergence and divergence of improper integrals that cannot be evaluated directly? Give examples of their use.

  6. What is a random variable? What is a continuous random variable? Give some specific examples.

  7. What is a probability density function? What is the probability that a continuous random variable has a value in the interval [𝑐,𝑑] ?

  8. What is an exponentially decreasing probability density function? What are some typical events that might be modeled by this distribution? What do we mean when we say such distributions are memoryless?

  9. What is the expected value of a continuous random variable? What is the expected value of an exponentially distributed random variable?

  10. What is the median of a continuous random variable? What is the median of an exponential distribution?

  11. What does the variance of a random variable measure? What is the standard deviation of a continuous random variable 𝑋 ?

  12. What probability density function describes the normal distribution? What are some examples typically modeled by a normal distribution? How do we usually calculate probabilities for a normal distribution?

  13. In a normal distribution, what percentage of the population lies within 1 standard deviation of the mean? Within 2 standard deviations?

CHAPTER 8 Practice Exercises

Integration by Parts

Evaluate the integrals in Exercises 1–8 using integration by parts.

  1. ∫ln⁡(𝑥 +1)𝑑𝑥

  2. ∫𝑥2ln⁡𝑥𝑑𝑥

  3. ∫arctan⁡3𝑥𝑑𝑥

  4. ∫cos−1⁡(𝑥2)𝑑𝑥

  5. ∫(𝑥 +1)2𝑒𝑥𝑑𝑥

  6. ∫𝑥2sin⁡(1 −𝑥)𝑑𝑥

  7. ∫𝑒𝑥cos⁡2𝑥𝑑𝑥

  8. ∫𝑥sin⁡𝑥cos⁡𝑥𝑑𝑥

Partial Fractions

Evaluate the integrals in Exercises 9–28. It may be necessary to use a substitution first.

  1. ∫𝑥𝑑𝑥𝑥2−3𝑥+2

  2. ∫𝑥𝑑𝑥𝑥2+4𝑥+3

  3. ∫𝑑𝑥𝑥(𝑥+1)2

  4. ∫𝑥+1𝑥2(𝑥−1)𝑑𝑥

  5. ∫sin⁡𝜃𝑑𝜃cos2⁡𝜃+cos⁡𝜃−2

  6. ∫cos⁡𝜃𝑑𝜃sin2⁡𝜃+sin⁡𝜃−6

  7. ∫3𝑥2+4𝑥+4𝑥3+𝑥𝑑𝑥

  8. ∫4𝑥𝑑𝑥𝑥3+4𝑥

  9. ∫𝑣+32𝑣3−8𝑣𝑑𝑣

  10. ∫(3𝑣−7)𝑑𝑣(𝑣−1)(𝑣−2)(𝑣−3)

  11. ∫𝑑𝑡𝑡4+4𝑡2+3

  12. ∫𝑡𝑑𝑡𝑡4−𝑡2−2

  13. ∫𝑥3+𝑥2𝑥2+𝑥−2𝑑𝑥

  14. ∫𝑥3+1𝑥3−𝑥𝑑𝑥

  15. ∫𝑥3+4𝑥2𝑥2+4𝑥+3𝑑𝑥

  16. ∫2𝑥3+𝑥2−21𝑥+24𝑥2+2𝑥−8𝑑𝑥

  17. ∫𝑑𝑥𝑥(3√𝑥+1)

  18. ∫𝑑𝑥𝑥(1+3√𝑥)

  19. ∫𝑑𝑠𝑒𝑠−1

  20. ∫𝑑𝑠√𝑒𝑠+1

Trigonometric Substitutions

Evaluate the integrals in Exercises 29–32 (a) without using a trigonometric substitution, (b) using a trigonometric substitution.

  1. ∫𝑦𝑑𝑦√16−𝑦2

  2. ∫𝑥𝑑𝑥√4+𝑥2

  3. ∫𝑥𝑑𝑥4−𝑥2

∫𝑡𝑑𝑡√4𝑡2−1

Evaluate the integrals in Exercises 33–36. 33. ∫𝑥𝑑𝑥9−𝑥2

  1. ∫𝑑𝑥𝑥(9−𝑥2)

  2. ∫𝑑𝑥9−𝑥2

  3. ∫𝑑𝑥√9−𝑥2

Trigonometric Integrals

Evaluate the integrals in Exercises 37–44.

  1. ∫sin3⁡𝑥cos4⁡𝑥𝑑𝑥

  2. ∫cos5⁡𝑥sin5⁡𝑥𝑑𝑥

  3. ∫tan4⁡𝑥sec2⁡𝑥𝑑𝑥

  4. ∫tan3⁡𝑥sec3⁡𝑥𝑑𝑥

  5. ∫sin⁡5𝜃cos⁡6𝜃𝑑𝜃

  6. ∫sec2⁡𝜃sin3⁡𝜃𝑑𝜃

  7. ∫√1+cos⁡(𝑡/2)𝑑𝑡

  8. ∫𝑒𝑡√tan2⁡𝑒𝑡+1𝑑𝑡

Numerical Integration

  1. According to the error-bound formula for Simpson’s Rule, how many subintervals should you use to be sure of estimating the value of
ln⁡3=∫311𝑥𝑑𝑥

by Simpson’s Rule with an error of no more than 10−4 in absolute value? (Remember that for Simpson’s Rule, the number of subintervals has to be even.)

  1. A brief calculation shows that if 0 ≤𝑥 ≤1 , then the second derivative of 𝑓(𝑥) =√1+𝑥4 lies between 0 and 8. Based on this, about how many subdivisions would you need to estimate the integral of 𝑓 from 0 to 1 with an error no greater than 10−3 in absolute value using the Trapezoidal Rule?

  2. A direct calculation shows that

∫𝜋02sin2⁡𝑥𝑑𝑥=𝜋.

How close do you come to this value by using the Trapezoidal Rule with 𝑛 =6 ? The Midpoint Rule with 𝑛 =6 ? Simpson’s Rule with 𝑛 =6 ? Try them and find out.

  1. You are planning to use Simpson’s Rule to estimate the value of the integral
∫21𝑓(𝑥)𝑑𝑥

with an error magnitude less than 10−5 . You have determined that |𝑓(4)(𝑥)| ≤3 throughout the interval of integration. How many subintervals should you use to ensure the required accuracy? (Remember that for Simpson’s Rule, the number has to be even.)

T 49. Mean temperature Compute the average value of the temperature function

𝑓(𝑥)=20sin⁡(2𝜋365(𝑥−101))−4

for a 365-day year. This is one way to estimate the annual mean air temperature in Fairbanks, Alaska. The National Weather Service’s official figure, a numerical average of the daily normal mean air temperatures for the year, is −3.5∘C , which is slightly higher than the average value of 𝑓(𝑥) .

  1. Heat capacity of a gas Heat capacity 𝐶𝑣 is the amount of heat required to raise the temperature of a given mass of gas with constant volume by 1 °C, measured in units of cal/deg-mol (calories per degree gram molecular weight). The heat capacity of oxygen depends on its temperature T and satisfies the formula
𝐶𝑣=8.27+10−5(26𝑇−1.87𝑇2).

Use Simpson’s Rule to find the average value of 𝐶𝑣 and the temperature at which it is attained for 20∘C ≤𝑇 ≤675∘C .

  1. Fuel efficiency An automobile computer gives a digital readout of fuel consumption in liters per hour. During a trip, a passenger recorded the fuel consumption every 5 min for a full hour of travel.
TimeL/hTimeL/h
02.5352.5
52.4402.4
102.3452.3
152.4502.4
202.4552.4
252.5602.3
302.6

a. Use the Trapezoidal Rule to approximate the total fuel consumption during the hour.

b. If the automobile covered 60 km in the hour, what was its fuel efficiency (in kilometers per liter) for that portion of the trip?

  1. A new parking lot To meet the demand for parking, your town has allocated the area shown here. As the town engineer, you have been asked by the town council to find out if the lot can be built for 11,000.𝑇ℎ𝑒𝑐𝑜𝑠𝑡𝑡𝑜𝑐𝑙𝑒𝑎𝑟𝑡ℎ𝑒𝑙𝑎𝑛𝑑𝑤𝑖𝑙𝑙𝑏𝑒1.00 a square meter, and the lot will cost 20.00𝑎𝑠𝑞𝑢𝑎𝑟𝑒𝑚𝑒𝑡𝑒𝑟𝑡𝑜𝑝𝑎𝑣𝑒.𝑈𝑠𝑒𝑆𝑖𝑚𝑝𝑠𝑜𝑛′𝑠𝑅𝑢𝑙𝑒𝑡𝑜𝑓𝑖𝑛𝑑𝑜𝑢𝑡𝑖𝑓𝑡ℎ𝑒𝑗𝑜𝑏𝑐𝑎𝑛𝑏𝑒𝑑𝑜𝑛𝑒𝑓𝑜𝑟11,000.

教材插图

Improper Integrals

Evaluate the improper integrals in Exercises 53–62.

  1. ∫30𝑑𝑥√9−𝑥2

  2. ∫10ln⁡𝑥𝑑𝑥

  3. ∫20𝑑𝑦(𝑦−1)2/3

  4. ∫0−2𝑑𝜃(𝜃+1)3/5

  5. ∫∞32𝑑𝑢𝑢2−2𝑢

  6. ∫∞13𝑣−14𝑣3−𝑣2𝑑𝑣

  7. ∫∞0𝑥2𝑒−𝑥𝑑𝑥

  8. ∫0−∞𝑥𝑒3𝑥𝑑𝑥

  9. ∫∞−∞𝑑𝑥4𝑥2+9

  10. ∫∞−∞4𝑑𝑥𝑥2+16

Which of the improper integrals in Exercises 63–68 converge and which diverge?

  1. ∫∞6𝑑𝜃√𝜃2+1

  2. ∫∞0𝑒−𝑢cos⁡𝑢𝑑𝑢

  3. ∫∞1ln⁡𝑧𝑧𝑑𝑧

  4. ∫∞1𝑒−𝑡√𝑡𝑑𝑡

  5. ∫∞−∞2𝑑𝑥𝑒𝑥+𝑒−𝑥

  6. ∫∞−∞𝑑𝑥𝑥2(1+𝑒𝑥)

Assorted Integrations

Evaluate the integrals in Exercises 69–134. The integrals are listed in random order so you need to decide which integration technique to use.

  1. ∫𝑥𝑒2𝑥𝑑𝑥

  2. ∫10𝑥2𝑒𝑥3𝑑𝑥

  3. ∫(tan2⁡𝑥 +sec2⁡𝑥)𝑑𝑥

  4. ∫𝜋/40cos2⁡2𝑥𝑑𝑥

  5. ∫𝑥sec2⁡𝑥𝑑𝑥

  6. ∫𝑥sec2⁡(𝑥2)𝑑𝑥

  7. ∫sin⁡𝑥cos2⁡𝑥𝑑𝑥

  8. ∫sin⁡2𝑥sin⁡(cos⁡2𝑥)𝑑𝑥

  9. ∫0−1𝑒𝑥𝑒𝑥+𝑒−𝑥𝑑𝑥

  10. ∫(𝑒2𝑥 +𝑒−𝑥)2𝑑𝑥

  11. ∫𝑥+1𝑥4−𝑥3𝑑𝑥

  12. ∫𝑒𝑥+1𝑒𝑥(𝑒2𝑥−4)𝑑𝑥

  13. ∫𝑒𝑥+𝑒3𝑥𝑒2𝑥𝑑𝑥

  14. ∫(𝑒𝑥 −𝑒−𝑥)(𝑒𝑥 +𝑒−𝑥)3𝑑𝑥

  15. ∫𝜋/30tan3⁡𝑥sec2⁡𝑥𝑑𝑥

  16. ∫tan4⁡𝑥sec4⁡𝑥𝑑𝑥

  17. ∫30(𝑥 +2)√𝑥+1𝑑𝑥

  18. ∫(𝑥 +1)√𝑥2+2𝑥𝑑𝑥

  19. ∫cot⁡𝑥csc3⁡𝑥𝑑𝑥

  20. ∫sin⁡𝑥(tan⁡𝑥 −cot⁡𝑥)2𝑑𝑥

  21. ∫𝑥𝑑𝑥1+√𝑥

  22. ∫𝑥3+24−𝑥2𝑑𝑥

  23. ∫√2𝑥−𝑥2𝑑𝑥

  24. ∫𝑑𝑥√−2𝑥−𝑥2

  25. ∫2−cos⁡𝑥+sin⁡𝑥sin2⁡𝑥𝑑𝑥

  26. ∫sin2⁡𝜃cos5⁡𝜃𝑑𝜃

  27. ∫9𝑑𝑣81−𝑣4

  28. ∫∞2𝑑𝑥(𝑥−1)2

  29. ∫𝜃cos⁡(2𝜃 +1)𝑑𝜃

  30. ∫𝑥3𝑑𝑥𝑥2−2𝑥+1

  31. ∫sin⁡2𝜃𝑑𝜃(1+cos⁡2𝜃)2

  32. ∫𝜋/2𝜋/4√1+cos⁡4𝑥𝑑𝑥

  33. ∫𝑥𝑑𝑥√2−𝑥

  34. ∫√1−𝑣2𝑣2𝑑𝑣

  35. ∫𝑑𝑦𝑦2−2𝑦+2

  36. ∫𝑥𝑑𝑥√8−2𝑥2−𝑥4

  37. ∫𝑧+1𝑧2(𝑧2+4)𝑑𝑧

  38. ∫𝑥2(𝑥 −1)1/3𝑑𝑥

  39. ∫𝑡𝑑𝑡√9−4𝑡2

  40. ∫arctan⁡𝑥𝑥2𝑑𝑥

  41. ∫𝑒𝑡𝑑𝑡𝑒2𝑡+3𝑒𝑡+2

  42. ∫tan3⁡𝑡𝑑𝑡

  43. ∫∞1ln⁡𝑦𝑦3𝑑𝑦

  44. ∫𝑦3/2(ln⁡𝑦)2𝑑𝑦

  45. ∫𝑒ln⁡√𝑥𝑑𝑥

  46. ∫𝑒𝜃√3+4𝑒𝜃𝑑𝜃

  47. ∫sin⁡5𝑡𝑑𝑡1+(cos⁡5𝑡)2

  48. ∫𝑑𝑣√𝑒2𝑣−1

  49. ∫𝑑𝑟1+√𝑟

  50. ∫4𝑥3−20𝑥𝑥4−10𝑥2+9𝑑𝑥

  51. ∫𝑥31+𝑥2𝑑𝑥

  52. ∫𝑥21+𝑥3𝑑𝑥

  53. ∫1+𝑥21+𝑥3𝑑𝑥

  54. ∫1+𝑥2(1+𝑥)3𝑑𝑥

  55. ∫√𝑥 ⋅√1+√𝑥𝑑𝑥

  56. ∫√1+√1+𝑥𝑑𝑥

  57. ∫1√𝑥⋅√1+𝑥𝑑𝑥

  58. ∫1/20√1+√1−𝑥2𝑑𝑥

  59. ∫ln⁡𝑥𝑥+𝑥ln⁡𝑥𝑑𝑥

  60. ∫1𝑥⋅ln⁡𝑥⋅ln⁡(ln⁡𝑥)𝑑𝑥

  61. ∫𝑥ln⁡𝑥ln⁡𝑥𝑥𝑑𝑥

  62. ∫(ln⁡𝑥)ln⁡𝑥[1𝑥+ln⁡(ln⁡𝑥)𝑥]𝑑𝑥

  63. ∫1𝑥√1−𝑥4𝑑𝑥

  64. ∫√1−𝑥𝑥𝑑𝑥

  65. ∫sin2⁡𝑥1+sin2⁡𝑥𝑑𝑥

  66. ∫1−cos⁡𝑥1+cos⁡𝑥𝑑𝑥

  67. Evaluate ∫𝜋/20sin⁡𝑥sin⁡𝑥+cos⁡𝑥𝑑𝑥 in two ways:

a. By evaluating ∫sin⁡𝑥sin⁡𝑥+cos⁡𝑥𝑑𝑥 , then using the Evaluation Theorem.

b. By showing that ∫𝑎0𝑓(𝑥)𝑑𝑥 =∫𝑎0𝑓(𝑎 −𝑥)𝑑𝑥 , then using this result.

CHAPTER 8 Additional and Advanced Exercises

Evaluating Integrals

Evaluate the integrals in Exercises 1–6.

  1. ∫(arcsin⁡𝑥)2𝑑𝑥

  2. ∫𝑑𝑥𝑥(𝑥+1)(𝑥+2)⋯(𝑥+𝑚)

  3. ∫𝑥arcsin⁡𝑥𝑑𝑥

  4. ∫sin−1⁡√𝑦𝑑𝑦

  5. ∫𝑑𝑡𝑡−√1−𝑡2

  6. ∫𝑑𝑥𝑥4+4

Evaluate the limits in Exercise 7 and 8.
7. lim𝑥→∞∫𝑥−𝑥sin⁡𝑡𝑑𝑡

  1. lim𝑥→0+𝑥∫1𝑥cos⁡𝑡𝑡2𝑑𝑡

Evaluate the limits in Exercise 9 and 10 by identifying them with definite integrals and evaluating the integrals.

  1. lim𝑛→∞∑𝑛𝑘=1ln⁡𝑛√1+𝑘𝑛

  2. lim𝑛→∞∑𝑛−1𝑘=01√𝑛2−𝑘2

Applications

  1. Finding arc length Find the length of the curve
𝑦=∫𝑥0√cos⁡2𝑡𝑑𝑡,0≤𝑥≤𝜋/4.
  1. Finding arc length Find the length of the graph of the function 𝑦 =ln⁡(1 −𝑥2) , 0 ≤𝑥 ≤1/2 .

  2. Finding volume The region in the first quadrant that is enclosed by the x-axis and the curve 𝑦 =3𝑥√1−𝑥 is revolved about the y-axis to generate a solid. Find the volume of the solid.

  3. Finding volume The region in the first quadrant that is enclosed by the x-axis, the curve 𝑦 =5/(𝑥√5−𝑥) , and the lines x = 1 and x = 4 is revolved about the x-axis to generate a solid. Find the volume of the solid.

  4. Finding volume The region in the first quadrant enclosed by the coordinate axes, the curve 𝑦 =𝑒𝑥 , and the line x = 1 is revolved about the y-axis to generate a solid. Find the volume of the solid.

  5. Finding volume The region in the first quadrant that is bounded above by the curve 𝑦 =𝑒𝑥 −1 , below by the x-axis, and on the right by the line 𝑥 =ln⁡2 is revolved about the line 𝑥 =ln⁡2 to generate a solid. Find the volume of the solid.

  6. Finding volume Let 𝑅 be the “triangular” region in the first quadrant that is bounded above by the line 𝑦 =1 , below by the curve 𝑦 =ln⁡𝑥 , and on the left by the line 𝑥 =1 . Find the volume of the solid generated by revolving 𝑅 about a. the 𝑥 -axis. b. the line 𝑦 =1 .

  7. Finding volume (Continuation of Exercise 17.) Find the volume of the solid generated by revolving the region R about a. the y-axis. b. the line x = 1.

  8. Finding volume The region between the x-axis and the curve

𝑦=𝑓(𝑥)={0,𝑥=0𝑥ln⁡𝑥,0<𝑥≤2

is revolved about the x-axis to generate the solid shown here.

a. Show that 𝑓 is continuous at 𝑥 =0 .

b. Find the volume of the solid.

教材插图

  1. Finding volume The infinite region bounded by the coordinate axes and the curve 𝑦 = −ln⁡𝑥 in the first quadrant is revolved about the x-axis to generate a solid. Find the volume of the solid.

  2. Centroid of a region Find the centroid of the region in the first quadrant that is bounded below by the x-axis, above by the curve 𝑦 =ln⁡𝑥 , and on the right by the line x = e.

  3. Centroid of a region Find the centroid of the region in the plane enclosed by the curves 𝑦 = ±(1 −𝑥2)−1/2 and the lines x = 0 and x = 1.

  4. Length of a curve Find the length of the curve 𝑦 =ln⁡𝑥 from 𝑥 =1 to 𝑥 =𝑒 .

  5. Finding surface area Find the area of the surface generated by revolving the curve in Exercise 23 about the y-axis.

  6. The surface generated by an astroid The graph of the equation 𝑥2/3 +𝑦2/3 =1 is an astroid (see accompanying figure). Find the area of the surface generated by revolving the curve about the 𝑥 -axis.

教材插图

  1. Length of a curve Find the length of the curve
𝑦=∫𝑥1√√𝑡−1𝑑𝑡,1≤𝑥≤16.
  1. For what value or values of a does
∫∞1(𝑎𝑥𝑥2+1−12𝑥)𝑑𝑥

converge? Evaluate the corresponding integral(s).

  1. For each 𝑥 >0 , let 𝐺(𝑥) =∫∞0𝑒−𝑥𝑡𝑑𝑡 . Prove that 𝑥𝐺(𝑥) =1 for each 𝑥 >0 .

  2. Infinite area and finite volume What values of 𝑝 have the following property? The area of the region between the curve 𝑦 =𝑥−𝑝 , 1 ≤𝑥 <∞ , and the 𝑥 -axis is infinite but the volume of the solid generated by revolving the region about the 𝑥 -axis is finite.

  3. Infinite area and finite volume What values of p have the following property? The area of the region in the first quadrant enclosed by the curve 𝑦 =𝑥−𝑝 , the y-axis, the line x = 1, and the interval [0, 1] on the x-axis is infinite, but the volume of the solid generated by revolving the region about one of the coordinate axes is finite.

  4. Integrating the square of the derivative If 𝑓 is continuously differentiable on [0,1] , and 𝑓(1) =𝑓(0) = −1/6 , prove that

∫10(𝑓′(𝑥))2𝑑𝑥≥2∫10𝑓(𝑥)𝑑𝑥+14.

Hint: Consider the inequality 0 ≤∫10(𝑓′(𝑥)+𝑥−12)2𝑑𝑥 .

Source: Mathematics Magazine, vol. 84, no. 4, Oct. 2011.

  1. (Continuation of Exercise 31.) If 𝑓 is continuously differentiable on [0,𝑎] for 𝑎 >0 , and 𝑓(𝑎) =𝑓(0) =𝑏 , prove that
∫𝑎0(𝑓′(𝑥))2𝑑𝑥≥2∫𝑎0𝑓(𝑥)𝑑𝑥−(2𝑎𝑏+𝑎312).

Hint: Consider the inequality 0 ≤∫𝑎0(𝑓′(𝑥)+𝑥−𝑎2)2𝑑𝑥 . Source: Mathematics Magazine, vol. 84, no. 4, Oct. 2011.

The Substitution 𝑧 =tan⁡(𝑥/2) The substitution

𝑧=tan⁡𝑥2(1)

reduces the problem of integrating a rational expression in sin⁡𝑥 and cos⁡𝑥 to a problem of integrating a rational function of z. This in turn can be integrated by partial fractions.

From the accompanying figure

教材插图

we can read the relation

tan⁡𝑥2=sin⁡𝑥1+cos⁡𝑥.

To see the effect of the substitution, we calculate

cos⁡𝑥=2cos2⁡(𝑥2)−1=2sec2⁡(𝑥/2)−1=21+tan2⁡(𝑥/2)−1=21+𝑧2−1cos⁡𝑥=1−𝑧21+𝑧2,(2)

and

sin⁡𝑥=2sin⁡𝑥2cos⁡𝑥2=2sin⁡(𝑥/2)cos⁡(𝑥/2)⋅cos2⁡(𝑥2)=2tan⁡𝑥2⋅1sec2⁡(𝑥/2)=2tan⁡(𝑥/2)1+tan2⁡(𝑥/2)sin⁡𝑥=2𝑧1+𝑧2.(3)

Finally, 𝑥 =2arctan⁡𝑧 , so

𝑑𝑥=2𝑑𝑧1+𝑧2.(4)

Examples

a.

∫11+cos⁡𝑥𝑑𝑥=∫1+𝑧222𝑑𝑧1+𝑧2=∫𝑑𝑧=𝑧+𝐶=tan⁡(𝑥2)+𝐶  b. ∫12+sin⁡𝑥𝑑𝑥=∫1+𝑧22+2𝑧+2𝑧22𝑑𝑧1+𝑧2=∫𝑑𝑧𝑧2+𝑧+1=∫𝑑𝑧(𝑧+(1/2))2+3/4=∫𝑑𝑢𝑢2+𝑎2=1𝑎arctan⁡(𝑢𝑎)+𝐶=2√3arctan⁡2𝑧+1√3+𝐶=2√3arctan⁡1+2tan⁡(𝑥/2)√3+𝐶

Use the substitutions in Equations (1)-(4) to evaluate the integrals in Exercises 33-40. Integrals like these arise in calculating the average angular velocity of the output shaft of a universal joint when the input and output shafts are not aligned.

  1. ∫𝑑𝑥1−sin⁡𝑥

  2. ∫𝑑𝑥1+sin⁡𝑥+cos⁡𝑥

  3. ∫𝜋/20𝑑𝑥1+sin⁡𝑥

  4. ∫𝜋/2𝜋/3𝑑𝑥1−cos⁡𝑥

  5. ∫𝜋/20𝑑𝜃2+cos⁡𝜃

  6. ∫2𝜋/3𝜋/2cos⁡𝜃𝑑𝜃sin⁡𝜃cos⁡𝜃+sin⁡𝜃

  7. ∫𝑑𝑡sin⁡𝑡−cos⁡𝑡

  8. ∫cos⁡𝑡𝑑𝑡1−cos⁡𝑡

Use the substitution 𝑧 =tan⁡(𝜃/2) to evaluate the integrals in Exercises 41 and 42.
41. ∫sec⁡𝜃𝑑𝜃

  1. ∫csc⁡𝜃𝑑𝜃

The Gamma Function and Stirling’s Formula

Euler’s gamma function Γ(𝑥) (“gamma of 𝑥 ”; Γ is a Greek capital 𝑔 ) uses an integral to extend the factorial function from the nonnegative integers to other real values. The formula is

Γ(𝑥)=∫∞0𝑡𝑥−1𝑒−𝑡𝑑𝑡,𝑥>0.

For each positive x, the number Γ(𝑥) is the integral of 𝑡𝑥−1𝑒−𝑡 with respect to t from 0 to ∞ . Figure 8.21 shows the graph of Γ near the origin. You will see how to calculate Γ(1/2) if you do Additional Exercise 23 in Chapter 14.

教材插图

FIGURE 8.21 Euler’s gamma function Γ(𝑥) is a continuous function of 𝑥 whose value at each positive integer 𝑛 +1 is 𝑛! . The defining integral formula for Γ is valid only for 𝑥 >0 , but we can extend Γ to negative noninteger values of 𝑥 with the formula Γ(𝑥) =(Γ(𝑥 +1))/𝑥 , which is the subject of Exercise 43.

  1. If 𝑛 is a nonnegative integer, then Γ(𝑛 +1) =𝑛!

a. Show that Γ(1) =1 .

b. Then apply integration by parts to the integral for Γ(𝑥 +1) to show that Γ(𝑥 +1) =𝑥Γ(𝑥) . This gives

Γ(2)=1Γ(1)=1Γ(3)=2Γ(2)=2Γ(4)=3Γ(3)=6⋮Γ(𝑛+1)=𝑛Γ(𝑛)=𝑛!(1)

c. Use mathematical induction to verify Equation (1) for every nonnegative integer 𝑛 .

  1. Stirling’s formula Scottish mathematician James Stirling (1692-1770) showed that
lim𝑥→∞(𝑒𝑥)𝑥√𝑥2𝜋Γ(𝑥)=1,

so, for large x,

Γ(𝑥)=(𝑥𝑒)𝑥√2𝜋𝑥(1+𝜀(𝑥))𝜀(𝑥)→0 as 𝑥→∞.(2)

Dropping 𝜀(𝑥) leads to the approximation

Γ(𝑥) ≈(𝑥𝑒)𝑥√2𝜋𝑥 (Stirling’s formula).

(3)

a. Stirling’s approximation for 𝑛! Use Equation (3) and the fact that 𝑛! =𝑛Γ(𝑛) to show that 𝑛! ≈(𝑛𝑒)𝑛√2𝑛𝜋 (Stirling’s approximation). (4)

As you will see if you do Exercise 114 in Section 9.1, Equation (4) leads to the approximation

𝑛√𝑛!≈𝑛𝑒.(5)

T b. Compare your calculator’s value for 𝑛! with the value given by Stirling’s approximation for 𝑛 =10,20,30,… , as far as your calculator can go.

T c. A refinement of Equation (2) gives

Γ(𝑥)=(𝑥𝑒)𝑥√2𝜋𝑥𝑒1/(12𝑥)(1+𝜀(𝑥))

or

Γ(𝑥)≈(𝑥𝑒)𝑥√2𝜋𝑥𝑒1/(12𝑥),

which tells us that

𝑛!≈(𝑛𝑒)𝑛√2𝑛𝜋𝑒1/(12𝑛).(6)

Compare the values given for 10! by your calculator, Stirling’s approximation, and Equation (6).

CHAPTER 8 Technology Application Projects

Mathematica/Maple Projects

Projects can be found within MyLab Math.

• Riemann, Trapezoidal, and Simpson Approximations

Part I: Visualize the error involved in using Riemann sums to approximate the area under a curve.

Part II: Build a table of values and compute the relative magnitude of the error as a function of the step size Δ𝑥 .

Part III: Investigate the effect of the derivative function on the error.

Parts IV and V: Trapezoidal Rule approximations.

Part VI: Simpson’s Rule approximations.

  • Games of Chance: Exploring the Monte Carlo Probabilistic Technique for Numerical Integration Graphically explore the Monte Carlo method for approximating definite integrals.

  • Computing Probabilities with Improper Integrals More explorations of the Monte Carlo method for approximating definite integrals.

教材插图

Infinite Sequences and Series

教材插图

OVERVIEW In this chapter we introduce the topic of infinite series. Such series give us precise ways to express many numbers and functions, both familiar and new, as arithmetic sums with infinitely many terms. For example, we will learn that

𝜋4=1−13+15−17+19−…

and

cos⁡𝑥=1−𝑥22+𝑥424−𝑥6720+𝑥840,320−….

We need to develop a method to make sense of such expressions. Everyone knows how to add two numbers together, or even several. But how do you add together infinitely many numbers? Or, when adding together functions, how do you add infinitely many powers of x? In this chapter we answer these questions, which are part of the theory of infinite sequences and series. As with the differential and integral calculus, limits play a major role in the development of infinite series.

One common and important application of series occurs in making computations with complicated functions. A hard-to-compute function is replaced by an expression that looks like an “infinite degree polynomial,” an infinite series in powers of x, as we see with the cosine function given above. Using the first few terms of this infinite series can allow for highly accurate approximations of functions by polynomials, enabling us to work with more general functions than those we have encountered before. These new functions are commonly obtained as solutions to differential equations arising in important applications of mathematics to science and engineering.

The terms “sequence” and “series” are sometimes used interchangeably in spoken language. In mathematics, however, each has a distinct meaning. A sequence is a type of infinite list, whereas a series is an infinite sum. To understand the infinite sums described by series, we first must understand infinite sequences.