I I

MULTIPLE INTEGRALS
I
A double integral of a positive function is a volume, which
is the limit of sums of volumes of rectangular columns.
In this chapter we extend the idea of a definite integral to double and triple integrals
of functions of two or three variables. These ideas are then used to compute volumes,
masses, and centroids of more general regions than we were able to consider in
Chapters 6 and 8. We also use double integrals to calculate probabilities when two
random variables are involved.
We will see that polar coordinates are useful in computing double integrals over some
types of regions. In a similar way, we will introduce two new coordinate systems in
three-dimensional space-cylindrical coordinates and spherical coordinates-that
greatly simplify the computation of triple integrals over certain commonly occurring
solid regions.
15.1 ' DOUBLE INTEGRALS OVER RECTANGLES
In much the same way that our attempt to solve the area problem led to the definition of a
definite integral, we now seek to find the volume of a solid and in the process we arrive at
the definition of a double integral.
REVIEW OF THE DEFINITE INTEGRAL
First let's recall the basic facts concerning definite integrals of functions of a single vari-
able. If f ( x ) is defined for a C x c b, we start by dividing the interval [a, b] inrn n sub-
intervals [x,-1, x, ] of equal width Ax = ( b - a)/n and we choose sample points .r" in these
subintervals. Then we form the Riemann sum
and take the limit of such sums as n + to obtain the definite integral of f from a to b:
In the special case where f ( x) 2 0, the Riemann sum can be interpreted as the sum of the
areas of the approximating rectangles in Figure 1, and Jf f ( x) dx represents the arsa under
the curve y = f ( x) from a to b.
FI GURE I
1
VOLUMES AND DOUBLE INTEGRALS
In a similar manner we consider a function f of two variables defined on a closed rectangle
R = [a, b] X [c, d ] = {(x, y) E R2 ( a c x c b, c
y d )
and we first suppose that f (x, y) 2 0. The graph off is a surface with equation i = - f ( x, y).
Let S be the solid that lies above R and under the graph of f , that is,
S = {(x, y, z ) E R3 10 c z G f ( x, y), (x, y) E R)
FI GURE 2
(See Figure 2.) Our goal is to find the volume of S.
The first step is to divide the rectangle R into subrectangles. We accomplish this by
dividing the interval [a, b] into m subintervals [xi-,, xi] of equal width A x = ( b - a) / m
and dividing [c, d ] into n subintervals [yj-1, yj] of equal width Ay = (d - c)/n. By draw-
952 11[1 CHAPTER 15 MULTIPLE INTEGRALS
FI GURE 3
Dividing R into subrectangles
ing lines parallel to the coordinate axes through the endpoints of these subintervals, as in
Figure 3, we hrm the subrectangles
Rv = [X,-1, X, ] X [)',-I, )',I = {(XI y) I Xl - I -x Xt, YJ-l y YJ}
each with area A A = A x Ay .
If we choose a sample point (xcT, y,;) in each R,,, then we can approximate the part of
S that lies above each R,, by a thin rectangular box (or '-column") with base R,, and height
f (x:, y,;) as shown in Figure 4. (Compare with Figure 1 .) The volume of this box is the
height of the box times the area of the base rectangle:
If we follow this procedure for all the rectangles and add the volumes of the corresponding
boxes, we get an approximation to the total volume of S:
(See Figure 5.) This double sum means that for each subrectangle we evaluate f at the cho-
sen point and multiply by the area of the subrectangle. and then we add the results.
FI GURE 4 FI GURE 5
SECT I ON 15.1 DOUBLE INTEGRALS OVER RECTANGLES 1111 953
1
Our intuition tells us that the approximat~on given In (3) becomes better as m and n
become larger and so we would expect that
B The meanlng of the double lhrn~t In Equatlon 4
1s t hat we can make t he double sum as close as
,Ti n
we like t o t he number V [for any choice of I3
V = lim 2 2 f (xi:, yiT) AA
m, n- - r x , =I , =I
(r,:. 1,:) In R, I by taklng m and n sufficiently
large
We use the expression in Equation 4 to define the volume of the solid S that lies under the
graph of f and above the rectangle R. (It can be shown that this definition is consistent with
our formula for volume in Section 6.2.)
Limits of the type that appear in Equation 4 occur frequently, not just in finding vol-
- >-even umes but in a variety of other situations as well-as we will see in Section
when f is not a positive function. So we make the following definition.
DEFINITION The double integral of f over the rectangle R is
~
The precise meaning of the limit in Definition 5 is that for every number E > 0 there is
E Although we have deflned t he double ~nt egr al
an integer N such that
t Notice t he similarity bet ween Def i ni t on 5
and t he def i r ~i t on of a single integral In
Equat on 2
by dividing R Into equal-sized subrectanges, w e
could have used subrectangles R,, of unequal
size But then w e woul d have t o ensure t hat all
of their dimensions approach 0 in t he l i mi t i ng
m n
[[fix, y) d~ = lim 2 2 f ( x z ~ , l ; : ) AA
* ,,I. ,,+I , = I
R
if this limit exists.
process
for all integers rn and n greater than N and for any choice of sample points (x;. T.") in R,,.
A function f is called integrable if the limit in Definition 5 exists. It is shown in courses
on advanced calculus that all continuous functions are integrable. In fact, the double inte-
gral of f exists provided that f is "not too discontinuous." In particular, if f is bounded
[that is, there is a constant M such that f(x, y) 6 M for all (x, y ) in R], and f is con-
tinuous there, except on a finite number of smooth curves, then f is integrable over R.
The sample point (xi:, y,;) can be chosen to be any point in the subrectangle R. . but if
we choose it to be the upper right-hand corner of R,j [namely (xi, yj), see Figure 31, then
the expression for the double integral looks simpler:
By comparing Definitions 4 and 5, we see that a volume can be written as a double
integral:
1
I f f (x, y ) 3 0. then the volume V of the solid that lies above the rectangle R and
below the surface z = f (x, y) is
9 5 4 ] I ( ( CHAPTER 15 MULTIPLE INTEGRALS
The sum in ~Lfinition 5,
P f ( x?. y.T) AA
, =I ,=I
is called a double Riemann sum and is used as an approximation to the value of the
double integral. [Notice how similar it is to the Riemann sum in (1) for a function of a
single variable.] If f happens to be a positzve function. then the double Riemann sum
represents the sum of volumes of columns, as in Figure 5. and is an approximation to the
volume under the graph of f and above the rectangle R.
m12.2i
PI EXAMPLE I Estimate the volume of the solid that lies above the square
R = [O, 21 X [O, 21 and below the elliptic parabolo~d I = 16 - x2 - 2y2 Divide R
into four equal squares and choose the sample point to be the upper right corner of
each square R,, Sketch the solid and the approximating rectangular boxes.
12, li
(1,11 SOLU1 I ON The squares are shown In Figure 6. The paraboloid is the graph of
R, , &I f (x, y) = 16 - x' - 2?' and the area of each square is 1. Approximating the volume
by the Riemann sum w~t h m = n = 2, we have
O I
1 2 k
2 2
FI GURE 6
FI GURE 7
FI GURE 8
The Rirmann sum approximations to
the volume under s = 16 - x' - 2y'
becwme more accurate as rn and
= 13(1) + 7(1) + lO(1) + 4(1) = 31
This is the volume of the approximating rectangular boxes shown in Figure 7.
We get better approximations to the volume in Example 1 if we increase the number of
squares. Figure 8 shows how the columns start to look more like the actual solid and the
corresponding approximations become more accurate when we use 16, 64, and 256
squares. In the next section we will be able to show that the exact volume is 48.
n increase. (a) m = n = 4, V = 41.5 ( b) m=n=S, V=44 875 (c) m = n = 16, V= 46.46875
PI EXAMPLE 2 If R = {( x, y) 1 - 1 S x < 1, -2
y
2) , evaluate the integral
Jj Jm dA
R
SECTI ON 15.1 DOUBLE INTEGRALS OVER RECTANGLES 1111 955
B
SOLUTION It would be very difficult to evaluate this integral directly from Definition 5 but,
because J i T 2 0, we can compute the integral by interpreting it as a volume. If
z - ,/-, then x2 + z2 = 1 and z 3 0, so the given double integral represents the
volume of the solid S that lies below the circular cylinder x2 + z2 = 1 and abol-e the
rectangle R. (See Figure 9.) The volume of S is the area of a semicircle with radius 1
times the length of the cylinder. Thus
/ [ J F F ~ A = ; n - ( l ) > x 4 =2n-
" 6
FIGURE 9 R
THE MIDPOINT RULE
FIGURE 10
The methods that we used for approximating single integrals (the Midpoint Rule. the
Trapezoidal Rule, Simpson's Rule) all have counterparts for double integrals. Here we
consider only the Midpoint Rule for double integrals. This means that we use a double Rie-
mann sum to approximate the double integral, where the sample point ( x: , j g) in R, is
-
chosen to be the center ( x, , yj) of R,,. In other words, T, is the midpoint of [ x , - ~ . .-i- ] and 4;
is the midpoint of [ y j - ] , y, ].
1 MIDPOINT RULE FOR DOUBLE INTEGRALS
1
where X, is the midpoint of [x.], x , ] and \, is the midpoint of [ y j - I . v,].
EXAMPLE 3 Use the Midpoint Rule with m = n = 2 to estimate the value of the
integral .(^.fR ( x - 3 y 2 ) dA, where R = { ( x, y ) 1 0 < x < 2, 1 6 y < 2).
SOLUTION In using the Midpoint Rule with m = n = 2, we evaluate f ( x , ?.) = x - 3!%t
the centers of the four subrectangles shown in Figure 10. So xl = $, X2 = ;, 7, = :. and
-
7
~2 = 4. The area of each subrectangle is AA = i. Thus
Thus we have l[ ( x - 3 y 2 ) dA = - 11.875
R
In the next section we will develop an efficient method f or com-
puting double integrals and then we will see that the exact value of the double integral in
Example 3 is - 12. (Remember that the interpretation of a double integral as a volume is
valid only when the integrand f is a positive function. The integrand in Example 3 is not
a positive function, so its integral is not a volume. In Examples 2 and 3 in Section 15.2 we
will discuss how to interpret integrals of functions that are not always positive in terms of
volumes.) If we keep dividing each subrectangle in Figure 10 into four smaller ones with
956 l l j l CHAPTER 15 MULTIPLE INTEGRALS
s~milar shape, we get the Midpoint Rule approximations displayed in the chart In the mar-
M~dpo~nt Rule
>uh~c.,-ln_rlej approxlrndtlon\
gin Notice how these approximdtlons dpproach the exact value ot the double Integral, - 12
AVERAGE VALUE
- 1 1 8750 ,
1 - I 1 9687 Recall from Section 6 5 that the average value of a funct~on f of one vanable defined on
L 1
- I I 9C)?7- an interval [a. b] i\
L
I
- -
--n I I ')980
i
1 ' h
l i 11 f * \ e =- \ f ( l ) d l B
b - u - u
In a \imilar fa\hion we define the average value of a functlon f of two variable\ defined
on a rectangle R to be
-
1 -,
,Live = - lj f(.x, .y) d'4
A( R) ;p
where A( R) is the area of R.
If f(x, y ) 3 0, the equation
*
A(R) X /,.. = 1 1 f (.r. i-i dA
".
R
says that the box with base R and height f,,, has the same volunle as the solid that lies
under the graph off: [If: = f(x. y ) describes a mountainous region and you chop off the
tops of the mountains at height fiLfe, then you can use them to fill in the valleys so that the
FI GURE I I region becomes completely flat. See Figure 11.1
EXAMPLE 4 The contour map in Figure 12 shows the snowfall. in inches, that fell on the
state of Colorado on December 20 and 21. 2006. (The state is in the shape of a rectangle
that measures 388 mi west to east and 276 mi south to north.) Use the contour map to
estimate the average snowfall for the entire state of Colorado on those days.
FI GURE I 2
FI GURE
SECTION 15.1 DOUBLE INTEGRALS OVER RECTANGLES 1111 957
SOLUTI ON Let's place the origin at the southwest comer of the state. Then 0 S x S 388,
0 S y S 276, and f ( x, y) is the snowfall, in inches, at a location x miles to the east and
y miles to the north of the origin. If R is the rectangle that represents Colorado. he n the
average snowfall for the state on December 20-21 was
1
hVe = - [[ f ( x, y) dA
A(R) ;;
where A(R) = 388 276. To estimate the value of this double integral, let's use the Mid-
point Rule with rn = n = 4. In other words, we divide R into 16 subrectangles of equal
size, as in Figure 13. The area of each subrectangle is
Using the contour map to estimate the value of f at the center of each subrectangle.
we get
Therefore
On December 20-2 1, 2006, Colorado received an average of approximately 13 inches of
snow.
7
958 1 ) ( 1 CHAPTER 15 MULTIPLE I NTEGRALS
a Double inr2;rals behave t h~s way because
the double s i r s mat define them behave
t h ~ s way.
PROPERTIES QF DOUBLE INTEGRALS
1
We list here three properties of doubl e integrals t hat can be proved in t he same manner as
i n Sect i on 5.2. We assume that all of the integrals exist. Properties 7 and 8 are referred t o
as t he linearity of t he integral.
c f (x. y) dA = c (j f (x, y) dA
whe r e c is a const ant
w -
R R
If f (x, y) 3 g(x, y) f or all (x, y ) i n R, then
, -
II.I (a) Estimate the volume of the solid that lies below (b) Estimate the double integral with m = n = 4 by choosing
the surface 3 = xy and above the rectangle the sample points to be the points farthest from the origin.
Use a Riemann sum with m = 3, n = 2, and take the sample
point t o be the upper right comer of each square.
(b) Ube the Midpoint Rule to estimate the volume of the solid
in ~ r v t (a).
2. If R = :- 1. 31 X [O, 21, use a Riemann sum with m = 4,
11 = 2 io estimate the value of !:jR (y' - 2.r') dA. Take the
samplc p i n t s to be the upper left comers of the squares.
3. ( a) Cw a Riemann sum with m = n = 2 to estimate the value
of :; sin(x + y) dA, where R = [0, rr] X [0, x]. Take the
sample points to be lower left corners.
(b) U= the Midpoint Rule to estimate the integral in part (a).
4. (a) Esumate the volume of the solid that lies below the surface
z = .r + 2y2 and above the rectangle R = [O, 21 X [O, 41.
Cse a Riemann sum with m = n = 2 and choose the
sample points to be lower right comers.
( b) U x the Midpoint Rule to estimate the volume in part (a).
5. A tablc of values is given for a function f (.r, y) defined on
R = [I. 31 X [O, 41.
(a) Estimate J'h f ( x , y) dA using the Midpoint Rule with
m = n = 2 .
6. A 20-ft-by-30-ft swimming pool is filled with water. The depth
is measured at 5-ft intemals. starting at one comer of the pool.
and the values are recorded in the table. Estimate the volume of
water in the pool.
Let V be the volume of the solid that lies under the graph of
f ( x, y) = J52 - x2 - y' and above the rectangle given by
2 S x S 4, 2 S y 5 6. We use the lines x = 3 and y = 4 to
SECTION 15.2 ITERATED INTEGRALS 1111 959
divide R into subrectangles. Let L and U be the Riemann sums
computed using lower left comers and upper right comers,
respectively. Without calculating the numbers V, L, and U,
arrange them in increasing order and explain your reasoning.
8. The figure shows level curves of a function f in the square
R = [O, 21 X [O, 21. Use the Midpoint Rule with rn = n = 2
to estimate I.[, f (x, y) dA. How could you improve your
estimate?
1 1 - 1 3 Evaluate the double integral by first identifying it a rhe
volume of a solid.
. .
A contour map is shown for a function f on the square
1 1 . ! J R3dA. R = {(x,~) / 2 S x G 2, 1 c y S 6)
R = [O,4] X [O,4]. 12. f i ( 5 - x ) dA, R = {(x,y) ( 0 S x < 5, 0 S y s 3)
(a) Use the Midpoint Rule with rn = 17 = 2 to estimate the
value of I(, f ( x, y) dA.
m. [ i R ( 4 - 2 y ) d A , R = [ O, I ] X [ O, l ]
(b) Estimate the average value of f .
14. The integral ((, J- dA, where R = [O, 41 X [O. 2;.
represents thevolume of a solid. Sketch the solid.
15. Use a programmable calculator or computer (or the sum
command on a CAS) to estimate
where R = [O, 1 1 X [0, I]. Use the Midpoint Rule with Lie
following numbers of squares of equal size: I , 4, 16, 64. 256.
and 1024.
16. Repeat Exercise 15 for the integral sin(* + A) dA
17. I f f is a constant function. f (x, y) = k , and
R = [u. b] X [c. dl , show that .[iR k dA = k(b - a) ( d - cl.
10. The contour map shows the temperature. in degrees Fahrenheit,
18. Use the result of Exercise 17 to show that
at 4:00 PM on February 26, 2007, in Colorado. (The state I
o s J"[ sin r x cos rjl d~ s -
measures 388 mi east to west and 276 mi north to south.) Use 32
R
the Midpoint Rule with m = n = 4 to estimate the average
temperature in Colorado at that time.
where R = 10, i] X [a, i].
Recall that it is usually difficult to evaluate single integrals directly from the definition of
an integral, but the Fundamental Theorem of Calculus provides a much easier method The
evaluation of double integrals from first principles is even more difficult, but in this sec-
15.2 ITERATED INTEGRALS
v
960 1 1 [ ( CHAPTER 15 MULTIPLE INTEGRALS
tion we see hoy to express a double integral as an iterated integral, which can then be eval-
uated by calculating two single integrals.
Suppose that f is a function of two variables that is integrable on the rectangle
R = [a, b] X [c, d l . We use the notation scd f ( x , Y ) d~ to mean that x is held fixed and
f i x, y) is integrated with respect to y from y = c to = d. This procedure is called par-
tial integration with respect to y. (Notice its similarity to partial differentiation.) Now
rd f ( x, y) dy is a number that depends on the value of .r. so it defines a function of x :
" C
I
A( x) = J y f i x, y) d~ 0
If we now integrate the function A with respect to x from x = a to x = b, we get
The integral on the right side of Equation 1 is called an iterated integral. Usually the
brackets are omitted. Thus
means that we first integrate with respect to y from c to d and then with respect to x from
a to b.
Similarly, the iterated integral
means that we first integrate with respect to x (holding y fixed) from x = a to x = b and
then we integrate the resulting function of y with respect to y from y = c to y = d. Notice
that in both Equations 2 and 3 we work from the inside out.
EXAMPLE I Evaluate the iterated integrals.
( a) lf j:x2y dy dx
SOLUTI ON
(a) Regarding x as a constant, we obtain
Thus the function A in the preceding discussion is given by A( x ) = 1x2 in this example.
We now integrate this function of x from 0 to 3:
s Theorem 4 IS named after t he Itallan rnathe-
mat cl an Gu~do Fub~r ~i (1879-19431, who provec
a very general version of t hi s theorem ~n 1907.
But t he version for continuous functions was
known t o t he French r nat hemat ~c~an Augustin-
LOUIS Cauchy almost a century earl er.
FI GURE I
Visual 15.2 illustrates Fubinis
Theorem by showing an animation of
Figures I and 2.
SECTION 15.2 ITERATED INTEGRALS jJll 961
(b) Here we first integrate with respect to x:
Notice that in Example 1 we obtained the same answer whether we inte=mred with
respect to y or x first. In general, it turns out (see Theorem 4) that the two iterated integrals
in Equations 2 and 3 are always equal; that is. the order of integration does nor matter.
(This is similar to Clairaut's Theorem on the equality of the mixed partial derivarives.)
The following theorem gives a practical method for evaluating a double inrsgral by
expressing it as an iterated integral (in either order).
FUBINI' S THEOREM I f f is continuous on the rectangle
R = {(x. x) ( a C x a b, c 6 y S d}, then
More generally. this is true if we assume that f is bounded on R, f is discontin-
uous only on a finite number of smooth curves, and the iterated integrals esisr-
The proof of Fubini's Theorem is too difficult to include in this book, but we can at least
give an intuitive indication of why it is true for the case where f(x, y) 2 0. R e d l that if
f is positive. then we can interpret the double integral jjR f (x. J*) dA as the \ olume 1' of
the solid S that lies above R and under the surface : = f (x, y). But we have another for-
mula that we used for volume in Chapter 6, namely,
where A(x) is the area of a cross-section of S in the plane through x- perpendicuiar to the
x-axis. From Figure 1 you can see that A(x) is the area under the curve C whose =quation
is z = f (x, y). where x is held constant and c a y d. Therefore
and we have
A similar argument, using cross-sections perpendicular to the y-axis as in Figure 2. shows
that
FI GURE 2
962 1 1 1 1 CHAPTER 15 MULTI PLE
a Notice t he n e 3 t i v e answer In Example 2;
nothing i s wrort; wi t h t hat . The I ~ ~ n c t ~ o n f In
t hat example is ~ 3 t a posltlve function, so i t s
integral doesn r -.present a volume. From
Figure 3 we SEE %at f I S always negative on
R, so the value :i t he integral is t he negative
of the volume i ~ e s above t he graph o f f
and below R
FI GURE 3
I NTEGRAL S
EXAMPLE 2 Evaluate the double integral ('& ( x - 3y2)d4. where
R = {(x. y) 0 S x S 2, 1 S y G 2). (compare with Exampl e 3 in Section 15.1.)
SOLUTI ON I Fubini's Theorem gives
SOLUTION 2 Again applying Fubini's Theorem. but this t i me integrating with respect to .r
first. we have
L!i EXAMPLE 3 Evaluate !:iR y sini.ry) dil. where R = [ I . 21 X [O. T I .
I
SOLUTION I If we first integrate with respect to .x, we pet
I
.. - -
( 1 y sin(xy) dA = (' i2 y sin(.ry) dx dy = I [-cosi.;y)]
dy
. . . o - 1 . I
R
=
(-cos 2y + cos yi 6.
-
- - sin 2y + sin j]," = O
SOLUTION 2 If we reverse the order of integration, we get
P For a il-~nc: c r f t hat takes on both posltlve
and negative ,.E ces, j:ik f (x. F) dA is a differ-
ence of volurr%s: I/, - V?, where V, is t he v o -
urne above R 273 bel ow t he graph o f f and V2 I S
. . .- .-
( 1 y sin(.~?) dA = ( - ( " Y sini .YJ -\) dy dx
" " " I . o '
R
t he vol l ~rne be -.v R anlJ above the graph The To evaluate the inner integral, we use integration by parts with
fact that the - q r a In Exarnpe 3 I S 0 means
t hat these nvc . ~i ur nes V, and V2 are equal
(See Figure .L
clv = sin(.xy) cly
COS(XY)
el11 = dy
21 = --
X
1
.- 0 and so I T y sin(.ry) dy = -
- 1
. 0
1
T
- -
n c o s TTTX 1
-
\=7
o x + 7 X- [sin(xy)];,=,
FI GURE 4 - -
n c o s n x sin 7. x
- +T
X 1-
SECTION 15.2 I T ERAT EDI NT EGRAL S 11() 963
m In Example 2 Sol ~~t l ons 1 and 2 are equallk
stralghtfonvard but i n Example 3 t he flrst s o u -
t l on IS much easier t han t he second one There-
fore when we evaluate double Integrals i t i s
wl se t o choose t he order of lntegratlon t hat gi ves
slmpler ~nt egr al s
If we now integrate the first term by parts with u = - l/x and dv = ncos ns dx. we get
du = dx/x2, u = sin nx, and
sinnx
dx=--- dx
x
( c x sin n-x sin n x
Therefore +-)dx= xu -- x
and so
-
sin 277
- --
+ sin .rr = O
2
EXAMPLE 4 Find the volume of the solid S that is bounded by the elliptic paraboloid
x2 + 2y2 + z = 16. the planes x = 2 and y = 2. and the three coordinate pl ans .
SOLUTI ON We first observe that S is the solid that lies under the surface z = 16 - .y2 - 2y2
and above the square R = [O. 21 X [O. 21. (See Figure 5. ) This solid was considxed in
Example 1 in Section 15.1. but we are now in a position to evaluate the double integral
using Fubini's Theorem. Therefore
FI GURE 5
In the special case where f (x. y) can be factored as the product of a function of x only
and a function of y only, the double integral of f can be written in a particulxl\ simple
form. To be specific. suppose that f ( x, g) = g(x)h(y) and R = [a. b] X [c. dl. Then
Fubini's Theorem gives
In the inner integral. y is a constant, so h(y) is a constant and we can write
since J: g(x) dx is a constant. Therefore, in this case, the double integral of f can be writ-
ten as the product of two single integrals:
JJ' q(x) h(g) dA = 6 : g(x) dx r i d C h( y ) dy
where R = [a, b] X [c, dl
R
I
964 ( 1 1 1 CHAPTER 15 MULTI PLE I NTEGRAL S
EXAMPLE 5 If R = [O. ~ / 2 ] X [0, 4 2 1 . then, by Equation 5,
"1
P The i n c t c - .'I I. yi = si n .Y cosy n
Examols 5 I S ccs :'be on R, so the Integral repre-
sents t?e vo,-5 2f the solid t hat lies above R
and be,c,v t he I-zch of j shown ~n Figure 6
FIGURE 6
-
15.2 EXERCI SES
1-2 Find ,' ' ,;I .I-. x) c1.r and 1; f (.r. g) dg.
-
- [ - c os x , 17:2 [ . sin!, 1.: 2 = 1 . 1 = 1
3- 14 Ca1cz;sts the iterated integral.
-- .- * I . I .-
3. 1 , I,, I ! - 4. y~) c1.u dg 4. !,! 1; 14.r' - 9.u2!') c f ~ t l i
. .
%T 2 - 5
5. 1' i-'.r sinyc1xti.r
6. j7
!II cosgc1.rdg
. , I . I ,
'I ' I
., .,
'
7. 1 / (I7 - ?.)"1.r d! 8. / - / - - dg ci*
. I .I1 . ' ) . I '
. - - . [ . 1 .
iig tb 10. jl, !; e' , -" r1.r clg
- ! .
. I . [ -
I I . 1,) jU ( [ L - ~1~ c1u c1u 12. .(" jn . Y ~ " . Y ~ + g2 dxdr
19. I'[ x sin(.u + g) tlA, R = [O. 7/61 X [O. n/3]
. .
R
" .Y
!!
clA. R = [O. !] X [O. 11
21. [ r e 1 R = [O. I ] X [ O 21
R
. .
2 . \ - 4 R = [ I . 21 * [O. 11
R
1
23-24 Sketch the solid whose \olurne 1s glven by the Iterated
f
lntesral
i
15-22 Calsxlate the double integral.
24. (2 - r - !' I i l l t1.r
15. 1.i(6.r'>-' - 5x4) d.4. R = {(.r.y) o S . r s 3. 0 c y s 1)
. .
R
25. Find the volume of the solid that lies under the plane
3x + 2~ + : = 12 and abo\-e the rectangle
R = { ( s . g ) ) O < . r S 1. -2 S g S 3).
26. Flnd the bolurne of the solid that lies under the hyperbolic
% I
.1
paraboloid: = 4 + .r2 - ?.' and above the square
R = [-1, 1] X [0, 21.
SECT I ON 15.3 DOUBLE INTEGRALS OVER GENERAL REGIONS 1111 965
Find the volume of the solid lying under the elliptic
-
L! L 34. Graph the solid that lies between the surfaces
paraboloid x2/4 + v2/9 + z = 1 and above the rectangle
z = e ~ x 2 c o s ( x 2 + y 2 ) a n d z = 2 - ~ 2 - Y 2 f o r ~ ~ ~ I ,
R = [ - I , 11 X [-2, 21. ( y I S 1. Use a computer algebra system to approximate the
volume of this solid correct to four decimal places.
28. Find the volume of the solid enclosed by the surface
I = 1 + exsi nyandt hepl anesr = ?1, y =O, y = v.
35-36 Find the average value off over the given rectangle.
and : = 0.
35. f ( x, y) =xZy, Rhasvert i ces(-1, 0), (-1, 5), (1. 5i . (1, 0)
29. Find the volume of the solid enclosed by the surface
z = x sec2! and the planes z = 0. x = 0, x = 2, y = 0, 36. f(x, y) = e 1 J Z , R = [O, 41 x [0, I ]
and y = n/4.
30. Find the volume of the solid in rhe first octant bounded by
i s , 37. Use your CAS to compute the iterated integrals
the cylinder i = 16 - x 2 and the plane y = 5.
x -.
3 1. Find the volume of the solid enclosed by the paraboloid
.o 1' .o [lady ( X + y)"
dx and
.o [ I .o [I- (x + !)> . &A- d~
z = 2 + x 2 + ( T- 2) 2a ndt he pl a ne s i = 1. x= 1 , x = - I ,
y = 0, and y = 4.
Do the answers contradict Fubini's Theorem? Explair. what
is happening.
32. Graph the solid that lies between the surface
,- = 2xy/(x2 + 1) and the plane 1 = x + 2y and 1s bounded
38. (a) In what way are the theorems of Fubini and Cl aki ut
by the planes x = 0, x = 2, y = 0. and y = 4. Then find its
similar?
volume.
(b) I f f (x, y) is continuous on [a, b] X [c, d ] and
33. Use a computer algebra system to find the exact value of the
integral j?^, xSv' eA' dA, where R = [0, 11 X [0, 11. Then use
g(x, y) = Jx o J' C f (.s, t ) dt ds
the CAS to draw the solid whose volume is given by the
integral. for a < x < b, c < v < d, show that g,, = g,, =.'(X,!I.
15.3 DOUBLE INTEGRALS OVER GENERAL REGIONS
For single integrals, the region over which we integrate is always an interval. But for
double integrals, we want to be able to integrate a function f not just over rectangles but
also over regions D of more general shape, such as the one illustrated in Figure 1. Rt sup-
pose that D is a bounded region, which means that D can be enclosed in a recTungular
region R as in Figure 2. Then we define a new function F with domain R by
{ x y if (x. y) is in D
F( x, y) =
if (x, y ) is in R but not in D
FI GURE I FI GURE 2
966 CHAPTER I S MULTIPLE INTEGRALS
If F is integrable over R, then we define the double integral of f over D by
graph of f
e
:
Y
Definition 2 makes sense because R is a rectangle and so J'S, F( x, y ) dA has been previ-
11f ( x, y ) dA = 11F( x, y) dA
where F is given by Equation 1
D R
ously defined in Section 15.1. The procedure that we have used is reasonable because the
FI GURE 3
values of F( x, y) are 0 when ( x, y ) lies outside D and s o they contribute nothing to the inte-
gral. This means that it doesn't matter what rectangle R we use as long as it contalns D.
In the case where f ( x, v) 3 0, we can still interpret fl, f ( x, y ) dA as the volume of the
A
solid that lies above D and under the surface z = f (x. y) (the graph of f ). You can see that
graph of
this is reasonable by companng the graphs of f and F in Flgures 3 and 4 and remember-
p/
ing that J'S, F( x, y ) dA is the volume under the graph of F.
Figure 4 also shows that F is likely to have discontinuities at the boundary points of
I '
D. Nonetheless, i f f is continuous on D and the boundary curve of D 1s "well behaved
r/\w
(in a sense outslde the scope of this book), then it can be shown that J'(k F( x, y) dA exists
and therefore ". ('(, f ( x , y ) dA exists. In particular, this is the case for the following types
of regions.
A plane region D 1s said to be of type I if it lies bemeen the graphs of two continuous
FI GURE 4 functions of x, that is,
where gj and g2 are continuous on [a, b] . Some examples of type I regions are shown in
Figure 5.
FI GURE 5 Some typeIregions
In order to evaluate Jk f ( x, y) dA when D is a re-on of type I, we choose a rectangle
)A
R = [a, b] x [c, d ] that contains D, as in Figure 6, and we let F be the function given by
Y = g2ix)
d ~ - - b
Equation 1; that is, F agrees with f on D and F is 0 outside D. Then, by Fubini's Theorem,
I
q
4
' b ' d
jj f ( x, Y ) dA = jj F( x. Y ) dA = 1 1 F( x. Y ) d y dx
1
s
3
.a .C
D R
G
7- 1 I Y = q, ( r ) I * Therefore
Observe that F( x, y ) = 0 if y < gl ( x) or y > g2( x) because ( x, y) then lies outside D.
O1
. b r
Jcd F( X. J ) d y = J 1 g 2 ( x ) F( x, V ) d y = jg (af(4 y) dy
FI GURE 6 gt(x) - sdx)
SECTION 15.3 DOUBLE INTEGRALS OVER GENERAL REGIONS 1111 967
FI GURE 7
Some type I1 regions
1 1 '
FI GURE 8
t'
6
because F( x, y) = f ( x , JJ) when gl ( x) 6 y 6 g2( x) . Thus we have the following formula
that enables us to evaluate the double integral as an iterated integral.
I f f is continuous on a type I region D such that
The integral on the right side of (3) is an iterated integral that is similar to ths ones we
considered in the preceding section, except that in the inner integral we regard .r as being
constant not only in f ( x, y) but also in the limits of integration, g ~ ( x ) and g2(x).
We also consider plane regions of type 11, which can be expressed as
rn
D = { ( x, y) I c y a d, S x a h?( y) }
where hl and hz are continuous. Two such regions are illustrated in Figure 7.
Using the same methods that were used in establishing (3), we can show that
1
where D is a type 11region given by Equation 4.
EXAMPLE I Evaluate [rD (x + 2y) dA, where D is the region bounded by the
. .
parabolas y = 2x2 and y = 1 + x'.
SOLUTI ON The parabolas intersect when 2x2 = 1 + I?, that is, x 2 = 1, so x = 1. We
note that the region D. sketched in Figure 8, is a type I region but not a type I 1 r s j o n
and we can write
Since the lower boundary is y = 2x2 and the upper boundary is y = 1 + .a2, Equation 3
gives
fi2.4,
m] When wh set up a double integral as in Example 1, it is essential to draw a
.:
' I
diagram. Often it is helpful to draw a vertical arrow as in Figure 8. Then the limits of
\ = -z
i
integration for the inner integral can be read from the diagram as follows: The arrow
I
starts at the lower boundary y = g,(x), which gives the lower limit in the integral, and
1 , the arrow ends at the upper boundary y = gz(x), whlch gives the upper limit of integration.
A L = t -
For a type I1 reg~on the arrow is drawn horizontally from the left boundary to the nght
i
I D
boundary.
p
EXAMPLE 2 Find the volume of the solid that lies under the paraboloid z = x' + y' and
1 0 : - 7 t
above the region D in the xy-plane bounded by the line y = 2.r and the parabola y = x2.
I
SOLUTION I From Figure 9 we see that D is a type I region and
FI GURE 9
\ A
1 - * 12. 41 Therefore the volume under 3 = .u' + > ' and above D is
I
FI GURE 1 0
D 3s a t! pe I1 region
SOLUTION 2 From Figure 10 we see that D can also be written as a type I1 region:
3 Figure " shc~,~,n -e solid whose volume
I S calculared in E: t -npi e 2. It lies above the D = { ( x , y ) 10 y 4, k y s .Y .,G}
. r ~- pl ane, below TE oaraboloid = x' + j',
'and bet wf en t h ~ : zne = 2x and the
parabolic c./linoe- .-.- = .r'.
Therefore another expression for V is
FI GURE I I
PI EXAMPLE 3 Evaluate jJL1xy dA, where D is the region bounded by the line y = x - 1
and the parabola y' = 2x + 6.
SOLUTION The region D is shown in Figure 12. Again D is both type I and type 11, but the
description of D as a type I region is more complicated because the lower boundary con-
sists of two parts. Therefore we prefer to express D as a type I1 region:
~ = { ( x , ~ ) 1 - 2 < y s 4, i Y2 - 3 ~ . x ~ y + I}
SECTI ON 15.3 DOUBLE INTEGRALS OVER GENERAL REGIONS 111) 969
/
J=-
FI GURE 1 2 (a) D as a type I region (b) D as a type I1 region
Then ( 5 ) gives
If we had expressed D as a type I region using Figure 12(a), then we would have
obtained
but this would have involved more work than the other method.
I
EXAMPLE 4 Find the volume of the tetrahedron bounded by the planes x + 2y - z = 2.
FI GURE 13
s = 2y, x = 0. and 2 = 0.
5OLUTl ON In a question such as this, it's wise to draw two diagrams: one of the three-
\ ' A dimensional solid and another of the plane region D over which it lies. Figure 13 shows
~ + 2 y = 2
the tetrahedron T bounded by the coordinate planes x = 0, z = 0, the vertical plan?
(or J = 1 - .x!2i
s = 2y, and the plane x + 2?; + z = 2. Since the plane x + 2y + z = 2 intersects the
I?;-plane (whose equation is z = 0) in the line x + 2y = 2, we see that T lies abo\-s the
' > (I. ; )
triangular region D in the xy-plane bounded by the lines x = 2y, x + 2y = 2, and -T = 0.
(See Figure 14.)
The plane x + 2y + z = 2 can be written as z = 2 - x - 2y, so the required volume
0
p
lies under the graph of the function z = 2 - x - 2y and above
1
Y
D = {(x,y) 0 Ss x S 1, x/ 2 s y 1 - X/2}
FI GURE 1 4
7
970 1111 CHAPTER I S MULTIPLE INTEGRALS
Therefore .
V = 11 (2 - x - 2y) dA = ( 2 - .r - 2y) d y d x
D
-
I " v=i -*/ ?
- j,, [2y - xY - y-].=y. d x
FI GURE 15
D as a t? pz I region
FI GURE 16
D as a t ype !J region
P1 EXAMPLE 5 Evaluate the iterated integral u u (': sin(?.') d y dx.
SOLUTI ON If we try to evaluate the integral as it stands. xve are faced with the task of first
evaluating /' sin(?" ddy. But it's impossible to do so in finite terms since 1 sin(?') d y is
not an eleGentary function. (See the end of Section 7.5.) So we must change the order
of integration. This is accomplished by first expressing the given iterated integral as a
double integral. Using (3) backward, we have
where D = { ( x , ~ ) 0 s .r a I , .r s .v < 1)
We sketch this region D in Figure 15. Then from Figure 16 we see that an alternative
description of D is
D = { ( x , y ) I O ~ y < 1. OG. r Cy }
This enables us to use (5) to express the double integral as an iterated integral in the
reverse order:
=; ( I - cos 1)
PROPERTIES OF DOUBLE INTEGRALS
We assume that all of the following integrals exist. The first three properties of double inte-
grals over a region D follow immediately from Definition 2 and Properties 7, 8, and 9 in
Section 15.1.
SECT I ON 15.3 DOUBLE INTEGRALS OVER GENERAL REGIONS 1111 971
I f f ( x , y) 2 g(x, y) for all ( x , y) in D, then
The next property of double integrals is similar to the property of single i nt e- ds given
by the equation jc: f ( x ) dx = [ : f ( x ) dx + -i(b f ( x ) dx.
If D = DL U D2, where D, and D2 don't overlap except perhaps on their boundaries
(see Figure 17), then
FI GURE 17
FI GURE 18
Property 9 can be used to evaluate double integrals over regions D that are neithcr type I
nor type I1 but can be expressed as a union of regions of type I or type 11. Figure I S illus- .
trates this procedure. (See Exercises 5 1 and 52.)
(a) D is neither type I nor type 11.
(b) D =Dl U D2. Dl is type I, D, is type ri.
The next property of integrals says that if we integrate the constant function f (.r. - I = 1
over a region D, we get the area of D:
El v, 1)
Y
Figure 19 illustrates why Equation 10 is true: A solid cylinder whose base is D and nhose
FI GURE 19
height is 1 has volume A(D) . 1 = A( D) , but we know that we can also write its volume
Cylinder with base D and height 1
as !:i, 1 dA.
Finally, we can combine Properties 7, 8, and 10 to prove the following propem-. (See
Exercise 57.)
/ If rn S ,f(x, y ) 6 M for all ( x, Y ) in D. then
i
F
972 1111 CHAPTER15 MULTIPLEINTEGRALS
EXAMPLE 6 Use Property 11 to estimate the inte-ml (jD e"""'"' dA, where D is the disk
with center the origin and radius 2.
SOLUTI ON Since - 1 S sin x S 1 and - 1 a cos J S 1, we have - 1 G sin x cos y S 1 and
therefore
Thus, using rn = e - ' = I/e, M = e, andA(D) = ~ ( 2 ) ' in Property 11, we obtain
-
15.3 EXERCISES
1-6 E.. aluate the iterated integral.
.. .
1. ) - 1 " s y ' d x d y
. . '
2. j' (? (X - y) dy dx
00 .?i
- ..
3. ( l ' ( 1 +2y) dyd. \ -
. ..
4. J; : 1: ' xy d,x d y
- 5, - rT: .COIH
- . ! o
eqn R dr do 6. j' i' JT-VL du dv
"0 "(I
- - - -- - - - - -- - - - -
7-1 8 E~~a l ua t e the double integral
13. I'(' r 20s ?. dA, D is bounded by y = 0, y = ,s2, ,s = 1
. .
n
, .
14. ( 1 I.r - y ) d A, Di s boundedbyy = , h a n d y =x'
. .
0
IS. [ ! ' d4.
D
D is d x triangular region with vertices ( 0, 2), (1, l), (3, 2)
16. jj .rvl o ~ . D is enclosed by x = o and x = Jm
-.
D
'm 1.i (2.w - y) dA,
I
I
d W
I> I
D is bounded by the circle with center the origin and radius 2
I
I
. .
I
18. J? 214. dA, D is the uiangular region with vertices (0. O),
1
I
I)
( I , 2) . and (0, 3)
19-28 Find the volume of the given solid.
19. Under the plane .r + 2~ - : = 0 and above the region
bounded by ji = x and = xJ
20. Under the surface r = 2 r + y' and above the region bounded
by, s =y' a nd. r =y:
Under the surface ; = r\ and above the tnangle with vertlcec
(1, I ) , (4, I ) , and ( 1, 2 )
22. Enclosed by the paraboloid 7 = x' + 3y' and the planes x = 0,
y = l , y =. r , z =O
23. Bounded by the coordinate planes and the plane
3 x + 2 v + z = 6
24. Bounded by the planes z = x, y = x, x + y = 2, and z = 0
25. Enclosed by the cylinders r = x2, y = x' and the planes
z = O. y = 4
26. Bounded by the cylinder x' + r' = 4 and the planes w = 7y,
x = O , - = 0 in the first octant
27. Bounded by the cylinder .I-' + y' = 1 and the planes y = 2 ,
x = 0 , z = 0 in the first octant
28. Bounded by the cylinden .s' + y2 = r' and y' + z' = r'
--
29. Use a graphing calculator or computer to estmate the
x-coordinates of the polnts of lntersectlon of the curves y = n4
and y = 3.r - x! If D 1s the region bounded by these curves.
estimate jjll . w dA.
. I
SECTI ON 15.3 DOUBLE INTEGRALS OVER GENERAL REGIONS 1111 973
30. Find the approximate volume of the solid in the first octant 5 1-52 Express D as a union of regions of type I or type Il and
that is bounded by the planes = x, 2 = 0, and z = x and evaluate the integral.
the cylinder J = cos x. (Use a graphing device to estimate
the points of intersection.)
". [[x2dA
52. ". /I ,. dA
I1 D
3 1-32 Flnd the volume of the solld by subtracting two volumes
1, \=I-,'
3 1. The solid enclosed by the pa r a bhc cylinders
I = 1 - 1'. \ = x' - I and the planes 1 + \ + : = 2,
\=11+1)y
2 x + 2 \ - ; + l O= O
- 1 , 0 , / I :
33-34 Sketch the solid whose volume is given by the iterated
I
integral.
32. The solid enclosed by the parabolic cyllnder 1 = .I' and the
planes ; = 3v. : = 2 + J
I
>
- 1
-
- 1 . I - *
53-54 Use Property 11 to estimate the value of the intfgral.
33. ( ( (1 - X - , . I ~ ~ > ~ X 34. (' [ ' - " ( I - ~ ) d ? , d ~ . .
.(I .O .O .(I
53. 1 1 e - ~ x 2 + , 2 ~ dA. Q is the quarter-circle with center the origin
/
/
,' I
1
and radius 4 in the first quadrant
-
tni 35-38 Use a computer algebra system to find the exact volume . .
of the solid.
54. ( 1 si nyx + q ) dA. T is the triangle enclosed by th- lines
. .
35. Under the surface r = x' j ' ' .~-y' and above the region
7
y = 0 , _ v = 21. and x = 1
bounded by the curves y = .x' - .Y and y = x 2 + .r for .u 2 0
36. Between the paraboloids z = 2.x' + y' and
55-56 Find the average value of,f over rezion D.
-
r - 8 - x' - 2y2 and inside the cylinder x' + !.' = 1
55. , f(a. y) = xq. D is the triangle with vertices (0, 0). 11. 0).
and (1. 3)
37. Enclosed by r = 1 - x' - j' and : = 0
56. f (s. y) = x sin y. D is enclosed by the curves y = 0.
38. Enclosed by r = x 2 + j' and r = 2). y = x2. and x- = 1
L---
57. Prove Property 11.
39-44 Sketch the region of integration and change the order of
~ntegration.
il
& In evaluating a double integral over a region D. a sum of
-1
' I - 4
40.
Ji f i r . j ) d! dx
iterated integrals was obtained as follows:
I )
Sketch the region D and express the double integral 3 an
iterated integral with reversed order of integration.
" - - 1111 . I -n,4
3: 1: l o ,f (x. y) dy dx
. .
44. jo Jarcun, f (1, ?.) d!. d x 59. Evaluate !:(,, ( x2 tan x + y' + 4) dA. where
D = { ( x , y ) 1 x' + y' G 2). [Hint: Exploit the fact that
D is symmetric with respect to both axes.]
45-50 Evaluate the integral by re\-ersing the order of integration.
60. Use symmetry to evaluate {I, ( 2 - 3x + 4?) d ~ , v-h-re D
-
is the region bounded by the square with vertices ( ~ 5 . 0)
46. i XT COS(X?) dx d y
.O . i
and (0, ?5) .
61. Compute J:j, J1 - x2 - y' dA, where D is the disk
x 2 + y' < I , by first identifying the integral as the \-olume
of a solid.
" I -,- 2
cos x J - dx d\
-
tnr 62. Graph the solid bounded by the plane x + 1 + r = 1 and
49. Ju Jarciln r
the paraboloid ; = 4 - x2 - J and find its exact \ ol ume
(Use your CAS to do the graphing. to find the equations of
50. [ el4 dr d )
the boundary curves of the region of integration, and to ekal-
uate the double ~ntegral.)