TRIPLE INTEGRALS

ELECTRONIC VERSION OF LECTURE

HoChiMinh City University of Technology

Faculty of Applied Science, Department of Applied Mathematics

(HCMUT-OISP) TRIPLE INTEGRALS 1 / 80

OUTLINE

1 DEFINITION OF TRIPLE INTEGRALS

2 TRIPLE INTEGRALS IN RECTANGULAR (CARTESIAN)
COORDINATE SYSTEM

3 TRIPLE INTEGRALS IN CYLINDRICAL COORDINATES

4 TRIPLE INTEGRALS IN SPHERICAL COORDINATES

5 APPLICATIONS OF TRIPLE INTEGRALS

6 MATL AB

(HCMUT-OISP) TRIPLE INTEGRALS 2 / 80

 Definition of Triple Integrals Evaluating mass of a solid

Suppose that the density function of a solid object

that occupies the region Ω is f (x, y, z), in units of
mass per unit volume, at any given point (x, y, z).

Find the mass m of this solid, if the solid Ω is the

box?
n o
Ω = (x, y, z) : a É x É b, c É y É d , r É z É s

(HCMUT-OISP) TRIPLE INTEGRALS 3 / 80

 Definition of Triple Integrals Evaluating mass of a solid

We divide Ω into sub-boxes by dividing [a, b] into m

subintervals [x i −1 , x i ] of equal width ∆x = b−a m
,
dividing [c, d ] into n subintervals [y j −1 , y j ] of equal
width ∆y = dn−c , and dividing [r, s] into p subintervals
[z k−1 , z k ] of equal width ∆z = s−r
p
·

(HCMUT-OISP) TRIPLE INTEGRALS 4 / 80

 Definition of Triple Integrals Evaluating mass of a solid

The volume of each sub-box Ωi j k is

∆V = ∆x.∆y.∆z.

The triple Riemann sum is approximately equal to

the mass of the solid Ω
m X
X p
n X
M≈ f (x i∗j k , y i∗j k , z i∗j k )∆x∆y.∆z.
i =1 j =1 k=1

where the sample point (x i∗j k , y i∗j k , z i∗j k ) ∈ Ωi j k .

(HCMUT-OISP) TRIPLE INTEGRALS 5 / 80

 Definition of Triple Integrals Definition

DEFINITION 1.1
The triple integral of f (x, y, z) over the box Ω is
Ñ Ñ
f (x, y, z)dV = f (x, y, z)d xd yd z =
Ω Ω

m X
X p
n X
= lim f (x i∗j k , y i∗j k , z i∗j k )∆x∆y.∆z.
m,n,p→∞
i =1 j =1 k=1

if this limit exists.

(HCMUT-OISP) TRIPLE INTEGRALS 6 / 80

Triple Integrals in rectangular (Cartesian) coordinate system Fubini’s Theorem

THEOREM 2.1
If f (x, y, z) is continuous on the rectangular box

Ω : a É x É b, c É y É d , r É z É s,

then
Ñ Zb Zd Z s
f (x, y, z)dV = f (x, y, z)d xd yd z
Ω a c r
Zb Zd Z s
   

=   f (x, y, z)d z  d y  d x.
a c r

(HCMUT-OISP) TRIPLE INTEGRALS 7 / 80

Triple Integrals in rectangular (Cartesian) coordinate system Fubini’s Theorem

EXAMPLE 2.1
Ñ
Evaluate I = x y z 2 dV, where Ω is the rectangular
Ω
box given by

Ω = 0 É x É 1, −1 É y É 2, 0 É z É 3
© ª

(HCMUT-OISP) TRIPLE INTEGRALS 8 / 80

Triple Integrals in rectangular (Cartesian) coordinate system Fubini’s Theorem

Z3 Z2 Z1 Z3 Z2 Z1
 

I= x y z 2 d xd yd z =  x y z 2 d x  d yd z
0 −1 0 0 −1 0
Z3 Z2 · ¸x=1 Z Z2
3
2x2 y z2
= yz · d yd z = d yd z
2 x=0 2
0 −1 0 −1
Z3 Z2 Z3 · Z3
 
¸ y=2
y z2  z2 y 2 3z 2
=  dy dz = dz = dz
2 4 y=−1 4
0 −1 0 0
3 ¸3
z
·
27
= = ·
4 0 4

(HCMUT-OISP) TRIPLE INTEGRALS 9 / 80

Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

Ñ
Evaluate I = f (x, y, z)dV
Ω
1
The projection of Ω onto the Ox y is D
2
The lower boundary of the solid Ω is the surface
with equation z = z 1 (x, y)
3
The upper boundary of the solid Ω is the surface
with equation z = z 2 (x, y)
 
Ñ Ï z 2Z(x,y)

I= f (x, y, z)dV = f (x, y, z)d z  d A

 

Ω D z 1 (x,y)

(HCMUT-OISP) TRIPLE INTEGRALS 10 / 80

Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

(HCMUT-OISP) TRIPLE INTEGRALS 11 / 80

Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

EXAMPLE 2.2
Evaluate Ñ
I= (x + 2y)dV,
Ω

where Ω is the region bounded by x 2 É y É x, 0 É z É x.

(HCMUT-OISP) TRIPLE INTEGRALS 12 / 80

Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

(HCMUT-OISP) TRIPLE INTEGRALS 13 / 80

Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

 x 
Ï Z Ï
£ ¤z=x
I=  (x + 2y)d z  d A = xz + 2y z z=0 d A
D 0 D
Ï
= (x 2 + 2y x)d A
D

where D = {(x, y) : 0 É x É 1, x 2 É y É x}.

Z1 Zx Z1
 
¡ 2 £ 2 ¤ y=x
x y + x y 2 y=x 2 d x
¢
I =  x + 2y x d y  d x =
0 x2 0
Z1
2
x3 − x4 + x3 − x5 d x = ·
¡ ¢
=
15
0
(HCMUT-OISP) TRIPLE INTEGRALS 14 / 80
Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

EXAMPLE 2.3
Evaluate Ñ p
I= x 2 + z 2 dV,
Ω

where Ω is the region bounded by x 2 + z 2 É y É 4.

(HCMUT-OISP) TRIPLE INTEGRALS 15 / 80

Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

(HCMUT-OISP) TRIPLE INTEGRALS 16 / 80

Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

The solid Ω is bounded by the right boundary y = 4,

and the left boundary y = x 2 + z 2 . To find the
projection of Ω onto Oxz, we can find the
intersection of y = x 2 + z 2 , y = 4.
y = x2 + z2 x2 + z2 = 4
½ ½
⇔
y =4 y =4

Therefore the intersection of y = x 2 + z 2 , y = 4 is the

intersection of x 2 + z 2 = 4 and y = 4. Thus the
projection D of solid Ω onto Oxz is the disk x 2 + z 2 É 4.

(HCMUT-OISP) TRIPLE INTEGRALS 17 / 80

Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

Z4 p
 
Ï
I=  x 2 + z 2d y  d A
D 2 2
Ï hxp+z i y=4
= 2 2
x + z .y dA
2 y=x +z 2
D
Ï p
2 2
= (4 − x − z ) x 2 + z 2 d A
D

where D = {(x, z) : x 2 + z 2 É 4}.

Convert to polar coordinates in xz−plane:
x = r cos ϕ, z = r sin ϕ. This gives
D = {(r, ϕ) : 0 É r É 2, 0 É ϕ É 2π}.
(HCMUT-OISP) TRIPLE INTEGRALS 18 / 80
Triple Integrals in rectangular (Cartesian) coordinate system The Triple Integral over a general bounded region

Z2π Z2
 

I =  4 − r 2 r.r d r  d ϕ
¡ ¢

0 0
2π
 2 
Z Z
¡ 2
4r − r 4 d r  d ϕ
¢
= 
0 0
¸2
4r 3 r 5
·
128π
= 2π. − = ·
3 5 0 15

(HCMUT-OISP) TRIPLE INTEGRALS 19 / 80

 Triple Integrals in Cylindrical Coordinates Cylindrical Coordinates

DEFINITION 3.1
In the cylindrical coordinate system, a point
M (x, y, z) in three-dimensional space is represented by
the ordered triple (r, ϕ, z), where r and ϕ are polar
coordinates of the projection of M onto the x y−plane
and z is the directed distance from the x y−plane to M .

(HCMUT-OISP) TRIPLE INTEGRALS 20 / 80

 Triple Integrals in Cylindrical Coordinates Cylindrical Coordinates

(HCMUT-OISP) TRIPLE INTEGRALS 21 / 80

 Triple Integrals in Cylindrical Coordinates Cylindrical Coordinates

TO CONVERT FROM CYLINDRICAL TO RECTANGULAR

COORDINATES , WE USE THE EQUATIONS

 x = r. cos ϕ

y = r. sin ϕ

 z = z

(HCMUT-OISP) TRIPLE INTEGRALS 22 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

Ω = {(r, ϕ, z) : α É ϕ É β, r 1 (ϕ) É r É r 2 (ϕ),

z 1 (x, y) É z É z 2 (x, y)}

(HCMUT-OISP) TRIPLE INTEGRALS 23 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

Volume element in cylindrical coordinates

dV = r d r d ϕ.d z.
(HCMUT-OISP) TRIPLE INTEGRALS 24 / 80
 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

Evaluating the triple integral with cylindrical

coordinates, we have to determine:
1
Lower boundary z = z 1 (r, ϕ)
2
Upper boundary z = z 2 (r, ϕ)
α Éϕ Éβ
½
3
Projection D :
r 1 (ϕ) É r É r 2 (ϕ)
4
Iterated integral is a triple integral
Z β Z r 2 (ϕ) Z z 2 (r,ϕ)
I= dϕ dr f (r cos ϕ, r sin ϕ, z)r d z
α r 1 (ϕ) z 1 (r,ϕ)

(HCMUT-OISP) TRIPLE INTEGRALS 25 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

EXAMPLE 3.1
Evaluate the triple integral with cylindrical
coordinates Ñ
I= 2xdV,
Ω

where Ω = {(x, y, z) : 0 É y É 2; 0 É x É
p
4 − y 2 ; 0 É z É x}

(HCMUT-OISP) TRIPLE INTEGRALS 26 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

Ω = {(x, y, z) : 0 É y É 2; 0 É x É
p
4 − y 2 ; 0 É z É x}

(HCMUT-OISP) TRIPLE INTEGRALS 27 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

1
The lower surface of Ω is z = 0.
2
The upper surface of Ω is z = x.
3
The projection of Ω onto x y−plane is
q
0 É y É 2; 0 É x É 4 − y 2 .

Therefore
 x 
Ï Z Ï h iz=x Ï
I=  2xd z  d A = 2xz dA= 2x 2 d A
z=0
D 0 D D
n p o
where D = (x, y) : 0 É y É 2; 0 É x É 4 − y2 .

(HCMUT-OISP) TRIPLE INTEGRALS 28 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

Convert to polar coordinates in x y−plane:

x = r cos ϕ, y = r sin ϕ, we get
n πo
D = (r, ϕ) : 0 É r É 2, 0 É ϕ É
2

Zπ/2 Z2
 

I =  2r 2 cos2 ϕr d r  d ϕ
0 0
π/2 ¸r =2
2r 4
Z ·
2
= cos ϕ. dϕ
4 r =0
0
Zπ/2 Zπ/2
= 8 cos2 ϕd ϕ = 4(1 + cos(2ϕ))d ϕ
0 0
(HCMUT-OISP)
· ¸π/2
TRIPLE INTEGRALS 29 / 80
 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

EXAMPLE 3.2
Evaluate the triple integral with cylindrical
coordinates Ñ
I= zdV,
Ω

where Ω is the solid that lies within the cylinder

x 2 + y 2 = 1, above the cone z = x 2 + y 2 and below the
cone z = 2x 2 + y 2 .

(HCMUT-OISP) TRIPLE INTEGRALS 30 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

(HCMUT-OISP) TRIPLE INTEGRALS 31 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

Ω is the solid bounded by

1
upper surface z = 2x 2 + y 2 ,
2
lower surface z = x 2 + y 2
3
projection of Ω onto x y−plane: x 2 + y 2 = 1.
Therefore,
 
Ï 2xZ2 +y 2

I= zd z  d A
 

D x 2 +y 2

(2x + y 2 )2 (x 2 + y 2 )2
2
Ï · ¸
= − dA
2 2
D

(HCMUT-OISP) TRIPLE INTEGRALS 32 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

Z2π Z1 2
r cos2 ϕ(r 2 cos2 ϕ + 2r 2 )
I = dϕ rdr
2
0 0
Z2π· ¸
1 2 2 11π
= cos ϕ(cos ϕ + 2) d ϕ = ·
12 48
0

(HCMUT-OISP) TRIPLE INTEGRALS 33 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

EXAMPLE 3.3
Evaluate the triple integral with cylindrical
coordinates
1
Ñ
I= dV,
x2 + y 2
Ω

where Ω is the solid bounded by z = 0, z = y, x 2 + y 2 = 1

and y Ê 0.

(HCMUT-OISP) TRIPLE INTEGRALS 34 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

(HCMUT-OISP) TRIPLE INTEGRALS 35 / 80

 Triple Integrals in Cylindrical Coordinates Evaluating Triple Integrals with Cylindrical Coordinates

The solid Ω is bounded by

1
upper surface z = y,
2
lower surface z = 0
3
projection of Ω onto x y−plane: x 2 + y 2 = 1, y Ê 0.
 y 
1
Ï Z
I=  d z d A
x2 + y 2
D 0
¸z=y
z y
Ï · Ï
= d A = dA
x 2 + y 2 z=0 x2 + y 2
D D
Z π Z π
r sin ϕ
Z 1
= dϕ rdr = sin ϕd ϕ = 2.
0 0 r2 0

(HCMUT-OISP) TRIPLE INTEGRALS 36 / 80

 Triple Integrals in Spherical Coordinates Spherical Coordinates

DEFINITION 4.1
In the spherical coordinate system, a point M (x, y, z)
in three-dimensional space is represented by the
ordered triple (θ, ϕ, ρ), where ρ is the distance from the
origin to M , ϕ is the same angle as in cylindrical
coordinates, and θ is the angle between the positive
z−axis and the line segment OM .

(HCMUT-OISP) TRIPLE INTEGRALS 37 / 80

 Triple Integrals in Spherical Coordinates Spherical Coordinates

The relationship between rectangular and spherical

coordinates 
 x = ρ. sin θ cos ϕ

y = ρ. sin θ sin ϕ
 z = ρ. cos θ


where 0 É θ É π, 0 É ϕ É 2π, ρ > 0.

(HCMUT-OISP) TRIPLE INTEGRALS 38 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

In the spherical coordinate system, the counterpart

of a rectangular box is a spherical wedge

Ω = {(ρ, ϕ, θ) : ρ 1 É ρ É ρ 2 , α É ϕ É β, c É θ É d },

where ρ 1 Ê 0, β − α É 2π, d − c É π.

(HCMUT-OISP) TRIPLE INTEGRALS 39 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

The Volume element dV = ρ 2 sin θd ρd ϕd θ.

(HCMUT-OISP) TRIPLE INTEGRALS 40 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

In order to evaluate the triple integral with spherical

coordinates, we have to determine:
 solid Ω in spherical coordinates
The
1


 θ1 É θ É θ2 , (0 É θ É π)
ϕ1 É ϕ É ϕ2 , (0 É ϕ É 2π)
 ρ (ϕ, θ) É ρ É ρ (ϕ, θ), (ρ Ê 0)

1 2
2
The formula for triple integration
Zθ2 Zϕ2 ρ 2Z(ϕ,θ)

I= dθ dϕ f (ρ cos ϕ sin θ, ρ sin ϕ sin θ, ρ cos θ)ρ 2 sin θd ρ

θ1 ϕ1 ρ 1 (ϕ,θ)

(HCMUT-OISP) TRIPLE INTEGRALS 41 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

3
We convert a triple integral from
rectangular coordinates to spherical
coordinates

 x = ρ. sin θ cos ϕ

y = ρ. sin θ sin ϕ
 z = ρ. cos θ


(
x 2 + y 2 + z 2 = ρ2
⇒
x 2 + y 2 = ρ. sin θ.
p

(HCMUT-OISP) TRIPLE INTEGRALS 42 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

4
Usually, spherical coordinates are used
in triple integrals when surfaces such
cones and spheres form the boundary of
the region of integration.

(HCMUT-OISP) TRIPLE INTEGRALS 43 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

EXAMPLE 4.1
Evaluate the triple integral with spherical coordinates
Ñ
I= (x + z)dV,
Ω

where Ω is the unit ball x 2 + y 2 + z 2 É 1.

(HCMUT-OISP) TRIPLE INTEGRALS 44 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

 solid Ω in spherical coordinate is bounded by

The

 0ÉθÉπ
0 É ϕ É 2π
 0 É ρ É 1.


(HCMUT-OISP) TRIPLE INTEGRALS 45 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

Zπ Z2π Z1
I= d θ d ϕ (ρ cos ϕ sin θ + ρ cos θ)ρ 2 sin θd ρ
0 0 0
Zπ Z2π
1
= dθ (cos ϕ sin θ + cos θ) sin θd ϕ
4
0 0
Zπ
1 ¡ 2 ¢ ¯¯ϕ=2π
= sin θ. sin ϕ + sin θ cos θ.ϕ ¯ dθ
4 ϕ=0
0
1 ¯π
= · π (− cos 2θ) ¯ = 0.
¯
8 0

(HCMUT-OISP) TRIPLE INTEGRALS 46 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

EXAMPLE 4.2
Evaluate the triple integral with spherical coordinates
Ñ
I= (x 2 + y 2 )dV,
Ω

where Ω is the solid hemisphere x 2 + y 2 + z 2 É 1, z Ê 0.

(HCMUT-OISP) TRIPLE INTEGRALS 47 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

The solid Ω in spherical coordinate is bounded by

π

 0 Éθ É


2
 0 É ϕ É 2π
 0 É ρ É 1.


(HCMUT-OISP) TRIPLE INTEGRALS 48 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

Zπ/2 Z2π Z1
I = d θ d ϕ [(ρ cos ϕ sin θ)2 + (ρ sin ϕ sin θ)2 ]ρ 2 sin θd ρ
0 0 0
Zπ/2
1 4π
= · 2π. sin3 θd θ = ·
5 15
0

(HCMUT-OISP) TRIPLE INTEGRALS 49 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

EXAMPLE 4.3
Evaluate the triple integral with spherical coordinates
Ñ
I= (x 2 + y 2 + z 2 )dV,
Ω

where Ω lies between the spheres x 2 + y 2 + z 2 = 1 and

x 2 + y 2 + z 2 = 4 in the first octant.

(HCMUT-OISP) TRIPLE INTEGRALS 50 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

(HCMUT-OISP) TRIPLE INTEGRALS 51 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

 solid Ω in spherical coordinate is bounded by

The
π

 0 É θ É
2

π

 0ÉϕÉ 2

 1 É ρ É 2.


x 2 +y 2 +z 2 = (ρ cos ϕ sin θ)2 +(ρ sin ϕ sin θ)2 +(ρ cos θ)2 = ρ 2
Zπ/2 Zπ/2 Z2 Zπ/2
31 π 31π
I = d θ d ϕ ρ 2 ρ 2 sin θd ρ = · · sin θd θ = ·
5 2 10
0 0 1 0

(HCMUT-OISP) TRIPLE INTEGRALS 52 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

EXAMPLE 4.4
Evaluate the triple integral with spherical coordinates
1
Ñ
I= dV,
x2 + y 2 + z2
Ω

where Ω is the solid that lies within the

psphere
2 2 2
x + y + z = 4 and above the cone z = x 2 + y 2

(HCMUT-OISP) TRIPLE INTEGRALS 53 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

n o
Ω = (x, y, z) : x + y + z É 4, z Ê x + y
2 2 2
p
2 2

(HCMUT-OISP) TRIPLE INTEGRALS 54 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

n o
Ω = (x, y, z) : x + y + z É 4, z Ê x + y
2 2 2
p
2 2

(HCMUT-OISP) TRIPLE INTEGRALS 55 / 80

 Triple Integrals in Spherical Coordinates Evaluating Triple Integrals with Spherical Coordinates

 solid Ω in spherical coordinate is bounded by

The
π
 0ÉθÉ


4
 0 É ϕ É 2π
 0 É ρ É 2.


Zπ/4 Z2π Z2
1 2
I = dθ dϕ ρ sin θd ρ =
ρ2
0 0 0

Zπ/4 p
= 2.2π. sin θd θ = (4 − 2 2)π.
0

(HCMUT-OISP) TRIPLE INTEGRALS 56 / 80

 Applications of Triple Integrals Evaluating the volume of a solid

THEOREM 5.1
If f (x, y, z) = 1 for all points in Ω, then the triple
integral represent the volume of Ω
Ñ
V= 1.dV (1)
Ω

(HCMUT-OISP) TRIPLE INTEGRALS 57 / 80

 Applications of Triple Integrals Evaluating the volume of a solid

EXAMPLE 5.1
Use spherical coordinates to findpthe volume of the
solid that lies above the cone z = x 2 + y 2 and below
the sphere x 2 + y 2 + z 2 = z.

(HCMUT-OISP) TRIPLE INTEGRALS 58 / 80

 Applications of Triple Integrals Evaluating the volume of a solid

n o
Ω = (x, y, z) : x + y + z É z, z Ê x + y
2 2 2
p
2 2

(HCMUT-OISP) TRIPLE INTEGRALS 59 / 80

 Applications of Triple Integrals Evaluating the volume of a solid

Substituting x = ρ. cos ϕ sin θ, y = ρ. sin ϕ sin θ,

z = ρ. cos θ into the equation x 2 + y 2 p
+ z 2 É z, we get
ρ É cos θ, and into the equation 2 2
h π i z Ê x + y , we get
ρ cos θ Ê ρ. sin θ, where θ ∈ 0, , we have tan θ É 1, or
2
π
0ÉθÉ
4
The solid Ω in spherical coordinate is bounded by
π

 0 Éθ É


4
 0 É ϕ É 2π
 0 É ρ É cos θ.


(HCMUT-OISP) TRIPLE INTEGRALS 60 / 80

 Applications of Triple Integrals Evaluating the volume of a solid

Zπ/4 Z2π Z θ
cos

V = dθ dϕ ρ 2 sin θd ρ
0 0 0
Z2π Zπ/4 ¸ρ=cos θ
ρ3
·
= d ϕ. sin θ. dθ
3 ρ=0
0 0
Zπ/4 ¸π/4
cos4 θ π
·
2π 2π
3
= · sin θ cos θd θ = − = ·
3 3 4 0 8
0

(HCMUT-OISP) TRIPLE INTEGRALS 61 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

SKETCH THE SOLID IN THE CARTESIAN COORDINATE

SYSTEM

[X , Y ] = meshg r i d (x, y)
x = [−1 0 1]; y = [−2 − 1 0 1 2];
[X , Y ] = meshg r i d (x, y);
−1 0 1 −2 −2 −2
−1 0 1 −1 −1 −1
X = −1 0 1 , Y = 0 0 0
−1 0 1 1 1 1
−1 0 1 2 2 2

(HCMUT-OISP) TRIPLE INTEGRALS 62 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

(HCMUT-OISP) TRIPLE INTEGRALS 63 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

x = [−1 0 1];
[X ] = meshg r i d (x);
−1 0 1
X = −1 0 1
−1 0 1

(HCMUT-OISP) TRIPLE INTEGRALS 64 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

THE SOLID IS THE BOX

1
Use [X , Y ] = meshg r i d (x, y) to receive projection
onto x y−plane
2
Sketch the upper and lower surfaces.
EXAMPLE 6.1
Sketch the solid bounded by
z = 0, z = 3, −1 É x É 1, −2 É y É 2

(HCMUT-OISP) TRIPLE INTEGRALS 65 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

Code
x = l i nspace(−1, 1, 50); y = l i nspace(−2, 2, 50);
z = l i nspace(0, 3, 50); [X , Y ] = meshg r i d (x, y);
sur f (X , Y , 3 + 0 ∗ X ) (or mesh(X , Y , Z ) or
sur f c(X , Y , Z ))
hold on
mesh(X , Y , 0 ∗ X )
[X , Z ] = meshg r i d (x, z);
sur f (X , 2 + 0 ∗ X , Z ,
′
F aceC ol or ′ ,′ g ′ ,′ E d g eC ol or ′ ,′ b ′ ,′ F ace Al pha ′ , 0.1);
sur f (X , −2 + 0 ∗ X , Z ,
′
F aceC ol or ′ ,′ g ′ ,′ E d g eC ol or ′ ,′ g ′ ,′ F ace Al pha ′ , 0.1);

(HCMUT-OISP) TRIPLE INTEGRALS 66 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

[Y , Z ] = meshg r i d (y, z);

sur f (−1 + 0 ∗ Y , Y , Z ,
′
F aceC ol or ′ ,′ b ′ ,′ E d g eC ol or ′ ,′ g ′ ,′ F ace Al pha ′ , 0.1);
sur f (1 + 0 ∗ Y , Y , Z ,
′
F aceC ol or ′ ,′ b ′ ,′ E d g eC ol or ′ ,′ g ′ ,′ F ace Al pha ′ , 0.1);
xl abel (′ x ′ ); yl abel (′ y ′ ); zl abel (′ z ′ );
axi s([−1 1 − 2 2 0 8])
view(120,12)
grid on
rotate3d on

(HCMUT-OISP) TRIPLE INTEGRALS 67 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

(HCMUT-OISP) TRIPLE INTEGRALS 68 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

THE SOLID LIES BELOW THE SURFACE z = z 2 , ABOVE THE

SURFACE z = z 1 WITH ITS PROJECTION D

EXAMPLE 6.2
Sketch the solid bounded by
z = 0, z = 1 − y, −1 É x É 1, x 2 É y É 1

Code
cl f
s = l i nspace(−1, 1, 30);
s1 = meshg r i d (s); t 1 = [];

(HCMUT-OISP) TRIPLE INTEGRALS 69 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

for i=1:length(s)
t am = l i nspace(s(i )b2, 1, 30);
t 1 = [t 1 t am ′ ];
end
x = s1; y = t 1; z = 1 − y; z1 = 0 ∗ x;
hold on
sur f (x, y, z,′ F aceC ol or ′ ,′ g ′ ,′ F ace Al pha ′ , 0.3);
sur f (x, y, z1,′ F aceC ol or ′ ,′ r ′ ,′ E d g eC ol or ′ ,′ none ′ );

(HCMUT-OISP) TRIPLE INTEGRALS 70 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

s = l i nspace(−1, 1, 30); t 1 = [];

for i=1:length(s)
t am = l i nspace(0, 1 − s(i )b2, 30);
t 1 = [t 1 t am ′ ];
end
x = s1; z = t 1; y = s1.b2;
sur f (x, y, z,′ F aceC ol or ′ ,′ b ′ ,′ F ace Al pha ′ , 0.5);
xl abel (′ x ′ ); yl abel (′ y ′ ); zl abel (′ z ′ );
vi ew(120, 12)
grid on
rotate3d on

(HCMUT-OISP) TRIPLE INTEGRALS 71 / 80

 MatLab Sketch the solid in the Cartesian Coordinate System

(HCMUT-OISP) TRIPLE INTEGRALS 72 / 80

 MatLab Sketch the solid in the Cylindrical Coordinate System

EXAMPLE 6.3
Sketch the solid bounded by
z = 0, z = 4 − x 2 − y 2 , x 2 + y 2 = 2

Code
clf
hold on
phi=linspace(0,2*pi,30); r=linspace(0,sqrt(2),30);
[r phi ] = meshg r i d (r, phi );
x=r.*cos(phi); y=r.*sin(phi);

(HCMUT-OISP) TRIPLE INTEGRALS 73 / 80

 MatLab Sketch the solid in the Cylindrical Coordinate System

z = 4 − x.b2 − y.b2; z1 = 0 ∗ x;
surf(x,y,z,’FaceColor’,’g’,’FaceAlpha’,0.3);
surf(x,y,z1,’FaceColor’,’b’,’FaceAlpha’,0.3);
phi=linspace(0,2*pi,30); z2=linspace(0,2,30);
[z2 phi ] = meshg r i d (z2, phi );
x1=sqrt(2).*cos(phi); y1=sqrt(2).*sin(phi);
surf(x1,y1,z2,’FaceColor’,’r’,’FaceAlpha’,0.3);
grid on; view (13,28); rotate3d on
xlabel(’x’); ylabel(’y’); zlabel(’z’);

(HCMUT-OISP) TRIPLE INTEGRALS 74 / 80

 MatLab Sketch the solid in the Cylindrical Coordinate System

(HCMUT-OISP) TRIPLE INTEGRALS 75 / 80

 MatLab Sketch the solid in the Spherical Coordinate System

EXAMPLE 6.4
Sketch the solid that lies between the spheres
x 2 + y 2 + z 2 = 1 and x 2 + y 2 + z 2 = 4 in the first octant.

Code
clf
phi=linspace(0,pi/2,30); theta=linspace(0,pi/2,30);
[p t]=meshgrid(phi,theta);
x=sin(t).*cos(p); y=sin(t).*sin(p); z=cos(t);
mesh(x,y,z,’FaceColor’,’r’,’FaceAlpha’,0.3);

(HCMUT-OISP) TRIPLE INTEGRALS 76 / 80

 MatLab Sketch the solid in the Spherical Coordinate System

hold on
x2=2*sin(t).*cos(p);
y2=2*sin(t).*sin(p);
z2=2*cos(t);
mesh(x2,y2,z2,’FaceColor’,’g’,’FaceAlpha’,0.3);
r=linspace(1,2,30)
[r th]=meshgrid(r,theta);
x3=sin(th).*r.*cos(0);
y3=sin(th).*r.*sin(0);
z3=r.*cos(th);
mesh(x3,y3,z3,’FaceColor’,’b’,’FaceAlpha’,0.3);

(HCMUT-OISP) TRIPLE INTEGRALS 77 / 80

 MatLab Sketch the solid in the Spherical Coordinate System

x4=sin(th).*r.*cos(pi/2);
y4=sin(th).*r.*sin(pi/2);
z4=r.*cos(th);
mesh(x4,y4,z4,’FaceColor’,’b’,’FaceAlpha’,0.3);
[r phi]=meshgrid(r,phi);
x5=sin(pi/2).*r.*cos(phi);
y5=sin(pi/2).*r.*sin(phi);
z5=r.*cos(pi/2);
mesh(x5,y5,z5,’FaceColor’,’r’,’FaceAlpha’,0.3);
rotate3d on
xlabel(’x’); ylabel(’y’); zlabel(’z’);

(HCMUT-OISP) TRIPLE INTEGRALS 78 / 80

 MatLab Sketch the solid in the Spherical Coordinate System

(HCMUT-OISP) TRIPLE INTEGRALS 79 / 80

 MatLab Sketch the solid in the Spherical Coordinate System

THANK YOU FOR YOUR ATTENTION

(HCMUT-OISP) TRIPLE INTEGRALS 80 / 80