书架/Thomas' Calculus

Chapter 14: Multiple Integrals

教材插图

OVERVIEW In this chapter we define the double integral of a function of two variables 𝑓(𝑥,𝑦) over a region in the plane as the limit of approximating Riemann sums. Just as a single integral can represent signed area, so can a double integral represent signed volume. Double integrals can be evaluated using the Fundamental Theorem of Calculus studied in Section 5.4, but now the evaluations are done twice by integrating with respect to each of the variables x and y in turn. Double integrals can be used to find areas of more general regions in the plane than those encountered in Chapter 5. Moreover, just as the Substitution Rule could simplify finding single integrals, we can sometimes use polar coordinates to simplify computing a double integral. We study more general substitutions for evaluating double integrals as well.

We also define the triple integral of a function of three variables 𝑓(𝑥,𝑦,𝑧) over a region in space. Triple integrals can be used to find volumes of still more general regions in space, and their evaluation is like that of double integrals with yet a third evaluation. Cylindrical or spherical coordinates can sometimes be used to simplify the calculation of a triple integral, and we investigate those techniques. Double and triple integrals have a number of applications, such as calculating the average value of a multivariable function, and finding moments and centers of mass.

14.1 Double and Iterated Integrals over Rectangles

教材插图

FIGURE 14.1 Rectangular grid partitioning the region R into small rectangles of area Δ𝐴𝑘 =Δ𝑥𝑘Δ𝑦𝑘 .

In Chapter 5 we defined the definite integral of a function 𝑓(𝑥) over an interval [𝑎,𝑏] as a limit of Riemann sums. In this section we extend this idea to define the double integral of a function of two variables 𝑓(𝑥,𝑦) over a bounded rectangle 𝑅 in the plane. The Riemann sums for the integral of a single-variable function 𝑓(𝑥) are obtained by partitioning a finite interval into thin subintervals, multiplying the width of each subinterval by the value of 𝑓 at a point 𝑐𝑘 inside that subinterval, and then adding together all the products. A similar method of partitioning, multiplying, and summing is used to construct double integrals as limits of approximating Riemann sums.

Double Integrals

We begin our investigation of double integrals by considering the simplest type of planar region, a rectangle. We consider a function 𝑓(𝑥,𝑦) defined on a rectangular region R,

𝑅:𝑎≤𝑥≤𝑏,𝑐≤𝑦≤𝑑.

We subdivide R into small rectangles using a network of lines parallel to the x- and y-axes (Figure 14.1). The lines divide R into n rectangular pieces, where the number of such pieces n gets large as the width and height of each piece gets small. These rectangles form a partition of R. A small rectangular piece of width Δ𝑥 and height Δ𝑦 has area Δ𝐴 =Δ𝑥Δ𝑦 . If we number the small pieces partitioning 𝑅 in some order, then their areas are given by numbers Δ𝐴1,Δ𝐴2,…,Δ𝐴𝑛 , where Δ𝐴𝑘 is the area of the 𝑘 th small rectangle.

FIGURE 14.2 Approximating solids with rectangular boxes leads us to define the volumes of more general solids as double integrals. The volume of the solid shown here is the double integral of 𝑓(𝑥,𝑦) over the base region R.

教材插图

To form a Riemann sum over R, we choose a point (𝑥𝑘,𝑦𝑘) in the kth small rectangle, multiply the value of f at that point by the area Δ𝐴𝑘 , and add together the products:

𝑆𝑛=𝑛∑𝑘=1𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘.

Depending on how we pick (𝑥𝑘,𝑦𝑘) in the 𝑘 th small rectangle, we may get different values for 𝑆𝑛 .

We are interested in what happens to these Riemann sums as the widths and heights of all the small rectangles in the partition of 𝑅 approach zero. The norm of a partition 𝑃 , written ||𝑃|| , is the largest width or height of any rectangle in the partition. If ||𝑃|| =0.1 , then all the rectangles in the partition of 𝑅 have width at most 0.1 and height at most 0.1. Sometimes the Riemann sums converge as the norm of 𝑃 goes to zero, which is written ||𝑃|| →0 . The resulting limit is then written as

lim||𝑃||→0𝑛∑𝑘=1𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘.

As ||𝑃|| →0 and the rectangles get narrow and short, their number 𝑛 increases, so we can also write this limit as

lim𝑛→∞𝑛∑𝑘=1𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘,

with the understanding that ||𝑃|| →0 , and hence Δ𝐴𝑘 →0 , as 𝑛 →∞ .

Many choices are involved in a limit of this kind. The collection of small rectangles is determined by the grid of vertical and horizontal lines that determine a rectangular partition of R. In each of the resulting small rectangles there is a choice of an arbitrary point (𝑥𝑘,𝑦𝑘) at which f is evaluated. These choices together determine a single Riemann sum. To form a limit, we repeat the whole process again and again, choosing partitions whose rectangle widths and heights both go to zero and whose number goes to infinity.

When a limit of the sums 𝑆𝑛 exists, giving the same limiting value no matter what choices are made, then the function f is said to be integrable and the limit is called the double integral of f over R, which is written as

∬𝑅𝑓(𝑥,𝑦)𝑑𝐴 or ∬𝑅𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.

It can be shown that if 𝑓(𝑥,𝑦) is a continuous function throughout R, then f is integrable, as in the single-variable case discussed in Chapter 5. Many discontinuous functions are also integrable, including functions that are discontinuous only on a finite number of points or smooth curves. We leave the proof of these facts to a more advanced text.

Double Integrals as Volumes

When 𝑓(𝑥,𝑦) is a positive function over a rectangular region R in the xy-plane, we may interpret the double integral of f over R as the volume of the three-dimensional solid region over the xy-plane bounded below by R and above by the surface 𝑧 =𝑓(𝑥,𝑦) (Figure 14.2). Each term 𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘 in the sum 𝑆𝑛 =∑𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘 is the volume of a vertical rectangular box that approximates the volume of the portion of the solid that stands directly above the base Δ𝐴𝑘 . The sum 𝑆𝑛 thus approximates what we want to call the total volume of the solid. We define this volume to be

 Volume =lim𝑛→∞𝑆𝑛=∬𝑅𝑓(𝑥,𝑦)𝑑𝐴,

where Δ𝐴𝑘 →0 as 𝑛 →∞ .

教材插图

FIGURE 14.4 To obtain the cross-sectional area 𝐴(𝑥) , we hold 𝑥 fixed and integrate with respect to 𝑦 .

As you might expect, this more general method of calculating volume agrees with the methods in Chapter 6, but we do not prove this here. Figure 14.3 shows Riemann sum approximations to the volume becoming more accurate as the number n of boxes increases.

教材插图

(a) n = 16

教材插图

(b) n = 64

教材插图

(c) n = 256

FIGURE 14.3 As n increases, the Riemann sum approximations approach the total volume of the solid shown in Figure 14.2.

Fubini’s Theorem for Calculating Double Integrals

Suppose that we wish to calculate the volume under the plane z = 4 - x - y over the rectangular region R: 0 ≤𝑥 ≤2 , 0 ≤𝑦 ≤1 in the xy-plane. If we apply the method of slicing from Section 6.1, with slices perpendicular to the x-axis (Figure 14.4), then the volume is

∫𝑥=2𝑥=0𝐴(𝑥)𝑑𝑥,(1)

where 𝐴(𝑥) is the cross-sectional area at x. For each value of x, we may calculate 𝐴(𝑥) as the integral

𝐴(𝑥)=∫𝑦=1𝑦=0(4−𝑥−𝑦)𝑑𝑦,(2)

which is the area under the curve z = 4 - x - y in the plane of the cross-section at x. In calculating 𝐴(𝑥) , x is held fixed and the integration takes place with respect to y. Combining Equations (1) and (2), we see that the volume of the entire solid is

 Volume =∫𝑥=2𝑥=0𝐴(𝑥)𝑑𝑥=∫𝑥=2𝑥=0(∫𝑦=1𝑦=0(4−𝑥−𝑦)𝑑𝑦)𝑑𝑥=∫𝑥=2𝑥=0[4𝑦−𝑥𝑦−𝑦22]𝑦=1𝑦=0𝑑𝑥=∫𝑥=2𝑥=0(72−𝑥)𝑑𝑥=[72𝑥−𝑥22]20=5.

We often omit parentheses separating the two integrals in the formula above and write

 Volume =∫20∫10(4−𝑥−𝑦)𝑑𝑦𝑑𝑥.(3)

The expression on the right, called an iterated or repeated integral, says that the volume is obtained by integrating 4 - x - y with respect to y from y = 0 to y = 1 while holding x fixed, and then integrating the resulting expression in x from x = 0 to x = 2. The limits of integration 0 and 1 are associated with y, so they are placed on the integral closest to dy. The other limits of integration, 0 and 2, are associated with the variable x, so they are placed on the outside integral symbol that is paired with dx.

教材插图

FIGURE 14.5 To obtain the cross-sectional area 𝐴(𝑦) , we hold y fixed and integrate with respect to x.

HISTORICAL BIOGRAPHY Guido Fubini

(1879–1943)

Fubini attended secondary school in Venice, Italy, where he showed that he was brilliant at mathematics. His advanced study was at the Scuola Normale Superiore di Pisa, where his doctoral thesis was in geometry. He then worked on harmonic functions in curved spaces. Fubini’s interests were wide ranging, from differential geometry to analysis and to the applications of differential equations.

To know more, visit the companion Website.

What would have happened if we had calculated the volume by slicing with planes perpendicular to the y-axis (Figure 14.5)? As a function of y, the typical cross-sectional area is

𝐴(𝑦)=∫𝑥=2𝑥=0(4−𝑥−𝑦)𝑑𝑥=[4𝑥−𝑥22−𝑥𝑦]𝑥=2𝑥=0=6−2𝑦.(4)

The volume of the entire solid is therefore

 Volume =∫𝑦=1𝑦=0𝐴(𝑦)𝑑𝑦=∫𝑦=1𝑦=0(6−2𝑦)𝑑𝑦=[6𝑦−𝑦2]10=5,

in agreement with our earlier calculation.

Again, we may give a formula for the volume as an iterated integral by writing

 Volume =∫10∫20(4−𝑥−𝑦)𝑑𝑥𝑑𝑦.

The expression on the right says we can find the volume by integrating 4 - x - y with respect to x from x = 0 to x = 2 as in Equation (4) and integrating the result with respect to y from y = 0 to y = 1. In this iterated integral, the order of integration is first x and then y, the reverse of the order in Equation (3).

What do these two volume calculations with iterated integrals have to do with the double integral

∬𝑅(4−𝑥−𝑦)𝑑𝐴

over the rectangle 𝑅 :0 ≤𝑥 ≤2,0 ≤𝑦 ≤1 ? The answer is that both iterated integrals give the value of the double integral. This is what we would reasonably expect, since the double integral measures the volume of the same region as the two iterated integrals. A theorem published in 1907 by Guido Fubini says that the double integral of any continuous function over a rectangle can be calculated as an iterated integral in either order of integration. (Fubini proved his theorem in greater generality, but this is what it says in our setting.)

THEOREM 1 – Fubini’s Theorem (First Form)

If 𝑓(𝑥,𝑦) is continuous throughout the rectangular region 𝑅 :𝑎 ≤𝑥 ≤𝑏 , 𝑐 ≤𝑦 ≤𝑑 , then

∬𝑅𝑓(𝑥,𝑦)𝑑𝐴=∫𝑑𝑐∫𝑏𝑎𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦=∫𝑏𝑎∫𝑑𝑐𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥.

Fubini’s Theorem says that double integrals over rectangles can be calculated as iterated integrals. Thus, we can evaluate a double integral by integrating with respect to one variable at a time using the Fundamental Theorem of Calculus.

Fubini’s Theorem also says that we may calculate the double integral by integrating in either order, a genuine convenience. When we calculate a volume by slicing, we may use either planes perpendicular to the 𝑥 -axis or planes perpendicular to the 𝑦 -axis.

EXAMPLE 1 Calculate ∬𝑅𝑓(𝑥,𝑦)𝑑𝐴 for

𝑓(𝑥,𝑦)=100−6𝑥2𝑦 and 𝑅:0≤𝑥≤2,−1≤𝑦≤1.

教材插图

∬𝑅𝑓(𝑥,𝑦)𝑑𝐴=∫1−1∫20(100−6𝑥2𝑦)𝑑𝑥𝑑𝑦=∫1−1[100𝑥−2𝑥3𝑦]𝑥=2𝑥=0𝑑𝑦=∫1−1(200−16𝑦)𝑑𝑦=[200𝑦−8𝑦2]1−1=400.

Reversing the order of integration gives the same answer:

Solution Figure 14.6 displays the volume beneath the surface. By Fubini’s Theorem,

FIGURE 14.6 The double integral ∬𝑅𝑓(𝑥,𝑦)𝑑𝐴 gives the volume under this surface over the rectangular region 𝑅 (Example 1).

∫20∫1−1(100−6𝑥2𝑦)𝑑𝑦𝑑𝑥=∫20[100𝑦−3𝑥2𝑦2]𝑦=1𝑦=−1𝑑𝑥=∫20[(100−3𝑥2)−(−100−3𝑥2)]𝑑𝑥=∫20200𝑑𝑥=400.

教材插图

EXAMPLE 2 Find the volume of the region bounded above by the elliptical paraboloid 𝑧 =10 +𝑥2 +3𝑦2 and below by the rectangle 𝑅 :0 ≤𝑥 ≤1,0 ≤𝑦 ≤2 .

FIGURE 14.7 The double integral ∬𝑅𝑓(𝑥,𝑦)𝑑𝐴 gives the volume under this surface over the rectangular region 𝑅 (Example 2).

Solution The surface and volume are shown in Figure 14.7. The volume is given by the double integral

𝑉=∬𝑅(10+𝑥2+3𝑦2)𝑑𝐴=∫10∫20(10+𝑥2+3𝑦2)𝑑𝑦𝑑𝑥=∫10[10𝑦+𝑥2𝑦+𝑦3]𝑦=2𝑦=0𝑑𝑥=∫10(28+2𝑥2)𝑑𝑥=[28𝑥+23𝑥3]10=863.

EXERCISES 14.1

Evaluating Iterated Integrals

In Exercises 1–14, evaluate the iterated integral.

  1. ∫21∫402𝑥𝑦𝑑𝑦𝑑𝑥

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

  3. ∫0−1∫1−1(𝑥 +𝑦 +1)𝑑𝑥𝑑𝑦

  4. ∫10∫10(1−𝑥2+𝑦22)𝑑𝑥𝑑𝑦

  5. ∫30∫20(4 −𝑦2)𝑑𝑦𝑑𝑥

  6. ∫30∫0−2(𝑥2𝑦 −2𝑥𝑦)𝑑𝑦𝑑𝑥

∫𝑐−1∫20(𝑥𝑦+1)𝑑𝑦𝑑𝑥=4+4𝑐.
  1. ∫10∫10𝑦1+𝑥𝑦𝑑𝑥𝑑𝑦

  2. ∫41∫40(𝑥2+√𝑦)𝑑𝑥𝑑𝑦

  3. ∫ln⁡20∫ln⁡51𝑒2𝑥+𝑦𝑑𝑦𝑑𝑥

  4. ∫10∫21𝑥𝑦𝑒𝑥𝑑𝑦𝑑𝑥

  5. ∫2−1∫𝜋/20𝑦sin⁡𝑥𝑑𝑥𝑑𝑦

  6. ∫2𝜋𝜋∫𝜋0(sin⁡𝑥 +cos⁡𝑦)𝑑𝑥𝑑𝑦

  7. ∫41∫𝑒1ln⁡𝑥𝑥𝑦𝑑𝑥𝑑𝑦

  8. ∫2−1∫21𝑥ln⁡𝑦𝑑𝑦𝑑𝑥

  9. Find all values of the constant 𝑐 so that ∫10∫𝑐0(2𝑥 +𝑦)𝑑𝑥𝑑𝑦 =3 .

  10. Find all values of the constant 𝑐 so that

Evaluating Double Integrals over Rectangles

In Exercises 17–24, evaluate the double integral over the given region R.

  1. ∬𝑅(6𝑦2 −2𝑥)𝑑𝐴, 𝑅 :0 ≤𝑥 ≤1, 0 ≤𝑦 ≤2

  2. ∬𝑅(√𝑥𝑦2)𝑑𝐴, 𝑅 :0 ≤𝑥 ≤4,1 ≤𝑦 ≤2

  3. ∬𝑅𝑥𝑦cos⁡𝑦𝑑𝐴, 𝑅 : −1 ≤𝑥 ≤1,0 ≤𝑦 ≤𝜋

  4. ∬𝑅𝑦sin⁡(𝑥 +𝑦)𝑑𝐴, 𝑅 : −𝜋 ≤𝑥 ≤0, 0 ≤𝑦 ≤𝜋

  5. ∬𝑅𝑒𝑥−𝑦𝑑𝐴, 𝑅 :0 ≤𝑥 ≤ln⁡2, 0 ≤𝑦 ≤ln⁡2

  6. ∬𝑅𝑥𝑦𝑒𝑥𝑦2𝑑𝐴, 𝑅 :0 ≤𝑥 ≤2, 0 ≤𝑦 ≤1

  7. ∬𝑅𝑥𝑦3𝑥2+1𝑑𝐴, 𝑅 :0 ≤𝑥 ≤1, 0 ≤𝑦 ≤2

  8. ∬𝑅𝑦𝑥2𝑦2+1𝑑𝐴, 𝑅 :0 ≤𝑥 ≤1, 0 ≤𝑦 ≤1

In Exercises 25 and 26, integrate f over the given region.

  1. Square 𝑓(𝑥,𝑦) =1/(𝑥𝑦) over the square 1 ≤𝑥 ≤2 , 1 ≤𝑦 ≤2

  2. Rectangle 𝑓(𝑥,𝑦) =𝑦cos⁡𝑥𝑦 over the rectangle 0 ≤𝑥 ≤𝜋 , 0 ≤𝑦 ≤1

In Exercises 27 and 28, sketch the solid whose volume is given by the specified integral.

  1. ∫10∫20(9 −𝑥2 −𝑦2)𝑑𝑦𝑑𝑥

  2. ∫30∫41(7 −𝑥 −𝑦)𝑑𝑥𝑑𝑦

  3. Find the volume of the region bounded above by the paraboloid 𝑧 =𝑥2 +𝑦2 and below by the square 𝑅 : −1 ≤𝑥 ≤1 , −1 ≤𝑦 ≤1 .

  4. Find the volume of the region bounded above by the elliptical paraboloid 𝑧 =16 −𝑥2 −𝑦2 and below by the square R: 0 ≤𝑥 ≤2 , 0 ≤𝑦 ≤2 .

  5. Find the volume of the region bounded above by the plane 𝑧 =2 −𝑥 −𝑦 and below by the square 𝑅 :0 ≤𝑥 ≤1,0 ≤𝑦 ≤1 .

  6. Find the volume of the region bounded above by the plane 𝑧 =𝑦/2 and below by the rectangle 𝑅 :0 ≤𝑥 ≤4,0 ≤𝑦 ≤2 .

  7. Find the volume of the region bounded above by the surface 𝑧 =2sin⁡𝑥cos⁡𝑦 and below by the rectangle 𝑅 :0 ≤𝑥 ≤𝜋/2,0 ≤𝑦 ≤𝜋/4.

  8. Find the volume of the region bounded above by the surface 𝑧 =4 −𝑦2 and below by the rectangle 𝑅 :0 ≤𝑥 ≤1,0 ≤𝑦 ≤2 .

  9. Find a value of the constant k so that ∫21∫30𝑘𝑥2𝑦𝑑𝑥𝑑𝑦 =1 .

  10. Evaluate ∫1−1∫𝜋/20𝑥sin⁡√𝑦𝑑𝑦𝑑𝑥.

  11. Use Fubini’s Theorem to evaluate

∫20∫10𝑥1+𝑥𝑦𝑑𝑥𝑑𝑦.38.$𝑈𝑠𝑒𝐹𝑢𝑏𝑖𝑛𝑖′𝑠𝑇ℎ𝑒𝑜𝑟𝑒𝑚𝑡𝑜𝑒𝑣𝑎𝑙𝑢𝑎𝑡𝑒$∫10∫30𝑥𝑒𝑥𝑦𝑑𝑥𝑑𝑦.

T 39. Use a software application to compute the integrals

𝐚.∫10∫20𝑦−𝑥(𝑥+𝑦)3𝑑𝑥𝑑𝑦𝐛.∫20∫10𝑦−𝑥(𝑥+𝑦)3𝑑𝑦𝑑𝑥

Explain why your results do not contradict Fubini’s Theorem.

  1. If 𝑓(𝑥,𝑦) is continuous over 𝑅 : 𝑎 ≤𝑥 ≤𝑏 , 𝑐 ≤𝑦 ≤𝑑 and
𝐹(𝑥,𝑦)=∫𝑥𝑎∫𝑦𝑐𝑓(𝑢,𝑣)𝑑𝑣𝑑𝑢

on the interior of 𝑅 , find the second partial derivatives 𝐹𝑥𝑦 and 𝐹𝑦𝑥 .

14.2 Double Integrals over General Regions

RΔ𝐴𝑘
Δ𝑦𝑘(𝑥𝑘,𝑦𝑘)
Δ𝑥𝑘

FIGURE 14.8 A rectangular grid partitioning a bounded, nonrectangular region into rectangular cells.

In this section we define and evaluate double integrals over bounded regions in the plane that are more general than rectangles. These double integrals are also evaluated as iterated integrals, with the main practical problem being that of determining the limits of integration. Since the region of integration may have boundaries other than line segments parallel to the coordinate axes, the limits of integration often involve variables, not just constants.

Double Integrals over Bounded, Nonrectangular Regions

To define the double integral of a function 𝑓(𝑥,𝑦) over a bounded, nonrectangular region R, such as the one in Figure 14.8, we again begin by covering R with a grid of small rectangular cells whose union contains all points of R. This time, however, we cannot exactly fill R with a finite number of rectangles lying inside R since its boundary is curved, and some of the small rectangles in the grid lie partly outside R. A partition of R is formed by taking the rectangles that lie completely inside it, not using any that are either partly or completely outside. For commonly arising regions, more and more of R is included as the norm of a partition (the largest width or height of any rectangle used) approaches zero.

教材插图

Volume = lim∑𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘 =∬𝑅𝑓(𝑥,𝑦) 𝑑𝐴

FIGURE 14.9 We define the volume of a solid with a curved base as a limit of the sums of volumes of approximating rectangular boxes.

教材插图

FIGURE 14.10 The area of the vertical slice shown here is 𝐴(𝑥) . To calculate the volume of the solid, we integrate this area from 𝑥 =𝑎 to 𝑥 =𝑏 :

∫𝑏𝑎𝐴(𝑥)𝑑𝑥=∫𝑏𝑎∫𝑔2(𝑥)𝑔1(𝑥)𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥.

Once we have a partition of R, we number the rectangles in some order from 1 to n and let Δ𝐴𝑘 be the area of the kth rectangle. We then choose a point (𝑥𝑘,𝑦𝑘) in the kth rectangle and form the Riemann sum

𝑆𝑛=𝑛∑𝑘=1𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘.

As the norm of the partition forming 𝑆𝑛 goes to zero, ‖𝑃‖ →0 , the width and height of each enclosed rectangle go to zero, their area Δ𝐴𝑘 goes to zero, and their number goes to infinity. If 𝑓(𝑥,𝑦) is a continuous function, then these Riemann sums converge to a limiting value that is not dependent on any of the choices we made. This limit is called the double integral of 𝑓(𝑥,𝑦) over R:

lim‖𝑃‖→0𝑛∑𝑘=1𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘=∬𝑅𝑓(𝑥,𝑦)𝑑𝐴.

The nature of the boundary of R introduces issues not found in integrals over an interval. When R has a curved boundary, the n rectangles of a partition lie inside R but do not cover all of R. In order for a partition to approximate R well, the parts of R covered by small rectangles lying partly outside R must become negligible as the norm of the partition approaches zero. This property of being nearly filled in by a partition of small norm is satisfied by all the regions that we will encounter. There is no problem with boundaries made from polygons, circles, and ellipses or from continuous graphs over an interval, joined end to end. A curve with a “fractal” type of shape would be problematic, but such curves arise rarely in most applications. A careful discussion of which types of regions R can be used for computing double integrals is left to a more advanced text.

Volumes

If 𝑓(𝑥,𝑦) is positive and continuous over R, we define the volume of the solid region between R and the surface 𝑧 =𝑓(𝑥,𝑦) to be ∬𝑅𝑓(𝑥,𝑦) 𝑑𝐴 , as before (Figure 14.9).

If R is a region like the one shown in the xy-plane in Figure 14.10, bounded “above” and “below” by the curves 𝑦 =𝑔2(𝑥) and 𝑦 =𝑔1(𝑥) and on the sides by the lines x = a, x = b, we may again calculate the volume by the method of slicing. We first calculate the cross-sectional area

𝐴(𝑥)=∫𝑦=𝑔2(𝑥)𝑦=𝑔1(𝑥)𝑓(𝑥,𝑦)𝑑𝑦

and then integrate 𝐴(𝑥) from 𝑥 =𝑎 to 𝑥 =𝑏 to get the volume as an iterated integral:

𝑉=∫𝑏𝑎𝐴(𝑥)𝑑𝑥=∫𝑏𝑎∫𝑔2(𝑥)𝑔1(𝑥)𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥.(1)

Similarly, if 𝑅 is a region like the one shown in Figure 14.11, bounded by the curves 𝑥 =ℎ2(𝑦) and 𝑥 =ℎ1(𝑦) and the lines 𝑦 =𝑐 and 𝑦 =𝑑 , then the volume calculated by slicing is given by the iterated integral

 Volume =∫𝑑𝑐∫ℎ2(𝑦)ℎ1(𝑦)𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.(2)

That the iterated integrals in Equations (1) and (2) both give the volume that we defined to be the double integral of 𝑓 over 𝑅 is a consequence of the following stronger form of Fubini’s Theorem.

教材插图

FIGURE 14.11 The volume of the solid shown here is

∫𝑑𝑐𝐴(𝑦)𝑑𝑦=∫𝑑𝑐∫ℎ2(𝑦)ℎ1(𝑦)𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.

For a given solid, Theorem 2 says we can calculate the volume as in Figure 14.10 or in the way shown here. Both calculations have the same result.

THEOREM 2—Fubini’s Theorem (Stronger Form)

Let 𝑓(𝑥,𝑦) be continuous on a region R.

  1. If 𝑅 is defined by 𝑎 ≤𝑥 ≤𝑏 , 𝑔1(𝑥) ≤𝑦 ≤𝑔2(𝑥) , with 𝑔1 and 𝑔2 continuous on [𝑎,𝑏] , then
∬𝑅𝑓(𝑥,𝑦)𝑑𝐴=∫𝑏𝑎∫𝑔2(𝑥)𝑔1(𝑥)𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥.
  1. If 𝑅 is defined by 𝑐 ≤𝑦 ≤𝑑 , ℎ1(𝑦) ≤𝑥 ≤ℎ2(𝑦) , with ℎ1 and ℎ2 continuous on [𝑐,𝑑] , then
∬𝑅𝑓(𝑥,𝑦)𝑑𝐴=∫𝑑𝑐∫ℎ2(𝑦)ℎ1(𝑦)𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.

Some iterated double integrals we will encounter later in this text will use variables of integration other than x and y. For instance, we may write

∫𝑢=𝑞𝑢=𝑝∫𝑣=𝐺2(𝑢)𝑣=𝐺1(𝑢)𝐹(𝑢,𝑣)𝑑𝑣𝑑𝑢=∫𝑞𝑝∫𝐺2(𝑢)𝐺1(𝑢)𝐹(𝑢,𝑣)𝑑𝑣𝑑𝑢.

Regardless of which specific variables of integration are used, the limits of an iterated double integral always satisfy these properties:

  • The limits of the outside integral are constants (they do not depend on either variable of integration), and

  • the limits of the inside integral are functions that may depend on the variable of the outside integral.

EXAMPLE 1 Find the volume of the right prism whose base is the triangle in the xy-plane bounded by the x-axis and the lines y = x and x = 1 and whose top lies in the plane

𝑧=𝑓(𝑥,𝑦)=3−𝑥−𝑦.

Solution See Figure 14.12a. For any x between 0 and 1, y may vary from y = 0 to y = x (Figure 14.12b). Hence,

𝑉=∫10∫𝑥0(3−𝑥−𝑦)𝑑𝑦𝑑𝑥=∫10[3𝑦−𝑥𝑦−𝑦22]𝑦=𝑥𝑦=0𝑑𝑥=∫10(3𝑥−3𝑥22)𝑑𝑥=[3𝑥22−𝑥32]𝑥=1𝑥=0=1.

When the order of integration is reversed (Figure 14.12c), the integral for the volume is

𝑉=∫10∫1𝑦(3−𝑥−𝑦)𝑑𝑥𝑑𝑦=∫10[3𝑥−𝑥22−𝑥𝑦]𝑥=1𝑥=𝑦𝑑𝑦=∫10(3−12−𝑦−3𝑦+𝑦22+𝑦2)𝑑𝑦=∫10(52−4𝑦+32𝑦2)𝑑𝑦=[52𝑦−2𝑦2+𝑦32]𝑦=1𝑦=0=1.

The two integrals are equal, as they should be.

Although Fubini’s Theorem assures us that a double integral may be calculated as an iterated integral in either order of integration, the value of one integral may be easier to find than the value of the other. The next example shows how this can happen.

教材插图

(a)

教材插图

(b)

教材插图

(c)

FIGURE 14.12 (a) Prism with a triangular base in the xy-plane. The volume of this prism is defined as a double integral over R. To evaluate it as an iterated integral, we may integrate first with respect to y and then with respect to x, or the other way around (Example 1). (b) Integration limits of

∫𝑥=1𝑥=0∫𝑦=𝑥𝑦=0𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥.

If we integrate first with respect to y, we integrate along a vertical line through R and then integrate from left to right to include all the vertical lines in R. (c) Integration limits of

∫𝑦=1𝑦=0∫𝑥=1𝑥=𝑦𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.

If we integrate first with respect to x, we integrate along a horizontal line through R and then integrate from bottom to top to include all the horizontal lines in R.

教材插图

EXAMPLE 2 Calculate

FIGURE 14.13 The region of integration in Example 2.

∬𝑅sin⁡𝑥𝑥𝑑𝐴,

where R is the triangle in the xy-plane bounded by the x-axis, the line y = x, and the line x = 1.

Solution The region of integration is shown in Figure 14.13. If we integrate first with respect to y and next with respect to x, then because x is held fixed in the first integration, we find

∫10(∫𝑥0sin⁡𝑥𝑥𝑑𝑦)𝑑𝑥=∫10[𝑦sin⁡𝑥𝑥]𝑦=𝑥𝑦=0𝑑𝑥=∫10sin⁡𝑥𝑑𝑥=−cos⁡(1)+1≈0.46.

教材插图

(a)

教材插图

(b)

教材插图

(c)

FIGURE 14.14 Finding the limits of integration when integrating first with respect to y and then with respect to x.

教材插图

FIGURE 14.15 Finding the limits of integration when integrating first with respect to x and then with respect to y.

If we reverse the order of integration and attempt to calculate

∫10∫1𝑦sin⁡𝑥𝑥𝑑𝑥𝑑𝑦,

we run into a problem because ∫((sin⁡𝑥)/𝑥)𝑑𝑥 cannot be expressed in terms of elementary functions (there is no simple antiderivative).

There is no general rule for predicting which order of integration will be the good one in circumstances like these. If the order you first choose doesn’t work, try the other. Sometimes neither order will work, and then we may need to use numerical approximations.

Finding Limits of Integration

We now give a procedure for finding limits of integration that applies for many regions in the plane. Regions that are more complicated, and for which this procedure fails, can often be split up into pieces on which the procedure works.

Using Vertical Cross-Sections When faced with evaluating ∬𝑅𝑓(𝑥,𝑦) 𝑑𝐴 , integrating first with respect to y and then with respect to x, do the following three steps:

  1. Sketch. Sketch the region of integration and label the bounding curves (Figure 14.14a).

  2. Find the y-limits of integration. Imagine a vertical line L cutting through R in the direction of increasing y. Mark the y-values where L enters and leaves. These are the y-limits of integration and are usually functions of x (instead of constants) (Figure 14.14b).

  3. Find the x-limits of integration. Choose x-limits that include all the vertical lines through R. These must be constants. The integral whose region of integration is shown in Figure 14.14c is

∬𝑅𝑓(𝑥,𝑦)𝑑𝐴=∫𝑥=1𝑥=0∫𝑦=√1−𝑥2𝑦=1−𝑥𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥.

Using Horizontal Cross-Sections To evaluate the same double integral as an iterated integral with the order of integration reversed, use horizontal lines instead of vertical lines in Steps 2 and 3 (see Figure 14.15). The integral is

∬𝑅𝑓(𝑥,𝑦)𝑑𝐴=∫10∫√1−𝑦21−𝑦𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.

EXAMPLE 3 Sketch the region of integration for the integral

∫20∫2𝑥𝑥2(4𝑥+2)𝑑𝑦𝑑𝑥

and write an equivalent integral with the order of integration reversed.

Solution The region of integration is given by the inequalities 𝑥2 ≤𝑦 ≤2𝑥 and 0 ≤𝑥 ≤2 . It is therefore the region bounded by the curves 𝑦 =𝑥2 and y = 2x between x = 0 and x = 2 (Figure 14.16a).

To find limits for integrating in the reverse order, we imagine a horizontal line passing from left to right through the region. It enters at x = y/2 and leaves at 𝑥 =√𝑦 . To include all such lines, we let y run from y = 0 to y = 4 (Figure 14.16b). The integral is

∫40∫√𝑦𝑦/2(4𝑥+2)𝑑𝑥𝑑𝑦.

The common value of these integrals is 8.

教材插图

教材插图

FIGURE 14.16 Region of integration for Example 3.

教材插图

FIGURE 14.17 The Additivity Property for rectangular regions holds for regions bounded by smooth curves.

Properties of Double Integrals

Like single integrals, double integrals of continuous functions have algebraic properties that are useful in computations and applications.

If 𝑓(𝑥,𝑦) and 𝑔(𝑥,𝑦) are continuous on the bounded region 𝑅 , then the following properties hold.

  1. Constant Multiple: ∬𝑅𝑐𝑓(𝑥,𝑦)𝑑𝐴 =𝑐∬𝑅𝑓(𝑥,𝑦)𝑑𝐴 (any number 𝑐 )

  2. Sum and Difference:

∬𝑅(𝑓(𝑥,𝑦)±𝑔(𝑥,𝑦))𝑑𝐴=∬𝑅𝑓(𝑥,𝑦)𝑑𝐴±∬𝑅𝑔(𝑥,𝑦)𝑑𝐴
  1. Domination:
∬𝑅𝑓(𝑥,𝑦)𝑑𝐴≥0 if 𝑓(𝑥,𝑦)≥0 on 𝑅 ∬𝑅𝑓(𝑥,𝑦)𝑑𝐴≥∬𝑅𝑔(𝑥,𝑦)𝑑𝐴 if 𝑓(𝑥,𝑦)≥𝑔(𝑥,𝑦) on 𝑅(b)
  1. Additivity: If R is the union of two nonoverlapping regions 𝑅1 and 𝑅2 , then
∬𝑅𝑓(𝑥,𝑦)𝑑𝐴=∬𝑅1𝑓(𝑥,𝑦)𝑑𝐴+∬𝑅2𝑓(𝑥,𝑦)𝑑𝐴

Property 4 assumes that the region of integration R is decomposed into nonoverlapping regions 𝑅1 and 𝑅2 with boundaries consisting of a finite number of line segments or smooth curves. Figure 14.17 illustrates an example of this property.

The idea behind these properties is that integrals behave like sums. If the function 𝑓(𝑥,𝑦) is replaced by its constant multiple 𝑐𝑓(𝑥,𝑦) , then a Riemann sum for f,

𝑆𝑛=𝑛∑𝑘=1𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘,

is replaced by a Riemann sum for cf:

𝑛∑𝑘=1𝑐𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘=𝑐𝑛∑𝑘=1𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘=𝑐𝑆𝑛.

Taking limits as 𝑛 →∞ shows that 𝑐lim𝑛→∞𝑆𝑛 =𝑐∬𝑅𝑓 𝑑𝐴 and lim𝑛→∞𝑐𝑆𝑛 =∬𝑅𝑐𝑓 𝑑𝐴 are equal. It follows that the Constant Multiple Property carries over from sums to double integrals.

The other properties are also easy to verify for Riemann sums, and carry over to double integrals for the same reason. While this discussion gives the idea, an actual proof that these properties hold requires a more careful analysis of how Riemann sums converge.

EXAMPLE 4 Find the volume of the wedgelike solid that lies beneath the surface 𝑧 =16 −𝑥2 −𝑦2 and above the region 𝑅 bounded by the curve 𝑦 =2√𝑥 , the line 𝑦 =4𝑥 −2 , and the 𝑥 -axis.

Solution Figure 14.18a shows the surface and the “wedgelike” solid whose volume we want to calculate. Figure 14.18b shows the region of integration in the xy-plane. If we integrate in the order dy dx (first with respect to y and then with respect to x), two integrations will be required because y varies from y = 0 to 𝑦 =2√𝑥 for 0 ≤𝑥 ≤0.5 , and then varies from y = 4x - 2 to 𝑦 =2√𝑥 for 0.5 ≤𝑥 ≤1 . So we choose to integrate in the order dx dy, which requires only one double integral whose limits of integration are indicated in Figure 14.18b. The volume is then calculated as the iterated integral:

教材插图

(a)

教材插图

(b)

FIGURE 14.18 (a) The solid “wedge-like” region whose volume is found in Example 4. (b) The region of integration R showing the order dx dy.

∬𝑅(16−𝑥2−𝑦2)𝑑𝐴=∫20∫(𝑦+2)/4𝑦2/4(16−𝑥2−𝑦2)𝑑𝑥𝑑𝑦=∫20[16𝑥−𝑥33−𝑥𝑦2]𝑥=(𝑦+2)/4𝑥=𝑦2/4𝑑𝑥=∫20[4(𝑦+2)−(𝑦+2)33⋅64−(𝑦+2)𝑦24−4𝑦2+𝑦63⋅64+𝑦44]𝑑𝑦=[191𝑦24+63𝑦232−145𝑦396−49𝑦4768+𝑦520+𝑦71344]20=208031680≈12.4.

Our development of the double integral has focused on its representation of the volume of the solid region between R and the surface 𝑧 =𝑓(𝑥,𝑦) of a positive continuous function. Just as we saw with signed area in the case of single integrals, when 𝑓(𝑥𝑘,𝑦𝑘) is negative, the product 𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘 is the negative of the volume of the rectangular box shown in Figure 14.9 that was used to form the approximating Riemann sum. So for an arbitrary continuous function f defined over R, the limit of any Riemann sum represents the signed volume (not the total volume) of the solid region between R and the surface. The double integral has other interpretations as well, and in the next section we will see how it is used to calculate the area of a general region in the plane.

EXERCISES 14.2

Sketching Regions of Integration

In Exercises 1–8, sketch the regions of integration associated with the given double integrals.

  1. ∫30∫2𝑥0𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥

  2. ∫2−1∫𝑥2𝑥−1𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥

  3. ∫2−2∫4𝑦2𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦

  4. ∫10∫2𝑦𝑦𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦

  5. ∫10∫𝑒𝑒𝑥𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥

  6. ∫𝑒21∫ln⁡𝑥0𝑓(𝑥,𝑦)𝑑𝑦𝑑𝑥

  7. ∫10∫arcsin⁡𝑦0𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦

  8. ∫80∫𝑦1/3𝑦/4𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦

Finding Limits of Integration

In Exercises 9–18, write an iterated integral for ∬𝑅𝑑𝐴 over the described region R using (a) vertical cross-sections, (b) horizontal cross-sections.

教材插图

教材插图

教材插图

教材插图

  1. Bounded by 𝑦 =√𝑥 , 𝑦 =0 , and 𝑥 =9

  2. Bounded by 𝑦 =tan⁡𝑥 , x = 0, and y = 1

  3. Bounded by 𝑦 =𝑒−𝑥 , 𝑦 =1 , and 𝑥 =ln⁡3

  4. Bounded by 𝑦 =0 , 𝑥 =0 , 𝑦 =1 , and 𝑦 =ln⁡𝑥

  5. Bounded by 𝑦 =3 −2𝑥 , 𝑦 =𝑥 , and 𝑥 =0

  6. Bounded by 𝑦 =𝑥2 and 𝑦 =𝑥 +2

Evaluating Iterated Integrals

In Exercises 19–26, evaluate the integral.

  1. ∫21∫2𝑥0𝑥𝑦3𝑑𝑦𝑑𝑥

  2. ∫31∫2𝑦𝑦𝑦𝑑𝑥𝑑𝑦

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

  4. ∫20∫𝑥30(𝑦2 −𝑥)𝑑𝑦𝑑𝑥

  5. ∫√𝜋0∫𝑥20𝑥sin⁡𝑦𝑑𝑦𝑑𝑥

  6. ∫10∫arctan⁡𝑦011+𝑦2𝑑𝑥𝑑𝑦

  7. ∫41∫𝑦2𝑦√𝑦𝑥𝑑𝑥𝑑𝑦

  8. ∫53∫𝑒𝑥11𝑥𝑦𝑑𝑦𝑑𝑥

Finding Regions of Integration and Double Integrals

In Exercises 27–32, sketch the region of integration and evaluate the integral.

  1. ∫𝜋0∫𝑥0𝑥sin⁡𝑦𝑑𝑦𝑑𝑥

  2. ∫𝜋0∫sin⁡𝑥0𝑦𝑑𝑦𝑑𝑥

  3. ∫ln⁡81∫ln⁡𝑦1𝑒𝑥+𝑦𝑑𝑥𝑑𝑦

  4. ∫21∫𝑦2𝑦𝑑𝑥𝑑𝑦

  5. ∫10∫𝑦203𝑦3𝑒𝑥𝑦𝑑𝑥𝑑𝑦

  6. ∫41∫√𝑥032𝑒𝑦/√𝑥𝑑𝑦𝑑𝑥

In Exercises 33–36, integrate f over the given region.

  1. Quadrilateral 𝑓(𝑥,𝑦) =𝑥/𝑦 over the region in the first quadrant bounded by the lines 𝑦 =𝑥,𝑦 =2𝑥,𝑥 =1 , and 𝑥 =2

  2. Triangle 𝑓(𝑥,𝑦) =𝑥2 +𝑦2 over the triangular region with vertices (0,0),(1,0) , and (0,1)

  3. Triangle 𝑓(𝑢,𝑣) =𝑣 −√𝑢 over the triangular region cut from the first quadrant of the uv-plane by the line 𝑢 +𝑣 =1

  4. Curved region 𝑓(𝑠,𝑡) =𝑒𝑠ln⁡𝑡 over the region in the first quadrant of the st-plane that lies above the curve 𝑠 =ln⁡𝑡 from t = 1 to t = 2

Each of Exercises 37–40 gives an integral over a region in a Cartesian coordinate plane. Sketch the region and evaluate the integral.

  1. ∫0−2∫−𝑣𝑣2𝑑𝑝𝑑𝑣 (the 𝑝𝑣 -plane)

  2. ∫10∫√1−𝑠208𝑡𝑑𝑡𝑑𝑠 (the 𝑠𝑡 -plane)

  3. ∫𝜋/3−𝜋/3∫sec⁡𝑡03cos⁡𝑡𝑑𝑢𝑑𝑡 (the 𝑡𝑢 -plane)

  4. ∫3/20∫4−2𝑢14−2𝑢𝑣2𝑑𝑣𝑑𝑢 (the uv-plane)

Reversing the Order of Integration

In Exercises 41–54, sketch the region of integration, and write an equivalent double integral with the order of integration reversed.

  1. ∫10∫4−2𝑥2𝑑𝑦𝑑𝑥

  2. ∫20∫0𝑦−2𝑑𝑥𝑑𝑦

  3. ∫10∫√𝑦𝑦𝑑𝑥𝑑𝑦

  4. ∫10∫1−𝑥21−𝑥𝑑𝑦𝑑𝑥

  5. ∫10∫𝑒𝑥1𝑑𝑦𝑑𝑥

  6. ∫ln⁡20∫2𝑒𝑦𝑑𝑥𝑑𝑦

  7. ∫3/20∫9−4𝑥2016𝑥𝑑𝑦𝑑𝑥

  8. ∫20∫4−𝑦20𝑦𝑑𝑥𝑑𝑦

  9. ∫10∫√1−𝑦2−√1−𝑦23𝑦𝑑𝑥𝑑𝑦

  10. ∫20∫√4−𝑥2−√4−𝑥26𝑥𝑑𝑦𝑑𝑥

  11. ∫𝑒1∫ln⁡𝑥0𝑥𝑦𝑑𝑦𝑑𝑥

  12. ∫𝜋/60∫1/2sin⁡𝑥𝑥𝑦2𝑑𝑦𝑑𝑥

  13. ∫30∫𝑒𝑦1(𝑥 +𝑦)𝑑𝑥𝑑𝑦

  14. ∫√30∫tan−1⁡𝑦0√𝑥𝑦𝑑𝑥𝑑𝑦

In Exercises 55–64, sketch the region of integration, reverse the order of integration, and evaluate the integral.

  1. ∫𝜋0∫𝜋𝑥sin⁡𝑦𝑦𝑑𝑦𝑑𝑥

  2. ∫20∫2𝑥2𝑦2sin⁡𝑥𝑦𝑑𝑦𝑑𝑥

  3. ∫10∫1𝑦𝑥2𝑒𝑥𝑦𝑑𝑥𝑑𝑦

  4. ∫20∫4−𝑥20𝑥𝑒2𝑦4−𝑦𝑑𝑦𝑑𝑥

  5. ∫2√ln⁡30∫√ln⁡3𝑦/2𝑒𝑥2𝑑𝑥𝑑𝑦

  6. ∫30∫1√𝑥/3𝑒𝑦3𝑑𝑦𝑑𝑥

  7. ∫1/160∫1/2𝑦1/4cos⁡(16𝜋𝑥5)𝑑𝑥𝑑𝑦

  8. ∫80∫23√𝑥𝑑𝑦𝑑𝑥𝑦4+1

  9. Square region ∬𝑅(𝑦 −2𝑥2)𝑑𝐴 where 𝑅 is the region bounded by the square |𝑥| +|𝑦| =1

  10. Triangular region ∬𝑅𝑥𝑦 𝑑𝐴 where 𝑅 is the region bounded by the lines 𝑦 =𝑥 , 𝑦 =2𝑥 , and 𝑥 +𝑦 =2

Volume Beneath a Surface 𝑧 =𝑓(𝑥,𝑦)

  1. Find the volume of the region bounded above by the paraboloid 𝑧 =𝑥2 +𝑦2 and below by the triangle enclosed by the lines y = x, x = 0, and 𝑥 +𝑦 =2 in the xy-plane.

  2. Find the volume of the solid that is bounded above by the cylinder 𝑧 =𝑥2 and below by the region enclosed by the parabola 𝑦 =2 −𝑥2 and the line y = x in the xy-plane.

  3. Find the volume of the solid whose base is the region in the xy-plane that is bounded by the parabola 𝑦 =4 −𝑥2 and the line y = 3x, while the top of the solid is bounded by the plane 𝑧 =𝑥 +4 .

  4. Find the volume of the solid in the first octant bounded by the coordinate planes, the cylinder 𝑥2 +𝑦2 =4 , and the plane 𝑧 +𝑦 =3 .

  5. Find the volume of the solid in the first octant bounded by the coordinate planes, the plane 𝑥 =3 , and the parabolic cylinder 𝑧 =4 −𝑦2 .

  6. Find the volume of the solid cut from the first octant by the surface 𝑧 =4 −𝑥2 −𝑦 .

  7. Find the volume of the wedge cut from the first octant by the cylinder 𝑧 =12 −3𝑦2 and the plane 𝑥 +𝑦 =2 .

  8. Find the volume of the solid cut from the square column |𝑥| +|𝑦| ≤1 by the planes 𝑧 =0 and 3𝑥 +𝑧 =3 .

  9. Find the volume of the solid that is bounded on the front and back by the planes x = 2 and x = 1, on the sides by the cylinders 𝑦 = ±1/𝑥 , and above and below by the planes 𝑧 =𝑥 +1 and z = 0.

  10. Find the volume of the solid bounded on the front and back by the planes 𝑥 = ±𝜋/3 , on the sides by the cylinders 𝑦 = ±sec⁡𝑥 , above by the cylinder 𝑧 =1 +𝑦2 , and below by the xy-plane.

In Exercises 75 and 76, sketch the region of integration and the solid whose volume is given by the double integral.

  1. ∫30∫2−2𝑥/30(1−13𝑥−12𝑦)𝑑𝑦𝑑𝑥

  2. ∫40∫√16−𝑦2−√16−𝑦2√25−𝑥2−𝑦2𝑑𝑥𝑑𝑦

Integrals over Unbounded Regions

Improper double integrals can often be computed similarly to improper integrals of one variable. The first iteration of the following improper integrals is conducted just as if they were proper integrals. One then evaluates an improper integral of a single variable by taking appropriate limits, as in Section 8.8. Evaluate the improper integrals in Exercises 77–80 as iterated integrals.

∫∞1∫1𝑒−𝑥1𝑥3𝑦𝑑𝑦𝑑𝑥77.∫1−1∫1/√1−𝑥2−1/√1−𝑥2(2𝑦+1)𝑑𝑦𝑑𝑥
  1. ∫∞−∞∫∞−∞1(𝑥2+1)(𝑦2+1)𝑑𝑥𝑑𝑦

  2. ∫∞0∫∞0𝑥𝑒−(𝑥+2𝑦)𝑑𝑥𝑑𝑦

Approximating Integrals with Finite Sums

In Exercises 81 and 82, approximate the double integral of 𝑓(𝑥,𝑦) over the region 𝑅 partitioned by the given vertical lines 𝑥 =𝑎 and horizontal lines 𝑦 =𝑐 . In each subrectangle, use (𝑥𝑘,𝑦𝑘) as indicated for your approximation.

∬𝑅𝑓(𝑥,𝑦)𝑑𝐴≈𝑛∑𝑘=1𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘
  1. 𝑓(𝑥,𝑦) =𝑥 +𝑦 over the region R bounded above by the semicircle 𝑦 =√1−𝑥2 and below by the x-axis, using the partition x = -1, -1/2, 0, 1/4, 1/2, 1 and y = 0, 1/2, 1 with (𝑥𝑘,𝑦𝑘) the lower left corner in the kth subrectangle (provided the subrectangle lies within R)

  2. 𝑓(𝑥,𝑦) =𝑥 +2𝑦 over the region R inside the circle (𝑥 −2)2 +(𝑦 −3)2 =1 using the partition x=1,3/2,2,5/2,3 and y=2,5/2,3,7/2,4 with (𝑥𝑘,𝑦𝑘) the center (centroid) in the kth subrectangle (provided the subrectangle lies within R)

Theory and Examples

  1. Circular sector Integrate 𝑓(𝑥,𝑦) =√4−𝑥2 over the smaller sector cut from the disk 𝑥2 +𝑦2 ≤4 by the rays 𝜃 =𝜋/6 and 𝜃 =𝜋/2 .

  2. Unbounded region Integrate 𝑓(𝑥,𝑦) =1/[(𝑥2−𝑥)(𝑦−1)2/3] over the infinite rectangle 2 ≤𝑥 <∞ , 0 ≤𝑦 ≤2 .

  3. Noncircular cylinder A solid right (noncircular) cylinder has its base 𝑅 in the 𝑥𝑦 -plane and is bounded above by the paraboloid 𝑧 =𝑥2 +𝑦2 . The cylinder’s volume is

𝑉=∫10∫𝑦0(𝑥2+𝑦2)𝑑𝑥𝑑𝑦+∫21∫2−𝑦0(𝑥2+𝑦2)𝑑𝑥𝑑𝑦.

Sketch the base region 𝑅 , and express the cylinder’s volume as a single iterated integral with the order of integration reversed. Then evaluate the integral to find the volume.

  1. Converting to a double integral Evaluate the integral
∫20(arctan⁡𝜋𝑥−arctan⁡𝑥)𝑑𝑥.

(Hint: Write the integrand as an integral.)

  1. Maximizing a double integral What region 𝑅 in the xy-plane maximizes the value of
∬𝑅(4−𝑥2−2𝑦2)𝑑𝐴?

Give reasons for your answer.

  1. Minimizing a double integral What region 𝑅 in the xy-plane minimizes the value of
∬𝑅(𝑥2+𝑦2−9)𝑑𝐴?

Give reasons for your answer.

  1. Is it possible to evaluate the integral of a continuous function 𝑓(𝑥,𝑦) over a rectangular region in the 𝑥𝑦 -plane and get different answers depending on the order of integration? Give reasons for your answer.

  2. How would you evaluate the double integral of a continuous function 𝑓(𝑥,𝑦) over the region R in the xy-plane enclosed by the triangle with vertices (0,1) , (2,0) , and (1,2) ? Give reasons for your answer.

  3. Unbounded region Prove that

∫∞−∞∫∞−∞𝑒−𝑥2−𝑦2𝑑𝑥𝑑𝑦=lim𝑏→∞∫𝑏−𝑏∫𝑏−𝑏𝑒−𝑥2−𝑦2𝑑𝑥𝑑𝑦=4(∫∞0𝑒−𝑥2𝑑𝑥)2.92.$𝐼𝑚𝑝𝑟𝑜𝑝𝑒𝑟𝑑𝑜𝑢𝑏𝑙𝑒𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙𝐸𝑣𝑎𝑙𝑢𝑎𝑡𝑒𝑡ℎ𝑒𝑖𝑚𝑝𝑟𝑜𝑝𝑒𝑟𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙$∫10∫30𝑥2(𝑦−1)2/3𝑑𝑦𝑑𝑥.

COMPUTER EXPLORATIONS

Use a CAS double-integral evaluator to estimate the values of the integrals in Exercises 93–96.

∫31∫𝑥11𝑥𝑦𝑑𝑦𝑑𝑥93.∫10∫10𝑒−(𝑥2+𝑦2)𝑑𝑦𝑑𝑥(94)
  1. ∫10∫10arctan⁡𝑥𝑦𝑑𝑦𝑑𝑥

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

Use a CAS double-integral evaluator to find the integrals in Exercises 97–102. Then reverse the order of integration and evaluate, again with a CAS.

  1. ∫10∫42𝑦𝑒𝑥2𝑑𝑥𝑑𝑦

  2. ∫30∫9𝑥2𝑥cos⁡(𝑦2)𝑑𝑦𝑑𝑥

  3. ∫20∫4√2𝑦𝑦3(𝑥2𝑦 −𝑥𝑦2)𝑑𝑥𝑑𝑦

  4. ∫20∫4−𝑦20𝑒𝑥𝑦𝑑𝑥𝑑𝑦

  5. ∫21∫𝑥201𝑥+𝑦𝑑𝑦𝑑𝑥

  6. ∫21∫8𝑦31√𝑥2+𝑦2𝑑𝑥𝑑𝑦

14.3 Area by Double Integration

In this section we show how to use double integrals to calculate the areas of bounded regions in the plane, and to find the average value of a function of two variables.

Areas of Bounded Regions in the Plane

If we take 𝑓(𝑥,𝑦) =1 in the definition of the double integral over a region R in the preceding section, the Riemann sums reduce to

𝑆𝑛=𝑛∑𝑘=1𝑓(𝑥𝑘,𝑦𝑘)Δ𝐴𝑘=𝑛∑𝑘=1Δ𝐴𝑘.(1)

This is simply the sum of the areas of the small rectangles in the partition of R, and it approximates what we would like to call the area of R. As the norm of a partition of R approaches zero, the height and width of all rectangles in the partition approach zero, and the coverage of R becomes increasingly complete (Figure 14.8). We define the area of R to be the limit

lim||𝑃||→0𝑛∑𝑘=1Δ𝐴𝑘=∬𝑅𝑑𝐴.(2)

DEFINITION The area of a closed, bounded plane region R is

𝐴=∬𝑅𝑑𝐴.

As with the other definitions in this chapter, the definition here applies to a greater variety of regions than does the earlier single-variable definition of area, but it agrees with the earlier definition on regions to which they both apply. To evaluate the integral in the definition of area, we integrate the constant function 𝑓(𝑥,𝑦) =1 over R.

教材插图

FIGURE 14.19 The region in Example 1.

EXAMPLE 1 Find the area of the region R bounded by y = x and 𝑦 =𝑥2 in the first quadrant.

Solution We sketch the region (Figure 14.19), noting where the two curves intersect at the origin and (1,1) , and calculate the area as

𝐴=∫10∫𝑥𝑥2𝑑𝑦𝑑𝑥=∫10[𝑦]𝑦=𝑥𝑦=𝑥2𝑑𝑥=∫10(𝑥−𝑥2)𝑑𝑥=[𝑥22−𝑥33]10=16.

Notice that the single-variable integral ∫10(𝑥 −𝑥2)𝑑𝑥 , obtained from evaluating the inside iterated integral, is the integral for the area between these two curves using the method of Section 5.6.

教材插图

教材插图

FIGURE 14.20 Calculating this area takes (a) two double integrals if the first integration is with respect to x, but (b) only one if the first integration is with respect to y (Example 2).

EXAMPLE 2 Find the area of the region R enclosed by the parabola 𝑦 =𝑥2 and the line 𝑦 =𝑥 +2 .

Solution If we divide R into the regions 𝑅1 and 𝑅2 shown in Figure 14.20a, we may calculate the area as

𝐴=∬𝑅1𝑑𝐴+∬𝑅2𝑑𝐴=∫10∫√𝑦−√𝑦𝑑𝑥𝑑𝑦+∫41∫√𝑦𝑦−2𝑑𝑥𝑑𝑦.

On the other hand, reversing the order of integration (Figure 14.20b) gives

𝐴=∫2−1∫𝑥+2𝑥2𝑑𝑦𝑑𝑥.

This second result, which requires only one integral, is simpler to evaluate, giving

𝐴=∫2−1[𝑦]𝑦=𝑥+2𝑦=𝑥2𝑑𝑥=∫2−1(𝑥+2−𝑥2)𝑑𝑥=[𝑥22+2𝑥−𝑥33]2−1=92.

EXAMPLE 3 Find the area of the playing field described by 𝑅 : −2 ≤𝑥 ≤2, −1 −√4−𝑥2 ≤𝑦 ≤1 +√4−𝑥2 , using

(a) Fubini’s Theorem

(b) simple geometry.

Solution The region R is shown in Figure 14.21a.

(a) From the symmetries observed in the figure, we see that the area of 𝑅 is 4 times its area in the first quadrant. As shown in Figure 14.21b, a vertical line at 𝑥 enters this part of the region at 𝑦 =0 and exits at 𝑦 =1 +√4−𝑥2 . Therefore, using Fubini’s Theorem, we have

𝐴=∬𝑅𝑑𝐴=4∫20∫1+√4−𝑥20𝑑𝑦𝑑𝑥=4∫20(1+√4−𝑥2)𝑑𝑥=4[𝑥+𝑥2√4−𝑥2+42sin−1⁡𝑥2]20=4(2+0+2⋅𝜋2−0)=8+4𝜋.

Integral Table Formula 45

(b) The region 𝑅 consists of a rectangle mounted on two sides by half disks of radius 2. The area can be computed by summing the area of the 4 ×2 rectangle and the area of a circle of radius 2, so

𝐴=8+𝜋22=8+4𝜋.

Average Value

The average value of an integrable function of one variable on a closed interval is the integral of the function over the interval divided by the length of the interval. For an integrable function of two variables defined on a bounded region in the plane, the average value is the integral over the region divided by the area of the region. This can be visualized by thinking of the region as being the base of a tank with vertical walls around the boundary of the region, and imagining that the tank is filled with water that is sloshing around. The value 𝑓(𝑥,𝑦) is then the height of the water that is directly above the point (𝑥,𝑦) . The average height of the water in the tank can be found by letting the water settle down to a constant height. This height is equal to the volume of water in the tank divided by the area of R. We therefore define the average value of an integrable function f over a region R as follows:

教材插图

(a)

教材插图

(b)

 Average value of 𝑓 over 𝑅=1 area of 𝑅∬𝑅𝑓𝑑𝐴.(3)

If 𝑓 is the temperature of a thin plate covering 𝑅 , then the double integral of 𝑓 over 𝑅 divided by the area of 𝑅 is the plate’s average temperature. If 𝑓(𝑥,𝑦) is the distance from the point (𝑥,𝑦) to a fixed point 𝑃 , then the average value of 𝑓 over 𝑅 is the average distance of points in 𝑅 from 𝑃 .

EXAMPLE 4 Find the average value of 𝑓(𝑥,𝑦) =𝑥cos⁡𝑥𝑦 over the rectangle 𝑅 :0 ≤𝑥 ≤𝜋,0 ≤𝑦 ≤1 .

Solution The value of the integral of f over R is

∫𝜋0∫10𝑥cos⁡𝑥𝑦𝑑𝑦𝑑𝑥=∫𝜋0[sin⁡𝑥𝑦]𝑦=1𝑦=0𝑑𝑥=∫𝜋0(sin⁡𝑥−0)𝑑𝑥=−cos⁡𝑥∣𝜋0=1+1=2.

FIGURE 14.21 (a) The playing field described by the region R in Example 3.

(b) First quadrant of the playing field.

The area of R is 𝜋 . The average value of f over R is 2/𝜋 .

EXERCISES 14.3

Area by Double Integrals

In Exercises 1–12, sketch the region bounded by the given lines and curves. Then express the region’s area as an iterated double integral and evaluate the integral.

  1. The coordinate axes and the line 𝑥 +𝑦 =2

  2. The lines 𝑥 =0 , 𝑦 =2𝑥 , and 𝑦 =4

  3. The parabola 𝑥 = −𝑦2 and the line 𝑦 =𝑥 +2

  4. The parabola 𝑥 =𝑦 −𝑦2 and the line 𝑦 = −𝑥

  5. The curve 𝑦 =𝑒𝑥 and the lines y = 0, x = 0, and 𝑥 =ln⁡2

  6. The curves 𝑦 =ln⁡𝑥 and 𝑦 =2ln⁡𝑥 and the line 𝑥 =𝑒 , in the first quadrant

  7. The parabolas 𝑥 =𝑦2 and 𝑥 =2𝑦 −𝑦2

  8. 𝑥 =𝑦2 −1 and 𝑥 =2𝑦2 −2

  9. The lines y = x, y = x/3, and y = 2

  10. The lines 𝑦 =1 −𝑥 and 𝑦 =2 and the curve 𝑦 =𝑒𝑥

  11. The lines 𝑦 =2𝑥 , 𝑦 =𝑥/2 , and 𝑦 =3 −𝑥

  12. The lines 𝑦 =𝑥 −2 and 𝑦 = −𝑥 and the curve 𝑦 =√𝑥

Identifying the Region of Integration

The integrals and sums of integrals in Exercises 13–18 give the areas of regions in the xy-plane. Sketch each region, label each bounding curve with its equation, and give the coordinates of the points where the curves intersect. Then find the area of the region.

  1. ∫60∫2𝑦𝑦2/3𝑑𝑥𝑑𝑦

  2. ∫30∫𝑥(2−𝑥)−𝑥𝑑𝑦𝑑𝑥

  3. ∫𝜋/40∫cos⁡𝑥sin⁡𝑥𝑑𝑦𝑑𝑥

  4. ∫2−1∫𝑦+2𝑦2𝑑𝑥𝑑𝑦

  5. ∫0−1∫1−𝑥−2𝑥𝑑𝑦𝑑𝑥 +∫20∫1−𝑥−𝑥/2𝑑𝑦𝑑𝑥

  6. ∫20∫0𝑥2−4𝑑𝑦𝑑𝑥 +∫40∫√𝑥0𝑑𝑦𝑑𝑥

Finding Average Values

  1. Find the average value of 𝑓(𝑥,𝑦) =sin⁡(𝑥 +𝑦) over

a. the rectangle 0 ≤𝑥 ≤𝜋 , 0 ≤𝑦 ≤𝜋 .

b. the rectangle 0 ≤𝑥 ≤𝜋 , 0 ≤𝑦 ≤𝜋/2 .

  1. Which do you think will be larger, the average value of 𝑓(𝑥,𝑦) =𝑥𝑦 over the square 0 ≤𝑥 ≤1 , 0 ≤𝑦 ≤1 , or the average value of f over the quarter circle 𝑥2 +𝑦2 ≤1 in the first quadrant? Calculate them to find out.

  2. Find the average height of the paraboloid 𝑧 =𝑥2 +𝑦2 over the square 0 ≤𝑥 ≤2 , 0 ≤𝑦 ≤2 .

  3. Find the average value of 𝑓(𝑥,𝑦) =1/(𝑥𝑦) over the square ln⁡2 ≤𝑥 ≤2ln⁡2 , ln⁡2 ≤𝑦 ≤2ln⁡2 .

Theory and Examples

  1. Geometric area Find the area of the region
𝑅:0≤𝑥≤2,2−𝑥≤𝑦≤√4−𝑥2,

using (a) Fubini’s Theorem, (b) simple geometry.

  1. Geometric area Find the area of the circular washer with outer radius 2 and inner radius 1, using (a) Fubini’s Theorem, (b) simple geometry.

  2. Bacterium population If 𝑓(𝑥,𝑦) =(10,000𝑒𝑦)/(1 +|𝑥|/2) represents the “population density” of a certain bacterium on the xy-plane, where x and y are measured in centimeters, find the total population of bacteria within the rectangle −5 ≤𝑥 ≤5 and −2 ≤𝑦 ≤0 .

  3. Regional population If 𝑓(𝑥,𝑦) =100(𝑦 +1) represents the population density of a planar region on Earth, where x and y are measured in kilometers, find the number of people in the region bounded by the curves 𝑥 =𝑦2 and 𝑥 =2𝑦 −𝑦2 .

  4. Average temperature in Texas According to the Texas Almanac, Texas has 254 counties and a National Weather Service station in each county. Assume that at time 𝑡0 , each of the 254 weather stations recorded the local temperature. Find a formula that would give a reasonable approximation of the average temperature in Texas at time 𝑡0 . Your answer should involve information that you would expect to be readily available in the Texas Almanac.

  5. If 𝑦 =𝑓(𝑥) is a nonnegative continuous function over the closed interval 𝑎 ≤𝑥 ≤𝑏 , show that the double integral definition of area for the closed plane region bounded by the graph of f, the vertical lines x = a and x = b, and the x-axis agrees with the definition for area beneath the curve in Section 5.3.

  6. Suppose 𝑓(𝑥,𝑦) is continuous over a region 𝑅 in the plane and that the area 𝐴(𝑅) of the region is defined. If there are constants 𝑚 and 𝑀 such that 𝑚 ≤𝑓(𝑥,𝑦) ≤𝑀 for all (𝑥,𝑦) ∈𝑅 , prove that

𝑚𝐴(𝑅)≤∬𝑅𝑓(𝑥,𝑦)𝑑𝐴≤𝑀𝐴(𝑅).
  1. Suppose 𝑓(𝑥,𝑦) is continuous and nonnegative over a region 𝑅 in the plane with a defined area 𝐴(𝑅) . If ∬𝑅𝑓(𝑥,𝑦)𝑑𝐴 =0 , prove that 𝑓(𝑥,𝑦) =0 at every point (𝑥,𝑦) ∈𝑅 .

14.4 Double Integrals in Polar Form

Double integrals are sometimes easier to evaluate if we change to polar coordinates. This section shows how to accomplish the change and how to evaluate double integrals over regions whose boundaries are given by polar equations.

Integrals in Polar Coordinates

When we defined the double integral of a function over a region R in the xy-plane, we began by cutting R into rectangles whose sides were parallel to the coordinate axes. These were the natural shapes to use because their sides have either constant x-values or constant y-values. In polar coordinates, the natural shape is a “polar rectangle” whose sides have constant r- and 𝜃 -values. To avoid ambiguities when describing the region of integration with polar coordinates, we use polar coordinate points (𝑟,𝜃) where 𝑟 ≥0 .

Suppose that a function 𝑓(𝑟,𝜃) is defined over a region R that is bounded by the rays 𝜃 =𝛼 and 𝜃 =𝛽 and by the continuous curves 𝑟 =𝑔1(𝜃) and 𝑟 =𝑔2(𝜃) . Suppose also that 0 ≤𝑔1(𝜃) ≤𝑔2(𝜃) ≤𝑎 for every value of 𝜃 between 𝛼 and 𝛽 . Then R lies in a fan-shaped region Q defined by the inequalities 0 ≤𝑟 ≤𝑎 and 𝛼 ≤𝜃 ≤𝛽 , where 0 ≤𝛽 −𝛼 ≤2𝜋 . See Figure 14.22.

We cover Q by a grid of circular arcs and rays. The arcs are cut from circles centered at the origin, with radii Δ𝑟 , 2Δ𝑟 , …, 𝑚Δ𝑟 , where Δ𝑟 =𝑎/𝑚 . The rays are given by

𝜃=𝛼,𝜃=𝛼+Δ𝜃,𝜃=𝛼+2Δ𝜃,…,𝜃=𝛼+𝑚′Δ𝜃=𝛽,

where Δ𝜃 =(𝛽 −𝛼)/𝑚′ . The arcs and rays partition Q into small patches called “polar rectangles.”

We number the polar rectangles that lie inside R (the order does not matter), calling their areas Δ𝐴1,Δ𝐴2,…,Δ𝐴𝑛 . We let (𝑟𝑘,𝜃𝑘) be any point in the polar rectangle whose area is Δ𝐴𝑘 . We then form the sum

𝑆𝑛=𝑛∑𝑘=1𝑓(𝑟𝑘,𝜃𝑘)Δ𝐴𝑘.

教材插图

教材插图

FIGURE 14.23 The observation that

Δ𝐴𝑘=( area of  large sector )−( area of  small sector )

leads to the formula Δ𝐴𝑘 =𝑟𝑘Δ𝑟Δ𝜃

FIGURE 14.22 The region 𝑅 : 𝑔1(𝜃) ≤𝑟 ≤𝑔2(𝜃) , 𝛼 ≤𝜃 ≤𝛽 , is contained in the fan-shaped region 𝑄 : 0 ≤𝑟 ≤𝑎 , 𝛼 ≤𝜃 ≤𝛽 , where 0 ≤𝛽 −𝛼 ≤2𝜋 . The partition of 𝑄 by circular arcs and rays induces a partition of 𝑅 .

If f is continuous throughout R, this sum will approach a limit as we refine the grid to make Δ𝑟 and Δ𝜃 go to zero. The limit is the double integral of f over R. In symbols,

lim𝑛→∞𝑆𝑛=∬𝑅𝑓(𝑟,𝜃)𝑑𝐴.

To evaluate this limit, we first have to write the sum 𝑆𝑛 in a way that expresses Δ𝐴𝑘 in terms of Δ𝑟 and Δ𝜃 . For convenience we choose 𝑟𝑘 to be the average of the radii of the inner and outer arcs bounding the kth polar rectangle Δ𝐴𝑘 . The radius of the inner arc bounding Δ𝐴𝑘 is then 𝑟𝑘 −(Δ𝑟/2) (Figure 14.23). The radius of the outer arc is 𝑟𝑘 +(Δ𝑟/2) .

The area of a wedge-shaped sector of a circle having radius r and central angle Δ𝜃 is

𝐴=12Δ𝜃⋅𝑟2,

as can be seen by multiplying 𝜋𝑟2 , the area of the circle, by Δ𝜃/2𝜋 , the fraction of the circle’s area contained in the wedge. So the areas of the circular sectors subtended by these arcs at the origin are

Area of small sector:

12(𝑟𝑘−Δ𝑟2)2Δ𝜃

Area of large sector:

12(𝑟𝑘+Δ𝑟2)2Δ𝜃.

Therefore,

Δ𝐴𝑘 = area of large sector - area of small sector

=Δ𝜃2[(𝑟𝑘+Δ𝑟2)2−(𝑟𝑘−Δ𝑟2)2]=Δ𝜃2(2𝑟𝑘Δ𝑟)=𝑟𝑘Δ𝑟Δ𝜃.

Combining this result with the sum defining 𝑆𝑛 gives

𝑆𝑛=𝑛∑𝑘=1𝑓(𝑟𝑘,𝜃𝑘)𝑟𝑘Δ𝑟Δ𝜃.

As 𝑛 →∞ and the values of Δ𝑟 and Δ𝜃 approach zero, these sums converge to the double integral

lim𝑛→∞𝑆𝑛=∬𝑅𝑓(𝑟,𝜃)𝑟𝑑𝑟𝑑𝜃.

教材插图

(a)

教材插图

(b)

教材插图

(c)

FIGURE 14.24 Finding the limits of integration in polar coordinates.

教材插图

FIGURE 14.25 Finding the limits of integration in polar coordinates for the region in Example 1.

Area Differential in Polar Coordinates

𝑑𝐴=𝑟𝑑𝑟𝑑𝜃

A version of Fubini’s Theorem says that the limit approached by these sums can be evaluated by repeated single integrations with respect to 𝑟 and 𝜃 as

∬𝑅𝑓(𝑟,𝜃)𝑑𝐴=∫𝜃=𝛽𝜃=𝛼∫𝑟=𝑔2(𝜃)𝑟=𝑔1(𝜃)𝑓(𝑟,𝜃)𝑟𝑑𝑟𝑑𝜃.

Finding Limits of Integration

The procedure for finding limits of integration in rectangular coordinates also works for polar coordinates. We illustrate this using the region R shown in Figure 14.24. To evaluate ∬𝑅𝑓(𝑟,𝜃) 𝑑𝐴 in polar coordinates, integrating first with respect to r and then with respect to 𝜃 , take the following steps.

  1. Sketch. Sketch the region and label the bounding curves (Figure 14.24a).

  2. Find the r-limits of integration. Imagine a ray L from the origin cutting through R in the direction of increasing r. Mark the r-values where L enters and leaves R. These are the r-limits of integration. They usually depend on the angle 𝜃 that L makes with the positive x-axis (Figure 14.24b).

  3. Find the 𝜃 -limits of integration. Find the smallest and largest 𝜃 -values that bound R. These are the 𝜃 -limits of integration (Figure 14.24c). The polar iterated integral is

∬𝑅𝑓(𝑟,𝜃)𝑑𝐴=∫𝜃=𝜋/2𝜃=𝜋/4∫𝑟=2𝑟=√2cos⁡𝜃𝑓(𝑟,𝜃)𝑟𝑑𝑟𝑑𝜃.

EXAMPLE 1 Find the limits of integration for integrating 𝑓(𝑟,𝜃) over the region R that lies inside the cardioid 𝑟 =1 +cos⁡𝜃 and outside the circle r = 1.

Solution

  1. We first sketch the region and label the bounding curves (Figure 14.25).

  2. Next we find the r-limits of integration. A typical ray from the origin enters R where r = 1 and leaves where 𝑟 =1 +cos⁡𝜃 .

  3. Finally, we find the 𝜃 -limits of integration. The rays from the origin that intersect 𝑅 run from 𝜃 = −𝜋/2 to 𝜃 =𝜋/2 . The integral is

∫𝜋/2−𝜋/2∫1+cos⁡𝜃1𝑓(𝑟,𝜃)𝑟𝑑𝑟𝑑𝜃.

If 𝑓(𝑟,𝜃) is the constant function whose value is 1, then the integral of f over R is the area of R.

Area in Polar Coordinates

The area of a closed and bounded region R in the polar coordinate plane is

𝐴=∬𝑅𝑟𝑑𝑟𝑑𝜃.

This formula for area is consistent with all earlier formulas.

教材插图

FIGURE 14.26 To integrate over the shaded region, we run r from 0 to √4cos⁡2𝜃 and 𝜃 from 0 to 𝜋/4 (Example 2).

教材插图

FIGURE 14.27 The semicircular region in Example 3 is the region

0≤𝑟≤1,0≤𝜃≤𝜋.

EXAMPLE 2 Find the area enclosed by the lemniscate 𝑟2 =4cos⁡2𝜃 .

Solution We graph the lemniscate to determine the limits of integration (Figure 14.26) and see from the symmetry of the region that the total area is 4 times the first-quadrant portion.

𝐴=4∫𝜋/40∫√4cos⁡2𝜃0𝑟𝑑𝑟𝑑𝜃=4∫𝜋/40[𝑟22]𝑟=√4cos⁡2𝜃𝑟=0𝑑𝜃=4∫𝜋/402cos⁡2𝜃𝑑𝜃=4sin⁡2𝜃]𝜋/40=4.

Changing Cartesian Integrals into Polar Integrals

The procedure for changing a Cartesian integral ∬𝑅𝑓(𝑥,𝑦) 𝑑𝑥 𝑑𝑦 into a polar integral has two steps. First substitute 𝑥 =𝑟cos⁡𝜃 and 𝑦 =𝑟sin⁡𝜃 , and replace dx dy by 𝑟 𝑑𝑟 𝑑𝜃 in the Cartesian integral. Then supply polar limits of integration for the boundary of R. The Cartesian integral then becomes

∬𝑅𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦=∬𝐺𝑓(𝑟cos⁡𝜃,𝑟sin⁡𝜃)𝑟𝑑𝑟𝑑𝜃,

where G denotes the same region of integration, but now described in polar coordinates. This is like the substitution method in Chapter 5 except that there are now two variables to substitute for instead of one. Notice that the area differential dx dy is replaced not by dr dθ but by r dr dθ. A more general discussion of changes of variables (substitutions) in multiple integrals is given in Section 14.8.

EXAMPLE 3 Evaluate

∬𝑅𝑒𝑥2+𝑦2𝑑𝑦𝑑𝑥,

where R is the semicircular region bounded by the x-axis and the curve 𝑦 =√1−𝑥2 (Figure 14.27).

Solution In Cartesian coordinates, the integral in question is a nonelementary integral and there is no direct way to integrate 𝑒𝑥2+𝑦2 with respect to either x or y. Yet this integral and others like it are important in mathematics—in statistics, for example—and we need to evaluate it. Polar coordinates make this possible. Substituting x = r cos 𝜃 and 𝑦 =𝑟sin⁡𝜃 and replacing dy dx by r dr d 𝜃 give

∬𝑅𝑒𝑥2+𝑦2𝑑𝑦𝑑𝑥=∫𝜋0∫10𝑒𝑟2𝑟𝑑𝑟𝑑𝜃=∫𝜋0[12𝑒𝑟2]𝑟=1𝑟=0𝑑𝜃=∫𝜋012(𝑒−1)𝑑𝜃=𝜋2(𝑒−1).

The r in the 𝑟 𝑑𝑟 𝑑𝜃 is what allowed us to integrate 𝑒𝑟2 . Without it, we would have been unable to find an antiderivative for the first (innermost) iterated integral.

EXAMPLE 4 Evaluate the integral

∫10∫√1−𝑥20(𝑥2+𝑦2)𝑑𝑦𝑑𝑥.

教材插图

FIGURE 14.28 The solid region in Example 5.

教材插图

FIGURE 14.29 The region R in Example 6.

Solution Integration with respect to y gives

∫10(𝑥2√1−𝑥2+(1−𝑥2)3/23)𝑑𝑥,

which is difficult to evaluate without tables. Things go better if we change the original integral to polar coordinates. The region of integration in Cartesian coordinates is given by the inequalities 0 ≤𝑦 ≤√1−𝑥2 and 0 ≤𝑥 ≤1 , which correspond to the interior of the unit quarter circle 𝑥2 +𝑦2 =1 in the first quadrant. (See Figure 14.27, first quadrant.) Substituting the polar coordinates 𝑥 =𝑟cos⁡𝜃 , 𝑦 =𝑟sin⁡𝜃 , 0 ≤𝜃 ≤𝜋/2 , and 0 ≤𝑟 ≤1 , and replacing dy dx by r dr dθ in the double integral, we get

∫10∫√1−𝑥20(𝑥2+𝑦2)𝑑𝑦𝑑𝑥=∫𝜋/20∫10(𝑟2)𝑟𝑑𝑟𝑑𝜃=∫𝜋/20[𝑟44]𝑟=1𝑟=0𝑑𝜃=∫𝜋/2014𝑑𝜃=𝜋8.

The polar coordinate transformation is effective here because 𝑥2 +𝑦2 simplifies to 𝑟2 and the limits of integration become constants.

EXAMPLE 5 Find the volume of the solid region bounded above by the paraboloid 𝑧 =9 −𝑥2 −𝑦2 and below by the unit circle in the xy-plane.

Solution The region of integration R is bounded by the unit circle 𝑥2 +𝑦2 =1 , which is described in polar coordinates by 𝑟 =1,0 ≤𝜃 ≤2𝜋 . The solid region is shown in Figure 14.28. The volume is given by the double integral

∬𝑅(9−𝑥2−𝑦2)𝑑𝐴=∫2𝜋0∫10(9−𝑟2)𝑟𝑑𝑟𝑑𝜃𝑟2=𝑥2+𝑦2,𝑑𝐴=𝑟𝑑𝑟𝑑𝜃.=∫2𝜋0∫10(9𝑟−𝑟3)𝑑𝑟𝑑𝜃=∫2𝜋0[92𝑟2−14𝑟4]𝑟=1𝑟=0𝑑𝜃=174∫2𝜋0𝑑𝜃=17𝜋2.

EXAMPLE 6 Using polar integration, find the area of the region R enclosed by the circle 𝑥2 +𝑦2 =4 , above the line y = 1, and below the line 𝑦 =√3𝑥 .

Solution A sketch of the region R is shown in Figure 14.29. First we note that the line 𝑦 =√3𝑥 has slope √3 =tan⁡𝜃 , so 𝜃 =𝜋/3 . Next we observe that the line y = 1 intersects the circle 𝑥2 +𝑦2 =4 when 𝑥2 +1 =4 , or 𝑥 =√3 . Moreover, the radial line from the origin through the point (√3,1) has slope 1/√3 =tan⁡𝜃 , giving its angle of inclination as 𝜃 =𝜋/6 . This information is shown in Figure 14.29.

Now, for the region R, as 𝜃 varies from 𝜋/6 to 𝜋/3 , the polar coordinate r varies from the horizontal line y = 1 to the circle 𝑥2 +𝑦2 =4 . Substituting r sin 𝜃 for y in the equation for the horizontal line, we have r sin 𝜃 =1 , or r = csc 𝜃 , which is the polar equation of the line. The polar equation for the circle is r = 2. So in polar coordinates, for 𝜋/6 ≤𝜃 ≤𝜋/3 , r varies from r = csc 𝜃 to r = 2. It follows that the iterated integral for the area is

∬𝑅𝑑𝐴=∫𝜋/3𝜋/6∫2csc⁡𝜃𝑟𝑑𝑟𝑑𝜃=∫𝜋/3𝜋/6[12𝑟2]𝑟=2𝑟=csc⁡𝜃𝑑𝜃=∫𝜋/3𝜋/612[4−csc2⁡𝜃]𝑑𝜃=12[4𝜃+cot⁡𝜃]𝜋/3𝜋/6=12(4𝜋3+1√3)−12(4𝜋6+√3)=𝜋−√33.

EXERCISES 14.4

Regions in Polar Coordinates

In Exercises 1–8, describe the given region in polar coordinates.

教材插图

教材插图

教材插图

教材插图

教材插图

教材插图

  1. The region enclosed by the circle 𝑥2 +𝑦2 =2𝑥

  2. The region enclosed by the semicircle 𝑥2 +𝑦2 =2𝑦 , 𝑦 ≥1

Evaluating Polar Integrals

In Exercises 9–22, change the Cartesian integral into an equivalent polar integral. Then evaluate the polar integral.

  1. ∫1−1∫√1−𝑥20𝑑𝑦𝑑𝑥

  2. ∫10∫√1−𝑦20(𝑥2 +𝑦2)𝑑𝑥𝑑𝑦

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

  4. ∫𝑎−𝑎∫√𝑎2−𝑥2−√𝑎2−𝑥2𝑑𝑦𝑑𝑥

  5. ∫60∫𝑦0𝑥𝑑𝑥𝑑𝑦

  6. ∫20∫𝑥0𝑦𝑑𝑦𝑑𝑥

  7. ∫√31∫𝑥1𝑑𝑦𝑑𝑥

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

  9. ∫0−1∫0−√1−𝑥221+√𝑥2+𝑦2𝑑𝑦𝑑𝑥

  10. ∫1−1∫√1−𝑥2−√1−𝑥22(1+𝑥2+𝑦2)2𝑑𝑦𝑑𝑥

  11. ∫ln⁡20∫√(ln⁡2)2−𝑦20𝑒√𝑥2+𝑦2𝑑𝑥𝑑𝑦

  12. ∫1−1∫√1−𝑦2−√1−𝑦2ln⁡(𝑥2 +𝑦2 +1)𝑑𝑥𝑑𝑦

  13. ∫10∫√2−𝑥2𝑥(𝑥 +2𝑦)𝑑𝑦𝑑𝑥

  14. ∫21∫√2𝑥−𝑥201(𝑥2+𝑦2)2𝑑𝑦𝑑𝑥

In Exercises 23–26, sketch the region of integration, and convert each polar integral or sum of integrals into a Cartesian integral or sum of integrals. Do not evaluate the integrals.

  1. ∫𝜋/20∫10𝑟3sin⁡𝜃cos⁡𝜃𝑑𝑟𝑑𝜃

  2. ∫𝜋/2𝜋/6∫csc⁡𝜃1𝑟2cos⁡𝜃𝑑𝑟𝑑𝜃

  3. ∫𝜋/40∫2sec⁡𝜃0𝑟5sin2⁡𝜃𝑑𝑟𝑑𝜃

  4. ∫arctan⁡430∫3sec⁡𝜃0𝑟7𝑑𝑟𝑑𝜃 +∫𝜋/2arctan⁡43∫4csc⁡𝜃0𝑟7𝑑𝑟𝑑𝜃

Area in Polar Coordinates

  1. Find the area of the region cut from the first quadrant by the curve 𝑟 =2(2 −sin⁡2𝜃)1/2 .

  2. Cardioid overlapping a circle Find the area of the region that lies inside the cardioid 𝑟 =1 +cos⁡𝜃 and outside the circle r = 1.

  3. One leaf of a rose Find the area enclosed by one leaf of the rose 𝑟 =12cos⁡3𝜃 .

  4. Snail shell Find the area of the region enclosed by the positive x-axis and spiral 𝑟 =4𝜃/3 , 0 ≤𝜃 ≤2𝜋 . The region looks like a snail shell.

  5. Cardioid in the first quadrant Find the area of the region cut from the first quadrant by the cardioid 𝑟 =1 +sin⁡𝜃 .

  6. Overlapping cardioids Find the area of the region common to the interiors of the cardioids 𝑟 =1 +cos⁡𝜃 and 𝑟 =1 −cos⁡𝜃 .

Average Values

In polar coordinates, the average value of a function over a region R (Section 14.3) is given by

1 Area (𝑅)∬𝑅𝑓(𝑟,𝜃)𝑟𝑑𝑟𝑑𝜃.
  1. Average height of a hemisphere Find the average height of the hemispherical surface 𝑧 =√𝑎2−𝑥2−𝑦2 above the disk 𝑥2 +𝑦2 ≤𝑎2 in the xy-plane.

  2. Average height of a cone Find the average height of the (single) cone 𝑧 =√𝑥2+𝑦2 above the disk 𝑥2 +𝑦2 ≤𝑎2 in the xy-plane.

  3. Average distance from interior of disk to center Find the average distance from a point 𝑃(𝑥,𝑦) in the disk 𝑥2 +𝑦2 ≤𝑎2 to the origin.

  4. Average distance squared from a point in a disk to a point in its boundary Find the average value of the square of the distance from the point 𝑃(𝑥,𝑦) in the disk 𝑥2 +𝑦2 ≤1 to the boundary point 𝐴(1,0) .

Theory and Examples

  1. Converting to a polar integral Integrate 𝑓(𝑥,𝑦) =[ln⁡(𝑥2 +𝑦2)]/√𝑥2+𝑦2 over the region 1 ≤𝑥2 +𝑦2 ≤𝑒 .

  2. Converting to a polar integral Integrate 𝑓(𝑥,𝑦) =[ln⁡(𝑥2 +𝑦2)]/(𝑥2 +𝑦2) over the region 1 ≤𝑥2 +𝑦2 ≤𝑒2 .

  3. Volume of noncircular right cylinder The region that lies inside the cardioid 𝑟 =1 +cos⁡𝜃 and outside the circle 𝑟 =1 is the base of a solid right cylinder. The top of the cylinder lies in the plane 𝑧 =𝑥 . Find the cylinder’s volume.

  4. Volume of noncircular right cylinder The region enclosed by the lemniscate 𝑟2 =2cos⁡2𝜃 is the base of a solid right cylinder whose top is bounded by the sphere 𝑧 =√2−𝑟2 . Find the cylinder’s volume.

  5. Converting to polar integrals

a. The usual way to evaluate the improper integral 𝐼 =∫∞0𝑒−𝑥2𝑑𝑥 is first to calculate its square:

𝐼2=(∫∞0𝑒−𝑥2𝑑𝑥)(∫∞0𝑒−𝑦2𝑑𝑦)=∫∞0∫∞0𝑒−(𝑥2+𝑦2)𝑑𝑥𝑑𝑦.

Evaluate the last integral using polar coordinates and solve the resulting equation for I.

b. Evaluate

lim𝑥→∞erf⁡(𝑥)=lim𝑥→∞∫𝑥02𝑒−𝑡2√𝜋𝑑𝑡.42.$𝐶𝑜𝑛𝑣𝑒𝑟𝑡𝑖𝑛𝑔𝑡𝑜𝑎𝑝𝑜𝑙𝑎𝑟𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙𝐸𝑣𝑎𝑙𝑢𝑎𝑡𝑒𝑡ℎ𝑒𝑖𝑛𝑡𝑒𝑔𝑟𝑎𝑙$∫∞0∫∞01(1+𝑥2+𝑦2)2𝑑𝑥𝑑𝑦.
  1. Existence Integrate the function 𝑓(𝑥,𝑦) =1/(1 −𝑥2 −𝑦2) over the disk 𝑥2 +𝑦2 ≤3/4 . Does the integral of 𝑓(𝑥,𝑦) over the disk 𝑥2 +𝑦2 ≤1 exist? Give reasons for your answer.

  2. Area formula in polar coordinates Use the double integral in polar coordinates to derive the formula

𝐴=∫𝛽𝛼12𝑟2𝑑𝜃

for the area of the fan-shaped region between the origin and the polar curve 𝑟 =𝑓(𝜃) , 𝛼 ≤𝜃 ≤𝛽 .

  1. Average distance to a given point inside a disk Let 𝑃0 be a point inside a circle of radius 𝑎 and let ℎ denote the distance from 𝑃0 to the center of the circle. Let 𝑑 denote the distance from an arbitrary point 𝑃 to 𝑃0 . Find the average value of 𝑑2 over the region enclosed by the circle. (Hint: Simplify your work by placing the center of the circle at the origin and 𝑃0 on the 𝑥 -axis.)

  2. Area Suppose that the area of a region in the polar coordinate plane is

𝐴=∫3𝜋/4𝜋/4∫2sin⁡𝜃csc⁡𝜃𝑟𝑑𝑟𝑑𝜃.

Sketch the region and find its area.

  1. Evaluate the integral ∬𝑅√𝑥2+𝑦2𝑑𝐴 , where 𝑅 is the region inside the upper semicircle of radius 2 centered at the origin, but outside the circle 𝑥2 +(𝑦 −1)2 =1 .

  2. Evaluate the integral ∬𝑅(𝑥2 +𝑦2)−2𝑑𝐴 , where 𝑅 is the region inside the circle 𝑥2 +𝑦2 =2 for 𝑥 ≤ −1 .

COMPUTER EXPLORATIONS

In Exercises 49–52, use a CAS to change the Cartesian integrals into an equivalent polar integral and evaluate the polar integral. Perform the following steps in each exercise.

a. Plot the Cartesian region of integration in the xy-plane.

b. Change each boundary curve of the Cartesian region in part (a) to its polar representation by solving its Cartesian equation for r and 𝜃 .

c. Using the results in part (b), plot the polar region of integration in the 𝑟𝜃 -plane.

d. Change the integrand from Cartesian to polar coordinates. Determine the limits of integration from your plot in part (c) and evaluate the polar integral using the CAS integration utility.

  1. ∫10∫1𝑥𝑦𝑥2+𝑦2𝑑𝑦𝑑𝑥 50.∫10∫𝑥/20𝑥𝑥2+𝑦2𝑑𝑦𝑑𝑥

  2. ∫10∫𝑦/3−𝑦/3𝑦√𝑥2+𝑦2𝑑𝑥𝑑𝑦

  3. ∫10∫2−𝑦𝑦√𝑥+𝑦𝑑𝑥𝑑𝑦

14.5 Triple Integrals in Rectangular Coordinates

教材插图

FIGURE 14.30 Partitioning a solid with rectangular cells of volume Δ𝑉𝑘 .

Just as double integrals allow us to deal with more general situations than could be handled by single integrals, triple integrals enable us to solve still more general problems. We use triple integrals to calculate the volumes of three-dimensional shapes and the average value of a function over a three-dimensional region. Triple integrals also arise in the study of vector fields and fluid flow in three dimensions, as we will see in Chapter 15.

Triple Integrals

If 𝐹(𝑥,𝑦,𝑧) is a function defined on a closed bounded solid region D in space, such as the region occupied by a solid ball or a lump of clay, then the integral of F over D may be defined in the following way. We partition a rectangular boxlike region containing D into rectangular cells by planes parallel to the coordinate axes (Figure 14.30). We number the cells that lie completely inside D from 1 to n in some order, the kth cell having dimensions Δ𝑥𝑘 by Δ𝑦𝑘 by Δ𝑧𝑘 and volume Δ𝑉𝑘 =Δ𝑥𝑘Δ𝑦𝑘Δ𝑧𝑘 . We choose a point (𝑥𝑘,𝑦𝑘,𝑧𝑘) in each cell and form the sum

𝑆𝑛=𝑛∑𝑘=1𝐹(𝑥𝑘,𝑦𝑘,𝑧𝑘)Δ𝑉𝑘.(1)

We are interested in what happens as D is partitioned by smaller and smaller cells, so that Δ𝑥𝑘,Δ𝑦𝑘,Δ𝑧𝑘 , and the norm of the partition ‖𝑃‖ , the largest value among Δ𝑥𝑘,Δ𝑦𝑘,Δ𝑧𝑘 , all approach zero. When a single limiting value is attained, no matter how the partitions and points (𝑥𝑘,𝑦𝑘,𝑧𝑘) are chosen, we say that F is integrable over D. As before, it can be shown that when F is continuous and the bounding surface of D is formed from finitely many smooth surfaces joined together along finitely many smooth curves, then F is integrable. In this case, as ‖𝑃‖ →0 and the number of cells n goes to ∞ , the sums 𝑆𝑛 approach a limit. We call this limit the triple integral of F over D and write

lim𝑛→∞𝑆𝑛=∭𝐷𝐹(𝑥,𝑦,𝑧)𝑑𝑉 or lim‖𝑃‖→0𝑆𝑛=∭𝐷𝐹(𝑥,𝑦,𝑧)𝑑𝑥𝑑𝑦𝑑𝑧.

The regions D over which continuous functions are integrable are those having “reasonably smooth” boundaries.

Volume of a Solid Region in Space

If 𝐹 is the constant function whose value is 1, then the sums in Equation (1) reduce to

𝑆𝑛=𝑛∑𝑘=1𝐹(𝑥𝑘,𝑦𝑘,𝑧𝑘)Δ𝑉𝑘=𝑛∑𝑘=11⋅Δ𝑉𝑘=𝑛∑𝑘=1Δ𝑉𝑘.

As Δ𝑥𝑘 , Δ𝑦𝑘 , and Δ𝑧𝑘 approach zero, the cells Δ𝑉𝑘 become smaller and more numerous and fill up more and more of D. We therefore define the volume of D to be the triple integral

lim𝑛→∞𝑛∑𝑘=1Δ𝑉𝑘=∭𝐷𝑑𝑉.

DEFINITION The volume of a closed and bounded solid region D in space is

𝑉=∭𝐷𝑑𝑉.

This definition is in agreement with our previous definitions of volume, although we omit the verification of this fact. As we will see in a moment, this integral enables us to calculate the volumes of solids enclosed by curved surfaces. These are more general solids than the ones encountered before (Chapter 6 and Section 14.2).

Iterated Integrals

We evaluate a triple integral by applying a three-dimensional version of Fubini’s Theorem (Section 14.2) to evaluate it by three repeated single integrations. As with double integrals, there is a geometric procedure for finding the limits of integration for these iterated integrals.

To evaluate

∭𝐷𝐹(𝑥,𝑦,𝑧)𝑑𝑉

over a solid region D, integrate first with respect to z, then with respect to y, and finally with respect to x. (You might choose a different order of integration, but the procedure is similar, as we illustrate in Example 2.)

  1. Sketch. Sketch the solid region D along with its “shadow” R (vertical projection) in the xy-plane. Label the upper and lower bounding surfaces of D and the upper and lower bounding curves of R.

教材插图

  1. Find the z-limits of integration. Draw a line M passing through a typical point (𝑥,𝑦) in R parallel to the z-axis. As z increases, M enters D at 𝑧 =𝑓1(𝑥,𝑦) and leaves at 𝑧 =𝑓2(𝑥,𝑦) . These are the z-limits of integration.

教材插图

  1. Find the y-limits of integration. Draw a line L through (𝑥,𝑦) parallel to the y-axis. As y increases, L enters R at 𝑦 =𝑔1(𝑥) and leaves at 𝑦 =𝑔2(𝑥) . These are the y-limits of integration.

教材插图

  1. Find the x-limits of integration. Choose x-limits that include all lines through R parallel to the y-axis (x = a and x = b in the preceding figure). These are the x-limits of integration. The integral is
∫𝑥=𝑏𝑥=𝑎∫𝑦=𝑔2(𝑥)𝑦=𝑔1(𝑥)∫𝑧=𝑓2(𝑥,𝑦)𝑧=𝑓1(𝑥,𝑦)𝐹(𝑥,𝑦,𝑧)𝑑𝑧𝑑𝑦𝑑𝑥.

Follow similar procedures if you change the order of integration. The “shadow” of the solid region D lies in the plane of the last two variables with respect to which the iterated integration takes place. The limits of an iterated triple integral satisfy these properties:

  • The limits of the outside integral are constants (they do not depend on any of the three variables of integration),

  • the limits of the middle integral are functions that may depend on the variable of the outside integral, and

  • the limits of the inside integral are functions that may depend on two variables: the middle integration variable and the outside integration variable.

The preceding procedure applies whenever a solid region D is bounded above and below by a surface, and when the “shadow” region R is bounded by a lower and upper curve. It does not apply to regions with more complicated shapes (such as regions containing holes); although, sometimes such regions can be subdivided into simpler regions for which the procedure does apply.

We illustrate this method of finding the limits of integration in our first example.

EXAMPLE 1 Let S be the sphere of radius 5 centered at the origin, and let D be the solid region under the sphere that lies above the plane z = 3. Set up the limits of integration for evaluating the triple integral of a function 𝐹(𝑥,𝑦,𝑧) over the region D.

Solution The solid region under the sphere that lies above the plane z = 3 is enclosed by the surfaces 𝑥2 +𝑦2 +𝑧2 =25 and z = 3.

To find the limits of integration, we first sketch the solid region, as shown in Figure 14.31. The “shadow region” R in the xy-plane is a circle of some radius centered at the origin. By considering a side view of the region D, we can determine that the radius of this circle is 4; see Figure 14.32a.

If we fix a point (𝑥,𝑦) in R and draw a vertical line M above (𝑥,𝑦) , then we see that this line enters the region D at the height z = 3 and leaves the region at the height 𝑧 =√25−𝑥2−𝑦2 ; see Figure 14.31. This gives us the z-limits of integration.

To find the 𝑦 -limits of integration, we consider a line 𝐿 that lies in the region 𝑅 , passes through the point (𝑥,𝑦) , and is parallel to the 𝑦 -axis. For clarity we have separately pictured the region 𝑅 and the line 𝐿 in Figure 14.32b. The line 𝐿 enters 𝑅 when 𝑦 = −√16−𝑥2 and exits when 𝑦 =√16−𝑥2 . This gives us the 𝑦 -limits of integration.

Finally, as L sweeps across R from left to right, the value of x varies from x = -4 to x = 4. This gives us the x-limits of integration. Therefore, the triple integral of F over the region D is given by

∭𝐷𝐹(𝑥,𝑦,𝑧)𝑑𝑧𝑑𝑦𝑑𝑥=∫4−4∫√16−𝑥2−√16−𝑥2∫√25−𝑥2−𝑦23𝐹(𝑥,𝑦,𝑧)𝑑𝑧𝑑𝑦𝑑𝑥.◼

教材插图

FIGURE 14.31 Finding the limits of integration for evaluating the triple integral of a function defined over the portion of the sphere of radius 5 that lies above the plane z = 3 (Example 1).

(b)

教材插图

教材插图

教材插图

FIGURE 14.33 (a) The tetrahedron in Example 2, showing how the limits of integration are found for the order dz dy dx. (b) The “shadow region” R shown face-on in the xy-plane.

FIGURE 14.32 (a) Side view of the solid region from Example 1, looking down the x-axis. The dashed right triangle has a hypotenuse of length 5 and sides of lengths 3 and 4. In this side view, the shadow region R lies between -4 and 4 on the y-axis. (b) The “shadow region” R shown face-on in the xy-plane.

The region D in Example 1 has a great deal of symmetry, which makes visualization easier. Even without symmetry, the steps in finding the limits of integration are the same, as shown in the next example.

EXAMPLE 2 Set up the limits of integration for evaluating the triple integral of a function 𝐹(𝑥,𝑦,𝑧) over the tetrahedron 𝐷 whose vertices are 𝑂(0,0,0) , 𝐴(1,1,0) , 𝐵(0,1,0) , and 𝐶(0,1,1) . Use the order of integration 𝑑𝑧 𝑑𝑦 𝑑𝑥 .

Solution The solid region D and its “shadow” R in the xy-plane are shown in Figure 14.33a. The “top” face is contained in the plane through the points O, A, and C. Following the procedure introduced in Example 7 of Section 11.5, we first form a normal vector to that plane:

𝐧=⟶𝑂𝐴×⟶𝑂𝐶=∣ ∣ ∣ ∣𝐢𝐣𝐤110011∣ ∣ ∣ ∣=𝐢−𝐣+𝐤,

and then use this vector and the coordinates of O to set up an equation for the plane:

1(𝑥−0)−1(𝑦−0)+1(𝑧−0)=0𝑥−𝑦+𝑧=0.

The “side” face of D is parallel to the xz-plane, the “back” face lies in the yz-plane, and the “bottom” face is contained in the xy-plane.

To find the z-limits of integration, fix a point (𝑥,𝑦) in the shadow region R, and consider the vertical line M that passes through (𝑥,𝑦) and is parallel to the z-axis. This line enters D at the height z = 0, and it exits at height z = y - x.

To find the y-limits of integration we again fix a point (𝑥,𝑦) in R, but now we consider a line L that lies in R, passes through (𝑥,𝑦) , and is parallel to the y-axis. This line is shown in Figure 14.33a and also in the face-on view of R that is pictured in Figure 14.33b. The line L enters R when y = x and exits when y = 1.

Finally, as L sweeps across R, the value of x varies from x = 0 to x = 1. Therefore, the triple integral of F over the region D is given by

∭𝐷𝐹(𝑥,𝑦,𝑧)𝑑𝑧𝑑𝑦𝑑𝑥=∫10∫1𝑥∫𝑦−𝑥0𝐹(𝑥,𝑦,𝑧)𝑑𝑧𝑑𝑦𝑑𝑥.

教材插图

FIGURE 14.34 Finding the limits of integration for evaluating the triple integral of a function defined over the tetrahedron D (Example 3).

In the next example we project the region D onto the xz-plane instead of the xy-plane, to show how to use a different order of integration.

EXAMPLE 3 Find the volume of the tetrahedron D from Example 2 by integrating 𝐹(𝑥,𝑦,𝑧) =1 over the region using the order dz dy dx. Then do the same calculation using the order dy dz dx.

Solution Using the limits of integration that we found in Example 2, we calculate the volume of the tetrahedron as follows:

𝑉=∫10∫1𝑥∫𝑦−𝑥0𝑑𝑧𝑑𝑦𝑑𝑥 Integrand is 1 when computing volume. =∫10∫1𝑥(𝑦−𝑥)𝑑𝑦𝑑𝑥 Integrate over 𝑧=∫10[12𝑦2−𝑥𝑦]𝑦=1𝑦=𝑥𝑑𝑥 and evaluate. =∫10(12−𝑥+12𝑥2)𝑑𝑥 Integrate over 𝑦.=[12𝑥−12𝑥2+16𝑥3]10 Evaluate. =16. Integrate over 𝑥. Evaluate. 

Now we will compute the volume using the order of integration dy dz dx. The procedure for finding the limits of integration is similar, except that we find the limits for y first, then for z, and then for x. The region D is the same tetrahedron as before, but now the “shadow region” R lies in the xz-plane, as shown in Figure 14.34.

To find the y-limits of integration, we fix a point (𝑥,𝑧) in the shadow R and consider the line M that passes through (𝑥,𝑧) and is parallel to the y-axis. As shown in Figure 14.34, this line enters D when 𝑦 =𝑥 +𝑧 , and it leaves when y = 1.

Next we find the z-limits of integration. The line L that passes through a point (𝑥,𝑧) in R and is parallel to the z-axis enters R when z = 0 and exits when z = 1 - x (see Figure 14.34).

Finally, as 𝐿 sweeps across 𝑅 , the value of 𝑥 varies from 𝑥 =0 to 𝑥 =1 . Therefore, the volume of the tetrahedron is

𝑉=∫10∫1−𝑥0∫1𝑥+𝑧𝑑𝑦𝑑𝑧𝑑𝑥=∫10∫1−𝑥0(1−𝑥−𝑧)𝑑𝑧𝑑𝑥=∫10[(1−𝑥)𝑧−12𝑧2]𝑧=1−𝑥𝑧=0𝑑𝑥=∫10[(1−𝑥)2−12(1−𝑥)2]𝑑𝑥=12∫10(1−𝑥)2𝑑𝑥=−16(1−𝑥)3]10=16.

Next we set up and evaluate a triple integral over a more complicated region.

EXAMPLE 4 Find the volume of the solid region D enclosed by the surfaces 𝑧 =𝑥2 +3𝑦2 and 𝑧 =8 −𝑥2 −𝑦2 .

Solution The volume is

𝑉=∭𝐷𝑑𝑧𝑑𝑦𝑑𝑥,

the integral of 𝐹(𝑥,𝑦,𝑧) =1 over D. To find the limits of integration for evaluating the integral, we first sketch the region. The surfaces (Figure 14.35) intersect on the elliptical cylinder 𝑥2 +3𝑦2 =8 −𝑥2 −𝑦2 or 𝑥2 +2𝑦2 =4 , z > 0. The boundary of the region R, the projection of D onto the xy-plane, is an ellipse with the same equation: 𝑥2 +2𝑦2 =4 . The “upper” boundary of R is the curve 𝑦 =√(4−𝑥2)/2 . The lower boundary is the curve 𝑦 = −√(4−𝑥2)/2 .

教材插图

FIGURE 14.35 The volume of the region enclosed by two paraboloids, calculated in Example 4.

Now we find the 𝑧 -limits of integration. The line 𝑀 passing through a typical point (𝑥,𝑦) in 𝑅 parallel to the 𝑧 -axis enters 𝐷 at 𝑧 =𝑥2 +3𝑦2 and leaves at 𝑧 =8 −𝑥2 −𝑦2 .

Next we find the 𝑦 -limits of integration. The line 𝐿 through (𝑥,𝑦) that lies parallel to the 𝑦 -axis enters the region 𝑅 when 𝑦 = −√(4−𝑥2)/2 and leaves when 𝑦 =√(4−𝑥2)/2 .

Finally, we find the x-limits of integration. As L sweeps across R, the value of x varies from 𝑥 = −2 at ( −2,0,0) to 𝑥 =2 at (2,0,0) . The volume of D is

𝑉=∭𝐷𝑑𝑧𝑑𝑦𝑑𝑥 Integrand is 1 when computing volume. =∫2−2∫√(4−𝑥2)/2−√(4−𝑥2)/2∫8−𝑥2−𝑦2𝑥2+3𝑦2𝑑𝑧𝑑𝑦𝑑𝑥 Form an iterated integral.  =∫2−2∫√(4−𝑥2)/2−√(4−𝑥2)/2(8−2𝑥2−4𝑦2)𝑑𝑦𝑑𝑥 Integrate over 𝑧 and evaluate. =∫2−2[(8−2𝑥2)𝑦−43𝑦3]𝑦=√(4−𝑥2)/2𝑦=−√(4−𝑥2)/2𝑑𝑥 Integrate over 𝑦.=∫2−2(2(8−2𝑥2)√4−𝑥22−83(4−𝑥22)3/2)𝑑𝑥 Evaluate. =∫2−2[8(4−𝑥22)3/2−83(4−𝑥22)3/2]𝑑𝑥=4√23∫2−2(4−𝑥2)3/2𝑑𝑥=8𝜋√2. After integration with the substitution 𝑥=2sin⁡𝜃

Average Value of a Function in Space

The average value of a function F over a solid region D in space is defined by the formula

Average value of

𝐹 over 𝐷=1 volume of 𝐷∭𝐷𝐹𝑑𝑉.(2)

For example, if 𝐹(𝑥,𝑦,𝑧) =√𝑥2+𝑦2+𝑧2 , then the average value of F over D is the average distance of points in D from the origin. If 𝐹(𝑥,𝑦,𝑧) is the temperature at (𝑥,𝑦,𝑧) on a solid that occupies a region D in space, then the average value of F over D is the average temperature of the solid.

EXAMPLE 5 Find the average value of 𝐹(𝑥,𝑦,𝑧) =𝑥𝑦𝑧 throughout the cubical region D bounded by the coordinate planes and the planes x = 2, y = 2, and z = 2 in the first octant.

教材插图

Solution We sketch the cube with enough detail to show the limits of integration (Figure 14.36). We then use Equation (2) to calculate the average value of F over the cube.

FIGURE 14.36 The region of integration in Example 5.

The volume of the region 𝐷 is (2)(2)(2) = 8. The value of the integral of 𝐹 over the cube is

∫20∫20∫20𝑥𝑦𝑧𝑑𝑥𝑑𝑦𝑑𝑧=∫20∫20[𝑥22𝑦𝑧]𝑥=2𝑥=0𝑑𝑦𝑑𝑧=∫20∫202𝑦𝑧𝑑𝑦𝑑𝑧=∫20[𝑦2𝑧]𝑦=2𝑦=0𝑑𝑧=∫204𝑧𝑑𝑧=[2𝑧2]20=8.

With these values, Equation (2) gives

 Average value of 𝑥𝑦𝑧 over the cube =1 volume ∭ cube 𝑥𝑦𝑧𝑑𝑉=(18)(8)=1.

In evaluating the integral, we chose the order dx dy dz, but any of the other five possible orders would have done as well.

Properties of Triple Integrals

Triple integrals have the same algebraic properties as double and single integrals. Simply replace the double integrals in the four properties given in Section 14.2, page 864, with triple integrals.

EXERCISES 14.5

Triple Integrals in Different Iteration Orders

  1. Evaluate the integral in Example 3, taking 𝐹(𝑥,𝑦,𝑧) =1 to find the volume of the tetrahedron in the order 𝑑𝑧 𝑑𝑥 𝑑𝑦 .

  2. Volume of rectangular solid Write six different iterated triple integrals for the volume of the rectangular solid in the first octant bounded by the coordinate planes and the planes x = 1, y = 2, and z = 3. Evaluate one of the integrals.

  3. Volume of tetrahedron Write six different iterated triple integrals for the volume of the tetrahedron cut from the first octant by the plane 6𝑥 +3𝑦 +2𝑧 =6 . Evaluate one of the integrals.

  4. Volume of solid Write six different iterated triple integrals for the volume of the solid region in the first octant enclosed by the cylinder 𝑥2 +𝑧2 =4 and the plane y = 3. Evaluate one of the integrals.

  5. Volume enclosed by paraboloids Let 𝐷 be the solid region bounded by the paraboloids 𝑧 =8 −𝑥2 −𝑦2 and 𝑧 =𝑥2 +𝑦2 . Write six different triple iterated integrals for the volume of 𝐷 . Evaluate one of the integrals.

  6. Volume inside paraboloid beneath a plane Let D be the solid region bounded by the paraboloid 𝑧 =𝑥2 +𝑦2 and the plane z = 2y. Write triple iterated integrals in the order dz dx dy and dz dy dx that give the volume of D. Do not evaluate either integral.

Evaluating Triple Iterated Integrals

Evaluate the integrals in Exercises 7–20.

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

  2. ∫√20∫3𝑦0∫8−𝑥2−𝑦2𝑥2+3𝑦2𝑑𝑧𝑑𝑥𝑑𝑦

  3. ∫𝑒1∫𝑒21∫𝑒311𝑥𝑦𝑧𝑑𝑥𝑑𝑦𝑑𝑧

  4. ∫10∫3−3𝑥0∫3−3𝑥−𝑦0𝑑𝑧𝑑𝑦𝑑𝑥

  5. ∫𝜋/60∫10∫3−2𝑦sin⁡𝑧𝑑𝑥𝑑𝑦𝑑𝑧

  6. ∫1−1∫10∫20(𝑥 +𝑦 +𝑧)𝑑𝑦𝑑𝑥𝑑𝑧

  7. ∫30∫√9−𝑥20∫√9−𝑥20𝑑𝑧𝑑𝑦𝑑𝑥

  8. ∫20∫√4−𝑦2−√4−𝑦2∫2𝑥+𝑦0𝑑𝑧𝑑𝑥𝑑𝑦

  9. ∫10∫2−𝑥0∫2−𝑥−𝑦0𝑑𝑧𝑑𝑦𝑑𝑥

  10. ∫10∫1−𝑥20∫4−𝑥2−𝑦3𝑥𝑑𝑧𝑑𝑦𝑑𝑥

  11. ∫𝜋0∫𝜋0∫𝜋0cos⁡(𝑢 +𝑣 +𝑤)𝑑𝑢𝑑𝑣𝑑𝑤 (uvw-space)

  12. ∫10∫√𝑒1∫𝑒1𝑠𝑒𝑠ln⁡𝑟(ln⁡𝑡)2𝑡𝑑𝑡𝑑𝑟𝑑𝑠 (rst-space)

  13. ∫𝜋/40∫ln⁡sec⁡𝑣0∫2𝑡−∞𝑒𝑥𝑑𝑥𝑑𝑡𝑑𝑣 (tvx-space)

  14. ∫70∫20∫√4−𝑞20𝑞𝑟+1𝑑𝑝𝑑𝑞𝑑𝑟 (pqr-space)

Finding Equivalent Iterated Integrals

  1. Here is the region of integration of the integral

教材插图

教材插图

Rewrite the integral as an equivalent iterated integral in the order

a. dy dz dx b. dy dx dz

c. 𝑑𝑥 dy dz d. 𝑑𝑥 dz dy

e. dz dx dy.

  1. Here is the region of integration of the integral

教材插图

教材插图

Rewrite the integral as an equivalent iterated integral in the order

a. dy dz dx b. dy dx dz

c. 𝑑𝑥 dy dz d. 𝑑𝑥 dz dy

e. dz dx dy.

Finding Volumes Using Triple Integrals

Find the volumes of the solid regions in Exercises 23–36.

  1. The region between the cylinder 𝑧 =𝑦2 and the 𝑥𝑦 -plane that is bounded by the planes 𝑥 =0,𝑥 =1,𝑦 = −1,𝑦 =1

教材插图

  1. The region in the first octant bounded by the coordinate planes and the planes 𝑥 +𝑧 =1 , 𝑦 +2𝑧 =2

教材插图

  1. The region in the first octant bounded by the coordinate planes, the plane 𝑦 +𝑧 =2 , and the cylinder 𝑥 =4 −𝑦2

教材插图

  1. The wedge cut from the cylinder 𝑥2 +𝑦2 =1 with 𝑧 ≥0 by the planes z = -y and z = 0

教材插图

  1. The tetrahedron in the first octant bounded by the coordinate planes and the plane passing through (1,0,0) , (0,2,0) , and (0,0,3)

教材插图

  1. The region in the first octant bounded by the coordinate planes, the plane 𝑦 =1 −𝑥 , and the surface 𝑧 =cos⁡(𝜋𝑥/2) , 0 ≤𝑥 ≤1

教材插图

  1. The region common to the interiors of the cylinders 𝑥2 +𝑦2 =1 and 𝑥2 +𝑧2 =1 , one-eighth of which is shown in the accompanying figure

教材插图

  1. The region in the first octant bounded by the coordinate planes and the surface 𝑧 =4 −𝑥2 −𝑦

教材插图

  1. The region in the first octant bounded by the coordinate planes, the plane 𝑥 +𝑦 =4 , and the cylinder 𝑦2 +4𝑧2 =16

教材插图

  1. The region cut from the cylinder 𝑥2 +𝑦2 =4 by the plane z = 0 and the plane 𝑥 +𝑧 =3

教材插图

  1. The region between the planes 𝑥 +𝑦 +2𝑧 =2 and 2𝑥 +2𝑦 +𝑧 =4 in the first octant

  2. The finite region bounded by the planes 𝑧 =𝑥 , 𝑥 +𝑧 =8 , 𝑧 =𝑦 , 𝑦 =8 , and 𝑧 =0

  3. The region cut from the solid elliptical cylinder 𝑥2 +4𝑦2 ≤4 by the 𝑥𝑦 -plane and the plane 𝑧 =𝑥 +2

  4. The region bounded in back by the plane x = 0, on the front and sides by the parabolic cylinder 𝑥 =1 −𝑦2 , on the top by the paraboloid 𝑧 =𝑥2 +𝑦2 , and on the bottom by the xy-plane

Average Values

In Exercises 37–40, find the average value of 𝐹(𝑥,𝑦,𝑧) over the given region.

  1. 𝐹(𝑥,𝑦,𝑧) =𝑥2 +9 over the cube in the first octant bounded by the coordinate planes and the planes x = 2, y = 2, and z = 2

  2. 𝐹(𝑥,𝑦,𝑧) =𝑥 +𝑦 −𝑧 over the rectangular box in the first octant bounded by the coordinate planes and the planes x = 1, y = 1, and z = 2

  3. 𝐹(𝑥,𝑦,𝑧) =𝑥2 +𝑦2 +𝑧2 over the cube in the first octant bounded by the coordinate planes and the planes 𝑥 =1 , 𝑦 =1 , and 𝑧 =1

  4. 𝐹(𝑥,𝑦,𝑧) =𝑥𝑦𝑧 over the cube in the first octant bounded by the coordinate planes and the planes x = 2, y = 2, and z = 2

Changing the Order of Integration

Evaluate the integrals in Exercises 41–44 by changing the order of integration in an appropriate way.

  1. ∫40∫10∫22𝑦4cos⁡(𝑥2)2√𝑧𝑑𝑥𝑑𝑦𝑑𝑧

  2. ∫10∫10∫1𝑥212𝑥𝑧𝑒𝑧𝑦2𝑑𝑦𝑑𝑥𝑑𝑧

  3. ∫10∫13√𝑧∫ln⁡30𝜋𝑒2𝑥sin⁡𝜋𝑦2𝑦2𝑑𝑥𝑑𝑦𝑑𝑧

  4. ∫20∫4−𝑥20∫𝑥0sin⁡2𝑧4−𝑧𝑑𝑦𝑑𝑧𝑑𝑥

Theory and Examples

  1. Finding an upper limit of an iterated integral Solve for a:
∫10∫4−𝑎−𝑥20∫4−𝑥2−𝑦𝑎𝑑𝑧𝑑𝑦𝑑𝑥=415.
  1. Ellipsoid For what value of 𝑐 is the volume of the ellipsoid 𝑥2 +(𝑦/2)2 +(𝑧/𝑐)2 =1 equal to 8𝜋 ?

  2. Minimizing a triple integral What domain D in space minimizes the value of the integral

∭𝐷(4𝑥2+4𝑦2+𝑧2−4)𝑑𝑉?

Give reasons for your answer.

  1. Maximizing a triple integral What domain D in space maximizes the value of the integral
∭𝐷(1−𝑥2−𝑦2−𝑧2)𝑑𝑉?

Give reasons for your answer.

COMPUTER EXPLORATIONS

In Exercises 49–52, use a CAS integration utility to evaluate the triple integral of the given function over the specified solid region.

  1. 𝐹(𝑥,𝑦,𝑧) =𝑥2𝑦2𝑧 over the solid cylinder bounded by 𝑥2 +𝑦2 =1 and the planes 𝑧 =0 and 𝑧 =1

  2. 𝐹(𝑥,𝑦,𝑧) =|𝑥𝑦𝑧| over the solid bounded below by the paraboloid 𝑧 =𝑥2 +𝑦2 and above by the plane z = 1

  3. 𝐹(𝑥,𝑦,𝑧) =𝑧(𝑥2+𝑦2+𝑧2)3/2 over the solid bounded below

by the cone 𝑧 =√𝑥2+𝑦2 and above by the plane z = 1

  1. 𝐹(𝑥,𝑦,𝑧) =𝑥4 +𝑦2 +𝑧2 over the solid sphere 𝑥2 +𝑦2 +𝑧2 ≤1

14.6 Applications

教材插图

FIGURE 14.37 To define an object’s mass, we first imagine it to be partitioned into a finite number of mass elements Δ𝑚𝑘 .

This section shows how to calculate the masses and moments of two- and three-dimensional objects in Cartesian coordinates. The definitions and ideas are similar to the single-variable case we studied in Section 6.6, but now we can consider more general situations.

Masses and First Moments

If 𝛿(𝑥,𝑦,𝑧) is the density (mass per unit volume) of an object occupying a solid region 𝐷 in space, the integral of 𝛿 over 𝐷 gives the mass of the object. To see why, imagine partitioning the object into 𝑛 mass elements like the one in Figure 14.37. The object’s mass is the limit

𝑀=lim𝑛→∞𝑛∑𝑘=1Δ𝑚𝑘=lim𝑛→∞𝑛∑𝑘=1𝛿(𝑥𝑘,𝑦𝑘,𝑧𝑘)Δ𝑉𝑘=∭𝐷𝛿(𝑥,𝑦,𝑧)𝑑𝑉.

The first moment of a solid region D about a coordinate plane is defined as the triple integral over D of the (signed) distance from a point (𝑥,𝑦,𝑧) in D to the plane multiplied by the density of the solid at that point. For instance, the first moment about the yz-plane is the integral

𝑀𝑦𝑧=∭𝐷𝑥𝛿(𝑥,𝑦,𝑧)𝑑𝑉.

The center of mass is found from the first moments. For instance, the 𝑥 -coordinate of the center of mass is ――𝑥 =𝑀𝑦𝑧/𝑀 .

For a two-dimensional object, such as a thin, flat plate, we calculate first moments about the coordinate axes by simply dropping the z-coordinate. So the first moment about the y-axis is the double integral over the region R forming the plate of the (signed) distance from the axis multiplied by the density, or

𝑀𝑦=∬𝑅𝑥𝛿(𝑥,𝑦)𝑑𝐴.

Table 14.1 summarizes the formulas.

TABLE 14.1 Mass and first moment formulas

THREE-DIMENSIONAL SOLID

Mass: 𝑀 =∭𝐷𝛿𝑑𝑉 𝛿 =𝛿(𝑥,𝑦,𝑧) is the density at (𝑥,𝑦,𝑧) .

First moments about the coordinate planes:

𝑀𝑦𝑧=∭𝐷𝑥𝛿𝑑𝑉,𝑀𝑥𝑧=∭𝐷𝑦𝛿𝑑𝑉,𝑀𝑥𝑦=∭𝐷𝑧𝛿𝑑𝑉

Center of mass: ――𝑥 =𝑀𝑦𝑧𝑀, ――𝑦 =𝑀𝑥𝑧𝑀, ――𝑧 =𝑀𝑥𝑦𝑀

TWO-DIMENSIONAL PLATE

Mass: 𝑀 =∬𝑅𝛿𝑑𝐴 𝛿 =𝛿(𝑥,𝑦) is the density at (𝑥,𝑦) .

First moments: 𝑀𝑦 =∬𝑅𝑥𝛿𝑑𝐴, 𝑀𝑥 =∬𝑅𝑦𝛿𝑑𝐴

Center of mass: ――𝑥 =𝑀𝑦𝑀, ――𝑦 =𝑀𝑥𝑀

教材插图

FIGURE 14.38 Finding the center of mass of a solid (Example 1).

EXAMPLE 1 Find the center of mass of a solid of constant density 𝛿 bounded below by the disk 𝑅 :𝑥2 +𝑦2 ≤4 in the plane z = 0 and above by the paraboloid 𝑧 =4 −𝑥2 −𝑦2 (Figure 14.38).

Solution By symmetry ――𝑥 =――𝑦 =0 . To find ――𝑧 , we first calculate

𝑀𝑥𝑦=∬𝑅∫𝑧=4−𝑥2−𝑦2𝑧=0𝑧𝛿𝑑𝑧𝑑𝑦𝑑𝑥=∬𝑅[𝑧22]𝑧=4−𝑥2−𝑦2𝑧=0𝛿𝑑𝑦𝑑𝑥=𝛿2∬𝑅(4−𝑥2−𝑦2)2𝑑𝑦𝑑𝑥 =𝛿2∫2𝜋0∫20(4−𝑟2)2𝑟𝑑𝑟𝑑𝜃 Polar coordinates simplify the integration. =𝛿2∫2𝜋0[−16(4−𝑟2)3]𝑟=2𝑟=0𝑑𝜃=16𝛿3∫2𝜋0𝑑𝜃=32𝜋𝛿3.

A similar calculation gives the mass:

FIGURE 14.40 To find an integral for the amount of energy stored in a rotating shaft, we first imagine the shaft to be partitioned into small blocks. Each block has its own kinetic energy. We add the contributions of the individual blocks to find the kinetic energy of the shaft.

𝑀=∬𝑅∫4−𝑥2−𝑦20𝛿𝑑𝑧𝑑𝑦𝑑𝑥=8𝜋𝛿.

教材插图

Therefore, ――𝑧 =(𝑀𝑥𝑦/𝑀) =4/3 and the center of mass is (――𝑥,――𝑦,――𝑧) =(0,0,4/3) .

教材插图

FIGURE 14.39 The centroid of this region is found in Example 2.

When the density of a solid object or plate is constant (as in Example 1), the center of mass is called the centroid of the object. To find a centroid, we set 𝛿 equal to 1 and proceed to find ――𝑥 , ――𝑦 , and ――𝑧 as before, by dividing first moments by masses. These calculations are also valid for two-dimensional objects.

EXAMPLE 2 Find the centroid of the region in the first quadrant that is bounded above by the line y = x and below by the parabola 𝑦 =𝑥2 .

Solution We sketch the region and include enough detail to determine the limits of integration (Figure 14.39). We then set 𝛿 equal to 1 and evaluate the appropriate formulas from Table 14.1:

𝑀=∫10∫𝑥𝑥21𝑑𝑦𝑑𝑥=∫10[𝑦]𝑦=𝑥𝑦=𝑥2𝑑𝑥=∫10(𝑥−𝑥2)𝑑𝑥=[𝑥22−𝑥33]10=16𝑀𝑥=∫10∫𝑥𝑥2𝑦𝑑𝑦𝑑𝑥=∫10[𝑦22]𝑦=𝑥𝑦=𝑥2𝑑𝑥=∫10(𝑥22−𝑥42)𝑑𝑥=[𝑥36−𝑥510]10=115𝑀𝑦=∫10∫𝑥𝑥2𝑥𝑑𝑦𝑑𝑥=∫10[𝑥𝑦]𝑦=𝑥𝑦=𝑥2𝑑𝑥=∫10(𝑥2−𝑥3)𝑑𝑥=[𝑥33−𝑥44]10=112.

From these values of 𝑀,𝑀𝑥 , and 𝑀𝑦 , we find

――𝑥=𝑀𝑦𝑀=1/121/6=12 and ――𝑦=𝑀𝑥𝑀=1/151/6=25.

The centroid is the point (1/2,2/5) .

Note that each coordinate of the centroid of a region is equal to the average value of the corresponding variable over the region.

Moments of Inertia

An object’s first moments (Table 14.1) give us information related to balance and to the torque the object experiences about different axes in a gravitational field. If the object is a rotating shaft, we are interested in how much energy is stored in the shaft and how much energy is generated by a shaft rotating at a particular angular velocity. This is captured by the second moment or moment of inertia.

Think of partitioning the shaft into small blocks of mass Δ𝑚𝑘 and let 𝑟𝑘 denote the distance from the 𝑘 th block’s center of mass to the axis of rotation (Figure 14.40). If the shaft rotates at a constant angular velocity of 𝜔 =𝑑𝜃/𝑑𝑡 radians per second, the block’s center of mass will trace its orbit at a linear speed of

𝑣𝑘=𝑑𝑑𝑡(𝑟𝑘𝜃)=𝑟𝑘𝑑𝜃𝑑𝑡=𝑟𝑘𝜔.

The block’s kinetic energy will be approximately

12Δ𝑚𝑘𝑣2𝑘=12Δ𝑚𝑘(𝑟𝑘𝜔)2=12𝜔2𝑟2𝑘Δ𝑚𝑘.

The kinetic energy of the shaft will be approximately

∑12𝜔2𝑟2𝑘Δ𝑚𝑘.

The integral approached by these sums as the shaft is partitioned into smaller and smaller blocks gives the shaft’s kinetic energy:

KE shaft =∫12𝜔2𝑟2𝑑𝑚=12𝜔2∫𝑟2𝑑𝑚.(1)

The factor

𝐼=∫𝑟2𝑑𝑚

is the moment of inertia of the shaft about its axis of rotation, and we see from Equation (1) that the shaft’s kinetic energy is

KE shaft =12𝐼𝜔2.

The moment of inertia of a shaft resembles in some ways the inertial mass of a locomotive. To start a locomotive with mass m moving at a linear velocity v, we need to provide a kinetic energy of KE =(1/2)𝑚𝑣2 . To stop the locomotive we have to remove this amount of energy. To start a shaft with moment of inertia I rotating at an angular velocity 𝜔 , we need to provide a kinetic energy of KE =(1/2)𝐼𝜔2 . To stop the shaft we have to take this amount of energy back out. The shaft’s moment of inertia is analogous to the locomotive’s mass. What makes the locomotive hard to start or stop is its mass. What makes the shaft hard to start or stop is its moment of inertia. The moment of inertia depends not only on the mass of the shaft but also on its distribution. Mass that is farther away from the axis of rotation contributes more to the moment of inertia.

教材插图

FIGURE 14.41 Distances from dV to the axes.

We now derive a formula for the moment of inertia for a solid in space. If 𝑟(𝑥,𝑦,𝑧) is the distance from the point (𝑥,𝑦,𝑧) in D to a line L, then the moment of inertia of the mass Δ𝑚𝑘 =𝛿(𝑥𝑘,𝑦𝑘,𝑧𝑘)Δ𝑉𝑘 about the line L (as in Figure 14.40) is approximately Δ𝐼𝑘 =𝑟2(𝑥𝑘,𝑦𝑘,𝑧𝑘)Δ𝑚𝑘 . The moment of inertia about L of the entire object is

𝐼𝐿=lim𝑛→∞𝑛∑𝑘=1Δ𝐼𝑘=lim𝑛→∞𝑛∑𝑘=1𝑟2(𝑥𝑘,𝑦𝑘,𝑧𝑘)𝛿(𝑥𝑘,𝑦𝑘,𝑧𝑘)Δ𝑉𝑘=∭𝐷𝑟2𝛿𝑑𝑉.

If L is the x-axis, then 𝑟2 =𝑦2 +𝑧2 (Figure 14.41) and

𝐼𝑥=∭𝐷(𝑦2+𝑧2)𝛿(𝑥,𝑦,𝑧)𝑑𝑉.

Similarly, if L is the y-axis or the z-axis, we have

𝐼𝑦=∭𝐷(𝑥2+𝑧2)𝛿(𝑥,𝑦,𝑧)𝑑𝑉 and 𝐼𝑧=∭𝐷(𝑥2+𝑦2)𝛿(𝑥,𝑦,𝑧)𝑑𝑉.

Table 14.2 summarizes the formulas for these moments of inertia (second moments because they invoke the squares of the distances). It shows the definition of the polar moment about the origin as well.

教材插图

FIGURE 14.42 Finding 𝐼𝑥 , 𝐼𝑦 , and 𝐼𝑧 for the block shown here. The origin lies at the center of the block (Example 3).

TABLE 14.2 Moments of inertia (second moments) formulas

THREE-DIMENSIONAL SOLID About the x-axis: 𝐼𝑥 =∭𝐷(𝑦2 +𝑧2)𝛿𝑑𝑉 𝛿 =𝛿(𝑥,𝑦,𝑧) About the y-axis: 𝐼𝑦 =∭𝐷(𝑥2 +𝑧2)𝛿𝑑𝑉 About the z-axis: 𝐼𝑧 =∭𝐷(𝑥2 +𝑦2)𝛿𝑑𝑉 About a line L: 𝐼𝐿 =∭𝐷𝑟2(𝑥,𝑦,𝑧)𝛿𝑑𝑉 𝑟(𝑥,𝑦,𝑧) =distance from the point(𝑥,𝑦,𝑧) to line 𝐿 TWO-DIMENSIONAL PLATE About the x-axis: 𝐼𝑥 =∬𝑅𝑦2𝛿𝑑𝐴 𝛿 =𝛿(𝑥,𝑦) About the y-axis: 𝐼𝑦 =∬𝑅𝑥2𝛿𝑑𝐴 About a line L: 𝐼𝐿 =∬𝑅𝑟2(𝑥,𝑦)𝛿𝑑𝐴 𝑟(𝑥,𝑦) =distance from(𝑥,𝑦) to 𝐿 About the origin (polar moment): 𝐼0 =∬𝑅(𝑥2 +𝑦2)𝛿𝑑𝐴 =𝐼𝑥 +𝐼𝑦

EXAMPLE 3 Find 𝐼𝑥 , 𝐼𝑦 , 𝐼𝑧 for the rectangular solid of constant density 𝛿 shown in Figure 14.42.

Solution The formula for 𝐼𝑥 gives

𝐼𝑥=∫𝑐/2−𝑐/2∫𝑏/2−𝑏/2∫𝑎/2−𝑎/2(𝑦2+𝑧2)𝛿𝑑𝑥𝑑𝑦𝑑𝑧.

We can avoid some of the work of integration by observing that (𝑦2 +𝑧2)𝛿 is an even function of x, y, and z since 𝛿 is constant. The rectangular solid consists of eight symmetric pieces, one in each octant. We can evaluate the integral on one of these pieces and then multiply by 8 to get the total value:

𝐼𝑥=8∫𝑐/20∫𝑏/20∫𝑎/20(𝑦2+𝑧2)𝛿𝑑𝑥𝑑𝑦𝑑𝑧=4𝑎𝛿∫𝑐/20∫𝑏/20(𝑦2+𝑧2)𝑑𝑦𝑑𝑧=4𝑎𝛿∫𝑐/20[𝑦33+𝑧2𝑦]𝑦=𝑏/2𝑦=0𝑑𝑧=4𝑎𝛿∫𝑐/20(𝑏324+𝑧2𝑏2)𝑑𝑧

教材插图

FIGURE 14.43 The triangular region covered by the plate in Example 4.

教材插图

教材插图

FIGURE 14.44 The greater the polar moment of inertia of the cross-section of a beam about the beam’s longitudinal axis, the stiffer the beam. Beams A and B have the same cross-sectional area, but A is stiffer.

=4𝑎𝛿(𝑏3𝑐48+𝑐3𝑏48)=𝑎𝑏𝑐𝛿12(𝑏2+𝑐2)=𝑀12(𝑏2+𝑐2).𝑀=𝑎𝑏𝑐𝛿

Similarly,

𝐼𝑦=𝑀12(𝑎2+𝑐2) and 𝐼𝑧=𝑀12(𝑎2+𝑏2).

EXAMPLE 4 A thin plate covers the triangular region bounded by the x-axis and the lines x = 1 and y = 2x in the first quadrant. The plate’s density at the point (𝑥,𝑦) is 𝛿(𝑥,𝑦) =6𝑥 +6𝑦 +6 . Find the plate’s moments of inertia about the coordinate axes and the origin.

Solution We sketch the plate and put in enough detail to determine the limits of integration for the integrals we have to evaluate (Figure 14.43). The moment of inertia about the x-axis is

𝐼𝑥=∫10∫2𝑥0𝑦2𝛿(𝑥,𝑦)𝑑𝑦𝑑𝑥=∫10∫2𝑥0(6𝑥𝑦2+6𝑦3+6𝑦2)𝑑𝑦𝑑𝑥=∫10[2𝑥𝑦3+32𝑦4+2𝑦3]𝑦=2𝑥𝑦=0𝑑𝑥=∫10(40𝑥4+16𝑥3)𝑑𝑥=[8𝑥5+4𝑥4]10=12.

Similarly, the moment of inertia about the y-axis is

𝐼𝑦=∫10∫2𝑥0𝑥2𝛿(𝑥,𝑦)𝑑𝑦𝑑𝑥=395.

Notice that we integrate 𝑦2 times density in calculating 𝐼𝑥 , and 𝑥2 times density to find 𝐼𝑦 . Since we know 𝐼𝑥 and 𝐼𝑦 , we do not need to evaluate an integral to find 𝐼0 ; we can use the equation 𝐼0 =𝐼𝑥 +𝐼𝑦 from Table 14.2 instead:

𝐼0=12+395=60+395=995.

The moment of inertia also plays a role in determining how much a horizontal metal beam will bend under a load. The stiffness of the beam is a constant times 𝐼 , the moment of inertia of a typical cross-section of the beam about the beam’s longitudinal axis. The greater the value of 𝐼 , the stiffer the beam and the less it will bend under a given load. That is why we use I-beams instead of beams whose cross-sections are square. The flanges at the top and bottom of the beam hold most of the beam’s mass away from the longitudinal axis to increase the value of 𝐼 (Figure 14.44).

Probability

The probability that a continuous random variable X takes values between a and b is found by integrating a probability density function f (Appendix A.8),

𝑃(𝑎≤𝑋≤𝑏)=∫𝑏𝑎𝑓(𝑥)𝑑𝑥.

A similar process applies to probabilities involving two continuous random variables. The probability that a pair of random variables (𝑋,𝑌) takes values lying within a particular region is determined by a joint probability density function f. Integrating the joint probability density function over a region R in the plane gives the probability that the pair of random variables take values in that region:

𝑃((𝑋,𝑌)∈𝑅)=∬𝑅𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.

If the region is a rectangle, then this expression has the simple form

𝑃(𝑎≤𝑋≤𝑏 and 𝑐≤𝑌≤𝑑)=∫𝑑𝑐∫𝑏𝑎𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.

A joint probability density function f is defined by three basic properties. The first property ensures that there are no negative probabilities, and the second implies that the total probability of all possible outcomes is one. The final property describes the connection of f to probabilities.

DEFINITION A joint probability density function 𝑓 is a function that satisfies three conditions:

  1. 𝑓(𝑥,𝑦) ≥0
∫∞−∞∫∞−∞𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦=1 3.𝑃((𝑋,𝑌)∈𝑅)=∬𝑅𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.

A pair of random variables has a uniform distribution on a region 𝑅 with finite area 𝐴 if 𝑓(𝑥,𝑦) =1/𝐴 for any (𝑥,𝑦) ∈𝑅 , and 𝑓(𝑥,𝑦) =0 otherwise.

EXAMPLE 5 A random number generator is used to generate two random real numbers X and Y in succession. The first number X is chosen between 0 and 10, and the second number Y is chosen between 0 and 5. The random number generation is done by a process that gives a uniform distribution. Find the joint probability density function f for the pair of numbers (𝑋,𝑌) and use it to compute the probability that X is larger than Y.

Solution The joint probability density function f is constant on the rectangle 0 ≤𝑥 ≤10 , 0 ≤𝑦 ≤5 , because (𝑋,𝑌) is uniformly distributed. The area of the rectangle is 50, so f takes the value 1/50 inside this rectangle:

教材插图

FIGURE 14.45 The pair of random variables X and Y take values anywhere in this rectangle with equal probability. In the shaded region we have X > Y.

𝑓(𝑥,𝑦)={1/50, if 0≤𝑥≤10 and 0≤𝑦≤5,0, otherwise. 

To compute the probability that X > Y, we integrate the joint probability density function f over the region in the rectangle where X > Y. This region is bounded on the left by the line x = y and on the right by the line x = 10. An integral over this region has limits of integration given by 𝑦 ≤𝑥 ≤10 , 0 ≤𝑦 ≤5 (see Figure 14.45). The probability is given by

𝑃(𝑋>𝑌)=∫50∫10𝑦150𝑑𝑥𝑑𝑦=34.

There is a 75% probability that the first number is larger than the second.

EXAMPLE 6 Using the joint probability density function

𝑓(𝑥,𝑦)={𝑒−(𝑥+𝑦), if 0<𝑥 and 0<𝑦0, otherwise 

find the probability that 1 < X < 2 and 2 < Y < 3.

Solution

𝑃(1<𝑋<2,2<𝑌<3)=∫32∫21𝑒−(𝑥+𝑦)𝑑𝑥𝑑𝑦=𝑒−5+𝑒−3−2𝑒−4≈0.01989.

There is slightly less than a 2% probability that X and Y fall within these bounds.

Means and Expected Values

The mean, or expected value, of a random variable is (Appendix A.8)

𝜇=∫∞−∞𝑥𝑓(𝑥)𝑑𝑥.

When X and Y have joint probability density function f, the expected value of X and the expected value of Y are

𝜇𝑋=∫∞−∞∫∞−∞𝑥𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦 and 𝜇𝑌=∫∞−∞∫∞−∞𝑦𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦.

These indicate the average value expected for each of X and Y. The expected values 𝜇𝑋 and 𝜇𝑌 are sometimes called the first moments of the distribution, because their defining formulas have the same form as those seen in Table 14.1 for the moments of a two-dimensional plate. The joint probability density function plays the role in computing 𝜇𝑋 that the mass density function plays in computing the x-coordinate ――𝑥 of the center of mass, and the same applies to 𝜇𝑌 and ――𝑦 . One can roughly think of the joint probability density function as measuring the probability concentration per unit area on the plane, just as density measures the mass per unit area for a plate.

EXAMPLE 7 Find the expected values 𝜇𝑋 and 𝜇𝑌 for the joint probability density function in Example 5.

Solution For the joint probability density function in Example 5, we compute

𝜇𝑋=∫50∫100𝑥(1/50)𝑑𝑥𝑑𝑦=5

and

𝜇𝑌=∫50∫100𝑦(1/50)𝑑𝑥𝑑𝑦=2.5.

The expected value of X is 5 and that of Y is 2.5.

EXERCISES 14.6

Plates of Constant Density

  1. Finding a center of mass Find the center of mass of a thin plate of density 𝛿 =3 bounded by the lines x = 0, y = x, and the parabola 𝑦 =2 −𝑥2 in the first quadrant.

  2. Finding moments of inertia Find the moments of inertia about the coordinate axes of a thin rectangular plate of constant density 𝛿gm/cm2 bounded by the lines 𝑥 =3 and 𝑦 =3 in the first quadrant.

  3. Finding a centroid Find the centroid of the region in the first quadrant bounded by the x-axis, the parabola 𝑦2 =2𝑥 , and the line 𝑥 +𝑦 =4 .

  4. Finding a centroid Find the centroid of the triangular region cut from the first quadrant by the line 𝑥 +𝑦 =3 .

  5. Finding a centroid Find the centroid of the region cut from the first quadrant by the circle 𝑥2 +𝑦2 =𝑎2 .

  6. Finding a centroid Find the centroid of the region between the x-axis and the arch 𝑦 =sin⁡𝑥,0 ≤𝑥 ≤𝜋 .

  7. Finding moments of inertia Find the moment of inertia about the 𝑥 -axis of a thin plate of density 𝛿 =1gm/cm2 bounded by the circle 𝑥2 +𝑦2 =4 . Then use your result to find 𝐼𝑦 and 𝐼0 for the plate.

  8. Finding a moment of inertia Find the moment of inertia with respect to the y-axis of a thin sheet of constant density 𝛿 =1 𝑔𝑚/𝑐𝑚2 bounded by the curve 𝑦 =(sin2⁡𝑥)/𝑥2 and the interval 𝜋 ≤𝑥 ≤2𝜋 of the x-axis.

  9. The centroid of an infinite region Find the centroid of the infinite region in the second quadrant enclosed by the coordinate axes and the curve 𝑦 =𝑒𝑥 . (Use improper integrals in the mass-moment formulas.)

  10. The first moment of an infinite plate Find the first moment about the y-axis of a thin plate of density 𝛿(𝑥,𝑦) =1 covering the infinite region under the curve 𝑦 =𝑒−𝑥2/2 in the first quadrant.

Plates with Varying Density

  1. Finding a moment of inertia Find the moment of inertia about the x-axis of a thin plate bounded by the parabola 𝑥 =𝑦 −𝑦2 and the line 𝑥 +𝑦 =0 if 𝛿(𝑥,𝑦) =𝑥 +𝑦 .

  2. Finding mass Find the mass of a thin plate occupying the smaller region cut from the ellipse 𝑥2 +4𝑦2 =12 by the parabola 𝑥 =4𝑦2 if 𝛿(𝑥,𝑦) =5𝑥kg/m2 .

  3. Finding a center of mass Find the center of mass of a thin triangular plate bounded by the 𝑦 -axis and the lines 𝑦 =𝑥 and 𝑦 =2 −𝑥 if 𝛿(𝑥,𝑦) =6𝑥 +3𝑦 +3 .

  4. Finding a center of mass and moment of inertia Find the center of mass and moment of inertia about the x-axis of a thin plate bounded by the curves 𝑥 =𝑦2 and 𝑥 =2𝑦 −𝑦2 if the density at the point (𝑥,𝑦) is 𝛿(𝑥,𝑦) =𝑦 +1 .

  5. Center of mass, moment of inertia Find the center of mass and the moment of inertia about the y-axis of a thin rectangular plate cut from the first quadrant by the lines x = 6 and y = 1 if 𝛿(𝑥,𝑦) =𝑥 +𝑦 +1 .

  6. Center of mass, moment of inertia Find the center of mass and the moment of inertia about the 𝑦 -axis of a thin plate bounded by the line 𝑦 =1 and the parabola 𝑦 =𝑥2 if the density is 𝛿(𝑥,𝑦) =𝑦 +1 .

  7. Center of mass, moment of inertia Find the center of mass and the moment of inertia about the 𝑦 -axis of a thin plate bounded by the 𝑥 -axis, the lines 𝑥 = ±1 , and the parabola 𝑦 =𝑥2 if 𝛿(𝑥,𝑦) =7𝑦 +1 .

  8. Center of mass, moments of inertia Find the center of mass and the moments of inertia about the 𝑥 -axis of a thin rectangular plate bounded by the lines 𝑥 =0 , 𝑥 =20 , 𝑦 = −1 , and 𝑦 =1 if 𝛿(𝑥,𝑦) =1 +(𝑥/20) .

  9. Center of mass, moments of inertia Find the center of mass, the moment of inertia about the coordinate axes, and the polar moment of inertia of a thin triangular plate bounded by the lines y = x, y = -x, and y = 1 if 𝛿(𝑥,𝑦) =𝑦 +1 kg/m2 .

  10. Center of mass, moments of inertia Repeat Exercise 19 for 𝛿(𝑥,𝑦) =3𝑥2 +1 𝑘𝑔/𝑚2 .

Solids with Constant Density

  1. Moments of inertia Find the moments of inertia of the rectangular box of constant density 𝛿(𝑥,𝑦,𝑧) =1 shown here with respect to its edges by calculating 𝐼𝑥 , 𝐼𝑦 , and 𝐼𝑧 .

教材插图

  1. Moments of inertia The coordinate axes in the figure run through the centroid of a solid wedge parallel to the labeled edges. Find 𝐼𝑥 , 𝐼𝑦 , and 𝐼𝑧 if a = b = 6, c = 4, and the density is 𝛿(𝑥,𝑦,𝑧) =1 .

教材插图

  1. Center of mass and moments of inertia A solid “trough” of constant density 𝛿(𝑥,𝑦,𝑧) =1 is bounded below by the surface 𝑧 =4𝑦2 , above by the plane z=4, and on the ends by the planes x=1 and x=-1. Find the center of mass and the moments of inertia with respect to the three axes.

  2. Center of mass A solid of constant density is bounded below by the plane z = 0, on the sides by the elliptical cylinder 𝑥2 +4𝑦2 =4 , and above by the plane z = 2 - x (see the accompanying figure).

a. Find ――𝑥 and ――𝑦 .

b. Evaluate the integral

𝑀𝑥𝑦=∫2−2∫(1/2)√4−𝑥2−(1/2)√4−𝑥2∫2−𝑥0𝑧𝑑𝑧𝑑𝑦𝑑𝑥,

using integral tables to carry out the final integration with respect to 𝑥 . Then divide 𝑀𝑥𝑦 by 𝑀 to verify that ――𝑧 =5/4 .

教材插图

  1. a. Center of mass Find the center of mass of a solid of constant density bounded below by the paraboloid 𝑧 =𝑥2 +𝑦2 and above by the plane z = 4.

b. Find the plane z = c that divides the solid into two parts of equal volume. This plane does not pass through the center of mass.

  1. Moments A solid cube of constant density 𝛿(𝑥,𝑦,𝑧) =1 , 2 units on a side, is bounded by the planes 𝑥 = ±1 , 𝑧 = ±1 , y=3, and y=5. Find the center of mass and the moments of inertia about the coordinate axes.

  2. Moment of inertia about a line A wedge like the one in Exercise 22 has 𝑎 =4 , 𝑏 =6 , 𝑐 =3 , and a constant density 𝛿(𝑥,𝑦,𝑧) =1 . Make a quick sketch to check for yourself that the square of the distance from a typical point (𝑥,𝑦,𝑧) of the wedge to the line 𝐿 : 𝑧 =0 , 𝑦 =6 is 𝑟2 =(𝑦 −6)2 +𝑧2 . Then calculate the moment of inertia of the wedge about 𝐿 .

  3. Moment of inertia about a line A wedge like the one in Exercise 22 has 𝑎 =4 , 𝑏 =6 , 𝑐 =3 , and a constant density 𝛿(𝑥,𝑦,𝑧) =1 . Make a quick sketch to check for yourself that the square of the distance from a typical point (𝑥,𝑦,𝑧) of the wedge to the line 𝐿 : 𝑥 =4 , 𝑦 =0 is 𝑟2 =(𝑥 −4)2 +𝑦2 . Then calculate the moment of inertia of the wedge about 𝐿 .

Solids with Varying Density

In Exercises 29 and 30, find

a. the mass of the solid. b. the center of mass.

  1. A solid region in the first octant is bounded by the coordinate planes and the plane 𝑥 +𝑦 +𝑧 =2 . The density of the solid is 𝛿(𝑥,𝑦,𝑧) =2𝑥 gm/cm3 .

  2. A solid in the first octant is bounded by the planes y = 0 and z = 0 and by the surfaces 𝑧 =4 −𝑥2 and 𝑥 =𝑦2 (see the accompanying figure). Its density function is 𝛿(𝑥,𝑦,𝑧) =𝑘𝑥𝑦 , k a constant.

教材插图

In Exercises 31 and 32, find

a. the mass of the solid.

b. the center of mass.

c. the moments of inertia about the coordinate axes.

  1. A solid cube in the first octant is bounded by the coordinate planes and by the planes x = 1, y = 1, and z = 1. The density of the cube is 𝛿(𝑥,𝑦,𝑧) =𝑥 +𝑦 +𝑧 +1 .

  2. A wedge like the one in Exercise 22 has dimensions 𝑎 =2 , 𝑏 =6 , and 𝑐 =3 . The density is 𝛿(𝑥,𝑦,𝑧) =𝑥 +1 . Notice that if the density is constant, the center of mass will be (0,0,0) .

  3. Mass Find the mass of the solid bounded by the planes 𝑥 +𝑧 =1 , x - z = -1, y = 0, and the surface 𝑦 =√𝑧 . The density of the solid is 𝛿(𝑥,𝑦,𝑧) =2𝑦 +5kg/m3 .

  4. Mass Find the mass of the solid region bounded by the parabolic surfaces 𝑧 =16 −2𝑥2 −2𝑦2 and 𝑧 =2𝑥2 +2𝑦2 if the density of the solid is 𝛿(𝑥,𝑦,𝑧) =√𝑥2+𝑦2 .

Theory and Examples

The Parallel Axis Theorem Let 𝐿𝑐.𝑚. be a line through the center of mass of a body of mass m and let L be a parallel line h units away from 𝐿𝑐.𝑚. . The Parallel Axis Theorem says that the moments of inertia 𝐼𝑐.𝑚. and 𝐼𝐿 of the body about 𝐿𝑐.𝑚. and L satisfy the equation

𝐼𝐿=𝐼c.m.+𝑚ℎ2.(2)

As in the two-dimensional case, the theorem gives a quick way to calculate one moment when the other moment and the mass are known.

  1. Proof of the Parallel Axis Theorem

a. Show that the first moment of a body in space about any plane through the body’s center of mass is zero. (Hint: Place the body’s center of mass at the origin and let the plane be the yz-plane. What does the formula ――𝑥 =𝑀𝑦𝑧/𝑀 then tell you?)

教材插图

b. To prove the Parallel Axis Theorem, place the body with its center of mass at the origin, with the line 𝐿𝑐.𝑚. along the z-axis and the line L perpendicular to the xy-plane at the point (ℎ,0,0) . Let D be the region of space occupied by the body. Then, in the notation of the figure,

𝐼𝐿=∭𝐷|𝐯−ℎ𝐢|2𝑑𝑚.

Expand the integrand in this integral and complete the proof.

  1. The moment of inertia about a diameter of a solid sphere of constant density and radius 𝑎 is (2/5)𝑚𝑎2 , where 𝑚 is the mass of the sphere. Find the moment of inertia about a line tangent to the sphere.

  2. The moment of inertia of the solid in Exercise 21 about the 𝑧 -axis is 𝐼𝑧 =𝑎𝑏𝑐(𝑎2 +𝑏2)/3 .

a. Use Equation (2) to find the moment of inertia of the solid about the line parallel to the 𝑧 -axis through the solid’s center of mass.

b. Use Equation (2) and the result in part (a) to find the moment of inertia of the solid about the line x = 0, y = 2b.

  1. If 𝑎 =𝑏 =6 and 𝑐 =4 , the moment of inertia of the solid wedge in Exercise 22 about the 𝑥 -axis is 𝐼𝑥 =208 . Find the moment of inertia of the wedge about the line 𝑦 =4 , 𝑧 = −4/3 (the edge of the wedge’s narrow end).

Joint Probability Density Functions

For Exercises 39–42, verify that f gives a joint probability density function. Then find the expected values 𝜇𝑋 and 𝜇𝑌 .

  1. 𝑓(𝑥,𝑦) ={𝑥+𝑦, if 0≤𝑥≤1 and 0≤𝑦≤1,0, otherwise .

  2. 𝑓(𝑥,𝑦) ={4𝑥𝑦, if 0≤𝑥≤1 and 0≤𝑦≤1,0, otherwise. 

  3. 𝑓(𝑥,𝑦) ={6𝑥2𝑦, if 0≤𝑥≤1 and 0≤𝑦≤1,0, otherwise. 

  4. 𝑓(𝑥,𝑦) ={32(𝑥2+𝑦2), if 0≤𝑥≤1 and 0≤𝑦≤1,0, otherwise. 

  5. Suppose that 𝑓 is a uniform joint probability density function on 0 ≤𝑥 <2 , 0 ≤𝑦 <3 . What is the formula for 𝑓 ? What is the probability that 𝑋 <𝑌 ?

  6. The following formula defines a joint probability density function. What is the value of C? What are the expected values 𝜇𝑋 and 𝜇𝑌 ?

𝑓(𝑥,𝑦)={𝐶𝑥𝑦, if 0≤𝑥≤2 and 0≤𝑦≤3,0, otherwise. 

14.7 Triple Integrals in Cylindrical and Spherical Coordinates

教材插图

FIGURE 14.46 The cylindrical coordinates of a point in space are r, 𝜃 , and z.

When a calculation in physics, engineering, or geometry involves a cylinder, cone, or sphere, we can often simplify our work by using cylindrical or spherical coordinates, which are introduced in this section. The procedure for transforming to these coordinates and evaluating the resulting triple integrals is similar to the transformation to polar coordinates in the plane discussed in Section 14.4.

教材插图

FIGURE 14.47 Constant-coordinate equations in cylindrical coordinates yield cylinders and planes.

Integration in Cylindrical Coordinates

We obtain cylindrical coordinates for space by combining polar coordinates in the xy-plane with the usual z-axis. This assigns to every point in space coordinate triples of the form (𝑟,𝜃,𝑧) , as shown in Figure 14.46. Here we require 𝑟 ≥0 .

DEFINITION Cylindrical coordinates represent a point P in space by ordered triples (𝑟,𝜃,𝑧) in which

  1. r and 𝜃 are polar coordinates for the vertical projection of P on the xy-plane, with 𝑟 ≥0 , and

  2. z is the rectangular vertical coordinate.

The values of 𝑥,𝑦,𝑟 , and 𝜃 in rectangular and cylindrical coordinates are related by the usual equations.

Equations Relating Rectangular (𝑥,𝑦,𝑧) and Cylindrical (𝑟,𝜃,𝑧) Coordinates

𝑥=𝑟cos⁡𝜃,𝑦=𝑟sin⁡𝜃,𝑧=𝑧,𝑟2=𝑥2+𝑦2,tan⁡𝜃=𝑦/𝑥

In cylindrical coordinates, the equation r = a describes not just a circle in the xy-plane but an entire cylinder about the z-axis (Figure 14.47). The z-axis is given by r = 0. The equation 𝜃 =𝜃0 describes the half-plane that contains the z-axis and makes an angle 𝜃0 with the positive x-axis. And, just as in rectangular coordinates, the equation 𝑧 =𝑧0 describes a plane perpendicular to the z-axis.

教材插图

FIGURE 14.48 In cylindrical coordinates the volume of the wedge is approximated by the product Δ𝑉 =𝑟Δ𝑧Δ𝑟Δ𝜃 .

Volume Differential in Cylindrical Coordinates

𝑑𝑉=𝑟𝑑𝑧𝑑𝑟𝑑𝜃

教材插图

FIGURE 14.49 Finding the limits of integration for evaluating an integral in cylindrical coordinates (Example 1).

Cylindrical coordinates are good for describing cylinders whose axes run along the z-axis and planes that either contain the z-axis or lie perpendicular to the z-axis. Surfaces like these have equations of constant coordinate value:

𝑟=4 Cylinder, radius 4, axis the z -axis𝜃=𝜋3 Half - plane containing the z -axis 𝑧=2. Plane perpendicular to the z -axis 

When computing triple integrals over a solid region D in cylindrical coordinates, we partition the region into n small cylindrical wedges, rather than into rectangular boxes. In the kth cylindrical wedge, 𝑟,𝜃 , and z change by Δ𝑟𝑘,Δ𝜃𝑘 , and Δ𝑧𝑘 , and the largest of these numbers among all the cylindrical wedges is called the norm of the partition. We express the triple integral as a limit of Riemann sums using these wedges. The volume of such a cylindrical wedge Δ𝑉𝑘 is obtained by taking the area Δ𝐴𝑘 of its base in the 𝑟𝜃 -plane and multiplying by the height Δ𝑧𝑘 (Figure 14.48).

For a point (𝑟𝑘,𝜃𝑘,𝑧𝑘) in the center of the kth wedge, we calculated in polar coordinates that Δ𝐴𝑘 =𝑟𝑘Δ𝑟𝑘Δ𝜃𝑘 . So Δ𝑉𝑘 =Δ𝑧𝑘𝑟𝑘Δ𝑟𝑘Δ𝜃𝑘 =𝑟𝑘Δ𝑧𝑘Δ𝑟𝑘Δ𝜃𝑘 , and a Riemann sum for f over D has the form

𝑆𝑛=𝑛∑𝑘=1𝑓(𝑟𝑘,𝜃𝑘,𝑧𝑘)𝑟𝑘Δ𝑧𝑘Δ𝑟𝑘Δ𝜃𝑘.

The triple integral of a function 𝑓 over 𝐷 is obtained by taking a limit of such Riemann sums with partitions whose norms approach zero:

lim𝑛→∞𝑆𝑛=∭𝐷𝑓𝑑𝑉=∭𝐷𝑓𝑟𝑑𝑧𝑑𝑟𝑑𝜃.

Triple integrals in cylindrical coordinates are then evaluated as iterated integrals, as in the following example. Although the definition of cylindrical coordinates makes sense without any restrictions on 𝜃 , in most situations when integrating, we will need to restrict 𝜃 to an interval of length 2𝜋 . So we impose the requirement that 𝛼 ≤𝜃 ≤𝛽 , where 0 ≤𝛽 −𝛼 ≤2𝜋 .

EXAMPLE 1 Find the limits of integration in cylindrical coordinates for integrating a function 𝑓(𝑟,𝜃,𝑧) over the solid region D bounded below by the plane z = 0, laterally by the circular cylinder 𝑥2 +(𝑦 −1)2 =1 , and above by the paraboloid 𝑧 =𝑥2 +𝑦2 .

Solution The base of 𝐷 is also the region’s projection 𝑅 on the 𝑥𝑦 -plane. The boundary of 𝑅 is the circle 𝑥2 +(𝑦 −1)2 =1 . Its polar coordinate equation is

𝑥2+(𝑦−1)2=1𝑥2+𝑦2−2𝑦+1=1𝑟2−2𝑟sin⁡𝜃=0𝑟=2sin⁡𝜃.

The region is sketched in Figure 14.49.

We find the limits of integration, starting with the 𝑧 -limits. A line 𝑀 through a typical point (𝑟,𝜃) in 𝑅 parallel to the 𝑧 -axis enters 𝐷 at 𝑧 =0 and leaves at 𝑧 =𝑥2 +𝑦2 =𝑟2 .

Next we find the r-limits of integration. A ray L through (𝑟,𝜃) from the origin enters R at r = 0 and leaves at 𝑟 =2sin⁡𝜃 .

Finally, we find the 𝜃 -limits of integration. As L sweeps across R, the angle 𝜃 it makes with the positive x-axis runs from 𝜃 =0 to 𝜃 =𝜋 . The integral is

∭𝐷𝑓(𝑟,𝜃,𝑧)𝑑𝑉=∫𝜋0∫2sin⁡𝜃0∫𝑟20𝑓(𝑟,𝜃,𝑧)𝑟𝑑𝑧𝑑𝑟𝑑𝜃.

Example 1 illustrates a good procedure for finding limits of integration in cylindrical coordinates. The procedure is summarized as follows.

How to Integrate in Cylindrical Coordinates

To evaluate

∭𝐷𝑓(𝑟,𝜃,𝑧)𝑑𝑉

over a solid region D in space in cylindrical coordinates, integrating first with respect to z, then with respect to r, and finally with respect to 𝜃 , take the following steps.

  1. Sketch. Sketch the solid region D along with its projection R on the xy-plane. Label the surfaces and curves that bound D and R.

教材插图

  1. Find the 𝑧 -limits of integration. Draw a line 𝑀 through a typical point (𝑟,𝜃) of 𝑅 parallel to the 𝑧 -axis. As 𝑧 increases, 𝑀 enters 𝐷 at 𝑧 =𝑔1(𝑟,𝜃) and leaves at 𝑧 =𝑔2(𝑟,𝜃) . These are the 𝑧 -limits of integration.

教材插图

教材插图

FIGURE 14.50 Example 2 shows how to find the centroid of this solid.

  1. Find the r-limits of integration. Draw a ray L through (𝑟,𝜃) from the origin. The ray enters R at 𝑟 =ℎ1(𝜃) and leaves at 𝑟 =ℎ2(𝜃) . These are the r-limits of integration.

教材插图

  1. Find the 𝜃 -limits of integration. As L sweeps across R, the angle 𝜃 it makes with the positive x-axis runs from 𝜃 =𝛼 to 𝜃 =𝛽 . These are the 𝜃 -limits of integration. The integral is
∭𝐷𝑓(𝑟,𝜃,𝑧)𝑑𝑉=∫𝜃=𝛽𝜃=𝛼∫𝑟=ℎ2(𝜃)𝑟=ℎ1(𝜃)∫𝑧=𝑔2(𝑟,𝜃)𝑧=𝑔1(𝑟,𝜃)𝑓(𝑟,𝜃,𝑧)𝑟𝑑𝑧𝑑𝑟𝑑𝜃.

EXAMPLE 2 Find the centroid ( 𝛿 =1 ) of the solid enclosed by the cylinder 𝑥2 +𝑦2 =4 , bounded above by the paraboloid 𝑧 =𝑥2 +𝑦2 , and bounded below by the xy-plane.

Solution We sketch the solid, bounded above by the paraboloid 𝑧 =𝑟2 and below by the plane z = 0 (Figure 14.50). Its base R is the disk 0 ≤𝑟 ≤2 in the xy-plane.

The solid’s centroid (――𝑥,――𝑦,――𝑧) lies on its axis of symmetry, here the 𝑧 -axis. This makes ――𝑥 =――𝑦 =0 . To find ――𝑧 , we divide the first moment 𝑀𝑥𝑦 by the mass 𝑀 .

To find the limits of integration for the mass and moment integrals, we continue with the four basic steps. We completed our initial sketch. The remaining steps give the limits of integration.

The z-limits. A line M through a typical point (𝑟,𝜃) in the base parallel to the z-axis enters the solid at z = 0 and leaves at 𝑧 =𝑟2 .

The 𝑟 -limits. A ray 𝐿 through (𝑟,𝜃) from the origin enters 𝑅 at 𝑟 =0 and leaves at 𝑟 =2 .

The 𝜃 -limits. As 𝐿 sweeps over the base like a clock hand, the angle 𝜃 it makes with the positive 𝑥 -axis runs from 𝜃 =0 to 𝜃 =2𝜋 . The value of 𝑀𝑥𝑦 is

𝑀𝑥𝑦=∫2𝜋0∫20∫𝑟20𝑧𝑟𝑑𝑧𝑑𝑟𝑑𝜃=∫2𝜋0∫20[𝑧22]𝑧=𝑟2𝑧=0𝑟𝑑𝑟𝑑𝜃=∫2𝜋0∫20𝑟52𝑑𝑟𝑑𝜃=∫2𝜋0[𝑟612]𝑟=2𝑟=0𝑑𝜃=∫2𝜋0163𝑑𝜃=32𝜋3.

The value of 𝑀 is

𝑀=∫2𝜋0∫20∫𝑟20𝑟𝑑𝑧𝑑𝑟𝑑𝜃=∫2𝜋0∫20[𝑧]𝑧=𝑟2𝑧=0𝑟𝑑𝑟𝑑𝜃=∫2𝜋0∫20𝑟3𝑑𝑟𝑑𝜃=∫2𝜋0[𝑟44]𝑟=2𝑟=0𝑑𝜃=∫2𝜋04𝑑𝜃=8𝜋.

𝜙 is the Greek letter phi, pronounced “fee.”

教材插图

FIGURE 14.51 The spherical coordinates 𝜌 , 𝜙 , and 𝜃 and their relation to x, y, z, and r.

教材插图

FIGURE 14.52 Constant-coordinate equations in spherical coordinates yield spheres, single cones, and half-planes.

Therefore,

――𝑧=𝑀𝑥𝑦𝑀=32𝜋318𝜋=43,

and the centroid is (0,0,4/3) . Notice that the centroid lies on the z-axis, outside the solid.

Spherical Coordinates and Integration

Spherical coordinates locate points in space with two angles and one distance, as shown in Figure 14.51. The first coordinate, 𝜌 =|⟶𝑂𝑃| , is the point’s distance from the origin and is never negative. The second coordinate, 𝜙 , is the angle ⟶𝑂𝑃 makes with the positive 𝑧 -axis. It is required to lie in the interval [0,𝜋] . The third coordinate is the angle 𝜃 as measured in cylindrical coordinates.

DEFINITION Spherical coordinates represent a point P in space by ordered triples (𝜌,𝜙,𝜃) in which

  1. 𝜌 is the distance from P to the origin ( 𝜌 ≥0 ).

  2. 𝜙 is the angle ⟶𝑂𝑃 makes with the positive z-axis ( 0 ≤𝜙 ≤𝜋 ).

  3. 𝜃 is the angle from cylindrical coordinates.

On maps of Earth, 𝜃 is related to the longitude of a point on the planet and 𝜙 to its latitude, while 𝜌 is related to elevation above Earth’s surface.

The equation 𝜌 =𝑎 describes the sphere of radius a centered at the origin (Figure 14.52). The equation 𝜙 =𝜙0 describes a single cone whose vertex lies at the origin and whose axis lies along the z-axis. (We broaden our interpretation to include the xy-plane as the cone 𝜙 =𝜋/2 .) If 𝜙0 is greater than 𝜋/2 , the cone 𝜙 =𝜙0 opens downward. The equation 𝜃 =𝜃0 describes the half-plane that contains the z-axis and makes an angle 𝜃0 with the positive x-axis.

Equations Relating Spherical Coordinates to Cartesian and Cylindrical Coordinates

𝑟=𝜌sin⁡𝜙,𝑥=𝑟cos⁡𝜃=𝜌sin⁡𝜙cos⁡𝜃,𝑧=𝜌cos⁡𝜙,𝑦=𝑟sin⁡𝜃=𝜌sin⁡𝜙sin⁡𝜃,𝜌=√𝑥2+𝑦2+𝑧2=√𝑟2+𝑧2.(1)

EXAMPLE 3 Find a spherical coordinate equation for the sphere 𝑥2 +𝑦2 +(𝑧 −1)2 =1 .

Solution We use Equations (1) to substitute for x, y, and z:

𝑥2+𝑦2+(𝑧−1)2=1𝜌2sin2⁡𝜙cos2⁡𝜃+𝜌2sin2⁡𝜙sin2⁡𝜃+(𝜌cos⁡𝜙−1)2=1𝜌2sin2⁡𝜙(cos2⁡𝜃+sin2⁡𝜃⏟___⏟___⏟1)+𝜌2cos2⁡𝜙−2𝜌cos⁡𝜙+1=1(Eqs. (1) 𝜌2(sin2⁡𝜙+cos2⁡𝜙)⏟___⏟___⏟1=2𝜌cos⁡𝜙 𝜌2=2𝜌cos⁡𝜙

教材插图

FIGURE 14.53 The sphere in Example 3.

教材插图

FIGURE 14.54 The cone in Example 4.

教材插图

FIGURE 14.55 In spherical coordinates we use the volume of a spherical wedge, which closely approximates that of a rectangular box.

Volume Differential in Spherical Coordinates

𝑑𝑉=𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃

The angle 𝜙 varies from 0 at the north pole of the sphere to 𝜋/2 at the south pole; the angle 𝜃 does not appear in the expression for 𝜌 , reflecting the symmetry about the z-axis (see Figure 14.53).

EXAMPLE 4 Find a spherical coordinate equation for the cone 𝑧 =√𝑥2+𝑦2 .

Solution 1 Use geometry. The cone is symmetric with respect to the z-axis and cuts the first quadrant of the yz-plane along the line z = y. The angle between the cone and the positive z-axis is therefore 𝜋/4 radians. The cone consists of the points whose spherical coordinates have 𝜙 equal to 𝜋/4 , so its equation is 𝜙 =𝜋/4 . (See Figure 14.54.)

Solution 2 Use algebra. If we use Equations (1) to substitute for x, y, and z, we obtain the same result:

𝑧=√𝑥2+𝑦2𝜌cos⁡𝜙=√𝜌2sin2⁡𝜙 Example 3 𝜌cos⁡𝜙=𝜌sin⁡𝜙𝜌≥0,sin⁡𝜙≥0cos⁡𝜙=sin⁡𝜙 Includes 𝜌=0( the origin )𝜙=𝜋4.0≤𝜙≤𝜋

Spherical coordinates are useful for describing spheres centered at the origin, half-planes hinged along the z-axis, and cones whose vertices lie at the origin and whose axes lie along the z-axis. Surfaces like these have equations of constant coordinate value:

𝜌=4 Sphere, radius 4, center at origin 𝜙=2𝜋3 Cone opening down from the origin, making an angle of 2𝜋/3 radians with the positive 𝑧-axis 𝜃=𝜋3. Half -plane, hinged along the 𝑧-axis, making an angle of 𝜋/3 radians with the positive 𝑥-axis 

When computing triple integrals over a solid region D in spherical coordinates, we partition the region into n spherical wedges. The size of the kth spherical wedge, which contains a point (𝜌𝑘,𝜙𝑘,𝜃𝑘) , is given by the changes Δ𝜌𝑘 , Δ𝜙𝑘 , and Δ𝜃𝑘 in 𝜌 , 𝜙 , and 𝜃 . Such a spherical wedge has one edge a circular arc of length 𝜌𝑘Δ𝜙𝑘 , another edge a circular arc of length 𝜌𝑘sin⁡𝜙𝑘Δ𝜃𝑘 , and thickness Δ𝜌𝑘 . The spherical wedge closely approximates a rectangular box of these dimensions when Δ𝜌𝑘 , Δ𝜙𝑘 , and Δ𝜃𝑘 are all small (Figure 14.55). It can be shown that the volume of this spherical wedge Δ𝑉𝑘 is Δ𝑉𝑘 =𝜌2𝑘sin⁡𝜙𝑘Δ𝜌𝑘Δ𝜙𝑘Δ𝜃𝑘 for (𝜌𝑘,𝜙𝑘,𝜃𝑘) , a point chosen inside the wedge.

The corresponding Riemann sum for a function 𝑓(𝜌,𝜙,𝜃) is

𝑆𝑛=𝑛∑𝑘=1𝑓(𝜌𝑘,𝜙𝑘,𝜃𝑘)𝜌2𝑘sin⁡𝜙𝑘Δ𝜌𝑘Δ𝜙𝑘Δ𝜃𝑘.

As the norm of a partition approaches zero, and the spherical wedges get smaller, the limit of the Riemann sums is the triple integral:

lim𝑛→∞𝑆𝑛=∭𝐷𝑓(𝜌,𝜙,𝜃)𝑑𝑉=∭𝐷𝑓(𝜌,𝜙,𝜃)𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃.

教材插图

FIGURE 14.56 The ice cream cone in Example 5.

To evaluate integrals in spherical coordinates, we usually integrate first with respect to 𝜌 . The procedure for finding the limits of integration is as follows. As with cylindrical coordinates, we restrict 𝜃 in the form 𝛼 ≤𝜃 ≤𝛽 and 0 ≤𝛽 −𝛼 ≤2𝜋 .

How to Integrate in Spherical Coordinates

To evaluate

∭𝐷𝑓(𝜌,𝜙,𝜃)𝑑𝑉

over a solid region D in space in spherical coordinates, integrating first with respect to 𝜌 , then with respect to 𝜙 , and finally with respect to 𝜃 , take the following steps.

  1. Sketch. Sketch the solid region D along with its projection R on the xy-plane. Label the surfaces that bound D.

教材插图

  1. Find the 𝜌 -limits of integration. Draw a ray M from the origin through D, making an angle 𝜙 with the positive z-axis. Also draw the projection of M on the xy-plane (call the projection L). The ray L makes an angle 𝜃 with the positive x-axis. As 𝜌 increases, M enters D at 𝜌 =𝑔1(𝜙,𝜃) and leaves at 𝜌 =𝑔2(𝜙,𝜃) . These are the 𝜌 -limits of integration shown in the above figure.

  2. Find the 𝜙 -limits of integration. For any given 𝜃 , the angle 𝜙 that M makes with the positive z-axis runs from 𝜙 =𝜙𝑚𝑖𝑛 to 𝜙 =𝜙𝑚𝑎𝑥 . The 𝜙 -limits of integration may depend on 𝜃 , but they are often constant.

  3. Find the 𝜃 -limits of integration. The ray L sweeps over R as 𝜃 runs from 𝛼 to 𝛽 . These are the 𝜃 -limits of integration. The integral is

∭𝐷𝑓(𝜌,𝜙,𝜃)𝑑𝑉=∫𝜃=𝛽𝜃=𝛼∫𝜙=𝜙max𝜙=𝜙min∫𝜌=𝑔2(𝜙,𝜃)𝜌=𝑔1(𝜙,𝜃)𝑓(𝜌,𝜙,𝜃)𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃.

EXAMPLE 5 Find the volume of the “ice cream cone” D bounded above by the sphere 𝜌 =1 and bounded below by the cone 𝜙 =𝜋/3 .

Solution The volume is 𝑉 =∭𝐷𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃 , the integral of 𝑓(𝜌,𝜙,𝜃) =1 over 𝐷 .

To find the limits of integration for evaluating the integral, we begin by sketching D and its projection R on the xy-plane (Figure 14.56).

The 𝜌 -limits of integration. We draw a ray M from the origin through D, making an angle 𝜙 with the positive z-axis. We also draw L, the projection of M on the xy-plane, along with the angle 𝜃 that L makes with the positive x-axis. Ray M enters D at 𝜌 =0 and leaves at 𝜌 =1 .

The 𝜙 -limits of integration. The cone 𝜙 =𝜋/3 makes an angle of 𝜋/3 with the positive z-axis. For any given 𝜃 , the angle 𝜙 can run from 𝜙 =0 to 𝜙 =𝜋/3 .

The 𝜃 -limits of integration. The ray L sweeps over R as 𝜃 runs from 0 to 2𝜋 . The volume is

𝑉=∭𝐷𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃=∫2𝜋0∫𝜋/30∫10𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃=∫2𝜋0∫𝜋/30[𝜌33]𝜌=1𝜌=0sin⁡𝜙𝑑𝜙𝑑𝜃=∫2𝜋0∫𝜋/3013sin⁡𝜙𝑑𝜙𝑑𝜃=∫2𝜋0[−13cos⁡𝜙]𝜙=𝜋/3𝜙=0𝑑𝜃=∫2𝜋0(−16+13)𝑑𝜃=16(2𝜋)=𝜋3.

EXAMPLE 6 A solid of constant density 𝛿 =1 occupies the solid region 𝐷 in Example 5. Find the solid’s moment of inertia about the 𝑧 -axis.

Solution In rectangular coordinates, the moment is

𝐼𝑧=∭𝐷(𝑥2+𝑦2)𝑑𝑉.

In spherical coordinates, 𝑥2 +𝑦2 =(𝜌sin⁡𝜙cos⁡𝜃)2 +(𝜌sin⁡𝜙sin⁡𝜃)2 =𝜌2sin2⁡𝜙 . Hence,

𝐼𝑧=∭𝐷(𝜌2sin2⁡𝜙)𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃=∭𝐷𝜌4sin3⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃.

For the region 𝐷 in Example 5, this becomes

𝐼𝑧=∫2𝜋0∫𝜋/30∫10𝜌4sin3⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃=∫2𝜋0∫𝜋/30[𝜌55]𝜌=1𝜌=0sin3⁡𝜙𝑑𝜙𝑑𝜃=15∫2𝜋0∫𝜋/30(1−cos2⁡𝜙)sin⁡𝜙𝑑𝜙𝑑𝜃=15∫2𝜋0[−cos⁡𝜙+cos3⁡𝜙3]𝜙=𝜋/3𝜙=0𝑑𝜃=15∫2𝜋0(−12+124+1−13)𝑑𝜃=15∫2𝜋0524𝑑𝜃=124(2𝜋)=𝜋12.

Coordinate Conversion Formulas

CYLINDRICAL TO

SPHERICAL TO

SPHERICAL TO

RECTANGULAR

RECTANGULAR

CYLINDRICAL

𝑥=𝑟cos⁡𝜃 𝑥=𝜌sin⁡𝜙cos⁡𝜃 𝑟=𝜌sin⁡𝜙 𝑦=𝑟sin⁡𝜃 𝑦=𝜌sin⁡𝜙sin⁡𝜃 𝑧=𝜌cos⁡𝜙 𝑧=𝑧 𝑧=𝜌cos⁡𝜙 𝜃=𝜃

Corresponding formulas for dV in triple integrals:

𝑑𝑉=𝑑𝑥𝑑𝑦𝑑𝑧=𝑟𝑑𝑧𝑑𝑟𝑑𝜃=𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃

In the next section we offer a more general procedure for determining dV in cylindrical and spherical coordinates. The results, of course, will be the same.

EXERCISES 14.7

In Exercises 1–12, sketch the region described by the following cylindrical coordinates in three-dimensional space.

  1. 𝑟 =2

  2. 𝜃 =𝜋4

  3. 𝑧 = −1

  4. 𝑧 =𝑟

  5. 𝑟 =𝜃

  6. 𝑧 =𝑟sin⁡𝜃

  7. 𝑟2 +𝑧2 =4

  8. 1 ≤𝑟 ≤2,0 ≤𝜃 ≤𝜋3

  9. 𝑟 ≤𝑧 ≤√9−𝑟2

  10. 0 ≤𝑟 ≤2sin⁡𝜃,1 ≤𝑧 ≤3

  11. 0 ≤𝑟 ≤4cos⁡𝜃,0 ≤𝜃 ≤𝜋2,0 ≤𝑧 ≤5

  12. 0 ≤𝑟 ≤3,−𝜋2 ≤𝜃 ≤𝜋2,0 ≤𝑧 ≤𝑟cos⁡𝜃

In Exercises 13–22, sketch the region described by the following spherical coordinates in three-dimensional space.

  1. 𝜌 =3

  2. 𝜙 =𝜋6

  3. 𝜃 =23𝜋

  4. 𝜌 =csc⁡𝜙

  5. 𝜌cos⁡𝜙 =4

  6. 1 ≤𝜌 ≤2sec⁡𝜙,0 ≤𝜙 ≤𝜋4

  7. 0 ≤𝜌 ≤3csc⁡𝜙

  8. 0 ≤𝜌 ≤1 , 𝜋2 ≤𝜙 ≤𝜋 , 0 ≤𝜃 ≤𝜋

  9. 0 ≤𝜌cos⁡𝜃sin⁡𝜙 ≤2,0 ≤𝜌sin⁡𝜃sin⁡𝜙 ≤3, 0 ≤𝜌cos⁡𝜙 ≤4

  10. 4sec⁡𝜙 ≤𝜌 ≤5,0 ≤𝜙 ≤𝜋2

Evaluating Integrals in Cylindrical Coordinates

Evaluate the cylindrical coordinate integrals in Exercises 23-28.

  1. ∫2𝜋0∫10∫√2−𝑟2𝑟𝑟𝑑𝑧𝑑𝑟𝑑𝜃

  2. ∫2𝜋0∫30∫√18−𝑟2𝑟2/3𝑟𝑑𝑧𝑑𝑟𝑑𝜃

  3. ∫2𝜋0∫𝜃/2𝜋0∫3+24𝑟20𝑟𝑑𝑧𝑑𝑟𝑑𝜃

  4. ∫𝜋0∫𝜃/𝜋0∫3√4−𝑟2−√4−𝑟2𝑧𝑟𝑑𝑧𝑑𝑟𝑑𝜃

  5. ∫2𝜋0∫10∫1/√2−𝑟2𝑟3𝑟𝑑𝑧𝑑𝑟𝑑𝜃

  6. ∫2𝜋0∫10∫1/2−1/2(𝑟2sin2⁡𝜃 +𝑧2)𝑟𝑑𝑧𝑑𝑟𝑑𝜃

The integrals we have seen so far suggest that there are preferred orders of integration for cylindrical coordinates, but other orders usually work well and are occasionally easier to evaluate. Evaluate the integrals in Exercises 29–32.

  1. ∫2𝜋0∫30∫𝑧/30𝑟3𝑑𝑟𝑑𝑧𝑑𝜃

  2. ∫1−1∫2𝜋0∫1+cos⁡𝜃04𝑟𝑑𝑟𝑑𝜃𝑑𝑧

  3. ∫10∫√𝑧0∫2𝜋0(𝑟2cos2⁡𝜃 +𝑧2)𝑟𝑑𝜃𝑑𝑟𝑑𝑧

  4. ∫20∫√4−𝑟2𝑟−2∫2𝜋0(𝑟sin⁡𝜃 +1)𝑟𝑑𝜃𝑑𝑧𝑑𝑟

  5. Let 𝐷 be the solid region bounded below by the plane 𝑧 =0 , above by the sphere 𝑥2 +𝑦2 +𝑧2 =4 , and on the sides by the cylinder 𝑥2 +𝑦2 =1 . Set up the triple integrals in cylindrical coordinates that give the volume of 𝐷 using the following orders of integration. a. 𝑑𝑧 𝑑𝑟 𝑑𝜃 b. 𝑑𝑟 𝑑𝑧 𝑑𝜃 c. 𝑑𝜃 𝑑𝑧 𝑑𝑟

  6. Let D be the solid region bounded below by the cone 𝑧 =√𝑥2+𝑦2 and above by the paraboloid 𝑧 =2 −𝑥2 −𝑦2 . Set up the triple integrals in cylindrical coordinates that give the volume of D using the following orders of integration. a. dz dr dθ b. dr dz dθ c. dθ dz dr

Finding Iterated Integrals in Cylindrical Coordinates

  1. Give the limits of integration for evaluating the integral
∭𝐷𝑓(𝑟,𝜃,𝑧)𝑟𝑑𝑧𝑑𝑟𝑑𝜃

as an iterated integral over the solid region D that is bounded below by the plane z = 0, on the side by the cylinder 𝑟 =cos⁡𝜃 , and on top by the paraboloid 𝑧 =3𝑟2 .

  1. Convert the integral
∫1−1∫√1−𝑦20∫𝑥0(𝑥2+𝑦2)𝑑𝑧𝑑𝑥𝑑𝑦

to an equivalent integral in cylindrical coordinates and evaluate the result.

In Exercises 37–42, set up the iterated integral for evaluating ∭𝐷𝑓(𝑟,𝜃,𝑧)𝑟𝑑𝑧𝑑𝑟𝑑𝜃 over the given solid region D.

  1. D is the right circular cylinder whose base is the circle 𝑟 =2sin⁡𝜃 in the xy-plane and whose top lies in the plane z = 4 - y.

教材插图

𝑟 =2cos⁡𝜃

  1. D is the right circular cylinder whose base is the circle 𝑟 =3cos⁡𝜃 and whose top lies in the plane z = 5 - x.

教材插图

  1. D is the solid right cylinder whose base is the region in the xy-plane that lies inside the cardioid 𝑟 =1 +cos⁡𝜃 and outside the circle r = 1 and whose top lies in the plane z = 4.

教材插图

  1. D is the solid right cylinder whose base is the region between the circles 𝑟 =cos⁡𝜃 and 𝑟 =2cos⁡𝜃 and whose top lies in the plane z = 3 - y.

教材插图

  1. D is the right prism whose base is the triangle in the xy-plane bounded by the x-axis and the lines y = x and x = 1 and whose top lies in the plane z = 2 - y.

教材插图

  1. D is the right prism whose base is the triangle in the xy-plane bounded by the y-axis and the lines y = x and y = 1 and whose top lies in the plane z = 2 - x.

教材插图

Evaluating Integrals in Spherical Coordinates

Evaluate the spherical coordinate integrals in Exercises 43-48.

  1. ∫𝜋0∫𝜋0∫2sin⁡𝜙0𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃

  2. ∫2𝜋0∫𝜋/40∫20(𝜌cos⁡𝜙)𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃

  3. ∫2𝜋0∫𝜋0∫(1−cos⁡𝜙)/20𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃

  4. ∫3𝜋/20∫𝜋0∫105𝜌3sin3⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃

  5. ∫2𝜋0∫𝜋/30∫2sec⁡𝜙3𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃

  6. ∫2𝜋0∫𝜋/40∫sec⁡𝜙0(𝜌cos⁡𝜙)𝜌2sin⁡𝜙𝑑𝜌𝑑𝜙𝑑𝜃

Changing the Order of Integration in Spherical Coordinates The previous integrals suggest there are preferred orders of integration for spherical coordinates, but other orders give the same value and are occasionally easier to evaluate. Evaluate the integrals in Exercises 49–52.

  1. ∫20∫0−𝜋∫𝜋/2𝜋/4𝜌3sin⁡2𝜙𝑑𝜙𝑑𝜃𝑑𝜌

  2. ∫𝜋/3𝜋/6∫2csc⁡𝜙csc⁡𝜙∫2𝜋0𝜌2sin⁡𝜙𝑑𝜃𝑑𝜌𝑑𝜙

  3. ∫10∫𝜋0∫𝜋/4012𝜌sin3⁡𝜙𝑑𝜙𝑑𝜃𝑑𝜌

  4. ∫𝜋/2𝜋/6∫𝜋/2−𝜋/2∫2csc⁡𝜙5𝜌4sin3⁡𝜙𝑑𝜌𝑑𝜃𝑑𝜙

  5. Let 𝐷 be the region in Exercise 33. Set up the triple integrals in spherical coordinates that give the volume of 𝐷 using the following orders of integration.

a. 𝑑𝜌𝑑𝜙𝑑𝜃

b. 𝑑𝜙 𝑑𝜌 𝑑𝜃

  1. Let D be the solid region bounded below by the cone 𝑧 =√𝑥2+𝑦2 and above by the plane z = 1. Set up the triple integrals in spherical coordinates that give the volume of D using the following orders of integration.

a. 𝑑𝜌 𝑑𝜙 𝑑𝜃

b. 𝑑𝜙 𝑑𝜌 𝑑𝜃

Finding Iterated Integrals in Spherical Coordinates

In Exercises 55–60, (a) find the spherical coordinate limits for the integral that calculates the volume of the given solid and then (b) evaluate the integral.

  1. The solid between the sphere 𝜌 =cos⁡𝜙 and the hemisphere 𝜌 =2,𝑧 ≥0

教材插图

  1. The solid bounded below by the hemisphere 𝜌 =1,𝑧 ≥0 , and above by the surface 𝜌 =1 +cos⁡𝜙

教材插图

  1. The solid enclosed by the surface 𝜌 =1 −cos⁡𝜙

  2. The upper portion cut from the solid in Exercise 57 by the xy-plane

  3. The solid bounded below by the sphere 𝜌 =2cos⁡𝜙 and above by the cone 𝑧 =√𝑥2+𝑦2

教材插图

  1. The solid bounded below by the xy-plane, on the sides by the sphere 𝜌 =2 , and above by the cone 𝜙 =𝜋/3

教材插图

Finding Triple Integrals

  1. Set up triple integrals for the volume of the sphere 𝜌 =2 in (a) spherical, (b) cylindrical, and (c) rectangular coordinates.

  2. Let D be the solid region in the first octant that is bounded below by the cone 𝜙 =𝜋/4 and above by the sphere 𝜌 =3 . Express the volume of D as an iterated triple integral in (a) cylindrical and (b) spherical coordinates. Then (c) find the volume.

  3. Let D be the smaller cap cut from a solid ball of radius 2 units by a plane 1 unit from the center of the sphere. Express the volume of D as an iterated triple integral in (a) spherical, (b) cylindrical, and (c) rectangular coordinates. Then (d) find the volume by evaluating one of the three triple integrals.

  4. Let 𝐷 be the solid hemisphere 𝑥2 +𝑦2 +𝑧2 ≤1 , 𝑧 ≥0 . If the density is 𝛿(𝑥,𝑦,𝑧) =1 , express the moment of intertia 𝐼𝑧 as an iterated integral in (a) cylindrical and (b) spherical coordinates. Then (c) find 𝐼𝑧 .

Volumes

Find the volumes of the solids in Exercises 65–70.

教材插图

教材插图

教材插图

教材插图

教材插图

教材插图

  1. Ball and cones Find the volume of the portion of the ball 𝜌 ≤𝑎 that lies between the cones 𝜙 =𝜋/3 and 𝜙 =2𝜋/3 .

  2. Ball and half-planes Find the volume of the region cut from the ball 𝜌 ≤𝑎 by the half-planes 𝜃 =0 and 𝜃 =𝜋/6 in the first octant.

  3. Ball and plane Find the volume of the smaller region cut from the ball 𝜌 ≤√2 by the plane z = 1.

  4. Cone and planes Find the volume of the solid enclosed by the cone 𝑧 =√𝑥2+𝑦2 between the planes z = 1 and z = 2.

  5. Cylinder and paraboloid Find the volume of the solid region bounded below by the plane z = 0, laterally by the cylinder 𝑥2 +𝑦2 =1 , and above by the paraboloid 𝑧 =𝑥2 +𝑦2 .

  6. Cylinder and paraboloids Find the volume of the solid region bounded below by the paraboloid 𝑧 =𝑥2 +𝑦2 , laterally by the cylinder 𝑥2 +𝑦2 =1 , and above by the paraboloid 𝑧 =𝑥2 +𝑦2 +1 .

  7. Cylinder and cones Find the volume of the solid cut from the thick-walled cylinder 1 ≤𝑥2 +𝑦2 ≤2 by the cones 𝑧 = ±√𝑥2+𝑦2 .

  8. Sphere and cylinder Find the volume of the solid region that lies inside the sphere 𝑥2 +𝑦2 +𝑧2 =2 and outside the cylinder 𝑥2 +𝑦2 =1 .

  9. Cylinder and planes Find the volume of the solid region enclosed by the cylinder 𝑥2 +𝑦2 =4 and the planes z = 0 and 𝑦 +𝑧 =4 .

  10. Cylinder and planes Find the volume of the solid region enclosed by the cylinder 𝑥2 +𝑦2 =4 and the planes z = 0 and 𝑥 +𝑦 +𝑧 =4 .

  11. Region trapped by paraboloids Find the volume of the solid region bounded above by the paraboloid 𝑧 =5 −𝑥2 −𝑦2 and below by the paraboloid 𝑧 =4𝑥2 +4𝑦2 .

  12. Paraboloid and cylinder Find the volume of the solid region bounded above by the paraboloid 𝑧 =9 −𝑥2 −𝑦2 , bounded below by the xy-plane, and lying outside the cylinder 𝑥2 +𝑦2 =1 .

  13. Cylinder and sphere Find the volume of the region cut from the solid cylinder 𝑥2 +𝑦2 ≤1 by the sphere 𝑥2 +𝑦2 +𝑧2 =4 .

  14. Sphere and paraboloid Find the volume of the solid region bounded above by the sphere 𝑥2 +𝑦2 +𝑧2 =2 and below by the paraboloid 𝑧 =𝑥2 +𝑦2 .

Average Values

  1. Find the average value of the function 𝑓(𝑟,𝜃,𝑧) =𝑟 over the solid region bounded by the cylinder 𝑟 =1 between the planes 𝑧 = −1 and 𝑧 =1 .

  2. Find the average value of the function 𝑓(𝑟,𝜃,𝑧) =𝑟 over the solid ball bounded by the sphere 𝑟2 +𝑧2 =1 . (This is the sphere 𝑥2 +𝑦2 +𝑧2 =1 .)

  3. Find the average value of the function 𝑓(𝜌,𝜙,𝜃) =𝜌 over the solid ball 𝜌 ≤1 .

  4. Find the average value of the function 𝑓(𝜌,𝜙,𝜃) =𝜌cos⁡𝜙 over the upper half of the solid ball 𝜌 ≤1,0 ≤𝜙 ≤𝜋/2 .

Masses, Moments, and Centroids

  1. Center of mass A solid of constant density is bounded below by the plane z = 0, above by the cone z = r, 𝑟 ≥0 , and on the sides by the cylinder r = 1. Find the center of mass.

  2. Centroid Find the centroid of the solid region in the first octant that is bounded above by the cone 𝑧 =√𝑥2+𝑦2 , below by the plane z = 0, and on the sides by the cylinder 𝑥2 +𝑦2 =4 and the planes x = 0 and y = 0.

  3. Centroid Find the centroid of the solid in Exercise 60.

  4. Centroid Find the centroid of the solid bounded above by the sphere 𝜌 =𝑎 and below by the cone 𝜙 =𝜋/4 .

  5. Centroid Find the centroid of the solid region that is bounded above by the surface 𝑧 =√𝑟 , on the sides by the cylinder r = 4, and below by the xy-plane.

  6. Centroid Find the centroid of the region cut from the solid ball 𝑟2 +𝑧2 ≤1 by the half-planes 𝜃 = −𝜋/3 , 𝑟 ≥0 , and 𝜃 =𝜋/3 , 𝑟 ≥0 .

  7. Moment of inertia of solid cone Find the moment of inertia of a solid right circular cone of base radius 1 and height 1 about an axis through the vertex parallel to the base if the density is 𝛿 =1 .

  8. Moment of inertia of ball Find the moment of inertia of a ball of radius 𝑎 about a diameter if the density is 𝛿 =1 .

  9. Moment of inertia of solid cone Find the moment of inertia of a solid right circular cone of base radius 𝑎 and height ℎ about its axis if the density is 𝛿 =1 . (Hint: Place the cone with its vertex at the origin and its axis along the 𝑧 -axis.)

  10. Variable density A solid is bounded on the top by the paraboloid 𝑧 =𝑟2 , on the bottom by the plane 𝑧 =0 , and on the sides by the cylinder 𝑟 =1 . Find the center of mass and the moment of inertia about the 𝑧 -axis if the density is a. 𝛿(𝑟,𝜃,𝑧) =𝑧 b. 𝛿(𝑟,𝜃,𝑧) =𝑟 .

  11. Variable density A solid is bounded below by the cone 𝑧 =√𝑥2+𝑦2 and above by the plane 𝑧 =1 . Find the center of mass and the moment of inertia about the 𝑧 -axis if the density is

𝛿(𝑟,𝜃,𝑧)=𝑧 𝛿(𝑟,𝜃,𝑧)=𝑧2.
  1. Variable density A solid ball is bounded by the sphere 𝜌 =𝑎 . Find the moment of inertia about the 𝑧 -axis if the density is
    a. 𝛿(𝜌,𝜙,𝜃) =𝜌2 b. 𝛿(𝜌,𝜙,𝜃) =𝑟 =𝜌sin⁡𝜙 .

  2. Centroid of solid semi-ellipsoid Show that the centroid of the solid semi-ellipsoid of revolution (𝑟2/𝑎2) +(𝑧2/ℎ2) ≤1 , 𝑧 ≥0 , lies on the z-axis three-eighths of the way from the base to the top. The special case h = a gives a solid hemisphere. Thus, the centroid of a solid hemisphere lies on the axis of symmetry three-eighths of the way from the base to the top.

  3. Centroid of solid cone Show that the centroid of a solid right circular cone is one-fourth of the way from the base to the vertex. (In general, the centroid of a solid cone or pyramid is one-fourth of the way from the centroid of the base to the vertex.)

  4. Density of center of a planet A planet is in the shape of a sphere of radius R and total mass M with spherically symmetric density distribution that increases linearly as one approaches its center. What is the density at the center of this planet if the density at its edge (surface) is taken to be zero?

  5. Mass of planet’s atmosphere A spherical planet of radius 𝑅 has an atmosphere whose density is 𝜇 =𝜇0𝑒−𝑐ℎ , where ℎ is the altitude above the surface of the planet, 𝜇0 is the density at sea level, and 𝑐 is a positive constant. Find the mass of the planet’s atmosphere.

Theory and Examples

  1. Vertical planes in cylindrical coordinates

a. Show that planes perpendicular to the x-axis have equations of the form 𝑟 =𝑎sec⁡𝜃 in cylindrical coordinates.

b. Show that planes perpendicular to the y-axis have equations of the form 𝑟 =𝑏csc⁡𝜃 .

  1. (Continuation of Exercise 105.) Find an equation of the form 𝑟 =𝑓(𝜃) in cylindrical coordinates for the plane 𝑎𝑥 +𝑏𝑦 =𝑐 , 𝑐 ≠0 .

  2. Symmetry What symmetry will you find in a surface that has an equation of the form 𝑟 =𝑓(𝑧) in cylindrical coordinates? Give reasons for your answer.

  3. Symmetry What symmetry will you find in a surface that has an equation of the form 𝜌 =𝑓(𝜙) in spherical coordinates? Give reasons for your answer.

14.8 Substitutions in Multiple Integrals

教材插图

Cartesian uv-plane

教材插图

Cartesian xy-plane

FIGURE 14.57 The equations

𝑥 =𝑔(𝑢,𝑣) and 𝑦 =ℎ(𝑢,𝑣) allow us to change an integral over a region R in the xy-plane into an integral over a region G in the uv-plane.

HISTORICAL BIOGRAPHY Carl Gustav Jacob Jacobi (1804–1851)

Jacobi, one of nineteenth-century Germany’s most accomplished scientists, developed the theory of determinants and transformations into a powerful tool for evaluating multiple integrals and solving differential equations. He also applied transformation methods to study integrals like the ones that arise in the calculation of arc length.

This section introduces the ideas involved in coordinate transformations to evaluate multiple integrals by substitution. The method replaces complicated integrals by ones that are easier to evaluate. Substitutions accomplish this by simplifying the integrand, the limits of integration, or both. A thorough discussion of multivariable transformations and substitutions is best left to a more advanced course, but our introduction here shows how the substitutions just studied reflect the general idea derived for single integral calculus.

To know more, visit the companion Website.

Substitutions in Double Integrals

The polar coordinate substitution of Section 14.4 is a special case of a more general substitution method for double integrals, a method that pictures changes in variables as transformations of regions.

Suppose that a region G in the uv-plane is transformed into the region R in the xy-plane by equations of the form

𝑥=𝑔(𝑢,𝑣),𝑦=ℎ(𝑢,𝑣),

as suggested in Figure 14.57. We assume the transformation is one-to-one on the interior of 𝐺 . We call 𝑅 the image of 𝐺 under the transformation, and 𝐺 the preimage of 𝑅 . Any function 𝑓(𝑥,𝑦) defined on 𝑅 can be thought of as a function 𝑓(𝑔(𝑢,𝑣),ℎ(𝑢,𝑣)) defined on 𝐺 as well. How is the integral of 𝑓(𝑥,𝑦) over 𝑅 related to the integral of 𝑓(𝑔(𝑢,𝑣),ℎ(𝑢,𝑣)) over 𝐺 ?

To gain some insight into the question, we look again at the single variable case. To be consistent with how we are using them now, we interchange the variables x and u used in the substitution method for single integrals in Chapter 5, so the equation is

∫𝑔(𝑏)𝑔(𝑎)𝑓(𝑥)𝑑𝑥=∫𝑏𝑎𝑓(𝑔(𝑢))𝑔′(𝑢)𝑑𝑢.𝑥=𝑔(𝑢),𝑑𝑥=𝑔′(𝑢)𝑑𝑢

To propose an analogue for substitution in a double integral ∬𝑅𝑓(𝑥,𝑦) 𝑑𝑥 𝑑𝑦 , we need a derivative factor like 𝑔′(𝑢) as a multiplier that transforms the area element du dv in the region G to its corresponding area element dx dy in the region R. We denote this factor by J. In continuing with our analogy, it is reasonable to assume that J is a function of both variables u and v, just as 𝑔′ is a function of the single variable u. Moreover, J should register instantaneous change, so partial derivatives are going to be involved in its expression. Since four partial derivatives are associated with the transforming equations 𝑥 =𝑔(𝑢,𝑣) and 𝑦 =ℎ(𝑢,𝑣) , it is also reasonable to assume that the factor 𝐽(𝑢,𝑣) we seek includes them all. These features are captured in the following definition, which is constructed from the partial derivatives and is named after the German mathematician Carl Jacobi.

Differential Area Change Substituting 𝑥 =𝑔(𝑢,𝑣),𝑦 =ℎ(𝑢,𝑣)

𝑑𝑥𝑑𝑦=∣𝜕(𝑥,𝑦)𝜕(𝑢,𝑣)∣𝑑𝑢𝑑𝑣

教材插图

↑ ↑ ↑ ↑ ↑ ↑↓𝑥=𝑟cos⁡𝜃↑ ↑ ↑ ↑ ↑ ↑↓𝑦=𝑟sin⁡𝜃

教材插图

𝑥 =𝑟cos⁡𝜃,𝑦 =𝑟sin⁡𝜃 transform G into R. The Jacobian factor r, calculated in Example 1, scales the differential rectangle dr dθ in G to match the differential area element dx dy in R.

DEFINITION The Jacobian determinant or Jacobian of the coordinate transformation 𝑥 =𝑔(𝑢,𝑣) , 𝑦 =ℎ(𝑢,𝑣) is

𝐽(𝑢,𝑣)=∣𝜕𝑥𝜕𝑢𝜕𝑥𝜕𝑣𝜕𝑦𝜕𝑢𝜕𝑦𝜕𝑣∣=𝜕𝑥𝜕𝑢𝜕𝑦𝜕𝑣−𝜕𝑦𝜕𝑢𝜕𝑥𝜕𝑣.(1)

The Jacobian can also be denoted by

𝐽(𝑢,𝑣)=𝜕(𝑥,𝑦)𝜕(𝑢,𝑣)

to help us remember how the determinant in Equation (1) is constructed from the partial derivatives of x and y. The array of partial derivatives in Equation (1) behaves just like the derivative 𝑔′ in the single variable situation. The Jacobian measures how much the transformation is expanding or contracting the area around the point (𝑢,𝑣) . Effectively, the factor |𝐽| converts the area of the differential rectangle du dv in G to match its corresponding differential area dx dy in R. We note that, in general, the value of the scaling factor |𝐽| depends on the point (𝑢,𝑣) in G; that is, the scaling changes as the point (𝑢,𝑣) varies through the region G. Our examples to follow will show how it scales the differential area du dv for specific transformations.

Now we can answer our original question concerning the relationship of the integral of 𝑓(𝑥,𝑦) over the region R to the integral of 𝑓(𝑔(𝑢,𝑣),ℎ(𝑢,𝑣)) over G.

THEOREM 3—Substitution for Double Integrals

Suppose that 𝑓(𝑥,𝑦) is continuous over the region R. Let G be the preimage of R under the transformation 𝑥 =𝑔(𝑢,𝑣) , 𝑦 =ℎ(𝑢,𝑣) , which is assumed to be one-to-one on the interior of G. If the functions g and h have continuous first partial derivatives within the interior of G, then

∬𝑅𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦=∬𝐺𝑓(𝑔(𝑢,𝑣),ℎ(𝑢,𝑣))∣𝜕(𝑥,𝑦)𝜕(𝑢,𝑣)∣𝑑𝑢𝑑𝑣.(2)

The derivation of Equation (2) is intricate and properly belongs to a course in advanced calculus, so we do not include it here. We now present examples illustrating the substitution method defined by the equation.

EXAMPLE 1 Find the Jacobian for the polar coordinate transformation 𝑥 =𝑟cos⁡𝜃 , 𝑦 =𝑟sin⁡𝜃 , and use Equation (2) to write the Cartesian integral ∬𝑅𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦 as a polar integral.

Solution Figure 14.58 shows how the equations 𝑥 =𝑟cos⁡𝜃 , 𝑦 =𝑟sin⁡𝜃 transform the rectangle 𝐺 :0 ≤𝑟 ≤1,0 ≤𝜃 ≤𝜋/2 , into the quarter of a circular disk 𝑅 bounded by 𝑥2 +𝑦2 =1 in the first quadrant of the 𝑥𝑦 -plane.

For polar coordinates, we have 𝑟 and 𝜃 in place of 𝑢 and 𝑣 . With 𝑥 =𝑟cos⁡𝜃 and 𝑦 =𝑟sin⁡𝜃 , the Jacobian is

𝐽(𝑟,𝜃)=∣𝜕𝑥𝜕𝑟𝜕𝑥𝜕𝜃𝜕𝑦𝜕𝑟𝜕𝑦𝜕𝜃∣=∣cos⁡𝜃−𝑟sin⁡𝜃sin⁡𝜃𝑟cos⁡𝜃∣=𝑟(cos2⁡𝜃+sin2⁡𝜃)=𝑟.

Since we assume 𝑟 ≥0 when integrating in polar coordinates, |𝐽(𝑟,𝜃)| =|𝑟| =𝑟 so that Equation (2) gives

∬𝑅𝑓(𝑥,𝑦)𝑑𝑥𝑑𝑦=∬𝐺𝑓(𝑟cos⁡𝜃,𝑟sin⁡𝜃)𝑟𝑑𝑟𝑑𝜃.(3)

This is the same formula we derived independently using a geometric argument for polar area in Section 14.4.

Here is an example of a substitution in which the image of a rectangle under the coordinate transformation is a trapezoid. Transformations like this one are called linear transformations, and their Jacobians are constant throughout G.

EXAMPLE 2 Evaluate

∫40∫𝑥=(𝑦/2)+1𝑥=𝑦/22𝑥−𝑦2𝑑𝑥𝑑𝑦

by applying the transformation

𝑢=2𝑥−𝑦2,𝑣=𝑦2(4)

and integrating over an appropriate region in the uv-plane.

Solution We sketch the region R of integration in the xy-plane and identify its boundaries (Figure 14.59).

教材插图

FIGURE 14.59 The equations 𝑥 =𝑢 +𝑣 and y = 2v transform G into R. Reversing the transformation by the equations 𝑢 =(2𝑥 −𝑦)/2 and v = y/2 transforms R into G (Example 2).

To apply Equation (2), we need to find the corresponding uv-region G and the Jacobian of the transformation. To find them, we first solve Equations (4) for x and y in terms of u and v. From those equations it is easy to find algebraically that

𝑥=𝑢+𝑣,𝑦=2𝑣.(5)

We then find the boundaries of G by substituting these expressions into the equations for the boundaries of R (Figure 14.59)

xy-equations for the boundary of RCorresponding uv-equations for the boundary of GSimplified uv-equations
x = y/2u + v = 2v/2 = vu = 0
x = (y/2) + 1u + v = (2v/2) + 1 = v + 1u = 1
y = 02v = 0v = 0
y = 42v = 4v = 2

教材插图

FIGURE 14.60 The equations 𝑥 =(𝑢/3) −(𝑣/3) and 𝑦 =(2𝑢/3) +(𝑣/3) transform G into R. Reversing the transformation by the equations 𝑢 =𝑥 +𝑦 and 𝑣 =𝑦 −2𝑥 transforms R into G (Example 3).

From Equations (5) the Jacobian of the transformation is

𝐽(𝑢,𝑣)=∣𝜕𝑥𝜕𝑢𝜕𝑥𝜕𝑣𝜕𝑦𝜕𝑢𝜕𝑦𝜕𝑣∣=∣𝜕𝜕𝑢(𝑢+𝑣)𝜕𝜕𝑣(𝑢+𝑣)𝜕𝜕𝑢(2𝑣)𝜕𝜕𝑣(2𝑣)∣=∣1102∣=2.

We now have everything we need to apply Equation (2):

∫40∫𝑥=(𝑦/2)+1𝑥=𝑦/22𝑥−𝑦2𝑑𝑥𝑑𝑦=∫𝑣=2𝑣=0∫𝑢=1𝑢=0𝑢|𝐽(𝑢,𝑣)|𝑑𝑢𝑑𝑣=∫20∫10(𝑢)(2)𝑑𝑢𝑑𝑣=∫20[𝑢2]𝑢=1𝑢=0𝑑𝑣=∫20𝑑𝑣=2.

EXAMPLE 3 Evaluate

∫10∫1−𝑥0√𝑥+𝑦(𝑦−2𝑥)2𝑑𝑦𝑑𝑥.

Solution We sketch the region R of integration in the xy-plane and identify its boundaries (Figure 14.60). The integrand suggests the transformation 𝑢 =𝑥 +𝑦 and 𝑣 =𝑦 −2𝑥 . Routine algebra produces x and y as functions of u and v:

𝑥=𝑢3−𝑣3,𝑦=2𝑢3+𝑣3.(6)

From Equations (6), we can find the boundaries of the uv-region G (Figure 14.60).

xy-equations for the boundary of RCorresponding uv-equations for the boundary of GSimplified uv-equations
x + y = 1(𝑢3−𝑣3) +(2𝑢3+𝑣3) =1u = 1
x = 0𝑢3 −𝑣3 =0v = u
y = 02𝑢3 +𝑣3 =0v = -2u

The Jacobian of the transformation in Equations (6) is

𝐽(𝑢,𝑣)=∣𝜕𝑥𝜕𝑢𝜕𝑥𝜕𝑣𝜕𝑦𝜕𝑢𝜕𝑦𝜕𝑣∣=∣13−132313∣=13.

Applying Equation (2), we evaluate the integral:

∫10∫1−𝑥0√𝑥+𝑦(𝑦−2𝑥)2𝑑𝑦𝑑𝑥=∫𝑢=1𝑢=0∫𝑣=𝑢𝑣=−2𝑢𝑢1/2𝑣2|𝐽(𝑢,𝑣)|𝑑𝑣𝑑𝑢=∫10∫𝑢−2𝑢𝑢1/2𝑣2(13)𝑑𝑣𝑑𝑢=13∫10𝑢1/2[13𝑣3]𝑣=𝑢𝑣=−2𝑢𝑑𝑢=19∫10𝑢1/2(𝑢3+8𝑢3)𝑑𝑢=∫10𝑢7/2𝑑𝑢=29𝑢9/2∣10=29.

In the next example we illustrate a nonlinear transformation of coordinates resulting from simplifying the form of the integrand. Like the polar coordinates’ transformation, nonlinear transformations can map a straight-line boundary of a region into a curved boundary (or vice versa with the inverse transformation). In general, nonlinear transformations are more complex to analyze than linear ones, and a complete treatment is left to a more advanced course.

教材插图

FIGURE 14.61 The region of integration R in Example 4.

教材插图

FIGURE 14.62 The boundaries of the region G correspond to those of region R in Figure 14.61. Notice that as we move counterclockwise around the region R, we move counterclockwise around the region G as well. The inverse transformation equations 𝑢 =√𝑥𝑦 , 𝑣 =√𝑦/𝑥 produce the region G from the region R.

EXAMPLE 4 Evaluate the integral

∫21∫𝑦1/𝑦√𝑦𝑥𝑒√𝑥𝑦𝑑𝑥𝑑𝑦.

Solution The square root terms in the integrand suggest that we might simplify the integration by substituting 𝑢 =√𝑥𝑦 and 𝑣 =√𝑦/𝑥 . Squaring these equations gives 𝑢2 =𝑥𝑦 and 𝑣2 =𝑦/𝑥 , which imply that 𝑢2𝑣2 =𝑦2 and 𝑢2/𝑣2 =𝑥2 . So we obtain the transformation (in the same ordering of the variables as discussed before)

𝑥=𝑢𝑣 and 𝑦=𝑢𝑣,

with 𝑢 >0 and 𝑣 >0 . Let’s first see what happens to the integrand itself under this transformation. The Jacobian of the transformation is not constant:

𝐽(𝑢,𝑣)=∣𝜕𝑥𝜕𝑢𝜕𝑥𝜕𝑣𝜕𝑦𝜕𝑢𝜕𝑦𝜕𝑣∣=∣1𝑣−𝑢𝑣2𝑣𝑢∣=2𝑢𝑣.

If G is the region of integration in the uv-plane, then by Equation (2) the transformed double integral under the substitution is

∬𝑅√𝑦𝑥𝑒√𝑥𝑦𝑑𝑥𝑑𝑦=∬𝐺𝑣𝑒𝑢|𝐽(𝑢,𝑣)|𝑑𝑢𝑑𝑣=∬𝐺𝑣𝑒𝑢2𝑢𝑣𝑑𝑢𝑑𝑣=∬𝐺2𝑢𝑒𝑢𝑑𝑢𝑑𝑣.

The transformed integrand function is easier to integrate than the original one, so we proceed to determine the limits of integration for the transformed integral.

The region of integration R of the original integral in the xy-plane is shown in Figure 14.61. From the substitution equations 𝑢 =√𝑥𝑦 and 𝑣 =√𝑦/𝑥 , we see that the image of the left-hand boundary xy = 1 for R is the vertical line segment u = 1, 2 ≥𝑣 ≥1 , in G (see Figure 14.62). Likewise, the right-hand boundary y = x of R maps to the horizontal line segment 𝑣 =1,1 ≤𝑢 ≤2 , in G. Finally, the horizontal top boundary y = 2 of R maps to 𝑢𝑣 =2,1 ≤𝑣 ≤2 , in G. As we move counterclockwise around the boundary of the region R, we also move counterclockwise around the boundary of G, as shown in Figure 14.62. Knowing the region of integration G in the uv-plane, we can now write equivalent iterated integrals:

∫21∫𝑦1/𝑦√𝑦𝑥𝑒√𝑥𝑦𝑑𝑥𝑑𝑦=∫21∫2/𝑢12𝑢𝑒𝑢𝑑𝑣𝑑𝑢.

Note the order of integration.

We now evaluate the transformed integral on the right-hand side:

∫21∫2/𝑢12𝑢𝑒𝑢𝑑𝑣𝑑𝑢=2∫21[𝑣𝑢𝑒𝑢]𝑣=2/𝑢𝑣=1𝑑𝑢=2∫21(2𝑒𝑢−𝑢𝑒𝑢)𝑑𝑢=2∫21(2−𝑢)𝑒𝑢𝑑𝑢=2[(2−𝑢)𝑒𝑢+𝑒𝑢]𝑢=2𝑢=1=2(𝑒2−(𝑒+𝑒))=2𝑒(𝑒−2).

Integrate by parts.

Determinants

2 ×2 and 3 ×3 determinants are evaluated as follows:

∣𝑎𝑏𝑐𝑑∣=𝑎𝑑−𝑏𝑐∣ ∣ ∣ ∣𝑎1𝑎2𝑎3𝑏1𝑏2𝑏3𝑐1𝑐2𝑐3∣ ∣ ∣ ∣=𝑎1∣𝑏2𝑏3𝑐2𝑐3∣−𝑎2∣𝑏1𝑏3𝑐1𝑐3∣+𝑎3∣𝑏1𝑏2𝑐1𝑐2∣

Substitutions in Triple Integrals

The cylindrical and spherical coordinate substitutions in Section 14.7 are special cases of a substitution method that pictures changes of variables in triple integrals as transformations of solid regions. The method is like the method for double integrals given by Equation (2) except that now we work in three dimensions instead of two.

Suppose that a solid region 𝐺 in uvw-space is transformed one-to-one into the solid region 𝐷 in xyz-space by differentiable equations of the form

𝑥=𝑔(𝑢,𝑣,𝑤),𝑦=ℎ(𝑢,𝑣,𝑤),𝑧=𝑘(𝑢,𝑣,𝑤),

as suggested in Figure 14.63. Then any function 𝐹(𝑥,𝑦,𝑧) defined on 𝐷 can be thought of as a function

𝐹(𝑔(𝑢,𝑣,𝑤),ℎ(𝑢,𝑣,𝑤),𝑘(𝑢,𝑣,𝑤))=𝐻(𝑢,𝑣,𝑤)

defined on G. If g, h, and k have continuous first partial derivatives, then the integral of 𝐹(𝑥,𝑦,𝑧) over D is related to the integral of 𝐻(𝑢,𝑣,𝑤) over G by the equation

∭𝐷𝐹(𝑥,𝑦,𝑧)𝑑𝑥𝑑𝑦𝑑𝑧=∭𝐺𝐻(𝑢,𝑣,𝑤)|𝐽(𝑢,𝑣,𝑤)|𝑑𝑢𝑑𝑣𝑑𝑤.(7)

教材插图

FIGURE 14.63 The equations 𝑥 =𝑔(𝑢,𝑣,𝑤) , 𝑦 =ℎ(𝑢,𝑣,𝑤) , and 𝑧 =𝑘(𝑢,𝑣,𝑤) allow us to change an integral over a region D in Cartesian xyz-space into an integral over a region G in Cartesian uvw-space using Equation (7).

The factor 𝐽(𝑢,𝑣,𝑤) , whose absolute value appears in this equation, is the Jacobian determinant

𝐽(𝑢,𝑣,𝑤)=∣ ∣ ∣ ∣𝜕𝑥𝜕𝑢𝜕𝑥𝜕𝑣𝜕𝑥𝜕𝑤𝜕𝑦𝜕𝑢𝜕𝑦𝜕𝑣𝜕𝑦𝜕𝑤𝜕𝑧𝜕𝑢𝜕𝑧𝜕𝑣𝜕𝑧𝜕𝑤∣ ∣ ∣ ∣=𝜕(𝑥,𝑦,𝑧)𝜕(𝑢,𝑣,𝑤).

This determinant measures how much the volume near a point in G is being expanded or contracted by the transformation from (𝑢,𝑣,𝑤) to (𝑥,𝑦,𝑧) coordinates. As in the two-dimensional case, the derivation of the change-of-variable formula in Equation (7) is omitted.

For cylindrical coordinates, 𝑟,𝜃 , and z take the place of u, v, and w. The transformation from Cartesian 𝑟𝜃𝑧 -space to Cartesian xyz-space is given by the equations

𝑥=𝑟cos⁡𝜃,𝑦=𝑟sin⁡𝜃,𝑧=𝑧

教材插图

FIGURE 14.64 The equations 𝑥 =𝑟cos⁡𝜃 , 𝑦 =𝑟sin⁡𝜃 , and z = z transform the rectangular box G into a cylindrical wedge D.

(Figure 14.64). The Jacobian of the transformation is

𝐽(𝑟,𝜃,𝑧)=∣ ∣ ∣ ∣𝜕𝑥𝜕𝑟𝜕𝑥𝜕𝜃𝜕𝑥𝜕𝑧𝜕𝑦𝜕𝑟𝜕𝑦𝜕𝜃𝜕𝑦𝜕𝑧𝜕𝑧𝜕𝑟𝜕𝑧𝜕𝜃𝜕𝑧𝜕𝑧∣ ∣ ∣ ∣=∣ ∣ ∣ ∣cos⁡𝜃−𝑟sin⁡𝜃0sin⁡𝜃𝑟cos⁡𝜃0001∣ ∣ ∣ ∣=𝑟cos2⁡𝜃+𝑟sin2⁡𝜃=𝑟.

The corresponding version of Equation (7) is

∭𝐷𝐹(𝑥,𝑦,𝑧)𝑑𝑥𝑑𝑦𝑑𝑧=∭𝐺𝐻(𝑟,𝜃,𝑧)|𝑟|𝑑𝑟𝑑𝜃𝑑𝑧.

We can drop the absolute value signs because 𝑟 ≥0 .

For spherical coordinates, 𝜌 , 𝜙 , and 𝜃 take the place of u, v, and w. The transformation from Cartesian 𝜌𝜙𝜃 -space to Cartesian xyz-space is given by

𝑥=𝜌sin⁡𝜙cos⁡𝜃,𝑦=𝜌sin⁡𝜙sin⁡𝜃,𝑧=𝜌cos⁡𝜙

(Figure 14.65). The Jacobian of the transformation (see Exercise 23) is

𝐽(𝜌,𝜙,𝜃)=∣ ∣ ∣ ∣ ∣𝜕𝑥𝜕𝜌𝜕𝑥𝜕𝜙𝜕𝑥𝜕𝜃𝜕𝑦𝜕𝜌𝜕𝑦𝜕𝜙𝜕𝑦𝜕𝜃𝜕𝑧𝜕𝜌𝜕𝑧𝜕𝜙𝜕𝑧𝜕𝜃∣ ∣ ∣ ∣ ∣=𝜌2sin⁡𝜙.

The corresponding version of Equation (7) is

∭𝐷𝐹(𝑥,𝑦,𝑧)𝑑𝑥𝑑𝑦𝑑𝑧=∭𝐺𝐻(𝜌,𝜙,𝜃)|𝜌2sin⁡𝜙|𝑑𝜌𝑑𝜙𝑑𝜃.

教材插图

FIGURE 14.65 The equations 𝑥 =𝜌sin⁡𝜙cos⁡𝜃 , 𝑦 =𝜌sin⁡𝜙sin⁡𝜃 , and 𝑧 =𝜌cos⁡𝜙 transform the rectangular box G into the spherical wedge D.

We can drop the absolute value signs because sin⁡𝜙 is never negative for 0 ≤𝜙 ≤𝜋 . Note that this is the same result we obtained in Section 14.7.

Here is an example of another substitution. Although we could evaluate the integral in this example directly, we have chosen it to illustrate the substitution method in a simple (and fairly intuitive) setting.

教材插图

教材插图

FIGURE 14.66 The equations 𝑥 =𝑢 +𝑣,𝑦 =2𝑣 , and 𝑧 =3𝑤 transform 𝐺 into 𝐷 . Reversing the transformation by the equations 𝑢 =(2𝑥 −𝑦)/2,𝑣 =𝑦/2 , and 𝑤 =𝑧/3 transforms 𝐷 into 𝐺 (Example 5).

EXAMPLE 5 Evaluate

∫30∫40∫𝑥=(𝑦/2)+1𝑥=𝑦/2(2𝑥−𝑦2+𝑧3)𝑑𝑥𝑑𝑦𝑑𝑧

by applying the transformation

𝑢=(2𝑥−𝑦)/2,𝑣=𝑦/2,𝑤=𝑧/3(8)

and integrating over an appropriate region in uvw-space.

Solution We sketch the solid region D of integration in xyz-space and identify its boundaries (Figure 14.66). In this case, the bounding surfaces are planes.

To apply Equation (7), we need to find the corresponding uvw-region G and the Jacobian of the transformation. To find them, we first solve Equations (8) for x, y, and z in terms of u, v, and w. Routine algebra gives

𝑥=𝑢+𝑣,𝑦=2𝑣,𝑧=3𝑤.(9)

We then find the boundaries of G by substituting these expressions into the equations for the boundaries of D.

xyz-equations for the boundary of DCorresponding uvw-equations for the boundary of GSimplified uvw-equations
x = y/2u + v = 2v/2 = vu = 0
x = (y/2) + 1u + v = (2v/2) + 1 = v + 1u = 1
y = 02v = 0v = 0
y = 42v = 4v = 2
z = 03w = 0w = 0
z = 33w = 3w = 1

The Jacobian of the transformation, again from Equations (9), is

𝐽(𝑢,𝑣,𝑤)=∣ ∣ ∣ ∣𝜕𝑥𝜕𝑢𝜕𝑥𝜕𝑣𝜕𝑥𝜕𝑤𝜕𝑦𝜕𝑢𝜕𝑦𝜕𝑣𝜕𝑦𝜕𝑤𝜕𝑧𝜕𝑢𝜕𝑧𝜕𝑣𝜕𝑧𝜕𝑤∣ ∣ ∣ ∣=∣ ∣ ∣ ∣110020003∣ ∣ ∣ ∣=6.

We now have everything we need to apply Equation (7):

∫30∫40∫𝑥=(𝑦/2)+1𝑥=𝑦/2(2𝑥−𝑦2+𝑧3)𝑑𝑥𝑑𝑦𝑑𝑧=∫10∫20∫10(𝑢+𝑤)|𝐽(𝑢,𝑣,𝑤)|𝑑𝑢𝑑𝑣𝑑𝑤=∫10∫20∫10(𝑢+𝑤)(6)𝑑𝑢𝑑𝑣𝑑𝑤=6∫10∫20[𝑢22+𝑢𝑤]𝑢=1𝑢=0𝑑𝑣𝑑𝑤=6∫10∫20(12+𝑤)𝑑𝑣𝑑𝑤=6∫10[𝑣2+𝑣𝑤]𝑣=2𝑣=0𝑑𝑤=6∫10(1+2𝑤)𝑑𝑤=6[𝑤+𝑤2]10=6(2)=12.

Jacobians and Transformed Regions in the Plane

  1. a. Solve the system
𝑢=𝑥−𝑦,𝑣=2𝑥+𝑦

for x and y in terms of u and v. Then find the value of the Jacobian 𝜕(𝑥,𝑦)/𝜕(𝑢,𝑣) .

b. Find the image under the transformation 𝑢 =𝑥 −𝑦 , 𝑣 =2𝑥 +𝑦 of the triangular region with vertices (0,0) , (1,1) , and (1, −2) in the 𝑥𝑦 -plane. Sketch the transformed region in the 𝑢𝑣 -plane.

  1. a. Solve the system
𝑢=𝑥+2𝑦,𝑣=𝑥−𝑦

for x and y in terms of u and v. Then find the value of the Jacobian 𝜕(𝑥,𝑦)/𝜕(𝑢,𝑣) .

b. Find the image under the transformation 𝑢 =𝑥 +2𝑦 , v = x - y of the triangular region in the xy-plane bounded by the lines y = 0, y = x, and 𝑥 +2𝑦 =2 . Sketch the transformed region in the uv-plane.

  1. a. Solve the system
𝑢=3𝑥+2𝑦,𝑣=𝑥+4𝑦

for x and y in terms of u and v. Then find the value of the Jacobian 𝜕(𝑥,𝑦)/𝜕(𝑢,𝑣) .

b. Find the image under the transformation 𝑢 =3𝑥 +2𝑦 , 𝑣 =𝑥 +4𝑦 of the triangular region in the xy-plane bounded by the x-axis, the y-axis, and the line 𝑥 +𝑦 =1 . Sketch the transformed region in the uv-plane.

  1. a. Solve the system
𝑢=2𝑥−3𝑦,𝑣=−𝑥+𝑦

for x and y in terms of u and v. Then find the value of the Jacobian 𝜕(𝑥,𝑦)/𝜕(𝑢,𝑣) .

b. Find the image under the transformation 𝑢 =2𝑥 −3𝑦 , 𝑣 = −𝑥 +𝑦 of the parallelogram 𝑅 in the 𝑥𝑦 -plane with boundaries 𝑥 = −3 , 𝑥 =0 , 𝑦 =𝑥 , and 𝑦 =𝑥 +1 . Sketch the transformed region in the 𝑢𝑣 -plane.

Substitutions in Double Integrals

  1. Evaluate the integral
∫40∫𝑥=(𝑦/2)+1𝑥=𝑦/22𝑥−𝑦2𝑑𝑥𝑑𝑦

from Example 1 directly by integration with respect to x and y to confirm that its value is 2.

  1. Use the transformation in Exercise 1 to evaluate the integral
∬𝑅(2𝑥2−𝑥𝑦−𝑦2)𝑑𝑥𝑑𝑦

for the region 𝑅 in the first quadrant bounded by the lines 𝑦 = −2𝑥 +4 , 𝑦 = −2𝑥 +7 , 𝑦 =𝑥 −2 , and 𝑦 =𝑥 +1 .

  1. Use the transformation in Exercise 3 to evaluate the integral
∬𝑅(3𝑥2+14𝑥𝑦+8𝑦2)𝑑𝑥𝑑𝑦

for the region 𝑅 bounded by the lines 𝑦 = −(3/2)𝑥 +1 , 𝑦 = −(3/2)𝑥 +3 , 𝑦 = −(1/4)𝑥 , and 𝑦 = −(1/4)𝑥 +1 .

  1. Use the transformation and parallelogram R in Exercise 4 to evaluate the integral
∬𝑅2(𝑥−𝑦)𝑑𝑥𝑑𝑦.
  1. Let R be the region in the first quadrant of the xy-plane bounded by the hyperbolas xy = 1, xy = 9 and the lines y = x, y = 4x. Use the transformation x = u/v, y = uv with u > 0 and v > 0 to rewrite
∬𝑅(√𝑦𝑥+√𝑥𝑦)𝑑𝑥𝑑𝑦

as an integral over an appropriate region G in the uv-plane. Then evaluate the uv-integral over G.

  1. a. Find the Jacobian of the transformation 𝑥 =𝑢 , 𝑦 =𝑢𝑣 and sketch the region 𝐺 :1 ≤𝑢 ≤2 , 1 ≤𝑢𝑣 ≤2 , in the 𝑢𝑣 -plane.

b. Then use Equation (2) to transform the integral

∫21∫21𝑦𝑥𝑑𝑦𝑑𝑥

into an integral over G, and evaluate both integrals.

  1. Polar moment of inertia of an elliptical plate A thin plate of constant density covers the region bounded by the ellipse 𝑥2/𝑎2 +𝑦2/𝑏2 =1 , a > 0, b > 0, in the xy-plane. Find the first moment of the plate about the origin. (Hint: Use the transformation 𝑥 =𝑎𝑟cos⁡𝜃 , 𝑦 =𝑏𝑟sin⁡𝜃 .)

  2. The area of an ellipse The area 𝜋𝑎𝑏 of the ellipse 𝑥2/𝑎2 +𝑦2/𝑏2 =1 can be found by integrating the function 𝑓(𝑥,𝑦) =1 over the region bounded by the ellipse in the xy-plane. Evaluating the integral directly requires a trigonometric substitution. An easier way to evaluate the integral is to use the transformation x = au, y = bv and evaluate the transformed integral over the disk 𝐺 :𝑢2 +𝑣2 ≤1 in the uv-plane. Find the area this way.

  3. Use the transformation in Exercise 2 to evaluate the integral

∫2/30∫2−2𝑦𝑦(𝑥+2𝑦)𝑒(𝑦−𝑥)𝑑𝑥𝑑𝑦

by first writing it as an integral over a region G in the uv-plane.

  1. Use the transformation 𝑥 =𝑢 +(1/2)𝑣 , 𝑦 =𝑣 to evaluate the integral
∫20∫(𝑦+4)/2𝑦/2𝑦3(2𝑥−𝑦)𝑒(2𝑥−𝑦)2𝑑𝑥𝑑𝑦

by first writing it as an integral over a region G in the uv-plane.

  1. Use the transformation 𝑥 =𝑢/𝑣 , 𝑦 =𝑢𝑣 to evaluate the integral sum
∫21∫𝑦1/𝑦(𝑥2+𝑦2)𝑑𝑥𝑑𝑦+∫42∫4/𝑦𝑦/4(𝑥2+𝑦2)𝑑𝑥𝑑𝑦.
  1. Use the transformation 𝑥 =𝑢2 −𝑣2 , y = 2uv to evaluate the integral
∫10∫2√1−𝑥0√𝑥2+𝑦2𝑑𝑦𝑑𝑥.

(Hint: Show that the image of the triangular region G with vertices (0,0) , (1,0) , (1,1) in the uv-plane is the region of integration R in the xy-plane defined by the limits of integration.)

Substitutions in Triple Integrals

  1. Evaluate the integral in Example 5 by integrating with respect to x, y, and z.

  2. Volume of a solid ellipsoid Find the volume of the solid ellipsoid

𝑥2𝑎2+𝑦2𝑏2+𝑧2𝑐2≤1.

(Hint: Let 𝑥 =𝑎𝑢 , 𝑦 =𝑏𝑣 , and 𝑧 =𝑐𝑤 . Then find the volume of an appropriate region in uvw-space.)

  1. Evaluate
∭𝐷|𝑥𝑦𝑧|𝑑𝑥𝑑𝑦𝑑𝑧

over the solid ellipsoid D,

𝑥2𝑎2+𝑦2𝑏2+𝑧2𝑐2≤1.

(Hint: Let 𝑥 =𝑎𝑢 , 𝑦 =𝑏𝑣 , and 𝑧 =𝑐𝑤 . Then integrate over an appropriate region in uvw-space.)

  1. Let D be the solid region in xyz-space defined by the inequalities
1≤𝑥≤2,0≤𝑥𝑦≤2,0≤𝑧≤1.

Evaluate

∭𝐷(𝑥2𝑦+3𝑥𝑦𝑧)𝑑𝑥𝑑𝑦𝑑𝑧

by applying the transformation

𝑢=𝑥,𝑣=𝑥𝑦,𝑤=3𝑧

and integrating over an appropriate region G in uvw-space.

Theory and Examples

  1. Find the Jacobian 𝜕(𝑥,𝑦)/𝜕(𝑢,𝑣) of the transformation

a. 𝑥 =𝑢cos⁡𝑣, 𝑦 =𝑢sin⁡𝑣

b. 𝑥 =𝑢sin⁡𝑣, 𝑦 =𝑢cos⁡𝑣.

  1. Find the Jacobian 𝜕(𝑥,𝑦,𝑧)/𝜕(𝑢,𝑣,𝑤) of the transformation

a. 𝑥 =𝑢cos⁡𝑣, 𝑦 =𝑢sin⁡𝑣, 𝑧 =𝑤

𝐛.𝑥=2𝑢−1,𝑦=3𝑣−4,𝑧=(1/2)(𝑤−4).
  1. How are double integrals used to calculate areas and average values. Give examples.

  2. How can you change a double integral in rectangular coordinates into a double integral in polar coordinates? Why might it be worthwhile to do so? Give an example.

  3. Define the triple integral of a function 𝑓(𝑥,𝑦,𝑧) over a bounded solid region in space.

  4. How are triple integrals in rectangular coordinates evaluated? How are the limits of integration determined? Give an example.

  5. Evaluate the appropriate determinant to show that the Jacobian of the transformation from Cartesian 𝜌𝜙𝜃 -space to Cartesian xyz-space is 𝜌2sin⁡𝜙 .

CHAPTER 14 Questions to Guide Your Review

  1. How are double integrals evaluated as iterated integrals? Does the order of integration matter? How are the limits of integration determined? Give examples.

  2. Define the double integral of a function of two variables over a bounded region in the coordinate plane.

  3. Substitutions in single integrals How can substitutions in single definite integrals be viewed as transformations of regions? What is the Jacobian in such a case? Illustrate with an example.

  4. Centroid of a solid semi-ellipsoid Assuming the result that the centroid of a solid hemisphere lies on the axis of symmetry three-eighths of the way from the base toward the top, show, by transforming the appropriate integrals, that the center of mass of a solid semi-ellipsoid (𝑥2/𝑎2) +(𝑦2/𝑏2) +(𝑧2/𝑐2) ≤1 , 𝑧 ≥0 , lies on the 𝑧 -axis three-eighths of the way from the base toward the top. (You can do this without evaluating any of the integrals.)

  5. Cylindrical shells In Section 6.2, we learned how to find the volume of a solid of revolution using the shell method. Specifically, if the region between the curve 𝑦 =𝑓(𝑥) and the 𝑥 -axis from 𝑎 to 𝑏 ( 0 <𝑎 <𝑏 ) is revolved about the 𝑦 -axis, the volume of the resulting solid is ∫𝑏𝑎2𝜋𝑥𝑓(𝑥)𝑑𝑥 . Prove that finding volumes by using triple integrals gives the same result. (Hint: Use cylindrical coordinates with the roles of 𝑦 and 𝑧 changed.)

  6. Inverse transform The equations 𝑥 =𝑔(𝑢,𝑣) , 𝑦 =ℎ(𝑢,𝑣) in Figure 14.57 transform the region G in the uv-plane into the region R in the xy-plane. Since the substitution transformation is one-to-one with continuous first partial derivatives, it has an inverse transformation, and there are equations 𝑢 =𝛼(𝑥,𝑦) , 𝑣 =𝛽(𝑥,𝑦) with continuous first partial derivatives transforming R back into G. Moreover, the Jacobian determinants of the transformations are related reciprocally by

𝜕(𝑥,𝑦)𝜕(𝑢,𝑣)=(𝜕(𝑢,𝑣)𝜕(𝑥,𝑦))−1.(10)

Equation (10) is proved in advanced calculus. Use it to find the area of the region R in the first quadrant of the xy-plane bounded by the lines y = 2x, 2y = x, and the curves xy = 2, 2xy = 1 for u = xy and v = y/x.

  1. (Continuation of Exercise 27.) For the region 𝑅 described in Exercise 27, evaluate the integral ∬𝑅𝑦2𝑑𝐴 .

  2. How are double and triple integrals in rectangular coordinates used to calculate volumes, average values, masses, moments, and centers of mass? Give examples.

  3. How are triple integrals defined in cylindrical and spherical coordinates? Why might one prefer working in one of these coordinate systems to working in rectangular coordinates?

  4. How are triple integrals in cylindrical and spherical coordinates evaluated? How are the limits of integration found? Give examples.

  5. How are substitutions in double integrals pictured as transformations of regions in the plane? Give a sample calculation.

  6. How are substitutions in triple integrals pictured as transformations of solid regions? Give a sample calculation.

CHAPTER 14 Practice Exercises

Evaluating Double Iterated Integrals

In Exercises 1–4, sketch the region of integration and evaluate the double integral.

  1. ∫101∫1/𝑦0𝑦𝑒𝑥𝑦𝑑𝑥𝑑𝑦

  2. ∫10∫𝑥30𝑒𝑦/𝑥𝑑𝑦𝑑𝑥

  3. ∫3/20∫√9−4𝑡2−√9−4𝑡2𝑡𝑑𝑠𝑑𝑡

  4. ∫10∫2−√𝑦√𝑦𝑥𝑦𝑑𝑥𝑑𝑦

In Exercises 5–8, sketch the region of integration and write an equivalent integral with the order of integration reversed. Then evaluate both integrals. 5. ∫40∫(𝑦−4)/2−√4−𝑦𝑑𝑥𝑑𝑦

  1. ∫10∫𝑥𝑥2√𝑥𝑑𝑦𝑑𝑥

  2. ∫3/20∫√9−4𝑦2−√9−4𝑦2𝑦𝑑𝑥𝑑𝑦

  3. ∫20∫4−𝑥202𝑥𝑑𝑦𝑑𝑥

Evaluate the integrals in Exercises 9–12.

  1. ∫10∫22𝑦4cos⁡(𝑥2)𝑑𝑥𝑑𝑦

  2. ∫20∫1𝑦/2𝑒𝑥2𝑑𝑥𝑑𝑦

  3. ∫80∫23√𝑥𝑑𝑦𝑑𝑥𝑦4+1

  4. ∫10∫13√𝑦2𝜋sin⁡𝜋𝑥2𝑥2𝑑𝑥𝑑𝑦

Areas and Volumes Using Double Integrals

  1. Area between line and parabola Find the area of the region enclosed by the line 𝑦 =2𝑥 +4 and the parabola 𝑦 =4 −𝑥2 in the xy-plane.

  2. Area bounded by lines and parabola Find the area of the “triangular” region in the xy-plane that is bounded on the right by the parabola 𝑦 =𝑥2 , on the left by the line 𝑥 +𝑦 =2 , and above by the line y = 4.

  3. Volume of the region under a paraboloid Find the volume under the paraboloid 𝑧 =𝑥2 +𝑦2 above the triangle enclosed by the lines 𝑦 =𝑥 , 𝑥 =0 , and 𝑥 +𝑦 =2 in the xy-plane.

  4. Volume of the region under a parabolic cylinder Find the volume under the parabolic cylinder 𝑧 =𝑥2 above the region enclosed by the parabola 𝑦 =6 −𝑥2 and the line y = x in the xy-plane.

Average Values

Find the average value of 𝑓(𝑥,𝑦) =𝑥𝑦 over the regions in Exercises 17 and 18.

  1. The square bounded by the lines 𝑥 =1,𝑦 =1 in the first quadrant

  2. The quarter circle 𝑥2 +𝑦2 ≤1 in the first quadrant

Polar Coordinates

Evaluate the integrals in Exercises 19 and 20 by changing to polar coordinates.

  1. ∫1−1∫√1−𝑥2−√1−𝑥22𝑑𝑦𝑑𝑥(1+𝑥2+𝑦2)2

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

  3. Integrating over a lemniscate Integrate the function 𝑓(𝑥,𝑦) =1/(1 +𝑥2 +𝑦2)2 over the region enclosed by one loop of the lemniscate (𝑥2 +𝑦2)2 −(𝑥2 −𝑦2) =0 .

  4. Integrate 𝑓(𝑥,𝑦) =1/(1 +𝑥2 +𝑦2)2 over

a. Triangular region The triangle with vertices (0,0),(1,0) , and (1,√3) .

b. First quadrant The first quadrant of the xy-plane.

Evaluating Triple Iterated Integrals Evaluate the integrals in Exercises 23–26. 23. ∫𝜋0∫𝜋0∫𝜋0cos⁡(𝑥 +𝑦 +𝑧)𝑑𝑥𝑑𝑦𝑑𝑧

  1. ∫ln⁡7ln⁡6∫ln⁡20∫ln⁡5ln⁡4𝑒(𝑥+𝑦+𝑧)𝑑𝑧𝑑𝑦𝑑𝑥

  2. ∫10∫𝑥20∫𝑥+𝑦0(2𝑥 −𝑦 −𝑧)𝑑𝑧𝑑𝑦𝑑𝑥

  3. ∫𝑒1∫𝑥1∫𝑧02𝑦𝑧3𝑑𝑦𝑑𝑧𝑑𝑥

Volumes and Average Values Using Triple Integrals

  1. Volume Find the volume of the wedge-shaped solid region enclosed on the side by the cylinder 𝑥 = −cos⁡𝑦 , −𝜋/2 ≤𝑦 ≤𝜋/2 , on the top by the plane 𝑧 = −2𝑥 , and below by the 𝑥𝑦 -plane.

教材插图

  1. Volume Find the volume of the solid that is bounded above by the cylinder 𝑧 =4 −𝑥2 , on the sides by the cylinder 𝑥2 +𝑦2 =4 , and below by the xy-plane.

教材插图

  1. Average value Find the average value of 𝑓(𝑥,𝑦,𝑧) =30𝑥𝑧√𝑥2+𝑦 over the rectangular solid in the first octant bounded by the coordinate planes and the planes x=1, y=3, z=1.

  2. Average value Find the average value of 𝜌 over the ball 𝜌 ≤𝑎 (spherical coordinates).

Cylindrical and Spherical Coordinates

  1. Cylindrical to rectangular coordinates Convert
∫2𝜋0∫√20∫√4−𝑟2𝑟3𝑟𝑑𝑧𝑑𝑟𝑑𝜃,𝑟≥0

to (a) rectangular coordinates with the order of integration dz dx dy and (b) spherical coordinates. Then (c) evaluate one of the integrals.

  1. Rectangular to cylindrical coordinates (a) Convert to cylindrical coordinates. Then (b) evaluate the new integral.
∫10∫√1−𝑥2−√1−𝑥2∫(𝑥2+𝑦2)−(𝑥2+𝑦2)21𝑥𝑦2𝑑𝑧𝑑𝑦𝑑𝑥
  1. Rectangular to spherical coordinates (a) Convert to spherical coordinates. Then (b) evaluate the new integral.
∫1−1∫√1−𝑥2−√1−𝑥2∫1√𝑥2+𝑦2𝑑𝑧𝑑𝑦𝑑𝑥
  1. Rectangular, cylindrical, and spherical coordinates Write an iterated triple integral for the integral of 𝑓(𝑥,𝑦,𝑧) =6 +4𝑦 over the region in the first octant bounded by the cone 𝑧 =√𝑥2+𝑦2 , the cylinder 𝑥2 +𝑦2 =1 , and the coordinate planes in (a) rectangular coordinates, (b) cylindrical coordinates, and (c) spherical coordinates. Then (d) find the integral of f by evaluating one of the triple integrals.

  2. Cylindrical to rectangular coordinates Set up an integral in rectangular coordinates equivalent to the integral

∫𝜋/20∫√31∫√4−𝑟21𝑟3(sin⁡𝜃cos⁡𝜃)𝑧2𝑑𝑧𝑑𝑟𝑑𝜃.

Arrange the order of integration to be z first, then y, then x.

  1. Rectangular to cylindrical coordinates The volume of a solid is
∫20∫√2𝑥−𝑥20∫√4−𝑥2−𝑦2−√4−𝑥2−𝑦2𝑑𝑧𝑑𝑦𝑑𝑥.

a. Describe the solid by giving equations for the surfaces that form its boundary.

b. Convert the integral to cylindrical coordinates, but do not evaluate the integral.

  1. Spherical versus cylindrical coordinates Triple integrals involving spherical shapes do not always require spherical coordinates for convenient evaluation. Some calculations may be accomplished more easily with cylindrical coordinates. As a case in point, find the volume of the solid region bounded above by the sphere 𝑥2 +𝑦2 +𝑧2 =8 and below by the plane z = 2 by using (a) cylindrical coordinates and (b) spherical coordinates.

Masses and Moments

  1. Finding 𝐼𝑧 in spherical coordinates Find the moment of inertia about the z-axis of a solid of constant density 𝛿 =1 that is bounded above by the sphere 𝜌 =2 and below by the cone 𝜙 =𝜋/3 (spherical coordinates).

  2. Moment of inertia of a “thick” sphere Find the moment of inertia of a solid of constant density 𝛿 bounded by two concentric spheres of radii a and 𝑏(𝑎 <𝑏) about a diameter.

  3. Moment of inertia of an apple Find the moment of inertia about the z-axis of a solid of density 𝛿 =1 enclosed by the spherical coordinate surface 𝜌 =1 −cos⁡𝜙 . The solid is the red curve rotated about the z-axis in the accompanying figure.

教材插图

  1. Centroid Find the centroid of the “triangular” region bounded by the lines x = 2, y = 2 and the hyperbola xy = 2 in the xy-plane.

  2. Centroid Find the centroid of the region between the parabola 𝑥 +𝑦2 −2𝑦 =0 and the line 𝑥 +2𝑦 =0 in the xy-plane.

  3. Polar moment Find the polar moment of inertia about the origin of a thin triangular plate of constant density 𝛿 =3 bounded by the y-axis and the lines y = 2x and y = 4 in the xy-plane.

  4. Polar moment Find the polar moment of inertia about the center of a thin rectangular sheet of constant density 𝛿 =1 bounded by the lines

a. 𝑥 = ±2, 𝑦 = ±1 in the xy-plane

b. 𝑥 = ±𝑎, 𝑦 = ±𝑏 in the xy-plane.

(Hint: Find 𝐼𝑥 . Then use the formula for 𝐼𝑥 to find 𝐼𝑦 , and add the two to find 𝐼0 .)

  1. Inertial moment Find the moment of inertia about the x-axis of a thin plate of constant density 𝛿 covering the triangle with vertices (0,0) , (3,0) , and (3,2) in the xy-plane.

  2. Plate with variable density Find the center of mass and the moments of inertia about the coordinate axes of a thin plate bounded by the line y = x and the parabola 𝑦 =𝑥2 in the xy-plane if the density is 𝛿(𝑥,𝑦) =𝑥 +1 .

  3. Plate with variable density Find the mass and first moments about the coordinate axes of a thin square plate bounded by the lines 𝑥 = ±1 , 𝑦 = ±1 in the 𝑥𝑦 -plane if the density is 𝛿(𝑥,𝑦) =𝑥2 +𝑦2 +1/3 .

  4. Triangles with same inertial moment Find the moment of inertia about the x-axis of a thin triangular plate of constant density 𝛿 whose base lies along the interval [0,𝑏] on the x-axis and whose vertex lies on the line y = h above the x-axis. As you will see, it does not matter where on the line this vertex lies. All such triangles have the same moment of inertia about the x-axis.

  5. Centroid Find the centroid of the region in the polar coordinate plane defined by the inequalities 0 ≤𝑟 ≤3 , −𝜋/3 ≤𝜃 ≤𝜋/3 .

  6. Centroid Find the centroid of the region in the first quadrant bounded by the rays 𝜃 =0 and 𝜃 =𝜋/2 and the circles r = 1 and r = 3.

  7. a. Centroid Find the centroid of the region in the polar coordinate plane that lies inside the cardioid 𝑟 =1 +cos⁡𝜃 and outside the circle r = 1.

b. Sketch the region and show the centroid in your sketch.

  1. a. Centroid Find the centroid of the plane region defined by the polar coordinate inequalities 0 ≤𝑟 ≤𝑎, −𝛼 ≤𝜃 ≤𝛼 ( 0 <𝛼 ≤𝜋 ). How does the centroid move as 𝛼 →𝜋− ?

b. Sketch the region for 𝛼 =5𝜋/6 and show the centroid in your sketch.

Substitutions

  1. Show that if u = x - y and v = y, then for any continuous f,
∫∞0∫𝑥0𝑒−𝑠𝑥𝑓(𝑥−𝑦,𝑦)𝑑𝑦𝑑𝑥=∫∞0∫∞0𝑒−𝑠(𝑢+𝑣)𝑓(𝑢,𝑣)𝑑𝑢𝑑𝑣.
  1. What relationship must hold between the constants 𝑎,𝑏 , and 𝑐 to make
∫∞−∞∫∞−∞𝑒−(𝑎𝑥2+2𝑏𝑥𝑦+𝑐𝑦2)𝑑𝑥𝑑𝑦=1?

(Hint: Let 𝑠 =𝛼𝑥 +𝛽𝑦 and 𝑡 =𝛾𝑥 +𝛿𝑦 , where (𝛼𝛿 −𝛽𝛾)2 =𝑎𝑐 −𝑏2 . Then 𝑎𝑥2 +2𝑏𝑥𝑦 +𝑐𝑦2 =𝑠2 +𝑡2 .)

CHAPTER 14 Additional and Advanced Exercises

Volumes

  1. Sand pile: double and triple integrals The base of a sand pile covers the region in the xy-plane that is bounded by the parabola 𝑥2 +𝑦 =6 and the line y = x. The height of the sand above the point (𝑥,𝑦) is 𝑥2 . Express the volume of sand as (a) a double integral and (b) a triple integral. Then (c) find the volume.

  2. Water in a hemispherical bowl A hemispherical bowl of radius 5 cm is filled with water to within 3 cm of the top. Find the volume of water in the bowl.

  3. Solid cylindrical region between two planes Find the volume of the portion of the solid cylinder 𝑥2 +𝑦2 ≤1 that lies between the planes z = 0 and 𝑥 +𝑦 +𝑧 =2 .

  4. Sphere and paraboloid Find the volume of the solid region bounded above by the sphere 𝑥2 +𝑦2 +𝑧2 =2 and below by the paraboloid 𝑧 =𝑥2 +𝑦2 .

  5. Two paraboloids Find the volume of the solid region bounded above by the paraboloid 𝑧 =3 −𝑥2 −𝑦2 and below by the paraboloid 𝑧 =2𝑥2 +2𝑦2 .

  6. Spherical coordinates Find the volume of the solid region enclosed by the spherical coordinate surface 𝜌 =2sin⁡𝜙 (see accompanying figure).

教材插图

  1. Hole in solid ball A circular cylindrical hole is bored through a ball, the axis of the hole being a diameter of the sphere. The volume of the remaining solid is
𝑉=2∫2𝜋0∫√30∫√4−𝑧21𝑟𝑑𝑟𝑑𝑧𝑑𝜃.

a. Find the radius of the hole and the radius of the sphere.

b. Evaluate the integral.

  1. Ball and cylinder Find the volume of material cut from the ball 𝑟2 +𝑧2 ≤9 by the cylinder 𝑟 =3sin⁡𝜃 .

  2. Two paraboloids Find the volume of the solid region enclosed by the surfaces 𝑧 =𝑥2 +𝑦2 and 𝑧 =(𝑥2 +𝑦2 +1)/2 .

  3. Cylinder and surface z = xy Find the volume of the solid region in the first octant that lies between the cylinders r = 1 and r = 2 and is bounded below by the xy-plane and above by the surface z = xy.

Changing the Order of Integration

  1. Evaluate the integral
∫∞0𝑒−𝑎𝑥−𝑒−𝑏𝑥𝑥𝑑𝑥.

(Hint: Use the relation

𝑒−𝑎𝑥−𝑒−𝑏𝑥𝑥=∫𝑏𝑎𝑒−𝑥𝑦𝑑𝑦

to form a double integral, and evaluate the integral by changing the order of integration.)

  1. a. Polar coordinates Show, by changing to polar coordinates, that
∫𝑎sin⁡𝛽0∫√𝑎2−𝑦2𝑦cot⁡𝛽ln⁡(𝑥2+𝑦2)𝑑𝑥𝑑𝑦=𝑎2𝛽(ln⁡𝑎−12),

where 𝑎 >0 and 0 <𝛽 <𝜋/2 .

b. Rewrite the Cartesian integral with the order of integration reversed.

  1. Reducing a double to a single integral By changing the order of integration, show that the following double integral can be reduced to a single integral:
∫𝑥0∫𝑢0𝑒𝑚(𝑥−𝑡)𝑓(𝑡)𝑑𝑡𝑑𝑢=∫𝑥0(𝑥−𝑡)𝑒𝑚(𝑥−𝑡)𝑓(𝑡)𝑑𝑡.

Similarly, it can be shown that

∫𝑥0∫𝑣0∫𝑢0𝑒𝑚(𝑥−𝑡)𝑓(𝑡)𝑑𝑡𝑑𝑢𝑑𝑣=∫𝑥0(𝑥−𝑡)22𝑒𝑚(𝑥−𝑡)𝑓(𝑡)𝑑𝑡.
  1. Transforming a double integral to obtain constant limits Sometimes a multiple integral with variable limits can be changed into one with constant limits. By changing the order of integration, show that
∫10𝑓(𝑥)(∫𝑥0𝑔(𝑥−𝑦)𝑓(𝑦)𝑑𝑦)𝑑𝑥=∫10𝑓(𝑦)(∫1𝑦𝑔(𝑥−𝑦)𝑓(𝑥)𝑑𝑥)𝑑𝑦=12∫10∫10𝑔(|𝑥−𝑦|)𝑓(𝑥)𝑓(𝑦)𝑑𝑥𝑑𝑦.

Masses and Moments

  1. Minimizing polar inertia A thin plate of constant density is to occupy the triangular region in the first quadrant of the xy-plane having vertices (0,0) , (𝑎,0) , and (𝑎,1/𝑎) . What value of a will minimize the plate’s polar moment of inertia about the origin?

  2. Polar inertia of triangular plate Find the polar moment of inertia about the origin of a thin triangular plate of constant density 𝛿 =3 bounded by the y-axis and the lines y = 2x and y = 4 in the xy-plane.

  3. Mass and polar inertia of a counterweight The counterweight of a flywheel of constant density 1 has the form of the smaller segment cut from a circle of radius 𝑎 by a chord at a distance 𝑏 from the center (𝑏 <𝑎) . Find the mass of the counterweight and its polar moment of inertia about the center of the wheel.

  4. Centroid of a boomerang Find the centroid of the boomerang-shaped region between the parabolas 𝑦2 = −4(𝑥 −1) and 𝑦2 = −2(𝑥 −2) in the xy-plane.

Theory and Examples

  1. Evaluate
∫𝑎0∫𝑏0𝑒max(𝑏2𝑥2,𝑎2𝑦2)𝑑𝑦𝑑𝑥,

where 𝑎 and 𝑏 are positive numbers and

max(𝑏2𝑥2,𝑎2𝑦2)={𝑏2𝑥2 if 𝑏2𝑥2≥𝑎2𝑦2𝑎2𝑦2 if 𝑏2𝑥2<𝑎2𝑦2.
  1. Show that
∬𝜕2𝐹(𝑥,𝑦)𝜕𝑥𝜕𝑦𝑑𝑥𝑑𝑦

over the rectangle 𝑥0 ≤𝑥 ≤𝑥1,𝑦0 ≤𝑦 ≤𝑦1 is

𝐹(𝑥1,𝑦1)−𝐹(𝑥0,𝑦1)−𝐹(𝑥1,𝑦0)+𝐹(𝑥0,𝑦0).
  1. Suppose that 𝑓(𝑥,𝑦) can be written as a product 𝑓(𝑥,𝑦) =𝐹(𝑥)𝐺(𝑦) of a function of x and a function of y. Then the integral of f over the rectangle 𝑅 :𝑎 ≤𝑥 ≤𝑏,𝑐 ≤𝑦 ≤𝑑 can be evaluated as a product as well, by the formula
∬𝑅𝑓(𝑥,𝑦)𝑑𝐴=(∫𝑏𝑎𝐹(𝑥)𝑑𝑥)(∫𝑑𝑐𝐺(𝑦)𝑑𝑦).(1)

The argument is that

∬𝑅𝑓(𝑥,𝑦)𝑑𝐴=∫𝑑𝑐(∫𝑏𝑎𝐹(𝑥)𝐺(𝑦)𝑑𝑥)𝑑𝑦(i) =∫𝑑𝑐(𝐺(𝑦)∫𝑏𝑎𝐹(𝑥)𝑑𝑥)𝑑𝑦(ii) =∫𝑑𝑐(∫𝑏𝑎𝐹(𝑥)𝑑𝑥)𝐺(𝑦)𝑑𝑦(iii) =(∫𝑏𝑎𝐹(𝑥)𝑑𝑥)∫𝑑𝑐𝐺(𝑦)𝑑𝑦.(iv)

a. Give reasons for Steps (i) through (iv).

When it applies, Equation (1) can be a time-saver. Use it to evaluate the following integrals.

𝐛.∫ln⁡20∫𝜋/20𝑒𝑥cos⁡𝑦𝑑𝑦𝑑𝑥𝐜.∫21∫1−1𝑥𝑦2𝑑𝑥𝑑𝑦
  1. Let 𝐷𝑢𝑓 denote the derivative of 𝑓(𝑥,𝑦) =(𝑥2 +𝑦2)/2 in the direction of the unit vector 𝑢 =𝑢1𝑖 +𝑢2𝑗 .

a. Finding average value Find the average value of 𝐷𝐮𝑓 over the triangular region cut from the first quadrant by the line 𝑥 +𝑦 =1 .

b. Average value and centroid Show in general that the average value of 𝐷𝑢𝑓 over a region in the xy-plane is the value of 𝐷𝑢𝑓 at the centroid of the region.

  1. The value of Γ(1/2) The gamma function,
Γ(𝑥)=∫∞0𝑡𝑥−1𝑒−𝑡𝑑𝑡,

extends the factorial function from the nonnegative integers to other real values. Of particular interest in the theory of differential equations is the number

Γ(12)=∫∞0𝑡(1/2)−1𝑒−𝑡𝑑𝑡=∫∞0𝑒−𝑡√𝑡𝑑𝑡.(2)

a. If you have not yet done Exercise 41 in Section 14.4, do it now to show that

𝐼=∫∞0𝑒−𝑦2𝑑𝑦=√𝜋2.

b. Substitute 𝑦 =√𝑡 in Equation (2) to show that Γ(1/2) =2𝐼 =√𝜋 .

  1. Total electrical charge over circular plate The electrical charge distribution on a circular plate of radius R meters is 𝜎(𝑟,𝜃) =𝑘𝑟(1 −sin⁡𝜃) coulomb/m 2 (k a constant). Integrate 𝜎 over the plate to find the total charge Q.

  2. A parabolic rain gauge A bowl is in the shape of the graph of 𝑧 =𝑥2 +𝑦2 from z = 0 to z = 30 cm. You plan to calibrate the bowl to make it into a rain gauge. What height in the bowl would correspond to 3 cm of rain? 9 cm of rain?

  3. Water in a satellite dish A parabolic satellite dish is 2 m wide and 1/2 m deep. Its axis of symmetry is tilted 30 degrees from the vertical.

a. Set up, but do not evaluate, a triple integral in rectangular coordinates that gives the amount of water the satellite dish will hold. (Hint: Put your coordinate system so that the satellite dish is in “standard position” and the plane of the water level is slanted.) (Caution: The limits of integration are not “nice.”)

b. What would be the smallest tilt of the satellite dish so that it holds no water?

  1. An infinite half-cylinder Let D be the interior of the infinite right circular half-cylinder of radius 1 with its single-end face suspended 1 unit above the origin and its axis the ray from (0,0,1) to ∞ . Use cylindrical coordinates to evaluate
∭𝐷𝑧(𝑟2+𝑧2)−5/2𝑑𝑉.
  1. Hypervolume We have learned that ∫𝑏𝑎1𝑑𝑥 is the length of the interval [𝑎,𝑏] on the number line (one-dimensional space), ∬𝑅1𝑑𝐴 is the area of region R in the xy-plane (two-dimensional space), and ∭𝐷1𝑑𝑉 is the volume of the region D in three-dimensional space (xyz-space). We could continue: If Q is a region in 4-space (xyzw-space), then ∭𝑄1𝑑𝑉 is the “hyper-volume” of Q. Use your generalizing abilities and a Cartesian coordinate system of 4-space to find the hypervolume inside the unit four-dimensional sphere 𝑥2 +𝑦2 +𝑧2 +𝑤2 =1 .

CHAPTER 14 Technology Application Projects

Mathematica/Maple Projects

Projects can be found within MyLab Math.

  • Take Your Chances: Try the Monte Carlo Technique for Numerical Integration in Three Dimensions Use the Monte Carlo technique to integrate numerically in three dimensions.

  • Means and Moments and Exploring New Plotting Techniques, Part II Use the method of moments in a form that makes use of geometric symmetry as well as multiple integration.