Active Calculus Multivariable - 2018 Edition
Active Calculus Multivariable - 2018 Edition
Active Calculus Multivariable - 2018 Edition
MULTIVARIABLE
2018 Edition – Updated
Steven Schlicker
David Austin Matthew Boelkins
Active Calculus - Multivariable
Active Calculus - Multivariable
Steve Schlicker
Grand Valley State University
Contributing Authors
David Austin
Grand Valley State University
Matt Boelkins
Grand Valley State University
Edition: 2018-updated
Website: https://activecalculus.org/
©2013–2019 Steven Schlicker
Permission is granted to copy, distribute and/or modify this document under
the terms of the Creative Commons Attribution-NonCommercial-ShareAlike
4.0 International License. The work may be used for free by any party so long
as attribution is given to the author(s), the work and its derivatives are used
in the spirit of “share and share alike”; no party may sell this work or any of
its derivatives for profit. All trademarks™ are the registered® marks of their
respective owners. The graphic
that may appear in other locations in the text shows that the work is li-
censed with the Creative Commons, that the work may be used for free by
any party so long as attribution is given to the author(s), that the work
and its derivatives are used in the spirit of “share and share alike,” and
that no party may sell this work or any of its derivatives for profit, with
the following exception: it is entirely acceptable for university bookstores to
sell bound photocopied copies of the activities workbook to students at their
standard markup above the copying expense. Full details may be found by
visiting creativecommons.org/licenses/by-nc-sa/4.0/ or sending a letter to
Creative Commons, 444 Castro Street, Suite 900, Mountain View, California,
94041, USA.
Features of the Text
Activities. Every section in the text contains several activities. These are
designed to engage students in an inquiry-based style that encourages
them to construct solutions to key examples on their own, working either
individually or in small groups.
Exercises. There are dozens of calculus texts with (collectively) tens of thou-
sands of exercises. Rather than repeat a large list of standard and rou-
tine exercises in this text, we recommend the use of WeBWorK with its
access to the National Problem Library and its many multivariable cal-
culus problems. In this text, each section begins with several anonymous
WeBWorK exercises, and follows with several challenging problems. The
WeBWorK exercises are best completed in the .html version of the
text. Almost every non-WeBWorK problem has multiple parts, requires
the student to connect several key ideas, and expects that the student
will do at least a modest amount of writing to answer the questions and
explain their findings. For instructors interested in a more conventional
source of exercises, consider the freely available APEX Calculus text by
Greg Hartmann et al.
v
vi
Summary of Key Ideas. Each section concludes with a summary of the key
ideas encountered in the preceding section; this summary normally re-
flects responses to the motivating questions that began the section.
Links to technological tools. Many of the ideas of multivariable calculus
are best understood dynamically, and we encourage readers to make
frequent use of technology to analyze graphs and data. Since tech-
nology changes so often, we refrain from indicating specific programs
to use in the text. However, aside from computer algebra systems like
Maple, Mathematics, or Sage, there are many free graphing tools avail-
able for drawing three-dimensional surfaces or curves. These programs
can be used by instructors and students to assist in the investigations and
demonstrations. The use of these freely available applets is in accord with
our philosophy that no one should be required to purchase materials to
learn calculus. We are indebted to everyone who allows their expertise
to be openly shared. Below is a list of a few of the technological tools
that are available (links active at the writing of this edition). Of course,
you can find your own by searching the web.
• Wolfram Alpha, useful for graphing surfaces in 2D and 3D, and for general
calculations
• Wolfram Alpha widgets, searchable site for simple to use programs using
Wolfram Alpha
• GeoGebra, all purpose graphing tool with some computer algebra capabil-
ities. Clicking on the magnifying glass icon allows you to search a large
database of GeoGebra applets.
• CalcPlot3D, good all-purpose 3D graphing tool
• A collection of Flash Mathlets for graphing surfaces, parametric
curves in 3D, spherical coordinates and other 3D tools. Requires Flash
Player.
Acknowledgments
This text is an extension of the single variable Active Calculus by Matt Boelkins.
The initial drafts of this multivariable edition were written by me; editing and
revisions were made by David Austin and Matt Boelkins. David Austin is
responsible for the beautiful full-color graphics in the text. Many of our col-
leagues at GVSU have shared their ideas and resources, which undoubtedly
had a significant influence on the product. We thank them for all of their
support. Most importantly, I want to thank thank the students who have used
this text and offered helpful advice and suggestions.
In advance, we also thank our colleagues throughout the mathematical
community who have or will read, edit, and use this book, and hence contribute
to its improvement through ongoing discussion. The following people have
used early drafts of this text and have generously offered suggestions that have
improved the text.
Feryâl Alayont Grand Valley State University
David Austin Grand Valley State University
Jon Barker and students St. Ignatius High School, Cleveland, OH
Matt Boelkins Grand Valley State University
Brian Drake Grand Valley State University
Brian Gleason Nevada State College
Mitch Keller Morningside College
The current .html version of the text is possible only because of the amazing
work of Rob Beezer and his development of PreTeXt. My ability to take
advantage of Rob’s work is largely due to the support of the American Institute
of Mathematics, which funded me for a weeklong workshop in Mathbook XML
in San Jose, CA, in April 2016. David Farmer also deserves credit for the
original conversion of the text from LATEX to PreTeXt.
I take full responsibility for all errors or inconsistencies in the text, and
welcome reader and user feedback to correct them, along with other suggestions
to improve the text.
David Austin, Matt Boelkins, Steven Schlicker, Allendale, MI, August,
2018.
vii
Active Calculus - Multivariable:
our goals
Several fundamental ideas in calculus are more than 2000 years old. As a for-
mal subdiscipline of mathematics, calculus was first introduced and developed
in the late 1600s, with key independent contributions from Sir Isaac Newton
and Gottfried Wilhelm Leibniz. Mathematicians agree that the subject has
been understood rigorously since the work of Augustin Louis Cauchy and Karl
Weierstrass in the mid 1800s when the field of modern analysis was devel-
oped, in part to make sense of the infinitely small quantities on which calculus
rests. As a body of knowledge, calculus has been completely understood for at
least 150 years. The discipline is one of our great human intellectual achieve-
ments: among many spectacular ideas, calculus models how objects fall under
the forces of gravity and wind resistance, explains how to compute areas and
volumes of interesting shapes, enables us to work rigorously with infinitely
small and infinitely large quantities, and connects the varying rates at which
quantities change to the total change in the quantities themselves.
While each author of a calculus textbook certainly offers their own creative
perspective on the subject, it is hardly the case that many of the ideas an
author presents are new. Indeed, the mathematics community broadly agrees
on what the main ideas of calculus are, as well as their justification and their
importance; the core parts of nearly all calculus textbooks are very similar.
As such, it is our opinion that in the 21st century—an age where the internet
permits seamless and immediate transmission of information—no one should
be required to purchase a calculus text to read, to use for a class, or to find a
coherent collection of problems to solve. Calculus belongs to humankind, not
any individual author or publishing company. Thus, a main purpose of this
work is to present a new multivariable calculus text that is free. In addition,
instructors who are looking for a calculus text should have the opportunity to
download the source files and make modifications that they see fit; thus this
text is open-source.
In Active Calculus - Multivariable, we endeavor to actively engage students
in learning the subject through an activity-driven approach in which the vast
majority of the examples are completed by students. Where many texts present
a general theory of calculus followed by substantial collections of worked ex-
amples, we instead pose problems or situations, consider possibilities, and then
ask students to investigate and explore. Following key activities or examples,
the presentation normally includes some overall perspective and a brief syn-
opsis of general trends or properties, followed by formal statements of rules or
theorems. While we often offer plausibility arguments for such results, rarely
do we include formal proofs. It is not the intent of this text for the instructor
or author to demonstrate to students that the ideas of calculus are coherent
viii
ix
and true, but rather for students to encounter these ideas in a supportive,
leading manner that enables them to begin to understand for themselves why
calculus is both coherent and true.
This approach is consistent with the following goals:
• To have students engage in an active, inquiry-driven approach, where
learners strive to construct solutions and approaches to ideas on their
own, with appropriate support through questions posed, hints, and guid-
ance from the instructor and text.
• To build in students intuition for why the main ideas in multivariable
calculus are natural and true. We strive to accomplish this by using spe-
cific cases to highlight the ideas for the general situation using contexts
that are common and familiar.
• To challenge students to acquire deep, personal understanding of multi-
variable calculus through reading the text and completing preview activ-
ities on their own, through working on activities in small groups in class,
and through doing substantial exercises outside of class time.
• To strengthen students’ written and oral communication skills by hav-
ing them write about and explain aloud the key ideas of multivariable
calculus.
How to Use this Text
Because the text is free, any professor or student may use the electronic version
of the text for no charge. For reading on laptops or mobile devices, the best
electronic version to use is the .html version of the text, but you can find
links to a pdf and hard copy of the text at https://activecalculus.org/.
Furthermore, because the text is open-source, any instructor may acquire the
full set of source files, which are available on GitHub.
This text may be used as a stand-alone textbook for a standard multivari-
able calculus course or as a supplement to a more traditional text. Chapter 9
introduces functions of several independent variables along with tools that will
be used to study these functions, namely vectors and vector-valued functions.
Chapter 10 studies differentiation of functions of several independent variables
in detail, addressing the typical topics including limits, partial derivatives, and
optimization, while Chapter 11 provides the standard topics of integration of
multivariable functions.
Electronic Edition Because students and instructors alike have access to the
book in electronic format, there are several advantages to the text over a
traditional print text. One is that the text may be projected on a screen
in the classroom (or even better, on a whiteboard) and the instructor may
reference ideas in the text directly, add comments or notation or features
to graphs, and indeed write right on the projected text itself. Students
can do the same when working at the board. In addition, students can
choose to print only whatever portions of the text are needed for them.
Also, the electronic versions of the text includes live .html links to on-
line programs, so student and instructor alike may follow those links to
additional resources that lie outside the text itself. Finally, students can
have access to a copy of the text anywhere they have a computer. The
.html version is far superior to the .pdf version; this is especially true for
viewing on a smartphone.
Note. In the .pdf version, there is not an obvious visual indicator of the
live .html links, so some availalable information is suppressed. If you are
using the text electronically in a setting with internet access, please know
that it is assumed you are using the .html version.
Activities Workbook Each section of the text has a preview activity and
at least three in-class activities embedded in the discussion. As it is the
expectation that students will complete all of these activities, it is ideal
for them to have room to work on them adjacent to the problem state-
ments themselves. A separate workbook of activities that includes only
the individual activity prompts, along with space provided for students
to write their responses, is in development.
x
xi
Acknowledgments vii
xii
CONTENTS xiii
Index 294
Chapter 9
1
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 2
a. Find the amount of money in the account after 7 years if you originally
invest 1000 dollars.
d. Now consider the situation where we want to know the amount of money
in the account after 10 years given various initial investments. Calculate
the amount of money in the account as indicated in Table 9.1.2. Round
payments to the nearest penny.
Table 9.1.2 Amount of money in an account after 10 years.
Initial investment (dollars) 500 1000 5000 7500 10000
Amount (dollars)
e. Describe as best you can the combinations of initial investments and time
that result in an account containing $10,000.
Activity 9.1.2 Identify the domain of each of the following functions. Draw
a picture of each domain in the xy-plane.
a. f (x, y) = x2 + y 2
p
b. f (x, y) = x2 + y 2
x+y
c. Q(x, y) = x2 −y 2
1
d. s(x, y) = √
1−xy 2
x2 sin(2y)
f (x, y) = ,
g
where g is the acceleration due to gravity. (Note that g is constant, 32 feet per
second squared. We will derive this equation in a later section.) To create a
table of values for f , we list the x-values down the first column and the y-values
across the first row. The value f (x, y) is then displayed in the location where
the x row intersects the y column, as shown in Table 9.1.5 (where we measure
x in feet per second and y in radians).
x2 sin(2y)
Table 9.1.5 Values of f (x, y) = g .
Activity 9.1.3 Complete Table 9.1.5 by filling in the missing values of the
function f . Round entries to the nearest tenth.
If f is a function of a single variable x, then we define the graph of f to be
the set of points of the form (x, f (x)), where x is in the domain of f . We then
plot these points using the coordinate axes in order to visualize the graph. We
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 4
can do a similar thing with functions of several variables. Table 9.1.5 identifies
points of the form (x, y, f (x, y)), and we define the graph of f to be the set of
these points.
Definition 9.1.6 The graph of a function f = f (x, y) is the set of points of
the form (x, y, f (x, y)), where the point (x, y) is in the domain of f . ♦
We also often refer to the graph of a function f of two variables as the sur-
face generated by f . Points in the form (x, y, f (x, y)) are in three dimensions,
so plotting these points takes a bit more work than graphs of functions in two
dimensions. To plot these three-dimensional points, we need to set up a coor-
dinate system with three mutually perpendicular axes — the x-axis, the y-axis,
and the z-axis (called the coordinate axes). There are essentially two different
ways we could set up a 3D coordinate system, as shown in Figure 9.1.7; thus,
before we can proceed, we need to establish a convention.
z z
x y
y x
Figure 9.1.7 Left: A left hand system. Right: A right hand system
The distinction between these two figures is subtle, but important. In
the coordinate system shown at left in Figure 9.1.7, imagine that you are
sitting on the positive z-axis next to the label “z.” Looking down at the x-
and y-axes, you see that the y-axis is obtained by rotating the x-axis by 90◦
in the counterclockwise direction. Again sitting on the positive z-axis in the
coordinate system at right in Figure 9.1.7, you see that the y-axis is obtained
by rotating the x-axis by 90◦ in the clockwise direction.
We call the coordinate system at right in Figure 9.1.7 a right-hand system;
if we point the index finger of our right hand along the positive x-axis and our
middle finger along the positive y-axis, then our thumb points in the direction
of the positive z-axis. Following mathematical conventions, we choose to use
a right-hand system throughout this book.
Now that we have established a convention for a right-hand system, we can
2
draw a graph of the distance function defined by f (x, y) = x sin(2y)
g . Note that
the function f is continuous in both variables, so when we plot these points in
the right hand coordinate system, we can connect them all to form a surface
in 3-space. The graph of the distance function f is shown in Figure 9.1.8.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 5
z
1500
1000
500
x
200
y 150
1.5 100
1.0 50
0.5
0
c. Consider the set of points (x, y, z) that satisfy the equation z = 0. De-
scribe this set as best you can.
Activity 9.1.4 shows that the equations where one independent variable is
constant lead to planes parallel to ones that result from a pair of the coordinate
axes. When we make the constant 0, we get the coordinate planes. The xy-
plane satisfies z = 0, the xz-plane satisfies y = 0, and the yz-plane satisfies
x = 0 (see Figure 9.1.9).
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 6
z z z
y y y
x x x
P = (x0 , y0 , z0 ) S
b. Now notice that the triangle P RQ whose hypotenuse is the blue segment
connecting the points P and Q with a leg as the hypotenuse P R of the
triangle found in part (a) lies entirely in a plane, so we can again use the
Pythagorean Theorem to find the length of its hypotenuse. Explain why
the length of this hypotenuse, which is the distance between the points
P and Q, is p
(x1 − x0 )2 + (y1 − y0 )2 + (z1 − z0 )2 .
The formula developed in Activity 9.1.5 is important to remember.
Equation (9.1.1) can be used to derive the formula for a sphere centered at
a point (x0 , y0 , z0 ) with radius r. Since the distance from any point (x, y, z)
on such a sphere to the point (x0 , y0 , z0 ) is r, the point (x, y, z) will satisfy the
equation p
(x − x0 )2 + (y − y0 )2 + (z − z0 )2 = r
Squaring both sides, we come to the standard equation for a sphere.
(x − x0 )2 + (y − y0 )2 + (z − z0 )2 = r2 .
(x − x0 )2 + (y − y0 )2 = r2 .
9.1.4 Traces
When we study functions of several variables we are often interested in how
each individual variable affects the function in and of itself. In Preview Ac-
tivity 9.1.1, we saw that the amount of money in an account depends on the
interest rate and the duration of the investment. However, if we fix the interest
rate, the amount of money in the account depends only on the duration of the
investment, and if we set the duration of the investment constant, then the
amount of money in the account depends only on the interest rate. This idea
of keeping one variable constant while we allow the other to change will be an
important tool for us when studying functions of several variables.
As another example, consider again the distance function f defined by
x2 sin(2y)
f (x, y) =
g
where x is the initial velocity of an object in feet per second, y is the launch
angle in radians, and g is the acceleration due to gravity (32 feet per second
squared). If we hold the launch angle constant at y = 0.6 radians, we can
consider f a function of the initial velocity alone. In this case we have
x2
f (x) = sin(2 · 0.6).
32
We can plot this curve on the surface by tracing out the points on the
surface when y = 0.6, as shown at left in Figure 9.1.11. The formula clearly
shows that f is quadratic in the x-direction. More descriptively, as we increase
the launch velocity while keeping the launch angle constant, the horizontal
distance the object travels increases proportional to the square of the initial
velocity.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 8
Similarly, if we fix the initial velocity at 150 feet per second, we can consider
the distance as a function of the launch angle only. In this case we have
1502 sin(2y)
f (y) = .
32
We can again plot this curve on the surface by tracing out the points on
the surface when x = 150, as shown at right in Figure 9.1.11. The formula
clearly show that f is sinusoidal in the y-direction. More descriptively, as
we increase the launch angle while keeping the initial velocity constant, the
horizontal distance traveled by the object is proportional to the sine of twice
the launch angle.
z z
1500 1500
1000 1000
500 500
x x
200 200
y 150 y 150
1.5 100 1.5 100
1.0 50 1.0 50
0.5 0.5
0 0
Figure 9.1.11 Left: The trace with y = 0.6. Right: The trace with x = 150.
The curves we define when we fix one of the independent variables in our
two variable function are called traces.
Definition 9.1.12 A trace of a function f of two independent variables x and
y in the x direction is a curve of the form z = f (x, c), where c is a constant.
Similarly, a trace of a function f of two independent variables x and y in the
y direction is a curve of the form z = f (c, y), where c is a constant. ♦
Understanding trends in the behavior of functions of two variables can be
challenging, as can sketching their graphs; traces help us with each of these
tasks.
Activity 9.1.6 In the following questions, we investigate the use of traces to
better understand a function through both tables and graphs.
a. Identify the y = 0.6 trace for the distance function f defined by f (x, y) =
x2 sin(2y)
g by highlighting or circling the appropriate cells in Table 9.1.5.
Write a sentence to describe the behavior of the function along this trace.
b. Identify the x = 150 trace for the distance function by highlighting or
circling the appropriate cells in Table 9.1.5. Write a sentence to describe
the behavior of the function along this trace.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 9
−4
−2 y
−4 −2 2 4
x 4
c. from that same initial point, determine the least steep path that leads to
the highest point.
Curves on a surface that describe points at the same height or level are
called level curves.
Definition 9.1.15 A level curve (or contour) of a function f of two in-
dependent variables x and y is a curve of the form k = f (x, y), where k is a
constant. ♦
Topographical maps can be used to create a three-dimensional surface from
the two-dimensional contours or level curves. For example, level curves of the
2
distance function defined by f (x, y) = x sin(2y)
32 plotted in the xy-plane are
shown at left in Figure 9.1.16. If we lift these contours and plot them at their
respective heights, then we get a picture of the surface itself, as illustrated at
right in Figure 9.1.16.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 11
z
1.5 y 500
250
1.0
x
0.5 200
y 150
100
x 1.5 1.0 50
0.5
50 100 150 200 0
Figure 9.1.16 Left: Level curves. Right: Level curves at appropriate heights.
The use of level curves and traces can help us construct the graph of a
function of two variables.
Activity 9.1.8
y y
x x
2 2
Figure 9.1.17p Left: Level curves for f (x, y) = x + y . Right: Level curves
for g(x, y) = x2 + y 2 .
c. Compare and contrast the graphs of f and g. How are they alike? How
are they different? Use traces for each function to help answer these
questions.
The traces and level curves of a function of two variables are curves in space.
In order to understand these traces and level curves better, we will first spend
some time learning about vectors and vector-valued functions in the next few
sections and return to our study of functions of several variables once we have
those more mathematical tools to support their study.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 12
Figure 9.1.18 z = x2 + y 2
p
Figure 9.1.20 z = x2 + y 2
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 13
x y
Figure 9.1.21 z = x2 − y 2
y
x
Figure 9.1.23 z = y 2 − x3 + x
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 14
2
−y 2
Figure 9.1.24 z = xye−x
9.1.7 Summary
• A function f of several variables is a rule that assigns a unique number to
an ordered collection of independent inputs. The domain of a function of
several variables is the set of all inputs for which the function is defined.
• In R3 , the distance between points P = (x0 , y0 , z0 ) and Q = (x1 , y1 , z1 )
(denoted as |P Q|) is given by the formula
p
|P Q| = (x1 − x0 )2 + (y1 − y0 )2 + (z1 − z0 )2 .
and thus the equation of a sphere with center (x0 , y0 , z0 ) and radius r is
(x − x0 )2 + (y − y0 )2 + (z − z0 )2 = r2 .
9.1.8 Exercises
1. Evaluate a function. Evaluate the function at the specified points.
f (x, y) = x + yx4 , (−1, 3) , (−4, 2) , (−1, −3)
At (−1, 3):
At (−4, 2):
At (−1, −3):
Answer 1. 2
Answer 2. 508
Answer 3. −4
2. Sketch a contour diagram of each function. Then, decide whether its
contours are predominantly lines, parabolas, ellipses, or hyperbolas.
(a) z = y − 2x2
(b) z = x2 + 3y 2
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 15
(c) z = x2 − 4y 2
(d) z = −4x2
3. Match the surfaces with the verbal description of the level curves by plac-
ing the letter of the verbal description to the left of the number of the
surface.
p
(a) z = (x2 + y 2 )
(b) z = 2x + 3y
(c) z = 2x2 + 3y 2
(d) z = x2 + y 2
(e) z = xy
1
(f) z = x−1
p
(g) z = (25 − x2 − y 2 )
(b) Find the equation for the set of points whose distance from the
x-axis equals the distance from the yz-plane.
x2 + y 2 = r 2
x2 + z 2 = r 2
y 2 = x2 + z 2
y2 + z 2 = r2
z 2 = x2 + y 2
x2 = y 2 + z 2
9. For each surface, decide whether it could be a bowl, a plate, or neither.
Consider a plate to be any fairly flat surface and a bowl to be anything
that could hold water, assuming the positive z-axis is up.
(a) x + y + z = 3
(b) z = 1 − x2 − y 2
(c) z = x2 + y 2
p
(d) z = − 3 − x2 − y 2
(e) z = 1
10. Consider the concentration, C, (in mg/liter) of a drug in the blood as a
function of the amount of drug given, x, and the time since injection, t.
For 0 ≤ x ≤ 5 mg and t ≥ 0 hours, we have
C = f (x, t) = 22te−(5−x)t
f (1, 4) =
Give a practical interpretation of your answer: f (1, 4) is
the change in concentration of a 4 mg dose in the blood 1 hours
after injection.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 17
circles
ellipses
parabolas
lines
hyperbolas
b. The plane parallel to the yz-plane that passes through the point
(7, −2, −3).
c. The sphere centered at the point (2, 1, 3) and has the point (−1, 0, −1)
on its surface.
d. The sphere whose diameter has endpoints (−3, 1, −5) and (7, 9, −1).
15. The Ideal Gas Law, P V = RT , relates the pressure (P , in pascals), tem-
perature (T , in Kelvin), and volume (V , in cubic meters) of 1 mole of a
gas (R = 8.314 molJ K is the universal gas constant), and describes the
behavior of gases that do not liquefy easily, such as oxygen and hydrogen.
We can solve the ideal gas law for the volume and hence treat the volume
as a function of the pressure and temperature:
8.314T
V (P, T ) = .
P
e. Based on all your work above, write a couple of sentences that de-
scribe the effects that temperature and pressure have on volume.
16. When people buy a large ticket item like a car or a house, they often
take out a loan to make the purchase. The loan is paid back in monthly
installments until the entire amount of the loan, plus interest, is paid. The
monthly payment that the borrower has to make depends on the amount
P of money borrowed (called the principal), the duration t of the loan in
years, and the interest rate r. For example, if we borrow $18,000 to buy
a car, the monthly payment M that we need to make to pay off the loan
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 19
c. Choose 4 different values from the range of h and plot the corre-
sponding level curves in the plane. What is the shape of a typical
level curve?
d. Choose 5 different values of x (including at least one negative value
and zero), and sketch the corresponding traces of the function h.
9.2 Vectors
Motivating Questions
• What is a vector?
• What does it mean for two vectors to be equal?
• How do we add two vectors together and multiply a vector by a scalar?
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 20
The key idea in these commands is that we start at the origin, then tell
Postscript that we want to start at the point (0, 0), draw a line from the
point (0, 0) to the point (1, 1) (this is what the lineto and stroke commands
do), then draw lines from (1, 1) to (1, −1) and (1, −1) back to the origin. Each
of these commands encodes two important pieces of information: a direction
in which to move and a distance to move. Mathematically, we can capture this
information succinctly in a vector. To do so, we record the movement on the
map in a pair hx, yi (this pair hx, yi is a vector), where x is the horizontal
displacement and y the vertical displacement from one point to another. So,
for example, the vector from the origin to the point (1, 1) is represented by the
vector h1, 1i.
a. What is the vector v1 = hx, yi that describes the displacement from the
point (1, 1) to the point (1, −1)? How can we use this vector to determine
the distance from the point (1, 1) to the point (1, −1)?
b. Suppose we want to draw the triangle with vertices A = (2, 3), B =
(−3, 1), and C = (4, −2). As a shorthand notation, we will denote the
−−→
vector from the point A to the point B as AB
−−→ −−→ −→
i. Determine the vectors AB, BC, and AC.
−−→ −−→ −→
ii. What relationship do you see among the vectors AB, BC, and AC?
Explain why this relationship should hold.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 21
y y
R(4, 5)
5
4
1
x Q(1, 1) x
3 1 4
Figure 9.2.3 Left: Standard position. Right: A vector between two points.
Alternatively, we may place the tail of the vector h3, 4i at another point,
such as Q(1, 1). After a displacement of three units to the right and four units
up, the tip of the vector is at the point R(4, 5) (see the vector at right in
Figure 9.2.3).
In this example, the vector led to the directed line segment from Q to R,
−−→
which we denote as QR. We may also turn the situation around: given the two
points Q and R, we obtain the vector h3, 4i because we move horizontally three
−−→
units and vertically four units to get from Q to R. In other words, QR = h3, 4i.
−−→
In general, the vector QR from the point Q = (q1 , q2 ) to R = (r1 , r2 ) is found
by taking the difference of coordinates, so that
−−→
QR = hr1 − q1 , r2 − q2 i.
y
P (3, 4)
−−→
OP = h3, 4i
v = hv1 , v2 , v3 , . . . , vn i.
The next activity will help us to become accustomed to vectors and oper-
ations on vectors in three dimensions.
Activity 9.2.2 An article by C.Kenneth Tanner of the University of Georgia
argues that, due to the concept of social distance, a secondary school classroom
for 20 students should have 1344 square feet of floor space. Suppose a classroom
is 32 feet by 42 feet by 8 feet. Set the origin O of the classroom to be its center.
In this classroom, a student is sitting on a chair whose seat is at location
A = (9, −6, −1.5), an overhead projector is located at position B = (0, 1, 7),
and the teacher is standing at point C = (−2, 20, −4), all distances measured
in feet. Determine the components of the indicated vectors and explain in
context what each represents.
−→ −−→ −−→ −−→ −→ −−→
a. OA b. OB c. OC d. AB e. AC f. BC
v − u = v + (−1)u.
The vectors i, j, and k are called the standard unit vectors (as we will learn
momentarily, unit vectors have length 1), and are important in the physical
sciences.
v + u = hv1 , v2 i + hu1 , u2 i
= hv1 + u1 , v2 + u2 i
= hu1 + v1 , u2 + v2 i
= hu1 , u2 i + hv1 , v2 i
= u + v.
2. (v + u) + w = v + (u + w)
3. The vector 0 = h0, 0, . . . , 0i has the property that v + 0 = v. The
vector 0 is called the zero vector.
4. (−1)v + v = 0. The vector (−1)v = −v is called the additive
inverse of the vector v.
5. (a + b)v = av + bv
6. a(v + u) = av + au
7. (ab)v = a(bv)
8. 1v = v.
u u u
v v
Figure 9.2.5 A vector sum (left), summing displacements (center), the paral-
lelogram law (right).
If we think of these vectors as displacements in the plane, we find a ge-
ometric way to envision vector addition. For instance, the vector u + v will
represent the displacement obtained by following the displacement u with the
displacement v. We may picture this by placing the tail of v at the tip of u,
as seen in the center of Figure 9.2.5.
Of course, vector addition is commutative so we obtain the same sum if
we place the tail of u at the tip of v. We therefore see that u + v appears as
the diagonal of the parallelogram determined by u and v, as shown at right in
Figure 9.2.5.
Vector subtraction has a similar interpretation. At left in Figure 9.2.6
we see vectors u, v, and w = u + v. If we rewrite v = w − u, we have
the arrangement shown at right in Figure 9.2.6. In other words, to form the
difference w − u, we draw a vector from the tip of u to the tip of w.
v w−u
u u
w =u+v w
2v
− 2v
u u
v v
b. What is 0v?
c. On the axes at right in Figure 9.2.8, sketch the vectors −3v, −2v, −1v,
2v, and 3v.
d. Give a geometric description of the set of terminal points of the vectors
tv where t is any scalar.
e. On the set of axes at right in Figure 9.2.8, sketch the vectors u − 3v,
u − 2v, u − v, u + v, and u + 2v.
f. Give a geometric description of the set of terminal points of the vectors
u + tv where t is any scalar.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 28
v v2
3
A
x x
2 4 v1
−−→
Figure 9.2.9 Left: AB. Right: An arbitrary vector, v.
a. Let A = (2, 3) and B = (4, 7), as shown at left in Figure 9.2.9. Compute
−−→
|AB|.
b. Let v = hv1 , v2 i be the vector in R2 with components v1 and v2 as shown
at right in Figure 9.2.9. Use the distance formula to find a general formula
for |v|.
c. Let v = hv1 , v2 , v3 i be a vector in R3 . Use the distance formula to find a
general formula for |v|.
d. Suppose that u = h2, 3i and v = h−1, 2i. Find |u|, |v|, and |u + v|. Is it
true that |u + v| = |u| + |v|?
e. Under what conditions will |u + v| = |u| + |v|? (Hint: Think about how
u, v, and u + v form the sides of a triangle.)
f. With the vector u = h2, 3i, find the lengths of 2u, 3u, and −2u, respec-
tively, and use proper notation to label your results.
g. If t is any scalar, how is |tu| related to |u|?
h. A unit vector is a vector whose magnitude is 1. Of the vectors i, j, and
i + j, which are unit vectors?
i. Find a unit vector v whose direction is the same as u = h2, 3i. (Hint:
Consider the result of part (g).)
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 29
9.2.7 Summary
• A vector is an object that possesses the attributes of magnitude and di-
rection. Examples of vector quantities are position, velocity, acceleration,
and force.
• Two vectors are equal if they have the same direction and magnitude.
Notice that position is not considered, so a vector is independent of its
location.
• If u = hu1 , u2 , . . . , un i and v = hv1 , v2 , . . . , vn i are two vectors in Rn ,
then their vector sum is the vector
u + v = hu1 + v1 , u2 + v2 , . . . , un + vn i.
9.2.8 Exercises
1. For each of the following, perform the indicated computation.
(a) (10 ĩ + 7 j̃ − 5 k̃) − (−6 ĩ + 4 j̃ + 7 k̃) =
(b) (10 ĩ + 6 j̃ − 3 k̃) − 2(−3 ĩ + 10 j̃ + 8 k̃) =
Answer 1. 16 ĩ + 3 j̃ − 12 k̃
Answer 2. 16 ĩ − 14 j̃ − 19 k̃
2. Find a vector a that has the same direction as h−6, 7, 6i but has length 5.
Answer: a =
Answer. h−2.72728, 3.18182, 2.72728i
3. Let a =< −3, −4, −4 > and b =< 2, 2, 4 >.
Show that there are scalars s and t so that sa + tb =< 20, 24, 32 >.
You might want to sketch the vectors to get some intuition.
s=
t=
Answer 1. −4
Answer 2. 4
4. Resolve the following vectors into components:
(a) The vector ~v in 2-space of length 5 pointing up at an angle of π/4
measured from the positive x-axis.
~v = ~i + ~j
(b) The vector w~ in 3-space of length 3 lying in the xz-plane pointing
upward at an angle of 2π/3 measured from the positive x-axis.
~v = ~i + ~j + ~k
Answer 1. 5 · 0.707107
Answer 2. 5 · 0.707107
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 30
Answer 3. 3 · (−0.5)
Answer 4. 0
Answer 5. 3 · 0.866025
5. Find all vectors ~v in 2 dimensions having ||~v || = 13 where the ĩ-component
of ~v is 5 ĩ.
vectors:
(If you find more than one vector, enter them in a comma-separated
list.)
Answer. 5 ĩ + 12 j̃, 5 ĩ − 12 j̃
6. Which is traveling faster, a car whose velocity vector is 26~i + 31~j, or a
car whose velocity vector is 40~i, assuming that the units are the same for
both directions?
( the first car the second car) is the faster car.
At what speed is the faster car traveling?
speed =
Answer 1. the first car
Answer 2. 40.4599
7. Let a = h2, 3, 3i and b = h0, 1, 4i.
Compute:
a+b = ( , , )
a−b = ( , , )
2a = ( , , )
3a + 4b = ( , , )
|a| =
Answer 1. 2
Answer 2. 4
Answer 3. 7
Answer 4. 2
Answer 5. 2
Answer 6. −1
Answer 7. 4
Answer 8. 6
Answer 9. 6
Answer 10. 6
Answer 11. 13
Answer 12. 25
Answer 13. 4.69041575982343
8. Find the length of the vectors
(a) 3 ĩ − j̃ − 3 k̃: length =
(b) −1.6 ĩ + 0.4 j̃ − 1.2 k̃: length =
Answer 1. 4.3589
Answer 2. 2.03961
9. For each of the following, perform the indicated operations on the vectors
~a = 5 j̃ + k̃, ~b = ĩ + 5 j̃ + k̃, ~z = ĩ + 4 j̃.
(a) 5~a + 4~b =
(b) 4~a + 5~b − 5~z =
Answer 1. 4 ĩ + 45 j̃ + 9 k̃
Answer 2. 25 j̃ + 9 k̃
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 31
~ = a2 ~i + 9 ~j.
10. Find the value(s) of a making ~v = 7a~i − 3 ~j parallel to w
a=
(If there is more than one value of a, enter the values as a comma-
separated list.)
Answer. 0, −21
11. (a) Find a unit vector from the point P = (1, 2) and toward the point
Q = (6, 14).
~u =
(b) Find a vector of length 26 pointing in the same direction.
~v =
Answer 1. 0.384615 ĩ + 0.923077 j̃
Answer 2. 10 ĩ + 24 j̃
12. A truck is traveling due north at 40 km/hr approaching a crossroad. On
a perpendicular road a police car is traveling west toward the intersection
at 35 km/hr. Both vehicles will reach the crossroad in exactly one hour.
Find the vector currently representing the displacement of the truck with
respect to the police car.
displacement d~ =
Answer. −35 ĩ − 40 j̃
13. Let v = h1, −2i, u = h0, 4i, and w = h−5, 7i.
a. Determine the components of the vector u − v.
b. Determine the components of the vector 2v − 3u.
c. Determine the components of the vector v + 2u − 7w.
d. Determine scalars a and b such that av + bu = w. Show all of your
work in finding a and b.
14. Let u = h2, 1i and v = h1, 2i.
a. Determine the components and draw geometric representations of
the vectors 2u, 21 u, (−1)u, and (−3)u on the same set of axes.
u · v = u1 v1 + u2 v2 .
d. On the axes in Figure 9.3.1, plot the vectors u = h1, 3i and v = h−3, 1i.
Then, find u · v. What is the angle between these vectors?
4 y
x
-4 -2 2 4
-2
-4
Figure 9.3.1 For part (d)
4 y
x
-4 -2 2 4
-2
-4
Figure 9.3.2 For part (e)
For each of the following vectors v, plot the vector on Figure 9.3.2 and
then compute the dot product u · v.
• v = h3, 2i.
• v = h3, 0i.
• v = h3, −1i.
• v = h3, −2i.
• v = h3, −4i.
f. Based upon the previous part of this activity, what do you think is the
sign of the dot product in the following three cases shown in Figure 9.3.3?
v
v v
u u u
u · v = u1 v1 + u2 v2 + . . . + un vn .
♦
(As we will see shortly, the dot product arises in physics to calculate the
work done by a vector force in a given direction. It might be more natural
to define the dot product in this context, but it is more convenient from a
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 35
v u u−v
θ θ
u
2π − θ
v
Figure 9.3.5 Left: The angle between u and v. Right: The triangle formed
by u, v, and u − v.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 36
To determine this angle, we may apply the Law of Cosines to the triangle
shown at right in Figure 9.3.5.
Using the fact that the dot product of a vector with itself gives us the
square of its length, together with the properties of the dot product, we find:
c. The angle between the vectors y = h1, 2, −3i and z = h−2, 1, 1i to the
nearest tenth of a degree.
d. If the angle between the vectors u and v is a right angle, what does the
expression u · v = |u||v| cos(θ) say about their dot product?
e. If the angle between the vectors u and v is acute—that is, less than
π/2—what does the expression u · v = |u||v| cos(θ) say about their dot
product?
f. If the angle between the vectors u and v is obtuse—that is, greater than
π/2—what does the expression u · v = |u||v| cos(θ) say about their dot
product?
More generally, the sign of the dot product gives us useful information
about the relative orientation of the vectors. If we remember that
v
v v
u u u
u·v >0 u·v =0 u·v <0
A θ B
|F| cos θ
Figure 9.3.7 A force F displacing an object.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 38
This means that the work is determined only by the magnitude of the force
applied parallel to the displacement. Consequently, if we are given two vectors
u and v, we would like to write u as a sum of two vectors, one of which is
parallel to v and one of which is orthogonal to v. We take up this task after
the next activity.
Activity 9.3.4 Determine the work done by a 25 pound force acting at a 30◦
angle to the direction of the object’s motion, if the object is pulled 10 feet. In
addition, is more work or less work done if the angle to the direction of the
object’s motion is 60◦ ?
9.3.5 Projections
u u
proj⊥v u proj⊥v u
θ
θ
projv u v projv u v
π
Figure 9.3.8 Left: projv u. Right: projv u when θ > 2.
Suppose we are given two vectors u and v as shown at left in Figure 9.3.8.
Motivated by our discussion of work, we would like to write u as a sum of two
vectors, one of which is parallel to v and one of which is orthogonal. That is,
we would like to write
u = projv u + proj⊥v u, (9.3.2)
where projv u is parallel to v and proj⊥v u is orthogonal to v. We call the
vector projv u the projection of u onto v. Note that, as the diagram at right
in Figure 9.3.8 illustrates, it is also possible to create a projection even if the
angle between the vectors u and v exceeds π2 .
To find the vector projv u, we will dot both sides of Equation (9.3.2) with
the vector v, to find that
u · v = (projv u + proj⊥v u) · v
= (projv u) · v + (proj⊥v u) · v
= (projv u) · v.
u · v = (projv u) · v = sv · v,
and hence
u·v u·v
projv u = v= v
v·v |v|2
It is sometimes useful to write projv u as a scalar times a unit vector in the
direction of v. We call this scalar the component of u along v and denote it
as compv u. We therefore have
u·v u·v v v
projv u = v= = compv u ,
|v|2 |v| |v| |v|
so that
u·v
compv u = .
|v|
The dot product and projections.
Moreover, since
u = projv u + proj⊥v u,
it follows that
proj⊥v u = u − projv u.
This shows that once we have computed projv u, we can find proj⊥v u simply
by calculating the difference of two known vectors.
Activity 9.3.5 Let u = h2, 6i.
a. Let v = h4, −8i. Find compv u, projv u and proj⊥v u, and draw a picture
to illustrate. Finally, express u as the sum of two vectors where one is
parallel to v and the other is perpendicular to v.
b. Now let v = h−2, 4i. Without doing any calculations, find projv u. Ex-
plain your reasoning. (Hint: Refer to the picture you drew in part (a).)
c. Find a vector w not parallel to z = h3, 4i such that projz w has length
10. Note that there are infinitely many different answers.
9.3.6 Summary
• The dot product of two vectors in Rn , u = hu1 , u2 , . . . , un i and v =
hv1 , v2 , . . . , vn i, is the scalar
u · v = u1 v1 + u2 v2 + · · · + un vn .
• Two vectors are orthogonal if the angle between them is π/2. In terms
of the dot product, the vectors u and v are orthogonal if and only if
u · v = 0.
• The projection of a vector u in Rn onto a vector v in Rn is the vector
u·v
projv u = v.
v·v
9.3.7 Exercises
1. Find a · b if
a = h−2, 2, −3i and b = h4, 0, 3i
a·b=
Is the angle between the vectors "acute", "obtuse" or "right"?
Answer 1. −17
Answer 2. obtuse
2. Determine if the pairs of vectors below are "parallel", "orthogonal", or
"neither".
a = h−1, −2, 2i and b = h4, 8, 10i are
a = h−1, −2, 2i and b = h4, 8, −8i are
a = h−1, −2, 2i and b = h2, 4, −5i are
Answer 1. orthogonal
Answer 2. parallel
Answer 3. neither
3. Perform the following operations on the vectors ~u = h0, 5, −4i, ~v =
h−2, 0, 3i, and w ~ = h−3, 0, 1i.
~u · w~=
(~u · ~v )~u =
((w~ · w)~
~ u) · ~u =
~u · ~v + ~v · w
~=
Answer 1. −4
Answer 2. h0, −60, 48i
Answer 3. 410
Answer 4. −3
π
4. Find a · b if |a| = 7, |b| = 7, and the angle between a and b is − 10
radians.
a·b =
Answer. 46.6017692984625
5. What is the angle in radians between the vectors
a = (6, -4, 9) and
b = (5, -1, 6)?
Angle:
(radians)
Answer. 0.249309305958578
6. Find a · b if |a| = 9, |b| = 10, and the angle between a and b is − π3
radians.
a·b =
Answer. 45
7. A constant force F = 9i + 3j + 6k moves an object along a straight line
from point (−5, −1, 5) to point (−4, 1, 3).
Find the work done if the distance is measured in meters and the
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 41
positive
10. Find the angle between the diagonal of a cube of side length 8 and the
diagonal of one of its faces. The angle should be measured in radians.
Answer. 0.615479708670387
11. Let v = h−2, 5i in R2 , and let y = h0, 3, −2i in R3 .
a. Is h2, −1i perpendicular to v? Why or why not?
f. θ
14. One of the properties of the dot product is that (u+v)·w = (u·w)+(v·w).
That is, the dot product distributes over vector addition on the right. Here
we investigate whether the dot product distributes over vector addition
on the left.
a. Let u = h1, 2, −1i, v = h4, −3, 6i, and v = h4, 7, 2i. Calculate
u · (v + w) and (u · v) + (u · w).
x · (y + z) = (x · y) + (x · z)
methane, we can place the center carbon atom at the point 12 , 12 , 12 and
the hydrogen atoms at the points (0, 0, 0), (1, 1, 0), (1, 0, 1), and (0, 1, 1).
Find the bond angle for methane to the nearest tenth of a degree.
d. Up to this point, the products you have seen, such as the product of
real numbers and the dot product of vectors, have been commutative,
meaning that the product does not depend on the order of the terms.
For instance, 2 · 5 = 5 · 2. The table above suggests, however, that
the cross product is anti-commutative: for any vectors u and v in R3 ,
u × v = −v × u. If we consider the case when u = v, this shows that
v × v = −(v × v). What does this tell us about v × v; in particular,
what vector is unchanged by scalar multiplication by −1?
e. It is not difficult to show that the cross product interacts with scalar
multiplication and vector addition as one would expect: that is
(cu) × v = c(u × v)
(u + v) × w = (u × w) + (v × w)
(2i + j) × k = (2i × k) + (j × k)
= 2(i × k) + (j × k)
= − 2j + i.
Using these properties along with Table 9.4.2, find the cross product u×v
if u = 2i + 3j and v = −i + k.
f. Verify that the cross product u×v you just found in part (e) is orthogonal
to both u and v.
g. Consider the vectors u and v in the xy-plane as shown below in Fig-
ure 9.4.3.
y
v
θ u x
Explain why u = |u|i and v = |v| cos(θ)i + |v| sin(θ)j. Then compute the
length of |u × v|.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 45
u × v = (u1 i + u2 j + u3 k) × (v1 i + v2 j + v3 k)
= u1 i × (v1 i + v2 j + v3 k) + u2 j × (v1 i + v2 j + v3 k)
+ u3 k × (v1 i + v2 j + v3 k)
= u1 v1 i × i + u1 v2 i × j + u1 v3 i × k + u2 v1 j × i + u2 v2 j × j
+ u2 v3 j × k + u3 v1 k × i + u3 v2 k × j + u3 v3 k × k
= u1 v2 k − u1 v3 j − u2 v1 k + u2 v3 i + u3 v1 j − u3 v2 i
= (u2 v3 − u3 v2 )i − (u1 v3 − u3 v1 )j + (u1 v2 − u2 v1 )k.
(Like the dot product, the cross product arises in physical applications, e.g.,
torque, but it is more convenient mathematically to begin from an algebraic
perspective.)
The previous calculations lead us to define the cross product of vectors in
R3 as follows.
Definition 9.4.4 The cross product u × v of vectors u = u1 i + u2 j + u3 k
and v = v1 i + v2 j + v3 k in R3 is the vector
♦
At first, this may look intimidating and difficult to remember. However, if
we rewrite the expression in Equation (9.4.1) using determinants, important
structure emerges. The determinant of a 2 × 2 matrix is
a b
= ad − bc.
c d
u2 u3 u1 u3 u1 u2
u×v= i− j+ k.
v2 v3 v1 v3 v1 v2
2. (u + v) × w = (u × w) + (v × w)
3. (cu) × w = c(u × w) = u × (cv)
4. u × v = 0 if u and v are parallel.
5. The cross product is not associative; that is, in general
(u × v) × w 6= u × (v × w).
Just as we found for the dot product, the cross product provides us with
useful geometric information. In particular, both the length and direction of
the cross product u × v encode information about the geometric relationship
between u and v.
Note that the third property stated above says that u × v = 0 if u and v
are parallel. This is reflected in Equation (9.4.2) since sin(θ) = 0 if u and v
are parallel, which implies that u × v = 0.
Equation (9.4.2) also has a geometric implication. Consider the parallelo-
gram formed by two vectors u and v, as shown in Figure 9.4.5.
v
|v| sin θ
θ u
The length, |u × v|, of the cross product of vectors u and v is the area
of the parallelogram determined by u and v.
b. Find the area of the parallelogram in R3 whose vertices are (1, 0, 1),
(0, 0, 1), (2, 1, 0), and (1, 1, 0). (Hint: It might be helpful to draw a
picture to see how the vertices are arranged so you can determine which
vectors you might use.)
=0
To summarize, we have u · (u × v) = 0, which implies that u is orthogonal
to u × v. In the same way, we can show that v is orthogonal to u × v. The
net effect is that u × v is a vector that is perpendicular to both u and v, and
hence u × v is perpendicular to the plane determined by u and v. Moreover,
the direction of u × v is determined by applying the right-hand rule to u and
v, as we saw in Preview Activity 9.4.1. In light of our earlier work that showed
|u||v| sin(θ) = |u × v|., we may now express u × v in the following different
way.
Suppose that u and v are not parallel and that n is the unit vector
perpendicular to the plane containing u and v determined by the right-
hand rule. Then
u × v = |u||v| sin(θ) n.
There is yet one more geometric implication we may draw from this result.
Suppose u, v, and w are vectors in R3 that are not coplanar and that form a
three-dimension parallelepiped as shown in Figure 9.4.6.
h
α v
d. Given the vectors u and v shown below in Figure 9.4.7, sketch the cross
products u × v and v × u.
z
v
u
F
θ
9.4.6 Summary
• The cross product is defined only for vectors in R3 . The cross product of
vectors u = u1 i + u2 j + u3 k and v = v1 i + v2 j + v3 k in R3 is the vector
9.4.7 Exercises
1. If a = i + j + 3k and b = i + j + 2k
Compute the cross product a × b.
a×b= i+ j+ k
Answer 1. −1
Answer 2. 1
Answer 3. 0
2. Suppose ~v · w
~ = 6 and ||~v × w||
~ = 4, and the angle between ~v and w
~ is θ.
Find
(a) tan θ =
(b) θ =
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 51
4
Answer 1. 6
Answer 2. tan−1 4
6
3. You are looking down at a map. A vector u with |u| = 7 points north and
a vector v with |v| = 6 points northeast. The crossproduct u × v points:
A) south
B) northwest
C) up
D) down
Please enter the letter of the correct answer: The magnitude
|u × v| =
Answer 1. D
Answer 2. 29.698484809835
4. If a = i + 8j + k and b = i + 10j + k, find a unit vector with positive first
coordinate orthogonal to both a and b.
i+ j+
k
Answer 1. 0.707106781186547
Answer 2. 0
Answer 3. −0.707106781186547
5. Sketch the triangle with vertices O, P = (0, 7, 6) and Q = (7, 0, 2) and
compute its area using cross products.
Area=
√
Answer. 12 4361
6. Let A = (5, 0, 0), B = (2, −2, −3), and P = (k, k, k). The vector from A
to B is perpendicular to the vector from A to P when k = .
Answer. 1.875
7. Find two unit vectors orthogonal to a = h−4, −4, 5i and b = h−1, −2, 3i
Enter your answer so that the first non-zero coordinate of the first
vector is positive.
First Vector: h , ,
i
Second Vector: h , ,
i
Answer 1. 0.240771706171538
Answer 2. −0.842700971600384
Answer 3. −0.481543412343077
Answer 4. −0.240771706171538
Answer 5. 0.842700971600384
Answer 6. 0.481543412343077
8. Use the geometric definition of the cross product and the properties of the
cross product to make the following calculations.
(a) ((~i + ~j) ×~i) × ~j =
(b) (~j + ~k) × (~j × ~k) =
(c) 5~i × (~i + ~j) =
(d) (~k + ~j) × (~k − ~j) =
Answer 1. h1, 0, 0i
Answer 2. h0, 1, −1i
Answer 3. h0, 0, 5i
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 52
Answer 4. h2, 0, 0i
9. Are the following statements true or false?
(a) For any scalar c and any vector ~v , we have ||c~v || = c||~v ||.
u × (v + w) and (u × v) + (u × w).
x × (y + z) = (x × y) + (x × z)
u1 u2 u3 v1 v2 v3
b. Show that v1 v2 v3 = − u1 u2 u3 . Conclude that
w1 w2 w3 w1 w2 w3
interchanging the first two rows in a 3 × 3 matrix changes the sign
of the determinant. In general (although we won’t show it here),
interchanging any two rows in a 3 × 3 matrix changes the sign of the
determinant.
c. Use the results of parts (a) and (b) to argue that
(u × v) · w = (w × u) · v = (v × w) · u.
2 y
1
x
-2 -1 1 2 3 4
-1
-2
-3
-4
Figure 9.5.1 The line through (2, −1) with slope 32 .
a. Suppose we increase x by 1 from the point (2, −1). How does the y-value
change? What is the point on the line with x-coordinate 3?
b. Suppose we decrease x by 3.25 from the point (2, −1). How does the
y-value change? What is the point on the line with x-coordinate −1.25?
c. Now, suppose we increase x by some arbitrary value 3t from the point
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 55
(2, −1). How does the y-value change? What is the point on the line
with x-coordinate 2 + 3t?
d. Observe that the slope of the line is related to any vector whose y-
component divided by the x-component is the slope of the line. For
the line in this exercise, we might use the vector h3, 2i, which describes
the direction of the line. Explain why the terminal points of the vectors
r(t), where
r(t) = h2, −1i + h3, 2it,
trace out the graph of the line through the point (2, −1) with slope 32 .
e. Now we extend this vector approach to R3 and consider a second example.
Let L be the line in R3 through the point (1, 0, 2) in the direction of the
vector h2, −1, 4i. Find the coordinates of three distinct points on line L.
Explain your thinking.
f. Find a vector in the form
whose terminal points trace out the line L that is described in (e). That
is, you should be able to locate any point on the line by determining a
corresponding value of t.
5 y 5 y
4 4
3 3
2 2 v
h1, mi
1 1
P
x x
-2 -1 1 2 3 4 -2 -1 1 2 3 4
-1 -1
Figure 9.5.2 A vector description of a line
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 56
Definition 9.5.3 A line in space is the set of terminal points of vectors ema-
nating from a given point P that are parallel to a fixed vector v. ♦
The fixed vector v in the definition is called a direction vector for the line.
As we saw in Preview Activity 9.5.1, to find an equation for a line through
point P in the direction of vector v, observe that any vector parallel to v will
have the form tv for some scalar t. So, any vector emanating from the point
P in a direction parallel to the vector v will be of the form
−−→
OP + vt (9.5.1)
4 4 4
3 3 3
tv
2 2 tv 2
tv
−−→ −−→ −−→
1 OP 1 OP 1 OP
-2 -1 1 2 3 4 -2 -1 1 2 3 4 -2 -1 1 2 3 4
-1 -1 -1
y y y
-3 -2 -1 1 2 3 -3 -2 -1 1 2 3 -3 -2 -1 1 2 3
1 -1 P2 1 -1 P2 1 -1 P2
L L L
2 2 2
-2 P1 -2 P1 -2 P1
x x x
3 3 3
-3 -3 -3
Activity 9.5.2 Let P1 = (1, 2, −1) and P2 = (−2, 1, −2). Let L be the line in
R3 through P1 and P2 , and note that three snapshots of this line are shown in
Figure 9.5.5.
a. Find a direction vector for the line L.
b. Find a vector equation of L in the form r(t) = r0 + tv.
c. Consider the vector equation s(t) = h−5, 0, −3i + th6, 2, 2i. What is the
direction of the line given by s(t)? Is this new line parallel to line L?
Notice that there are many different parametric equations for the same line.
For example, choosing another point P on the line or another direction vector
v produces another set of parametric equations. It is sometimes useful to think
of t as a time parameter and the parametric equations as telling us where we
are on the line at each time. In this way, the parametric equations describe a
particular walk taken along the line; there are, of course, many possible ways
to walk along a line.
Activity 9.5.3 Let P1 = (1, 2, −1) and P2 = (−2, 1, −2), and let L be the line
in R3 through P1 and P2 , which is the same line as in Activity 9.5.2.
a. Find parametric equations of the line L.
b. Does the point (1, 2, 1) lie on L? If so, what value of t results in this
point?
c. Consider another line, K, whose parametric equations are
x(s) = 11 + 4s, y(s) = 1 − 3s, z(s) = 3 + 2s.
What is the direction of the line K?
d. Do the lines L and K intersect? If so, provide the point of intersection
and the t and s values, respectively, that result in the point. If not,
explain why.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 58
p
n = ha, b, ci
P = (x, y, z)
P0 = (x0 , y0 , z0 )
Equations of a plane.
We may take the scalar equation of a plane a little further and note that
since
a(x − x0 ) + b(y − y0 ) + c(z − z0 ) = 0,
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 59
P2 p
n
P
P0 P1
Activity 9.5.5 Let P0 = (1, 2, −1), P1 = (1, 0, −1), and P2 = (0, 1, 3) and let
p be the plane containing P0 , P1 , and P2 .
−−−→ −−−→
a. Determine the components of the vectors P0 P1 and P0 P2 .
b. Find a normal vector n to the plane p.
c. Find a scalar equation of the plane p.
d. Consider a second plane, q, with scalar equation −3(x − 1) + 4(y + 3) +
2(z − 5) = 0. Find two different points on plane q, as well as a vector m
that is normal to q.
e. The angle between two planes is the acute angle between their respective
normal vectors. What is the angle between planes p and q?
9.5.4 Summary
• While lines in R3 do not have a slope, like lines in R2 they can be char-
acterized by a point and a direction vector. Indeed, we define a line in
space to be the set of terminal points of vectors emanating from a given
point that are parallel to a fixed vector.
• Vectors play a critical role in representing the equation of a line. In
particular, the terminal points of the vector r(t) = r0 + tv define a linear
function r in space through the terminal point of the vector r0 in the
direction of the vector v, tracing out a line in space.
• A plane in space is the set of all terminal points of vectors emanating
from a given point perpendicular to a fixed vector.
−−−→
• If P1 , P2 , and P3 are non-collinear points in space, the vectors P1 P2 and
−−−→ −−−→ −−−→
and P1 P3 are vectors in the plane and the vector n = P1 P2 × P1 P3 is
a normal vector to the plane. So any point P in the plane satisfies the
−−→
equation P P1 · n = 0. If we let P = (x, y, z), n = ha, b, ci be the normal
vector, and P1 = (x0 , y0 , z0 ), we can also represent the plane with the
equation
a(x − x0 ) + b(y − y0 ) + c(z − z0 ) = 0.
9.5.5 Exercises
1. Rewrite the vector equation r(t) = (−2t)i + (3 − 3t)j + (1 + 3t)k as the
corresponding parametric equations for the line.
x(t) =
y(t) =
z(t) =
Answer 1. 0 + (−2) t
Answer 2. 3 + (−3) t
Answer 3. 1 + 3t
2. Find the vector and parametric equations for the line through the point
P(2, -1, -5) and parallel to the vector −5i − 3j − 3k.
Vector Form: r = h , , -5 i + th
, , -3 i
Parametric form (parameter t, and passing through P when t = 0):
x = x(t) =
y = y(t) =
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 61
z = z(t) =
Answer 1. 2
Answer 2. −1
Answer 3. −5
Answer 4. −3
Answer 5. 2 + t(−5)
Answer 6. −1 + t(−3)
Answer 7. −5 + t(−3)
3. Consider the line which passes through the point P(3, -5, -1), and which
is parallel to the line x = 1 + 6t, y = 2 + 2t, z = 3 + 6t
Find the point of intersection of this new line with each of the coordi-
nate planes:
xy-plane: ( , , )
xz-plane: ( , , )
yz-plane: ( , , )
Answer 1. 4
Answer 2. −4.66666666666667
Answer 3. 0
Answer 4. 18
Answer 5. 0
Answer 6. 14
Answer 7. 0
Answer 8. −6
Answer 9. −4
4. Find the point at which the line h4, 2, 4i + th−3, −3, −4i intersects the
plane −5x + 5y − 3z = 2.
( , ,
)
Answer 1. −2
Answer 2. −4
Answer 3. −4
5. Find an equation of a plane containing the three points (5, -2, -2), (2, -5,
1), (2, -4, 3) in which the coefficient of x is -9.
= 0.
Answer. −9(x − 5) + 6(y − (−2)) + (−3)(z − (−2))
6. Find an equation for the plane containing the line in the xy-plane where
y = 1, and the line in the xz-plane where z = 2.
equation:
Answer. z = −2y + 2
7. Find the angle in radians between the planes 4x + z = 1 and 5y + z = 1.
Answer. 1.52321322351792
8. A store sells CDs at one price and DVDs at another price. The figure
below shows the revenue (in dollars) of the music store as a function of
the number, c, of CDs and the number, d, of DVDs that it sells. The
values of the revenue are shown on each line.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 62
(Hint: for this problem there are many possible ways to estimate the
requisite values; you should be able to find information from the figure that
allows you to give an answer that is essentially exact.)
(a) What is the price of a CD?
dollars
(b) What is the price of a DVD?
dollars
Answer 1. 10
Answer 2. 14
9. The table below gives the number of calories burned per minute for some-
one roller-blading, as a function of the person’s weight in pounds and speed
in miles per hour [from the August 28,1994, issue of Parade Magazine].
calories burned per minute
weight\speed 8 9 10 11
120 4.2 5.8 7.4 8.9
140 5.1 6.7 8.3 9.9
160 6.1 7.7 9.2 10.8
180 7 8.6 10.2 11.7
200 7.9 9.5 11.1 12.6
(a) Suppose that a 180 lb person and a 200 person both go 8 miles,
the first at 9 mph and the second at 8 mph.
How many calories does the 180 lb person burn?
How many calories does the 200 lb person burn?
(b) We might also be interested in the number of calories each person
burns per pound of their weight.
How many calories per pound does the 180 lb person burn?
How many calories per pound does the 200 lb person burn?
480·8.6
Answer 1. 9
480·7.9
Answer 2. 8
458.667
Answer 3. 180
474
Answer 4. 200
10. The vector and parametric forms of a line allow us to easily describe line
segments in space.
Let P1 = (1, 2, −1) and P2 = (−2, 1, −2), and let L be the line in R3
through P1 and P2 as in Activity 9.5.2.
a. What value of the parameter t makes (x(t), y(t), z(t)) = P1 ? What
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 63
c. What about the line segment (along the same line) from (7, 4, 1) to
(−8, −1, −4)?
d. Now, consider a segment that lies on a different line: parameterize
the segment that connects point R = (4, −2, 7) to Q = (−11, 4, 27)
in such a way that t = 0 corresponds to point Q, while t = 2
corresponds to R.
11. This exercise explores key relationships between a pair of lines. Consider
the following two lines: one with parametric equations x(s) = 4 − 2s,
y(s) = −2 + s, z(s) = 1 + 3s, and the other being the line through
(−4, 2, 17) in the direction v = h−2, 1, 5i.
a. Find a direction vector for the first line, which is given in parametric
form.
b. Find parametric equations for the second line, written in terms of
the parameter t.
c. Show that the two lines intersect at a single point by finding the
values of s and t that result in the same point. Then find the point
of intersection.
d. Find the acute angle formed where the two lines intersect, noting
that this angle will be given by the acute angle between their re-
spective direction vectors.
e. Find an equation for the plane that contains both of the lines de-
scribed in this problem.
12. This exercise explores key relationships between a pair of planes. Consider
the following two planes: one with scalar equation 4x − 5y + z = −2,
and the other which passes through the points (1, 1, 1), (0, 1, −1), and
(4, 2, −1).
a. Find a vector normal to the first plane.
b. Find a scalar equation for the second plane.
c. Find the angle between the planes, where the angle between them
is defined by the angle between their respective normal vectors.
• hcos(0), sin(0)i,
• cos π2 , sin π2 ,
• cos 3π 3π
4 , sin 4 ,
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 65
5π 5π
• cos 4 , sin 4 , and
7π 7π
• cos 4 , sin 4
and thus we have three coordinate functions that enable us to represent the
curve. The variable t is called a parameter and the equations x = x(t), y = y(t),
and z = z(t) are called parametric equations (or a parameterization of the
curve). The function r whose output is the vector from the origin to a point
on the curve is defined by
Note that the input of r is the real-valued parameter t and the correspond-
ing output is vector hx(t), y(t), z(t)i. Such a function is called a vector-valued
function because each real number input generates a vector output. More
formally, we state the following definition.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 66
d. r(t) = ht2 sin(t) cos(t), 0.9t cos(t2 ), sin(t)i (Note that this defines a curve
in 3-space.)
e. Experiment with different formulas for x(t) and y(t) and ranges for t
to see what other interesting curves you can generate. Share your best
results with peers.
Recall from our earlier work that the traces and level curves of a function
are themselves curves in space. Thus, we may determine parameterizations for
them. For example, if z = f (x, y) = cos(x2 +y 2 ), the y = 1 trace of the function
is given by setting y = 1 and letting x be parameterized by the variable t; then,
the trace is the curve whose parameterization is ht, 1, cos(t2 + 1)i.
Activity 9.6.4 Consider the paraboloid defined by f (x, y) = x2 + y 2 .
a. Find a parameterization for the x = 2 trace of f . What type of curve
does this trace describe?
b. Find a parameterization for the y = −1 trace of f . What type of curve
does this trace describe?
c. Find a parameterization for the level curve f (x, y) = 25. What type of
curve does this trace describe?
d. How do your responses change to all three of the preceding questions if
you instead consider the function g defined by g(x, y) = x2 − y 2 ? (Hint
for generating one of the parameterizations: sec2 (t) − tan2 (t) = 1.)
9.6.2 Summary
• A vector-valued function is a function whose input is a real parameter t
and whose output is a vector that depends on t. The graph of a vector-
valued function is the set of all terminal points of the output vectors with
their initial points at the origin.
9.6.3 Exercises
1. Find the domain of the vector function
√
1
r(t) = ln(11t), t + 10, √
12 − t
using interval notation.
Domain:
Answer. (0, 12)
2. Find a parametrization of the circle of radius 6 in the xy-plane, centered
at the origin, oriented clockwise. The point (6, 0) should correspond to
t = 0. Use t as the parameter for all of your answers.
x(t) =
y(t) =
Answer. 6 cos(t) ; −6 sin(t)
3. Find a vector parametrization of the circle of radius 7 in the xy-plane, cen-
tered at the origin, oriented clockwise so that the point (7, 0) corresponds
to t = 0 and the point (0, −7) corresponds to t = 1.
~r(t) =
Answer. 7 cos πt πt
2 , (−7) sin 2
4. Find a vector parametric equation ~r(t) for the line through the points
P = (4, 1, 1) and Q = (9, −3, 3) for each of the given conditions on the
parameter t.
(a) If ~r(0) = h4, 1, 1i and ~r(8) = h9, −3, 3i, then
~r(t) =
(b) If ~r(6) = P and ~r(10) = Q, then
~r(t) =
(c) If the points P and Q correspond to the parameter values t = 0
and t = −2, respectively, then
~r(t) =
Answer 1. (4, 1, 1) + 8t h5, −4, 2i
t−6
Answer 2. (4, 1, 1) + 4 h5, −4, 2i
− 2t h5, −4, 2i
Answer 3. (4, 1, 1) +
5. Suppose parametric equations for the line segment between (9, −6) and
(−2, 5) have the form:
x = a + bt
y = c + dt
If the parametric curve starts at (9, −6) when t = 0 and ends at (−2, 5)
at t = 1, then find a, b, c, and d.
a= ,b = ,c = ,d =
.
Answer 1. 9
Answer 2. −11
Answer 3. −6
Answer 4. 11
6. Find a parametrization of the curve x = −5z 2 in the xz-plane. Use t as
the parameter for all of your answers.
x(t) =
y(t) =
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 69
z(t) =
Answer. −5t2 ; 0; t
7. Find parametric equations for the quarter-ellipse from (2, 0, 9) to (0, −3, 9)
centered at (0, 0, 9) in the plane z = 9. Use the interval 0 ≤ t ≤ π/2.
x(t) =
y(t) =
z(t) =
Answer. 2 cos(t) ; −3 sin(t) ; 9
8. Are the following statements true or false?
(a) The line parametrized by x = 7, y = 5t, z = 6 + t is parallel to the
x-axis.
(b) A parametrization of the graph of y = ln(x) for x > 0 is given by
x = et , y = t for −∞ < t < ∞.
(c) The parametric curve x = (3t + 4)2 , y = 5(3t + 4)2 − 9, for 0 ≤ t ≤ 3
is a line segment.
9. Find a vector function that represents the curve of intersection of the
paraboloid z = 5x2 + 5y 2 and the cylinder y = 5x2 . Use the variable t for
the parameter.
r(t) = ht, , i
Answer 1. 5tt
Answer 2. 5tt + 125t4
10. A bicycle wheel has radius R. Let P be a point on the spoke of a wheel
at a distance d from the center of the wheel. The wheel begins to roll to
the right along the the x-axis. The curve traced out by P is given by the
following parametric equations:
x = 15θ − 10 sin(θ)
y = 15 − 10 cos(θ)
What must we have for R and d?
R= d=
Answer 1. 15
Answer 2. 10
11. A standard parameterization for the unit circle is hcos(t), sin(t)i, for 0 ≤
t ≤ 2π.
a. Find a vector-valued function r that describes a point traveling
√ √ along
the unit circle so that at time t = 0 the point is at 22 , 22 and
travels clockwise along the circle as t increases.
b. Find a vector-valued function r that describes a point traveling
√ √ along
the unit circle so that at time t = 0 the point is at 22 , 22 and
travels counter-clockwise along the circle as t increases.
c. Find a vector-valued function r that describes a pointtraveling
√ √
along
2 2
the unit circle so that at time t = 0 the point is at − 2 , 2 and
travels clockwise along the circle as t increases.
d. Find a vector-valued function r that describes a point traveling along
the unit circle so that at time t = 0 the point is at (0, 1) and makes
one complete revolution around the circle in the counter-clockwise
direction on the interval [0, π].
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 70
12. Let a and b be positive real numbers. You have probably seen the equation
(x−h)2 2
a2 + (y−k)
b2 = 1 that generates an ellipse, centered at (h, k), with a
horizontal axis of length 2a and a vertical axis of length 2b.
a. Explain why the vector function r defined by r(t) = ha cos(t), b sin(t)i,
2 2
0 ≤ t ≤ 2π is one parameterization of the ellipse xa2 + yb2 = 1.
x2 y2
b. Find a parameterization of the ellipse 4 + 16 = 1 that is traversed
counterclockwise.
(x+3)2 (y−2)2
c. Find a parameterization of the ellipse 4 + 9 = 1.
r(t + h)
r(t)
b. Is r(t+h)−r(t)
h a vector or a scalar? Sketch a representative vector r(t+h)−r(t)
h
with h < 1 in Figure 9.7.2.
c. Think of r(t) as providing the position of an object moving along the
curve these vectors trace out. What do you think that the vector r(t+h)−r(t)
h
measures? Why? (Hint: You might think analogously about difference
quotients such as f (x+h)−f
h
(x)
or s(t+h)−s(t)
h from calculus I.)
r(t + h) − r(t)
lim .
h→0 h
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 73
r(t + h)
r(t + h)
r(t + h)
r(t + h) − r(t)
r0 (t) = lim
h→0 h
is the instantaneous rate of change of r(t) at time t (for those values of t for
which the limit exists), so r0 (t) = v(t) is the instantaneous velocity of the
object at time t. Furthermore, we can interpret the derivative r0 (t) as the
direction vector of the line tangent to the graph of r at the value t.
Similarly,
v(t + h) − v(t)
v0 (t) = r00 (t) = lim
h→0 h
is the instantaneous rate of change of the velocity of the object at time t,
for those values of t for which the limits exists, and thus v0 (t) = a(t) is the
acceleration of the moving object.
Note well: Both the velocity and acceleration are vector quantities: they
have magnitude and direction. By contrast, the magnitude of the velocity
vector, |v(t)|, which is the speed of the object at time t, is a scalar quantity.
r(t + h) − r(t)
r0 (t) = lim
h→0 h
[x(t + h) − x(t)]i + [y(t + h) − y(t)]j + [z(t + h) − z(t)]k
= lim
h→0 h
x(t + h) − x(t) y(t + h) − y(t)
= lim i + lim j
h→0 h h→0 h
z(t + h) − z(t)
+ lim k
h→0 h
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 74
d
r(t) = x0 (t)i + y 0 (t)j + z 0 (t)k
dt
for those values of t at which x, y, and z are differentiable.
Activity 9.7.3 For each of the following vector-valued functions, find r0 (t).
a. r(t) = hcos(t), t sin(t), ln(t)i.
b. r(t) = ht2 + 3t, e−2t , t2 +1
t
i.
2. d
dt [f (t)r(t)] = f (t)r0 (t) + f 0 (t)r(t)
3. d
dt [r(t) · s(t)] = r0 (t) · s(t) + r(t) · s0 (t)
4. d
dt [r(t) × s(t)] = r0 (t) × s(t) + r(t) × s0 (t)
5. d
dt [r(f (t))] = f 0 (t)r0 (f (t)).
Note well. When applying these properties, use care to interpret the quan-
tities involved as either scalars or vectors. For example, r(t) · s(t) defines a
scalar function because we have taken the dot product of two vector-valued
functions. However, r(t) × s(t) defines a vector-valued function since we have
taken the cross product of two vector-valued functions.
Activity 9.7.4 The left side of Figure 9.7.4 shows the curve described by the
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 75
In the same way, we expect that a smooth curve in 3-space will be locally
linear. In the following activity, we investigate how to find the tangent line to
such a curve. Recall from our work in Section 9.5 that the vector equation of
a line that passes through the point at the tip of the vector L0 = hx0 , y0 , z0 i
in the direction of the vector u = ha, b, ci can be written as
L(t) = L0 + tu.
c. Find the parametric equations of the line tangent to the graph of r when
t = π.
d. Sketch a plot of the curve r(t) and its tangent line near the point where
t = π. In addition, include a sketch of r0 (π). What is the important role
of r0 (π) in this activity?
We see that our work in Activity 9.7.5 can be generalized. Given a differ-
entiable vector-valued function r, the tangent line to the curve at the input
value a is given by
L(t) = r(a) + tr0 (a). (9.7.1)
Here we see that because the tangent line is determined entirely by a given
point and direction, the point is provided by the function r, evaluated at t = a,
while the direction is provided by the derivative, r0 , again evaluated at t = a.
Note how analogous the formula for L(t) is to the tangent line approximation
from single-variable calculus: in that context, for a given function y = f (x) at
a value x = a, we found that the tangent line can be expressed by the linear
function y = L(x) whose formula is
Equation (9.7.1) for the tangent line L(t) to the vector-valued function r(t)
is nearly identical. Indeed, because there are multiple parameterizations for a
single line, it is even possible to write the parameterization as
R0 (t) = r(t).
R
The indefinite integral r(t) dt of a vector-valued function r is the gen-
eral antiderivative of r and represents the collection of all antiderivatives of r.
♦
The same reasoning that allows us to differentiate a vector-valued function
componentwise applies to integrating as well. Recall that the integral of a sum
is the sum of the integrals and also that we can remove constant factors from
integrals. So, given r(t) = x(t)i + y(t)j + z(t)k, it follows that we can integrate
componentwise. Expressed more formally,
2 z
y
-2 -1 1
1 -1
x
2
Figure 9.7.6 The position graph for the function in Activity 9.7.6.
v(0)
θ
(x0 , y0 )
Figure 9.7.7 Projectile motion.
Assume we fire a projectile from a launcher and the only force acting on
the fired object is the force of gravity pulling down on the object. That is, we
assume no effect due to spin, wind, or air resistance. With these assumptions,
the motion of the object will be planar, so we can also assume that the motion
occurs in two-dimensional space. Suppose we launch the object from an initial
position (x0 , y0 ) at an angle θ with the positive x-axis as illustrated in Fig-
ure 9.7.7, and that we fire the object with an initial speed of v0 = |v(0)|, where
v(t) is the velocity vector of the object at time t. Assume g is the positive
constant acceleration force due to gravity, which acts to pull the fired object
toward the ground (in the negative y direction). Note particularly that there
is no external force acting on the object to move it in the x direction.
We first observe that since gravity only acts in the downward direction
and that the acceleration due to gravity is constant, the acceleration vector
is h0, −gi. That is, a(t) = h0, −gi. We may use this fact about acceleration,
together with the initial position and initial velocity in order to fully determine
the position r(t) of the object at time t. In Exercise 9.7.7.17, you can work
through the details to show that the following general formula holds.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 79
This assumes that the only force acting on the object is the acceleration g
due to gravity.
9.7.6 Summary
• If r is a vector-valued function, then the derivative of r is defined by
r(t + h) − r(t)
r0 (t) = lim
h→0 h
for those values of t at which the limit exists, and is computed compo-
nentwise by the formula
9.7.7 Exercises
The WeBWorK problems are written by many different authors. Some authors
use parentheses when writing vectors, e.g., (x(t), y(t), z(t)) instead of angle
brackets hx(t), y(t), z(t)i. Please keep this in mind when working WeBWorK
exercises.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 80
Answer 1. 5 sin(−5t)
Answer 2. −5 cos(−5t)
Answer 3. 6
Answer 4. (−0.2) sin(−5t)
Answer 5. 0.2 cos(−5t)
Answer 6. 3tt
6. For the given position vectors r(t),
compute the (tangent) velocity vector r0 (t) for the given value of t .
A) Let r(t) = (cos t, sin t).
Then r0 ( π4 )= ( , )?
B) Let r(t) = (t2 , t3 ).
Then r0 (2)= ( , )?
C) Let r(t) = et i + e−2t j + tk.
Then r0 (1)= i+ j+ k?
Answer 1. −0.707106781186547
Answer 2. 0.707106781186548
Answer 3. 4
Answer 4. 12
Answer 5. 2.71828182845905
Answer 6. −0.270670566473225
Answer 7. 1
7. Suppose ~r(t) = cos(πt) i + sin(πt) j + 3tk represents the position of a
particle on a helix, where z is the height of the particle.
(a) What is t when the particle has height 6?
t=
(b) What is the velocity of the particle when its height is 6?
~v =
(c) When the particle has height 6, it leaves the helix and moves along
the tangent line at the constant velocity found in part (b). Find a vector
parametric equation for the position of the particle (in terms of the original
parameter t) as it moves along this tangent line.
L(t) =
Answer 1. 2
Answer 2. 3.14159j + 3k
Answer 3. i + 6k + (t − 2) · 3.14159j + 3k
8. Suppose the displacement of a particle in motion at time t is given by the
parametric equations
2
x(t) = (3t − 1) , y(t) = 7, z(t) = 54t3 − 27t2 .
(a) Find the speed of the particle when t = 3.
Speed =
(b) Find t when the particle is stationary.
t=
Answer 1. 1296.89
Answer 2. 0.333333
9. Find the derivative of the vector function
r(t) = ta × (b + tc), where
a = h−2, 4, 4i, b = h4, −4, −3i, and c = h3, −3, −5i.
r0 (t) = h , ,
i
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 82
Answer 1. 4 + 2t(−8)
Answer 2. 10 + 2t · 2
Answer 3. −8 + 2t(−6)
10. Let c1 (t) = (e2t , sin(5t), −2t3 ), and c2 (t) = (et , cos(4t), 2t3 )
d
[c1 (t) · c2 (t)] =
dt
d
[c1 (t) × c2 (t)] = i+
dt
j+
k
Answer 1. 3 exp(3t) + 5 cos(5t) cos(4t) − 4 sin(5t) sin(4t) + (−24)t5
Answer 2. 6t2 sin(5t) + 10t3 cos(5t) + (−8)t3 sin(4t) − (−6)t2 cos(4t)
Answer 3. (−6)t2 exp(1t) + (−2)t3 exp(1t) − 4t3 exp(2t) − 6t2 exp(2t)
Answer 4. −4 exp(2t) sin(4t) + 2 exp(2t) cos(4t) − 5 exp(1t) cos(5t) −
1 exp(1t) sin(5t)
11. A gun has a muzzle speed of 90 meters per second. What angle of elevation
should be used to hit an object 180 meters away? Neglect air resistance
and use g = 9.8 m/sec2 as the acceleration of gravity.
Answer:
radians
Answer. No correct answer specified
12. A child wanders slowly down a circular staircase from the top of a tower.
With x, y, z in feet and the origin at the base of the tower, her position t
minutes from the start is given by
and D π E
w(s) = s2 , sin s ,s .
2
a. Determine the point of intersection of the curves generated by r(t)
and w(s). To do so, you will have to find values of a and b that
result in r(a) and w(b) being the same vector.
b. Use the value of a you determined in (a) to find a vector form of the
tangent line to r(t) at the point where t = a.
c. Use the value of b you determined in (a) to find a vector form of the
tangent line to w(s) at the point where s = b.
d. Suppose that z = f (x, y) is a function that generates a surface in
three-dimensional space, and that the curves generated by r(t) and
w(s) both lie on this surface. Note particularly that the point of
intersection you found in (a) lies on this surface. In addition, observe
that the two tangent lines found in (b) and (c) both lie in the tangent
plane to the surface at the point of intersection. Use your preceding
work to determine the equation of this tangent plane.
15. In this exercise, we determinep the equation of a plane tangent to the
surface defined by f (x, y) = x2 + y 2 at the point (3, 4, 5).
a. Find a parameterization for the x = 3 trace of f . What is a direction
vector for the line tangent to this trace at the point (3, 4, 5)?
b. Find a parameterization for the y = 4 trace of f . What is a direction
vector for the line tangent to this trace at the point (3, 4, 5)?
c. The direction vectors in parts (a) and (b) form a plane containing
the point (3, 4, 5). What is a normal vector for this plane?
d. Use your work in parts (a), (b), and (c) to determine an equation
for the tangent plane. Then, use appropriate technology to draw
the graph of f and the plane you determined on the same set of
axes. What do you observe? (We will discuss tangent planes in
more detail in Chapter 10.)
R
16. For each given function r, determine r(t) dt. In addition, recalling the
Fundamental
R 1 Theorem of Calculus for functions of a single variable, also
evaluate 0 r(t) dt for each given function r. Is the resulting quantity a
scalar or a vector?
D What does
E it measure?
1
a. r(t) = cos(t), t+1 , tet
17. In this exercise, we develop the formula for the position function of a
projectile that has been launched at an initial speed of |v0 | and a launch
angle of θ. Recall that a(t) = h0, −gi is the constant acceleration of the
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 84
c. Next, find the specific velocity vector function v for the projectile.
That is, combine your work in (a) and (b) in order to determine
expressions in terms of |v0 | and θ for the constants that arose when
integrating.
d. Find all possible position vectors for the velocity vector v(t) you
determined in (c).
e. Let r(t) denote the position vector function for the given projectile.
Use the fact that the object is fired from the position (x0 , y0 ) to
show it follows that
D g E
r(t) = |v0 | cos(θ)t + x0 , − t2 + |v0 | sin(θ)t + y0 .
2
18. A central force is one that acts on an object so that the force F is parallel to
the object’s position r. Since Newton’s Second Law says that an object’s
acceleration is proportional to the force exerted on it, the acceleration a
of an object moving under a central force will be parallel to its position r.
For instance, the Earth’s acceleration due to the gravitational force that
the sun exerts on the Earth is parallel to the Earth’s position vector as
shown in Figure 9.7.8.
Earth
r
Sun
4 4 4
2 2 2
-1 -1 -1
y y y
-1 1 -1 1 -1 1
x x x
1 Estimate: 8.15 1 Estimate: 8.69 1 Estimate: 8.80
a. Write a formula for the length of the line segment that connects the
endpoints of the curve on the ith subinterval [ti−1 , ti ]. (This length is
our approximation of the length of the curve on this interval.)
b. Use your formula in part (a) to write a sum that adds all of the approx-
imations to the lengths on each subinterval.
c. What do we need to do with the sum in part (b) in order to obtain the
exact value of the length of the graph of r(t) on the interval [0, 2π]?
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 86
We now want to take the limit of this sum as n goes to infinity, but in
its present form it might be difficult to see how. We first introduce ∆t by
multiplying by ∆t
∆t , and see that
n p
X
L≈ (x(ti ) − x(ti−1 ))2 + (y(ti ) − y(ti−1 ))2 + (z(ti ) − z(ti−1 ))2
i=1
n p
X ∆t
= (x(ti ) − x(ti−1 ))2 + (y(ti ) − y(ti−1 ))2 + (z(ti ) − z(ti−1 ))2
i=1
∆t
n p
X ∆t
= (x(ti ) − x(ti−1 ))2 + (y(ti ) − y(ti−1 ))2 + (z(ti ) − z(ti−1 ))2 p
(∆t)2
i=1
To get the difference quotients under the radical, we use properties of the
square root function to see further that
n
s
X 1
L≈ [(x(ti ) − x(ti−1 ))2 + (y(ti ) − y(ti−1 ))2 + (z(ti ) − z(ti−1 )2 ] 2
∆t
i=1
(∆t)
s
n 2 2 2
X x(ti ) − x(ti−1 ) y(ti ) − y(ti−1 ) z(ti ) − z(ti−1 )
= + + ∆t.
i=1
∆t ∆t ∆t
is equal to
Z b p
(x0 (t))2 + (y 0 (t))2 + (z 0 (t))2 dt.
a
Noting further that
p
|r0 (t)| = (x0 (t))2 + (y 0 (t))2 + (z 0 (t))2 ,
If r(t) defines a smooth curve C on an interval [a, b], then the length L
of C is given by
Z b
L= |r0 (t)| dt. (9.8.1)
a
Note that formula (9.8.1) applies to curves in any dimensional space. More-
over, this formula has a natural interpretation: if r(t) records the position of a
moving object, then r0 (t) is the object’s velocity and |r0 (t)| its speed. Formula
(9.8.1) says that we simply integrate the speed of an object traveling over the
curve to find the distance traveled by the object, which is the same as the
length of the curve, just as in one-variable calculus.
Activity 9.8.2 Here we calculate the arc length of two familiar curves.
a. Use Equation (9.8.1) to calculate the circumference of a circle of radius
r.
b. Find the exact length of the spiral defined by r(t) = hcos(t), sin(t), ti on
the interval [0, 2π].
We can adapt the arc length formula to curves in 2-space that define y as
a function of x as the following activity shows.
Activity 9.8.3 Let y = f (x) define a smooth curve in 2-space. Parameterize
this curve and use Equation (9.8.1) to show that the length of the curve defined
by f on an interval [a, b] is
Z b p
1 + [f 0 (t)]2 dt.
a
2 y 2 y
t =2.0
s =2.5
1 t =1.5 1 s =2.0
s =1.5
t =1.0
s =1.0
t =0.5 x s =0.5 x
t =0.0 1 2 s =0.0 1 2
which measures the length along the curve from r(a) to r(t).
The Fundamental Theorem of Calculus shows us that
ds p
= L0 (t) = (x0 (t))2 + (y 0 (t))2 + (z 0 (t))2 = |r0 (t)| (9.8.2)
dt
and so Z t
d
L(t) = r(w) dw.
a dw
If we assume that r0 (t) is never 0, then L0 (t) > 0 for all t and s = L(t) is
always increasing. This should seem reasonable: unless we stop, the distance
traveled along the curve increases as we move along the curve.
Since s = L(t) is an increasing function, it is invertible, which means we
may view the time t as a function of the distance traveled; that is, we have
the relationship t = L−1 (s). We then obtain the arc length parametrization
by composing r(t) with t = L−1 (s) to obtain r(s). Let’s illustrate this with an
example.
Example 9.8.3 Consider a circle of radius 5 in 2-space centered at the origin.
We know that we can parameterize this circle as
where t runs from 0 to 2π. We see that r0 (t) = h−5 sin(t), 5 cos(t)i, and hence
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 89
dr dr ds dr 0
= = |r (t)|,
dt ds dt ds
so
dr 1 dr
= 0 = 1,
ds |r (t)| dt
which means that we move along the curve with unit speed when we parame-
terize by arc length. This is clearly seen in Example 9.8.3 where |r0 (s)| = 1. It
follows that the parameter s is the distance traveled along the curve, as shown
by: Z s Z s
d
L(s) = r(w) dw = 1 dw = s.
0 ds 0
Activity 9.8.4 In this activity we parameterize a line in 2-space in terms of
arc length. Consider the line with parametric equations
b. Use the formula from (a) for s in terms of t to write t in terms of s. Then
explain why a parameterization of the line in terms of arc length is
a b
x(s) = x0 + √ s and y(s) = y0 + √ s. (9.8.4)
a2
+b 2 a + b2
2
t:
Z tp
s(t) = (x0 (w))2 + (y 0 (w))2 + (z 0 (w))2 dw
0
Z tq
= (2w)2 + (4w1/2 )2 + (4)2 dw
0
Z tp
= 4w2 + 16w + 16 dw
0
Z tp
=2 (w + 2)2 dw
0
Z t
=2 w + 2 dw
0
t
= w2 + 4w
0
= t2 + 4t.
These examples illustrate a general method. Of course, evaluating an arc
length integral and finding a formula for the inverse of a function can be diffi-
cult, so while this process is theoretically possible, it is not always practical to
parameterize a curve in terms of arc length. However, we can guarantee that
such a parameterization exists, and this observation plays an important role
in the next section.
9.8.3 Curvature
For a smooth space curve, the curvature measures how fast the curve is bending
or changing direction at a given point. For example, we expect that a line
should have zero curvature everywhere, while a circle (which is bending the
same at every point) should have constant curvature. Circles with larger radii
should have smaller curvatures.
To measure the curvature, we first need to describe the direction of the
curve at a point. We may do this using a continuously varying tangent vector
to the curve, as shown at left in Figure 9.8.5. The direction of the curve is
then determined by the angle φ each tangent vector makes with a horizontal
vector, as shown at right in Figure 9.8.5.
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 91
φ
1 1
φ
-2 -1 1 2 -2 -1 1 2
φ
-1 -1
dT
κ = κ(s) = .
ds
♦
Note that κ is the Greek lowercase letter “kappa”.
Activity 9.8.5
a. We should expect that the curvature of a line is 0 everywhere. To show
that our definition of curvature measures this correctly in 2-space, recall
that (9.8.4) gives us the arc length parameterization
a b
x(s) = x0 + √ s and y(s) = y0 + √ s
a2 + b2 a2 + b2
of a line. Use this information to explain why the curvature of a line is 0
everywhere.
b. Recall that an arc length parameterization of a circle in 2-space of radius
a centered at the origin is, from (9.8.3),
D s s E
r(s) = a cos , a sin .
a a
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 92
Show that the curvature of this circle is the constant a1 . What can you
say about the relationship between the size of the radius of a circle and
the value of its curvature? Why does this make sense?
The definition of curvature relies on our ability to parameterize curves in
terms of arc length. Since we have seen that finding an arc length parametriza-
tion can be difficult, we would like to be able to express the curvature in terms
of a more general parametrization r(t).
To begin, we need to describe the vector T, which is a vector tangent to
the curve having unit length. Of course, the velocity vector r0 (t) is tangent to
the curve; we simply need to normalize its length to be one. This means that
we may take
r0 (t)
T(t) = 0 . (9.8.5)
|r (t)|
Then the curvature of the curve defined by r is
dT
κ=
ds
dT dt
=
dt ds
dT
dt
= ds
dt
0
|T (t)|
= .
|r0 (t)|
Activity 9.8.6 Use one of the two formulas for κ in terms of t to help you
answer the following questions.
x2 y2
a. The ellipse a2 + b2 = 1 has parameterization
Find the curvature of the ellipse. Assuming 0 < b < a, at what points
is the curvature the greatest and at what points is the curvature the
smallest? Does this agree with your intuition?
b. The standard helix has parameterization r(t) = cos(t)i+sin(t)j+tk. Find
the curvature of the helix. Does the result agree with your intuition?
The curvature has another interpretation. Recall that the tangent line to
a curve at a point is the line that best approximates the curve at that point.
The curvature at a point on a curve describes the circle that best approximates
the curve at that point. Remembering that a circle of radius a has curvature
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 93
1/a, then the circle that best approximates the curve near a point on a curve
whose curvature is κ has radius 1/κ and will be tangent to the tangent line
at that point and has its center on the concave side of the curve. This circle,
called the osculating circle of the curve at the point, is shown in Figure 9.8.7
for a portion of a parabola.
4 y
x
-2 -1 1 2
9.8.4 Summary
• The integration process shows that the length L of a smooth curve defined
by r(t) on an interval [a, b] is
Z b
L= |r0 (t)| dt.
a
9.8.5 Exercises
The WeBWorK problems are written by many different authors. Some authors
use parentheses when writing vectors, e.g., (x(t), y(t), z(t)) instead of angle
brackets hx(t), y(t), z(t)i. Please keep this in mind when working WeBWorK
exercises.
1. Find the length of the curve
x = 5 + 5t, y = 3 + 4t, z = t − 4,
for 2 ≤ t ≤ 3.
length =
(Think of second way that you could calculate this length, too, and see
that you get the same result.)
CHAPTER 9. MULTIVARIABLE AND VECTOR FUNCTIONS 94
√
Answer. 5 · 5 + 4 · 4 + 1 · 1(3 − 2)
2. Consider the curve r = (e−5t cos(2t), e−5t sin(2t), e−5t ).
Compute the arclength function s(t): (with initial point t = 0).
7.34847(exp(−5t)−1)
Answer. −5
3. Find the length of the given curve:
where −4 ≤ t ≤ 1.
Answer. 21.2132034355964
π
4. Find the curvature of y = sin (−2x) at x = 4.
Answer. 4
5. Consider the path r(t) = (10t, 5t2 , 5 ln t) defined for t > 0.
Find the length of the curve between the points (10, 5, 0) and (40, 80, 5 ln(4)).
Answer. 75 + 5 ln(4)
6. Find the curvature κ(t) of the curve r(t) = (−3 sin t) i + (−3 sin t) j +
(1 cos t) k
4.24264
Answer. 1.5
((1 sin(t))2 +2(−3 cos(t))2 )
7. A factory has a machine which bends wire at a rate of 9 unit(s) of curvature
per second. How long does it take to bend a straight wire into a circle of
radius 7?
seconds
Answer. 0.0158730158730159
8. Find the unit tangent vector at the indicated point of the vector function
Answer 7. 0
Answer 8. 56t6
Answer 9. 6t1
Answer 10. (−120)t8
Answer 11. (−6)t1
Answer 12. 56t6
10. Starting from the point (5, 4, −2), reparametrize the curve
x(t) = (5 − 2t, 4 + 2t, −2 + t) in terms of arclength.
y(s) = ( , ,
)
s(−2)
Answer 1. 5 + 3
Answer 2. 4 + s·2
3
Answer 3. −2 + s·13
11. Consider the moving particle whose position at time t in seconds is given
by the vector-valued function r defined by r(t) = 5ti+4 sin(3t)j+4 cos(3t)k.
Use this function to answer each of the following questions.
a. Find the unit tangent vector, T(t), to the space curve traced by r(t)
at time t. Write one sentence that explains what T(t) tells us about
the particle’s motion.
b. Determine the speed of the particle moving along the space curve
with the given parameterization.
c. Find the exact distance traveled by the particle on the time interval
[0, π/3].
d. Find the average velocity of the particle on the time interval [0, π/3].
e. Determine the parameterization of the given curve with respect to
arc length.
12. Let y = f (x) define a curve in the plane. We can consider this curve as a
curve in three-space with z-coordinate 0.
a. Find a parameterization of the form r(t) = hx(t), y(t), z(t)i of the
curve y = f (x) in three-space.
c. For the given function r with unit tangent vector T(t) (from (a)),
determine N(t) = |T01(t)| T0 (t).
d. What geometric properties does N(t) have? That is, how long is
this vector, and how is it situated in comparison to T(t)?
e. Let B(t) = T(t) × N(t), and compute B(t) in terms of your results
in (a) and (c).
f. What geometric properties does B(t) have? That is, how long is
this vector, and how is it situated in comparison to T(t) and N(t)?
g. Sketch a plot of the given helix, and compute and sketch T(π/2),
N(π/2), and B(π/2).
15. In this exercise we verify the curvature formula
a. Explain why
ds
|r0 (t)| = .
dt
r0 (t)
b. Use the fact that T(t) = |r0 (t)| and |r0 (t)| = ds
dt to explain why
ds
r0 (t) = T(t).
dt
|T0 (t)|
e. Finally, use the fact that κ = |r0 (t)| to verify that
γ = r + ρN,
−−→
where r = OP .
e. Finally, use the previous work to find the center of the osculating
circle for f at the point (1, 1). Draw pictures of the curve and the
osculating circle to verify your work.
Chapter 10
Derivatives of Multivariable
Functions
10.1 Limits
Motivating Questions
• What do we mean by the limit of a function f of two variables at a point
(a, b)?
• What techniques can we use to show that a function of two variables does
not have a limit at a point (a, b)?
• What does it mean for a function f of two variables to be continuous at
a point (a, b)?
lim f (x) = L.
x→a
f (x) = 3 − x.
98
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 99
5 y
1
x
-3 -2 -1 1 2 3
-1
Figure 10.1.2 The graph of f (x) = 3 − x.
c. Next, consider
x
g(x) = .
|x|
Complete Table 10.1.3 with values near x = 0, the point at which g is
not defined.
x
Table 10.1.3 Values of g(x) = |x| .
-1 1
-1
x
Figure 10.1.4 The graph of g(x) = |x| .
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 100
f (x, y) = 3 − x − 2y.
1.0
1.5 y
0.5
3
x
x
0.5 1.0 1.5
g. Finally, consider
2xy
g(x, y) = ,
x2 + y 2
which is not defined at (0, 0). Complete Table 10.1.7. Round to three
decinal places.
2xy
Table 10.1.7 Values of g(x, y) = x2 +y 2 .
h. Explain how your results are reflected in Figure 10.1.8. Compare this
limit to the limit in part (b). How are the results similar and how are
they different?
z
x
0
-1
x
0
-1
-1 0 y 1
2xy
Figure 10.1.8 Left: The graph of g(x, y) = x2 +y 2 . Right: A contour
plot.
lim f (x, y) = L.
(x,y)→(a,b)
♦
To investigate the limit of a single variable function, limx→a f (x), we often
consider the behavior of f as x approaches a from the right and from the left.
Similarly, we may investigate limits of two-variable functions, lim(x,y)→(a,b) f (x, y)
by considering the behavior of f as (x, y) approaches (a, b) from various direc-
tions. This situation is more complicated because there are infinitely many
ways in which (x, y) may approach (a, b). In the next activity, we see how it is
important to consider a variety of those paths in investigating whether or not
a limit exists.
Activity 10.1.2 Consider the function f , defined by
y
f (x, y) = p ,
x + y2
2
-1
x
0
-1
-1 0 y 1
y
Figure 10.1.10 The graph of f (x, y) = √ .
x2 +y 2
a. Is f defined at the point (0, 0)? What, if anything, does this say about
whether f has a limit at the point (0, 0)?
b. Values of f (to three decimal places) at several points close to (0, 0) are
shown in Table 10.1.11.
Table 10.1.11 Values of a function f .
x\y −1.000 −0.100 0.000 0.100 1.000
−1.000 −0.707 — 0.000 — 0.707
−0.100 — −0.707 0.000 0.707 —
0.000 −1.000 −1.000 — 1.000 1.000
0.100 — −0.707 0.000 0.707 —
1.000 −0.707 — 0.000 — 0.707
Based on these calculations, state whether f has a limit at (0, 0) and give
an argument supporting your statement. (Hint: The blank spaces in the
table are there to help you see the patterns.)
c. Now we formalize the conjecture from the previous part by considering
what happens if we restrict our attention to different paths. First, we
look at f for points in the domain along the x-axis; that is, we consider
what happens when y = 0. What is the behavior of f (x, 0) as x → 0? If
we approach (0, 0) by moving along the x-axis, what value do we find as
the limit?
d. What is the behavior of f along the line y = x when x > 0; that is, what
is the value of f (x, x) when x > 0? If we approach (0, 0) by moving along
the line y = x in the first quadrant (thus considering f (x, x) as x → 0+ ,
what value do we find as the limit?
e. In general, if lim(x,y)→(0,0) f (x, y) = L, then f (x, y) approaches L as (x, y)
approaches (0, 0), regardless of the path we take in letting (x, y) → (0, 0).
Explain what the last two parts of this activity imply about the existence
of lim(x,y)→(0,0) f (x, y).
f. Shown below in Figure 10.1.12 is a set of contour lines of the function
f . What is the behavior of f (x, y) as (x, y) approaches (0, 0) along any
straight line? How does this observation reinforce your conclusion about
the existence of lim(x,y)→(0,0) f (x, y) from the previous part of this activ-
ity? (Hint: Use the fact that a non-vertical line has equation y = mx for
some constant m.)
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 103
y
Figure 10.1.12 Contour lines of f (x, y) = √ .
x2 +y 2
-1 1
-1
x
Figure 10.1.13 The graph of g(x) = |x| .
If f (x, y) has two different limits as (x, y) approaches (a, b) along two
different paths, then lim(x,y)→(a,b) f (x, y) does not exist.
x2 y
g(x, y) =
x4 + y 2
lim f (x, y) = 0,
(x,y)→(0,0)
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 105
In spite of the fact that these two most recent examples illustrate some
of the complications that arise when studying limits of two-variable functions,
many of the properties that are familiar from our study of single variable
functions hold in precisely the same way.
Properties of Limits.
lim f (x, y)
f (x, y) (x,y)→(a,b)
5. lim = if lim g(x, y) 6= 0.
(x,y)→(a,b) g(x, y) lim g(x, y) (x,y)→(a,b)
(x,y)→(a,b)
We can use these properties and results from single variable calculus to
verify that many limits exist. For example, these properties show that the
function f defined by
The reason for this is that polynomial functions of a single variable have
limits at every point.
10.1.2 Continuity
Recall that a function f of a single variable x is said to be continuous at x = a
provided that the following three conditions are satisfied:
1. f (a) exists,
2. limx→a f (x) exists, and
3. limx→a f (x) = f (a).
Using our understanding of limits of multivariable functions, we can define
continuity in the same way.
Definition 10.1.15 A function f = f (x, y) is continuous at the point (a, b)
provided that
1. f is defined at the point (a, b),
2. lim(x,y)→(a,b) f (x, y) exists, and
♦
For instance, we have seen that the function f defined by f (x, y) = 3x2 y 3 +
2xy 2 − 3x + 1 is continous at every point. And just as with single variable
functions, continuity has certain properties that are based on the properties of
limits.
Properties of continuity.
3. f − g is continuous at (a, b)
4. f g is continuous at (a, b)
f
5. g is continuous at (a, b) if g(a, b) 6= 0
Using these properties, we can apply results from single variable calculus to
decide about continuity of multivariable functions. For example, the coordinate
functions f and g defined by f (x, y) = x and g(x, y) = y are continuous
at every point. We can then use properties of continuity listed to conclude
that every polynomial function in x and y is continuous at every point. For
example, g(x, y) = x2 and h(x, y) = y 3 are continuous functions, so their
product f (x, y) = x2 y 3 is a continuous multivariable function.
10.1.3 Summary
• A function f = f (x, y) has a limit L at a point (a, b) provided that we
can make f (x, y) as close to L as we like by taking (x, y) sufficiently close
(but not equal) to (a, b).
• If (x, y) has two different limits as (x, y) approaches (a, b) along two
different paths, we can conclude that lim(x,y)→(a,b) f (x, y) does not exist.
• Properties similar to those for one-variable functions allow us to conclude
that many limits exist and to evaluate them.
• A function f = f (x, y) is continuous at a point (a, b) in its domain if f
has a limit at (a, b) and
f (a, b) = lim f (x, y).
(x,y)→(a,b)
10.1.4 Exercises
1. Find the limits, if they exist, or type DNE for any which do not exist.
1x2
lim
(x,y)→(0,0) 3x2 + 4y 2
1
Answer 3. 3+4m2
Answer 4. dne
2. Determining the limit of a funtion. In this problem we show that
the function
6x − y
f (x, y) =
x+y
does not have a limit as (x, y) → (0, 0).
(a) Suppose that we consider (x, y) → (0, 0) along the curve y = 4x.
Find the limit in this case:
6x−y
lim x+y =
(x,4x)→(0,0)
(b) Now consider (x, y) → (0, 0) along the curve y = 5x. Find the
limit in this case:
6x−y
lim x+y =
(x,5x)→(0,0)
(c) Note that the results from (a) and (b) indicate that f has no limit
as (x, y) → (0, 0) (be sure you can explain why!).
To show this more generally, consider (x, y) → (0, 0) along the curve
y = mx, for arbitrary m. Find the limit in this case:
6x−y
lim x+y =
(x,mx)→(0,0)
(Be sure that you can explain how this result also indicates that f has
no limit as (x, y) → (0, 0).
6−41
Answer 1. 1+41
6−51
Answer 2. 1+51
6−m
Answer 3. 1+m
3. Show that the function
x3 y
f (x, y) = .
x6+ y3
(Be sure that you are able to explain why the results in (a) and (b)
indicate that f does not have a limit at (0,0)!
Answer 1. 0
Answer 2. 1
4. Find the limit, if it exists, or type N if it does not exist.
3x2
lim =
(x,y)→(0,0) 3x + 5y 2
2
Answer. N
5. Find the limit, if it exists, or type N if it does not exist.
(x + 13y)2
lim 2 =
(x,y)→(0,0) x2 + 13 y 2
Answer. N
6. Find the limit,√if it exists, or type ’DNE’ if it does not exist.
2 2
lim e 4x +3y =
(x,y)→(1,3)
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 108
Answer. 261.848047157352
7. Find the limit, if it exists, or type N if it does not exist.
4xy + 2yz + 2xz
lim =
(x,y,z)→(0,0,0) 16x2 + 4y 2 + 4z 2
Answer. N
8. Find the limit, if it exists, or type N if it does not exist.
2 2
2zex +y
lim =
(x,y,z)→(5,1,4) 5x2 + y 2 + 4z 2
Answer. 8241246712.79321
9. Find the limit (enter ’DNE’ if the limit does not exist)
Hint: rationalize the denominator.
(9x2 + 2y 2 )
lim p
(x,y)→(0,0) (9x2 + 2y 2 + 1) − 1
Answer. 2
10. The largest set on which the function f (x, y) = 1/(3−x2 −y 2 ) is continuous
is
A. All of the xy-plane
B. The interior of the circle x2 + y 2 = 3
C. The exterior of the circle x2 + y 2 = 3
Answer. E
xy
11. Consider the function f defined by f (x, y) = x2 +y 2 +1 .
a. What is the domain of f ?
2x2 y
13. Consider the function h defined by h(x, y) = x4 +y 2 .
a. What is the domain of h?
b. Evaluate the limit of h at (0, 0) along all linear paths the contain the
origin. What does this tell us about lim(x,y)→(0,0) h(x, y)? (Hint: A
non-vertical line throught the origin has the form y = mx for some
constant m.)
c. Does lim(x,y)→(0,0) h(x, y) exist? Verify your answer. Check by using
appropriate technology to sketch both surface and contour plots of
h near (0, 0). Write several sentences to say how your plots affirm
your findings about lim(x,y)→(0,0) h(x, y).
14. For each of the following prompts, provide an example of a function of
two variables with the desired properties (with justification), or explain
why such a function does not exist.
a. A function p that is defined at (0, 0), but lim(x,y)→(0,0) p(x, y) does
not exist.
b. A function q that does not have a limit at (0, 0), but that has the
same limiting value along any line y = mx as x → 0.
c. A function r that is continuous at (0, 0), but lim(x,y)→(0,0) r(x, y)
does not exist.
d. A function s such that
sin(x)
b. The function g defined by g(x, y) = 1+ey
750 750
500 500
250 250
t r
1 2 3 4 5 6 7 8 9 10 0.02 0.04 0.06 0.08 0.10
Figure 10.2.1 Left: Graphs f (t) = M (0.03, t). Right: Graph g(r) =
M (r, 4).
c. Find the instantaneous rate of change f 0 (4) and state the units on this
quantity. What information does f 0 (4) tell us about our car loan? What
information does f 0 (4) tell us about the graph you sketched in (b)?
d. Express M as a function of r alone, using a fixed time of t = 4. That is,
let g(r) = M (r, 4). Sketch the graph of g on the right of Figure 10.2.1.
Explain the meaning of the function g.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 111
e. Find the instantaneous rate of change g 0 (0.03) and state the units on
this quantity. What information does g 0 (0.03) tell us about our car loan?
What information does g 0 (0.03) tell us about the graph you sketched in
(d)?
x2 sin(2y)
f (x, y) = ,
32
which measures the range, or horizontal distance, in feet, traveled by a projec-
tile launched with an initial speed of x feet per second at an angle y radians
to the horizontal. The graph of this function is given again on the left in
Figure 10.2.2. Moreover, if we fix the angle y = 0.6, we may view the trace
f (x, 0.6) as a function of x alone, as seen at right in Figure 10.2.2.
z 1000
f (x, 0.6)
1500
1000 800
500
600
x
200
400
y 150
1.5 100
1.0
0.5
50 200
0
x
50 100 150 200
x2 sin(2y)
Figure 10.2.2 Left: The trace of z = 32 with y = 0.6.
Since the trace is a one-variable function, we may consider its derivative
just as we did in the first semester of calculus. With y = 0.6, we have
sin(1.2) 2
f (x, 0.6) = x ,
32
and therefore
d sin(1.2)
[f (x, 0.6)] = x.
dx 16
When x = 150, this gives
d sin(1.2)
[f (x, 0.6)]|x=150 = 150 ≈ 8.74 feet per feet per second,
dx 16
which gives the slope of the tangent line shown on the right of Figure 10.2.2.
Thinking of this derivative as an instantaneous rate of change implies that if
we increase the initial speed of the projectile by one foot per second, we expect
the horizontal distance traveled to increase by approximately 8.74 feet if we
hold the launch angle constant at 0.6 radians.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 112
1000 800
500
600
x
200
400
y 150
1.5 100
1.0
0.5
50 200
0
y
0.25 0.50 0.75 1.00 1.25
x2 sin(2y)
Figure 10.2.3 The trace of z = 32 with x = 150.
This gives
1502
f (150, y) = sin(2y),
32
and therefore
d 1502
[f (150, y)] = cos(2y).
dy 16
If we evaluate this quantity at y = 0.6, we have
d 1502
[f (150, y)]|y=0.6 = cos(1.2) ≈ 509.5 feet per radian.
dy 16
Once again, the derivative gives the slope of the tangent line shown on the
right in Figure 10.2.3. Thinking of the derivative as an instantaneous rate of
change, we expect that the range of the projectile increases by 509.5 feet for
every radian we increase the launch angle y if we keep the initial speed of the
projectile constant at 150 feet per second.
By holding x fixed and differentiating with respect to y, we obtain the first-
order partial derivative of f with respect to y. As before, we denote this partial
derivative as fy and write
d f (150, 0.6 + h) − f (150, 0.6)
fy (150, 0.6) = f (150, y)|y=0.6 = lim .
dy h→0 h
As with the partial derivative with respect to x, we may express this quan-
tity more generally at an arbitrary point (a, b). To recap, we have now arrived
at the formal definition of the first-order partial derivatives of a function of
two variables.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 113
f (a + h, b) − f (a, b)
fx (a, b) = lim , and
h→0 h
f (a, b + h) − f (a, b)
fy (a, b) = lim ,
h→0 h
provided the limits exist. ♦
Activity 10.2.2 Consider the function f defined by
xy 2
f (x, y) =
x+1
at the point (1, 2).
a. Write the trace f (x, 2) at the fixed value y = 2. On the left side of
Figure 10.2.5, draw the graph of the trace with y = 2 around the point
where x = 1, indicating the scale and labels on the axes. Also, sketch
the tangent line at the point x = 1.
xy 2
Figure 10.2.5 Traces of f (x, y) = x+1 .
b. Find the partial derivative fx (1, 2) and relate its value to the sketch you
just made.
c. Write the trace f (1, y) at the fixed value x = 1. On the right side of
Figure 10.2.5, draw the graph of the trace with x = 1 indicating the scale
and labels on the axes. Also, sketch the tangent line at the point y = 2.
d. Find the partial derivative fy (1, 2) and relate its value to the sketch you
just made.
As these examples show, each partial derivative at a point arises as the
derivative of a one-variable function defined by fixing one of the coordinates. In
addition, we may consider each partial derivative as defining a new function of
the point (x, y), just as the derivative f 0 (x) defines a new function of x in single-
variable calculus. Due to the connection between one-variable derivatives and
partial derivatives, we will often use Leibniz-style notation to denote partial
derivatives by writing
∂f ∂f
(a, b) = fx (a, b), and (a, b) = fy (a, b).
∂x ∂y
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 114
xy 2
b. If f (x, y) = , find the partial derivatives fx and fy .
x+1
c. If g(r, s) = rs cos(r), find the partial derivatives gr and gs .
d. Assuming f (w, x, y) = (6w + 1) cos(3x2 + 4xy 3 + y), find the partial
derivatives fw , fx , and fy .
x2t z 3
e. Find all possible first-order partial derivatives of q(x, t, z) = .
1 + x2
1000 1000
500 500
x x
200 200
y 150 y 150
1.5 100 1.5 100
1.0 50 1.0 50
0.5 0.5
0 0
a. Estimate the partial derivative wv (20, −10). What are the units on this
quantity and what does it mean? (Recall that we can estimate a partial
derivative of a single variable function f using the symmetric difference
quotient f (x+h)−f
2h
(x−h)
for small values of h. A partial derivative is a
derivative of an appropriate trace.)
b. Estimate the partial derivative wT (20, −10). What are the units on this
quantity and what does it mean?
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 116
c. Use your results to estimate the wind chill w(18, −10). (Recall from single
variable calculus that for a function f of x, f (x + h) ≈ f (x) + hf 0 (x).)
2
-1
1
0
x
1
-3 -2 -1 1 2 3
2
3 -1
4
5 -2
6
-3
a. Estimate the partial derivative fx (−2, −1). (Hint: How can you find
values of f that are of the form f (−2 + h) and f (−2 − h) so that you can
use a symmetric difference quotient?)
4 y 4 y
3 3
2 2
1 1
x x
1 2 3 4 1 2 3 4
Figure 10.2.9 Plots for contours of g and h.
h. Suppose you have yet another function h, and you know that h(2, 2) = 4,
hx (2, 2) < 0, and hy (2, 2) > 0. Using this information, sketch a possi-
ble contour h(x, y) = 4 passing through (2, 2) on the right side of Fig-
ure 10.2.9. Then include possible contours h(x, y) = 3 and h(x, y) = 5.
10.2.4 Summary
• If f = f (x, y) is a function of two variables, there are two first order
partial derivatives of f : the partial derivative of f with respect to x,
∂f f (x + h, y) − f (x, y)
(x, y) = fx (x, y) = lim ,
∂x h→0 h
and the partial derivative of f with respect to y,
∂f f (x, y + h) − f (x, y)
(x, y) = fy (x, y) = lim ,
∂y h→0 h
where each partial derivative exists only at those points (x, y) for which
the limit exists.
• The partial derivative fx (a, b) tells us the instantaneous rate of change of
f with respect to x at (x, y) = (a, b) when y is fixed at b. Geometrically,
the partial derivative fx (a, b) tells us the slope of the line tangent to the
y = b trace of the function f at the point (a, b, f (a, b)).
• The partial derivative fy (a, b) tells us the instantaneous rate of change of
f with respect to y at (x, y) = (a, b) when x is fixed at a. Geometrically,
the partial derivative fy (a, b) tells us the slope of the line tangent to the
x = a trace of the function f at the point (a, b, f (a, b)).
10.2.5 Exercises
1. Find the first partial derivatives of
x − 4y
f (x, y) = at the point (x, y) = (4, 1).
x + 4y
∂f
(4, 1) =
∂x
∂f
(4, 1) =
∂y
Answer 1. 0.125
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 118
Answer 2. −0.5
2. Find the first partial derivatives of f (x, y) = sin(x − y) at the point (-4,
-4).
A. fx (−4, −4) =
B. fy (−4, −4) =
Answer 1. 1
Answer 2. −1
3. Find the partial derivatives of the function
p
w = 4r2 + 3s2 + 7t2
∂w
∂r =
∂w
∂s =
∂w
∂t =
4r
Answer 1. √
4r 2 +3s2 +7t2
3s
Answer 2. √
4r 2 +3s2 +7t2
7t
Answer 3. √
4r 2 +3s2 +7t2
4. Suppose that f (x, y) is a smooth function and that its partial derivatives
have the values, fx (0, 9) = −4 and fy (0, 9) = −2. Given that f (0, 9) =
1, use this information to estimate the value of f (1, 10). Note this is
analogous to finding the tangent line approximation to a function of one
variable. In fancy terms, it is the first Taylor approximation.
Estimate of (integer value) f (0, 10)
Estimate of (integer value) f (1, 9)
Estimate of (integer value) f (1, 10)
Answer 1. −1
Answer 2. −3
Answer 3. −5
5. The gas law for a fixed mass m of an ideal gas at absolute temperature T ,
pressure P , and volume V is P V = mRT , where R is the gas constant.
Find the partial derivatives
∂P
=
∂V
∂V
=
∂T
∂T
=
∂P
∂P ∂V ∂T
= (an integer)
∂V ∂T ∂P
Answer 1. −mRT V2
mR
Answer 2. P
V
Answer 3. mR
Answer 4. −1
6. Find the first partial derivatives of f (x, y, z) = z arctan( xy ) at the point
(4, 4, -3).
A. ∂f
∂x (4, 4, −3) =
∂f
B. ∂y (4, 4, −3) =
∂f
C. ∂z (4, 4, −3) =
Answer 1. 0.375
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 119
Answer 2. −0.375
Answer 3. 0.785398163397448
7. Find the partial derivatives of the function
Z x
f (x, y) = cos(7t2 + 8t − 2) dt
y
fx (x, y) =
fy (x, y) =
Answer 1. cos 7x2 + 8x + (−2)
Answer 2. − cos 7y 2 + 8y + (−2)
8. Let f (x, y) = e−2x sin(4y).
(a) Using difference quotients with ∆x = 0.1 and ∆y = 0.1, we esti-
mate
fx (2, −2) ≈
fy (2, −2) ≈
(b) Using difference quotients with ∆x = 0.01 and ∆y = 0.01, we find
better estimates:
fx (2, −2) ≈
fy (2, −2) ≈
Answer 1. 0.0328473
Answer 2. 0.00392661
Answer 3. 0.0358814
Answer 4. −0.0092074
9. Determine the sign of fx and fy at each indicated point using the contour
diagram of f shown below. (The point P is that in the first quadrant, at
a positive x and y value; Q through T are located clockwise from P , so
that Q is at a positive x value and negative y, etc.)
(a) At point P ,
fx is ( positive negative) and
fy is ( positive negative) .
(b) At point Q,
fx is ( positive negative) and
fy is ( positive negative) .
(c) At point S,
fx is ( positive negative) and
fy is ( positive negative) .
Answer 1. negative
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 120
Answer 2. negative
Answer 3. negative
Answer 4. positive
Answer 5. positive
Answer 6. negative
10. Your monthly car payment in dollars is P = f (P0 , t, r), where $P0 is the
amount you borrowed, t is the number of months it takes to pay off the
loan, and r percent is the interest rate.
(a) Is ∂P/∂t positive or negative? ( positive negative)
Suppose that your bank tells you that the magnitude of ∂P/∂t is 15.
What are the units of this value?
(For this problem, write our your units in full, writing dollars for
$, months for months, percent for %, etc. Note that fractional units
generally have a plural numerator and singular denominator.)
(b) Is ∂P/∂r positive or negative? ( positive negative)
Suppose that your bank tells you that the magnitude of ∂P/∂r is 25.
What are the units of this value?
(For this problem, write our your units in full, writing dollars for
$, months for months, percent for %, etc. Note that fractional units
generally have a plural numerator and singular denominator.)
For both parts of this problem, be sure you can explain what the prac-
tical meanings of the partial derivatives are.
Answer 1. negative
Answer 2. dollars/month
Answer 3. positive
Answer 4. dollars/percent
11. An experiment to measure the toxicity of formaldehyde yielded the data
in the table below. The values show the percent, P = f (t, c), of rats
surviving an exposure to formaldehyde at a concentration of c (in parts
per million, ppm) after t months.
t = 14 t = 16 t = 18 t = 20 t = 22 t = 24
c=0 100 100 100 99 97 95
c=2 100 99 98 97 95 92
c=6 96 95 93 90 86 80
c = 15 96 93 82 70 58 36
(a) Estimate ft (18, 6):
ft (18, 6) ≈
(b) Estimate fc (18, 6):
fc (18, 6) ≈
(Be sure that you can give the practical meaning of these two values
in terms of formaldehyde toxicity.)
90−95
Answer 1. 20−16
82−98
Answer 2. 15−2
12. An airport can be cleared of fog by heating the air. The amount of heat
required depends on the air temperature and the wetness of the fog. The
figure below shows the heat H(T, w) required (in calories per cubic meter
of fog) as a function of the temperature T (in degrees Celsius) and the
water content w (in grams per cubic meter of fog). Note that this figure
is not a contour diagram, but shows cross-sections of H with w fixed at
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 121
a. State the limit definition of the value IT (94, 75). Then, estimate
IT (94, 75), and write one complete sentence that carefully explains
the meaning of this value, including its units.
b. State the limit definition of the value IH (94, 75). Then, estimate
IH (94, 75), and write one complete sentence that carefully explains
the meaning of this value, including its units.
c. Suppose you are given that IT (92, 80) = 3.75 and IH (92, 80) = 0.8.
Estimate the values of I(91, 80) and I(92, 78). Explain how the
partial derivatives are relevant to your thinking.
d. On a certain day, at 1 p.m. the temperature is 92 degrees and the
relative humidity is 85%. At 3 p.m., the temperature is 96 degrees
and the relative humidity 75%. What is the average rate of change
of the heat index over this time period, and what are the units on
your answer? Write a sentence to explain your thinking.
14. Let f (x, y) = 12 xy 2 represent the kinetic energy in Joules of an object of
mass x in kilograms with velocity y in meters per second. Let (a, b) be
the point (4, 5) in the domain of f .
a. Calculate fx (a, b).
b. Explain as best you can in the context of kinetic energy what the
partial derivative
f (a + h, b) − f (a, b)
fx (a, b) = lim
h→0 h
tells us about kinetic energy.
f (a, b + h) − f (a, b)
fy (a, b) = lim
h→0 h
tells us about kinetic energy.
e. Often we are given certain graphical information about a function
instead of a rule. We can use that information to approximate par-
tial derivatives. For example, suppose that we are given a contour
plot of the kinetic energy function (as in Figure 10.2.11) instead of a
formula. Use this contour plot to approximate fx (4, 5) and fy (4, 5)
as best you can. Compare to your calculations from earlier parts of
this exercise.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 123
8 y
7
6
5
70 80
4 50 60
3 40
30
20
2
10
1
x
1 2 3 4 5 6 7 8
b. If an ant is on the metal plate, standing at the point (2, 3), and starts
walking in the direction parallel to the positive y axis, at what rate
will the temperature the ant is experiencing change? Explain, and
include appropriate units.
c. If an ant is walking along the line y = 3 in the positive x direction, at
what instantaneous rate will the temperature the ant is experiencing
change when the ant passes the point (1, 3)?
d. Now suppose the ant is stationed at the point (6, 3) and walks in a
straight line towards the point (2, 0). Determine the average rate of
change in temperature (per unit distance traveled) the ant encoun-
ters in moving between these two points. Explain your reasoning
carefully. What are the units on your answer?
16. Consider the function f defined by f (x, y) = 8 − x2 − 3y 2 .
a. Determine fx (x, y) and fy (x, y).
b. Find parametric equations in R3 for the tangent line to the trace
f (x, 1) at x = 2.
e. Determine the equation of the plane that passes through the point
(2, 1, f (2, 1)) whose normal vector is orthogonal to the direction vec-
tors of the two lines found in (b) and (c).
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 124
f 0 (x + h) − f 0 (x)
f 00 (x) = lim .
h→0 h
In what follows, we begin exploring the four different second-order partial
derivatives of a function of two variables and seek to understand what these
various derivatives tell us about the function’s behavior.
Preview Activity 10.3.1 Once again, let’s consider the function f defined
2
by f (x, y) = x sin(2y)
32 that measures a projectile’s range as a function of its
initial speed x and launch angle y. The graph of this function, including traces
with x = 150 and y = 0.6, is shown in Figure 10.3.1.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 125
z z
1500 1500
1000 1000
500 500
x x
200 200
y 150 y 150
1.5 100 1.5 100
1.0 50 1.0 50
0.5 0.5
0 0
Figure 10.3.1 The distance function with traces x = 150 and y = 0.6.
800
600
400
200
x
50 100 150 200
c. Determine the partial derivative fy , and then find the partial derivative
fyy = (fy )y . Evaluate fyy (150, 0.6).
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 126
1000
f (150, y)
800
600
400
200
y
0.25 0.50 0.75 1.00 1.25
d. Figure 10.3.3 shows the trace f (150, y) and includes three tangent lines.
Explain how the value of fyy (150, 0.6) is reflected in this figure.
e. Because fx and fy are each functions of both x and y, they each have
two partial derivatives. Not only can we compute fxx = (fx )x , but also
fxy = (fx )y ; likewise, in addition to fyy = (fy )y , but also fyx = (fy )x .
2
For the range function f (x, y) = x sin(2y)
32 , use your earlier computations
of fx and fy to now determine fxy and fyx . Write one sentence to explain
how you calculated these “mixed” partial derivatives.
The first two are called unmixed second-order partial derivatives while the
last two are called the mixed second-order partial derivatives.
One aspect of this notation can be a little confusing. The notation
∂2f
∂ ∂f
=
∂y∂x ∂y ∂x
means that we first differentiate with respect to x and then with respect to y;
this can be expressed in the alternate notation fxy = (fx )y . However, to find
the second partial derivative
fyx = (fy )x
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 127
d. How many second order partial derivatives does the function h defined
by h(x, y, z) = 9x9 z − xyz 9 + 9 have? Find hxz and hzx (you do not need
to find the other second order partial derivatives).
In Preview Activity 10.3.1 and Activity 10.3.2, you may have noticed that
the mixed second-order partial derivatives are equal. This observation holds
generally and is known as Clairaut’s Theorem.
Clairaut’s Theorem.
Let f be a function of several variables for which the partial derivatives
fxy and fyx are continuous near the point (a, b). Then
z z z
6 6 6
4 4 4
2 2 2
-2 1 -2 1 -2 1
y -1 2 x y -1 2 x y -1 2 x
3 3 3
Consider, for example, f (x, y) = sin(x)e−y . Figure 10.3.4 shows the graph
of this function along with the trace given by y = −1.5. Also shown are three
tangent lines to this trace, with increasing x-values from left to right among
the three plots in Figure 10.3.4.
That the slope of the tangent line is decreasing as x increases is reflected,
as it is in one-variable calculus, in the fact that the trace is concave down.
Indeed, we see that fx (x, y) = cos(x)e−y and so fxx (x, y) = − sin(x)e−y < 0,
since e−y > 0 for all values of y, including y = −1.5.
In the following activity, we further explore what second-order partial
derivatives tell us about the geometric behavior of a surface.
Activity 10.3.3 We continue to consider the function f defined by f (x, y) =
sin(x)e−y .
a. In Figure 10.3.5, we see the trace of f (x, y) = sin(x)e−y that has x held
constant with x = 1.75. We also see three different lines that are tangent
to the trace of f in the x direction at values of y that are increasing
from left to right in the figure. Write a couple of sentences that describe
whether the slope of the tangent lines to this curve increase or decrease as
y increases, and, after computing fyy (x, y), explain how this observation
is related to the value of fyy (1.75, y). Be sure to address the notion of
concavity in your response.(You need to be careful about the directions
in which x and y are increasing.)
z z z
6 6 6
4 4 4
2 2 2
-2 1 -2 1 -2 1
y -1 2 x y -1 2 x y -1 2 x
3 3 3
6 6 6
4 4 4
2 2 2
-2 1 -2 1 -2 1
y -1 2 x y -1 2 x y -1 2 x
3 3 3
d. We know that fxx (1.75, −1.5) measures the concavity of the y = −1.5
trace, and that fyy (1.75, −1.5) measures the concavity of the x = 1.75
trace. What do you think the quantity fxy (1.75, −1.5) measures?
e. On Figure 10.3.6, sketch the trace with y = −1.5, and sketch three tan-
gent lines whose slopes correspond to the value of fyx (x, −1.5) for three
different values of x, the middle of which is x = −1.5. Is fyx (1.75, −1.5)
positive or negative? Why? What does fyx (1.75, −1.5) measure?
Just as with the first-order partial derivatives, we can approximate second-
order partial derivatives in the situation where we have only partial information
about the function.
Activity 10.3.4 As we saw in Activity 10.2.5, the wind chill w(v, T ), in degrees
Fahrenheit, is a function of the wind speed, in miles per hour, and the air
temperature, in degrees Fahrenheit. Some values of the wind chill are recorded
in Table 10.3.7.
Table 10.3.7 Wind chill as a function of wind speed and temperature.
v\T -30 -25 -20 -15 -10 -5 0 5 10 15 20
5 -46 -40 -34 -28 -22 -16 -11 -5 1 7 13
10 -53 -47 -41 -35 -28 -22 -16 -10 -4 3 9
15 -58 -51 -45 -39 -32 -26 -19 -13 -7 0 6
20 -61 -55 -48 -42 -35 -29 -22 -15 -9 -2 4
25 -64 -58 -51 -44 -37 -31 -24 -17 -11 -4 3
30 -67 -60 -53 -46 -39 -33 -26 -19 -12 -5 1
35 -69 -62 -55 -48 -41 -34 -27 -21 -14 -7 0
40 -71 -64 -57 -50 -43 -36 -29 -22 -15 -8 -1
a. Estimate the partial derivatives wT (20, −15), wT (20, −10), and wT (20, −5).
Use these results to estimate the second-order partial wT T (20, −10).
b. In a similar way, estimate the second-order partial wvv (20, −10).
c. Estimate the partial derivatives wT (20, −10), wT (25, −10), and wT (15, −10),
and use your results to estimate the partial wT v (20, −10).
10.3.3 Summary
• There are four second-order partial derivatives of a function f of two
independent variables x and y:
fxx = (fx )x , fxy = (fx )y , fyx = (fy )x , and fyy = (fy )y .
• The unmixed second-order partial derivatives, fxx and fyy , tell us about
the concavity of the traces. The mixed second-order partial derivatives,
fxy and fyx , tell us how the graph of f twists.
10.3.4 Exercises
1. Calculate all four second-order partial derivatives of f (x, y) = 4x2 y+8xy 3 .
fxx (x, y) =
fxy (x, y) =
fyx (x, y) =
fyy (x, y) =
Answer 1. 8y
Answer 2. 4 · 2x + 8 · 3y 2
Answer 3. 4 · 2x + 8 · 3y 2
Answer 4. 8x · 3 · 2y
2. Find all the first and second order partial derivatives of f (x, y) = 3 sin(2x+
y) − 4 cos(x − y).
A. ∂f
∂x = fx =
∂f
B. ∂y = fy =
∂2f
C. ∂x2 = fxx =
∂2f
D. ∂y 2 = fy y =
∂2f
E. ∂x∂y = fy x =
∂2f
F. ∂y∂x = fxy =
Answer 1. 6 cos(2x + y) − (−4) sin(x − y)
Answer 2. 3 cos(2x + y) + (−4) sin(x − y)
Answer 3. −4 · 3 sin(2x + y) − (−4) cos(x − y)
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 131
f (x, y) = xye4y
fx (x, y) =
fy (x, y) =
fxy (x, y) =
fyx (x, y) =
Answer 1. y exp(4y)
Answer 2. x(4y exp(4y) + exp(4y))
Answer 3. 4y exp(4y) + exp(4y)
Answer 4. 4y exp(4y) + exp(4y)
5x
4. Calculate all four second-order partial derivatives of f (x, y) = sin .
y
fxx (x, y) =
fxy (x, y) =
fyx (x, y) =
fyy (x, y) =
Answer 1. − 5y 5y
y 2 y 2 sin y
5x
2
Answer 2. 5y (y−5y·2y2) 2 cos 5x
y + 5y 5x
2
y y 2 sin 5x
y
2
5y 5y
Answer 3. − (y2 )2 cos 5x y − 5x
y 2 y 2 sin y
5x
Answer 4. − 5x 5x
y 2 y 2 sin y
5x
− 5x·2y
(y 2 )2
cos 5x y
4 5
5. Given F (r, s, t) = r 9s − t , compute:
Frst =
Answer. 0
6. Calculate all four second-order partial derivatives and check that fxy =
fyx . Assume the variables are restricted to a domain on which the function
is defined.
f (x, y) = e2xy
fxx =
fyy =
fxy =
fyx =
Answer 1. 4y 2 e2xy
Answer 2. 4x2 e2xy
Answer 3. 2e2xy + 4xye2xy
Answer 4. 2e2xy + 4xye2xy
7. Calculate all four second-order partial derivatives of f (x, y) = (2x + 4y) ey .
fxx (x, y) =
fxy (x, y) =
fyx (x, y) =
fyy (x, y) =
Answer 1. 0
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 132
Answer 2. 2ey
Answer 3. 2ey
Answer 4. 4ey + 4ey + (2x + 4y) ey
6
8. Let f (x, y) = (− (2x + y)) . Then
∂ 2f
∂x∂y =
∂ 3f
∂x∂y∂x =
∂ 3f
∂x2 ∂y =
4
Answer 1. 60(− (2x + y))
3
Answer 2. −480(− (2x + y))
3
Answer 3. −480(− (2x + y))
9. If zxy = 5y and all of the second order partial derivatives of z are contin-
uous, then
(a) zyx =
(b) zxyx =
(c) zxyy =
Answer 1. 5y
Answer 2. 0
Answer 3. 5
10. If z = f (x) + yg(x), what can we say about zyy ?
zyy = 0
zyy = y
zyy = zxx
zyy = g(x)
y
5
4
4
3
2
2
1
x
2 4
4 y 4 y
3 3
2 2
1 1
x x
1 2 3 4 1 2 3 4
Figure 10.3.10 Plots for contours of g and h.
a. State the limit definition of the value IT T (94, 75). Then, estimate
IT T (94, 75), and write one complete sentence that carefully explains
the meaning of this value, including units.
b. State the limit definition of the value IHH (94, 75). Then, estimate
IHH (94, 75), and write one complete sentence that carefully explains
the meaning of this value, including units.
c. Finally, do likewise to estimate IHT (94, 75), and write a sentence to
explain the meaning of the value you found.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 134
13. The temperature on a heated metal plate positioned in the first quadrant
of the xy-plane is given by
2
−(y−1)3
C(x, y) = 25e−(x−1) .
c. Calculate Cyy (1.1, 1.2). Suppose instead that an ant is walking past
the point (1.1, 1.2) along the line x = 1.1. Write a sentence to
explain the meaning of the value of Cyy (1.1, 1.2), including units.
d. Determine Cxy (x, y) and hence compute Cxy (1.1, 1.2). What is the
meaning of this value? Explain, in terms of an ant walking on the
heated metal plate.
14. Let f (x, y) = 8 − x2 − y 2 and g(x, y) = 8 − x2 + 4xy − y 2 .
a. Determine fx , fy , fxx , fyy , fxy , and fyx .
b. Evaluate each of the partial derivatives in (a) at the point (0, 0).
c. What do the values in (b) suggest about the behavior of f near
(0, 0)? Plot a graph of f and compare what you see visually to what
the values suggest.
d. Determine gx , gy , gxx , gyy , gxy , and gyx .
e. Evaluate each of the partial derivatives in (d) at the point (0, 0).
x2
f (x, y) = 6 − − y2 ,
2
whose graph is shown in Figure 10.4.1.
6
4
z
2
-2
x
-2 y 2
6 6
4 4
z z
2 2
0 0.5
x x
0 y 2 0.5 y 1.5
Figure 10.4.2 The graph of f (x, y) = 6 − x2 /2 − y 2 .
Just as the graph of a differentiable single-variable function looks like a
line when viewed on a small scale, we see that the graph of this particular two-
variable function looks like a plane, as seen in Figure 10.4.3. In the following
preview activity, we explore how to find the equation of this plane.
6
4
z
2
0.5
x
0.5 y 1.5
graph of the function f contain the point (x0 , y0 , z0 ). Use this observation
to determine z0 in the expression z = z0 + a(x − x0 ) + b(y − y0 ).
x2
c. Sketch the traces of f (x, y) = 6 − 2 − y 2 for y = y0 = 1 and x = x0 = 1
below in Figure 10.4.4.
5.0 5.0
z = f (x, 1) z = f (1, y)
4.5 4.5
x y
4.0 4.0
0.5 1.0 1.5 0.5 1.0 1.5
d. Determine the equation of the tangent line of the trace that you sketched
in the previous part with y = 1 (in the x direction) at the point x0 = 1.
6 6
4 4
z z
2 2
0.5 0.5
x x
0.5 y 1.5 0.5 y 1.5
Figure 10.4.5 The traces of f (x, y) and the tangent plane.
e. Figure 10.4.5 shows the traces of the function and the traces of the tan-
gent plane. Explain how the tangent line of the trace of f , whose equation
you found in the last part of this activity, is related to the tangent plane.
How does this observation help you determine the constant a in the equa-
tion for the tangent plane z = z0 + a(x − x0 ) + b(y − y0 )? (Hint: How do
you think fx (x0 , y0 ) should be related to zx (x0 , y0 )?)
f. In a similar way to what you did in (d), determine the equation of the
tangent line of the trace with x = 1 at the point y0 = 1. Explain how
this tangent line is related to the tangent plane, and use this observation
to determine the constant b in the equation for the tangent plane z =
z0 + a(x − x0 ) + b(y − y0 ). (Hint: How do you think fy (x0 , y0 ) should be
related to zy (x0 , y0 )?)
g. Finally, write the equation z = z0 + a(x − x0 ) + b(y − y0 ) of the tangent
plane to the graph of f (x, y) = 6−x2 /2−y 2 at the point (x0 , y0 ) = (1, 1).
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 138
Differentiablity.
locally linear at (x0 , y0 ) (because the graph of f does not look linear when we
zoom in around the point (x0 , y0 )). In such a case this plane is not tangent to
the graph. Differentiability for a function of two variables implies the existence
of a tangent plane, but the existence of the two first order partial derivatives
of a function at a point does not imply differentiaility. This is quite different
than what happens in single variable calculus.
Finally, one important note about the form of the equation for the tangent
plane, z = f (x0 , y0 ) + fx (x0 , y0 )(x − x0 ) + fy (x0 , y0 )(y − y0 ). Say, for exam-
ple, that we have the particular tangent plane z = 7 − 2(x − 3) + 4(y + 1).
Observe that we can immediately read from this form that fx (3, −1) = −2
and fy (3, −1) = 4; furthermore, fx (3, −1) = −2 is the slope of the trace to
both f and the tangent plane in the x-direction at (−3, 1). In the same way,
fy (3, −1) = 4 is the slope of the trace of both f and the tangent plane in the
y-direction at (3, −1).
Activity 10.4.2
a. Find the equation of the tangent plane to f (x, y) = 2 + 4x − 3y at the
point (1, 2). Simplify as much as possible. Does the result surprise you?
Explain.
10.4.2 Linearization
In single variable calculus, an important use of the tangent line is to approx-
imate the value of a differentiable function. Near the point x0 , the tangent
line to the graph of f at x0 is close to the graph of f near x0 , as shown in
Figure 10.4.6.
y y
y = L(x) y = L(x)
y = f (x) y = f (x)
x x
x0 x0
4
z
2
z = f (x, y)
0.5
x
0.5 y 1.5
Use the data to first estimate the appropriate partial derivatives, and
then find the linearization L(v, T ) at the point (20, −10). Finally, use
the linearization to estimate w(10, −10), w(20, −12), and w(18, −12).
Compare your results to what you obtained in Activity 10.2.5
c. Figure 10.4.9 gives a contour plot of a continuously differentiable function
f.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 141
y
5
4
4
3
2
2
1
x
2 4
10.4.3 Differentials
As we have seen, the linearization L(x, y) enables us to estimate the value of
f (x, y) for points (x, y) near the base point (x0 , y0 ). Sometimes, however, we
are more interested in the change in f as we move from the base point (x0 , y0 )
to another point (x, y).
z
df
∆f
y
(x0 , y0 )
dx = ∆x
dy = ∆y (x, y)
x
A(x, y) = xy,
so that the area of a perfectly manufactured rectangle is A(20, 10) = 200 square
centimeters. Since the machine isn’t perfect, we would like to know how much
the area of a given manufactured rectangle could differ from the perfect rect-
angle. We will estimate the uncertainty in the area using (10.4.2), and find
that
∆A ≈ dA = Ax (20, 10) dx + Ay (20, 10) dy.
Since Ax = y and Ay = x, we have
That is, we estimate that the area in our rectangles could be off by as much
as 10 square centimeters.
Activity 10.4.4 The questions in this activity explore the differential in several
different contexts.
a. Suppose that the elevation of a landscape is given by the function h, where
we additionally know that h(3, 1) = 4.35, hx (3, 1) = 0.27, and hy (3, 1) =
−0.19. Assume that x and y are measured in miles in the easterly and
northerly directions, respectively, from some base point (0, 0). Your GPS
device says that you are currently at the point (3, 1). However, you
know that the coordinates are only accurate to within 0.2 units; that is,
dx = ∆x = 0.2 and dy = ∆y = 0.2. Estimate the uncertainty in your
elevation using differentials.
b. The pressure, volume, and temperature of an ideal gas are related by the
equation
P = P (T, V ) = 8.31T /V,
where P is measured in kilopascals, V in liters, and T in kelvin. Find
the pressure when the volume is 12 liters and the temperature is 310 K.
Use differentials to estimate the change in the pressure when the volume
increases to 12.3 liters and the temperature decreases to 305 K.
c. Refer to Table 10.4.8, the table of values of the wind chill w(v, T ), in
degrees Fahrenheit, as a function of temperature, also in degrees Fahren-
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 143
heit, and wind speed, in miles per hour. Suppose your anemometer says
the wind is blowing at 25 miles per hour and your thermometer shows a
reading of −15◦ degrees. However, you know your thermometer is only
accurate to within 2◦ degrees and your anemometer is only accurate to
within 3 miles per hour. What is the wind chill based on your mea-
surements? Estimate the uncertainty in your measurement of the wind
chill.
10.4.4 Summary
• A function f of two independent variables is locally linear at a point
(x0 , y0 ) if the graph of f looks like a plane as we zoom in on the graph
around the point (x0 , y0 ). In this case, the equation of the tangent plane
is given by
• The tangent plane L(x, y) = f (x0 , y0 )+fx (x0 , y0 )(x−x0 )+fy (x0 , y0 )(y −
y0 ), when considered as a function, is called the linearization of a differ-
entiable function f at (x0 , y0 ) and may be used to estimate values of
f (x, y); that is, f (x, y) ≈ L(x, y) for points (x, y) near (x0 , y0 ).
• A function f of two independent variables is differentiable at (x0 , y0 )
provided that both fx and fy exist and are continuous in an open disk
containing the point (x0 , y0 ).
• The differential df of a function f = f (x, y) is related to the differentials
dx and dy by
df = fx (x0 , y0 )dx + fy (x0 , y0 )dy.
We can use this relationship to approximate small changes in f that
result from small changes in x and y.
10.4.5 Exercises
p
1. Find the linearization L (x, y) of the function f (x, y) = 161 − 9x2 − 4y 2
at (−4, −2).
L (x, y) =
Note: Your answer should be an expression in x and y; e.g. “3x - 5y
+ 9”
Answer. 36(x − (−4)) + 8(y − (−2)) + 1
2. Find the equation of the tangent plane to the surface z = e2x/17 ln (1y) at
the point (3, 2, 0.9865).
z=
Note: Your answer should be an expression of x and y; e.g. “5x + 2y
- 3”
Answer. 0.116061x + 0.711624y + (−0.784911)
3. A student was asked to find the equation of the tangent plane to the
surface z = x4 − y 5 at the point (x, y) = (4, 5). The student’s answer was
z = −2869 + 4x3 (x − 4) − 5y 4 (y − 5).
(a) At a glance, how do you know this is wrong. What mistakes did
the student make? Select all that apply.
The answer is not a linear function. The (x - 4) and (y -
5) should be x and y. The partial derivatives were not evaluated
a the point. The -2869 should not be in the answer. All
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 144
of the above
(b) Find the correct equation for the tangent plane.
z=
Answer. AC
4. (a) Check the local linearity of f (x, y) = e−x cos(y) near x = 1, y = 1.5
by filling in the following table of values of f for x = 0.9, 1, 1.1 and
y = 1.4, 1.5, 1.6. Express values of f with 4 digits after the decimal
point.
x= 0.9 1 1.1
y = 1.4
y = 1.5
y = 1.6
(b) Next, fill in the table for the values x = 0.99, 1, 1.01 and y =
1.49, 1.5, 1.51, again showing 4 digits after the decimal point.
x= 0.99 1 1.01
y = 1.49
y = 1.5
y = 1.51
Notice if the two tables look nearly linear, and whether the second
looks more linear than the first (in particular, think about how you would
decide if they were linear, or if the one were more closely linear than the
other).
(c) Give the local linearization of f (x, y) = e−x cos(y) at (1, 1.5):
Using the second of your tables:
f (x, y) ≈
Using the fact that fx (x, y) = −e−x cos(y) and fy (x, y) = −e−x sin(y):
f (x, y) ≈
Answer 1. 0.0691
Answer 2. 0.0625
Answer 3. 0.0566
Answer 4. 0.0288
Answer 5. 0.026
Answer 6. 0.0235
Answer 7. −0.0119
Answer 8. −0.0107
Answer 9. −0.0097
Answer 10. 0.03
Answer 11. 0.0297
Answer 12. 0.0294
Answer 13. 0.0263
Answer 14. 0.026
Answer 15. 0.0258
Answer 16. 0.0226
Answer 17. 0.0224
Answer 18. 0.0221
Answer 19. 0.026 + (0.0258−0.0263)(x−1)
0.02 + (0.0224−0.0297)(y−1.5)
0.02
Answer 20. 0.0260228 + −0.0260228(x − 1) + −0.366958(y − 1.5)
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 145
it increases slightly
it does not change
it decreases slightly
(b) How does the energy change if the temperature is held constant
and the volume is increased slightly?
it decreases slightly
0
-1
x
0
-1
-1 0 y 1
a. Determine
f (0 + h, 0) − f (0, 0)
lim .
h→0 h
What does this limit tell us about fx (0, 0)?
b. Note that f (x, y) = f (y, x), and this symmetry implies that fx (0, 0) =
fy (0, 0). So both partial derivatives of f exist at (0, 0). A picture of
the surface defined by f near (0, 0) is shown in Figure 10.4.12. Based
on this picture, do you think f is locally linear at (0, 0)? Why?
Suppose you are given that IT (94, 75) = 3.75 and IH (94, 75) = 0.9.
Use this given information and one other value from the table to es-
timate the value of I(93.1, 77) using the linearization at (94, 75). Us-
ing proper terminology and notation, explain your work and think-
ing.
c. Just as we can find a local linearization for a differentiable function
of two variables, we can do so for functions of three or more variables.
By extending the concept of the local linearization from two to three
variables, find the linearization of the function h(x, y, z) = e2x (y +
z 2 ) at the point (x0 , y0 , z0 ) = (0, 1, −2). Then, use the linearization
to estimate the value of h(−0.1, 0.9, −1.8).
14. In the following questions, we investigate two different applied settings
using the differential.
a. Let f represent the vertical displacement in centimeters from the rest
position of a string (like a guitar string) as a function of the distance
x in centimeters from the fixed left end of the string and y the time
in seconds after the string has been plucked. (An interesting video of
this can be seen at https://www.youtube.com/watch?v=TKF6nFzpHBUA.)
A simple model for f could be
f (x0 + h) − f (x0 )
f 0 (x0 ) = lim
h→0 h
exists. We saw in single variable calculus that the existence of f 0 (x0 )
means that the graph of f is locally linear at x = x0 . In other
words, the graph of f looks like its linearization L(x) = f (x0 ) +
f 0 (x0 )(x − x0 ) for x close to x0 . That is, the values of f (x) can be
closely approximated by L(x) as long as x is close to x0 . We can
measure how good the approximation of L(x) is to f (x) with the
error function
at (0, 0) has a limit of 0 as (x, y) approaches (0, 0). (Use the lin-
earization you found in part (b).) This shows that just because an
error term goes to 0 as (x, y) approaches (x0 , y0 ), we cannot conclude
that a function is locally linear at (x0 , y0 ).
d. As the previous part illustrates, having the error term go to 0 does
not ensure that a function of two variables is locally linear. Instead,
we need a notation of a relative error. To see how this works, let
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 149
E(x0 + h)
.
h
Show that, for a function f = f (x) of a single variable, the limit of
the relative error is 0 as h approaches 0.
e. Even though the error term for a function of two variables might
have a limit of 0 at a point, our example shows that the function
may not be locally linear at that point. So we use the concept of
relative error to define differentiability of a function of two variables.
When we consider differentiability of a function f = f (x, y) at a
point (x0 , y0 ), then √if x = x0 + h and y = y0 + k, the distance from
(x, y) to (x0 , y0 ) is h2 + k 2 .
Definition 10.4.14 A function f = f (x, y) is differentiable at a
point (x0 , y0 ) if there is a linear function L = L(x, y) = f (x0 , y0 ) +
m(x − x0 ) + n(y − y0 ) such that the relative error
E(x0 + h, y0 + k)
√ ,
h2 + k 2
has at limit of 0 at (h, k) = (0, 0), where E(x, y) = f (x, y) − L(x, y),
h = x − x0 , and k = y − y0 . ♦
A function f is differentiable if it is differentiable at every point in
its domain. The function L in thepdefinition is the linearization of f
at (x0 , y0 ). Verify that g(x, y) = |xy| is not differentiable at (0, 0)
by showing that the relative error at (0, 0) does not have a limit at
(0, 0). Conclude that the existence of partial derivatives at a point is
not enough to ensure differentiability at that point. (Hint: Consider
the limit along different paths.)
16. Suppose that a function f = f (x, y) is differentiable at a point (x0 , y0 ).
Let L = L(x, y) = f (x0 , y0 ) + m(x − x0 ) + n(y − y0 ) as in the conditions of
Definition 10.4.14. Show that m = fx (x0 , y0 ) and n = fy (x0 , y0 ). (Hint:
Calculate the limits of the relative errors when h = 0 and k = 0.)
17. We know that if a function of a single variable is differentiable at a point,
then that function is also continuous at that point. In this exercise we
determine that the same property holds for functions of two variables. A
function f of the two variables x and y is continuous at a point (x0 , y0 )
in its domain if
lim f (x, y) = f (x0 , y0 )
(x,y)→(x0 ,y0 )
or (letting x = x0 + h and y = y0 + k,
1
y
x
x
1 2 3
Figure 10.5.1 Left: Your position in the plane. Right: The corresponding
temperature.
Suppose, furthermore, that the temperature at a point in the plane is given
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 151
by
1 1
T (x, y) = 10 − x2 − y 2 ,
2 5
and note that the surface generated by T is shown on the right of Figure 10.5.1.
Therefore, as time passes, your position (x(t), y(t)) changes, and, as your po-
sition changes, the temperature T (x, y) also changes.
a. The position function r provides a parameterization x = x(t) and y = y(t)
of the position at time t. By substituting x(t) for x and y(t) for y in the
formula for T , we can write T = T (x(t), y(t)) as a function of t. Make
these substitutions to write T as a function of t and then use the Chain
Rule from single variable calculus to find dT
dt . (Do not do any algebra to
simplify the derivative, either before taking the derivative, nor after.)
b. Now we want to understand how the result from part (a) can be obtained
from T as a multivariable function. Recall from the previous section that
small changes in x and y produce a change in T that is approximated by
∆T ≈ Tx ∆x + Ty ∆y.
The Chain Rule tells us about the instantaneous rate of change of T , and
this can be found as
∆T Tx ∆x + Ty ∆y
lim = lim . (10.5.1)
∆t→0 ∆t ∆t→0 ∆t
Use equation (10.5.1) to explain why the instantaneous rate of change of
T that results from a change in t is
dT ∂T dx ∂T dy
= + . (10.5.2)
dt ∂x dt ∂y dt
Since z is a function of the two variables x and y, the derivatives in the Chain
Rule for z with respect to x and y are partial derivatives. However, since
x = x(t) and y = y(t) are functions of the single variable t, their derivatives
are the standard derivatives of functions of one variable. When we compose z
with x(t) and y(t), we then have z as a function of the single variable t, making
the derivative of z with respect to t a standard derivative from single variable
calculus as well.
To understand why this Chain Rule works in general, suppose that some
quantity z depends on x and y so that
∂z ∂z
dz = dx + dy. (10.5.4)
∂x ∂y
Next, suppose that x and y each depend on another quantity t, so that
dx dy
dx = dt and dy = dt. (10.5.5)
dt dt
Combining Equations (10.5.4) and (10.5.5), we find that
∂z dx ∂z dy dz
dz = dt + dt = dt,
∂x dt ∂y dt dt
which is the Chain Rule in this particular context, as expressed in Equa-
tion (10.5.3).
Activity 10.5.2 In the following questions, we apply the Chain Rule in several
different contexts.
a. Suppose that we have a function z defined by z(x, y) = x2 + xy 3 . In
addition, suppose that x and y are restricted to points that move around
the plane by following a circle of radius 2 centered at the origin that is
parameterized by
i. Use the Chain Rule to find the resulting instantaneous rate of change
dz
dt .
ii. Substitute x(t) for x and y(t) for y in the rule for z to write z in
terms of t and calculate dz
dt directly. Compare to the result of part
(i.).
c. Suppose that you are walking along a surface whose elevation is given
by a function f . Furthermore, suppose that if you consider how your
location corresponds to points in the xy-plane, you know that when you
pass the point (2, 1), your velocity vector is v = h−1, 2i. If some contours
of f are as shown in Figure 10.5.2, estimate the rate of change df /dt when
you pass through (2, 1).
y
5
4
4
3
2
2
1
x
2 4
z z
∂z ∂z
∂x ∂y
x y x y
dx dy
dt dt
t t t t
Figure 10.5.3 A tree diagram illustrating dependencies.
To represent the Chain Rule, we label every edge of the diagram with the
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 154
w x y
u v u v u v
t t t t t t
x = r cos(θ)
y = r sin(θ).
10.5.3 Summary
• The Chain Rule is a tool for differentiating a composite for functions. In
its simplest form, it says that if f (x, y) is a function of two variables and
x(t) and y(t) depend on t, then
df ∂f dx ∂f dy
= + .
dt ∂x dt ∂y dt
• A tree diagram can be used to represent the dependence of variables on
other variables. By following the links in the tree diagram, we can form
chains of partial derivatives or derivatives that can be combined to give
a desired partial derivative.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 155
10.5.4 Exercises
dz
1. Use the chain rule to find dt , where
z = x2 y + xy 2 , x = 1 − t5 , y = 3 − t2
3s
z = exy tan y, x = 4s + 5t, y =
3t
First the pieces:
∂z ∂z
∂x = ∂y =
∂x ∂x
∂s = ∂t =
∂y ∂y
∂s = ∂t =
Answer 2. 2.75
4. If z = (x + y) ey and x = u2 + v 2 and y = u2 − v 2 , find the following
partial derivatives using the chain rule. Enter your answers as functions
of u and v.
∂z
=
∂u
∂z
=
∂v
2 2 2 2
Answer 1. 2 · 2ueu −v + 2u2 eu −v · 2u
2 2
Answer 2. −2u2 eu −v · 2v
5. If
z = sin x2 + y 2 ,
x = v cos(u) , y = v sin(u) ,
find ∂z/∂u and ∂z/∂v. The variables are restricted to domains on which
the functions are defined.
∂z/∂u =
∂z/∂v =
Answer 1. 0
Answer 2. cos v 2 · 2v
6. Let z = g(u, v) and u(r, s), v(r, s). How many terms are there in the
expression for ∂z/∂r?
terms
Answer. 2
7. Let W (s, t) = F (u(s, t), v(s, t)) where
∂z ∂z
B. Find the numerical values of ∂s and ∂t when (s, t) = (−3, 3).
∂z
∂s (−3, 3) =
∂z
∂t (−3, 3) =
Answer 1. 8sx sin(y) + (−6)tx2 cos(y)
Answer 2. 4tx sin(y) + (−6)sx2 cos(y)
Answer 3. 11244.2989304861
Answer 4. −11063.2512786582
11. Find the indicated derivative. In each case, state the version of the Chain
Rule that you are using.
df
a. dt , if f (x, y) = 2x2 y, x = cos(t), and y = ln(t).
∂f 2z+1
b. ∂w , if f (x, y) = 2x2 y, x = w + z 2 , and y = w
∂f
c. if f (x, y, z) = 2x2 y + z 3 , x = u − v + 2w, y = w2v − u3 , and
∂v ,
z = u2 − v
12. Let z = u2 − v 2 and suppose that
u = ex cos(y)
v = ex sin(y)
x = ρ sin(φ) cos(θ)
y = ρ sin(φ) sin(θ)
z = ρ cos(φ)
6 y
1
x
1 2 3 4 5 6
a. Suppose that a person is walking due east, and thus parallel to the x-axis.
At what instantaneous rate is the temperature changing with respect to
x at the moment the walker passes the point (2, 1)? What are the units
on this rate of change?
b. Next, determine the instantaneous rate of change of temperature with
respect to distance at the point (2, 1) if the person is instead walking due
north. Again, include units on your result.
c. Now, rather than walking due east or due north, let’s suppose that the
person is walking with velocity given by the vector v = h3, 4i, where time
is measured in seconds. Note that the person’s speed is thus |v| = 5 feet
per second. Find parametric equations for the person’s path; that is,
parameterize the line through (2, 1) using the direction vector v = h3, 4i.
Let x(t) denote the x-coordinate of the line, and y(t) its y-coordinate.
Make sure your parameterization places the walker at the point (2, 1)
when t = 0.
d. With the parameterization in (c), we can now view the temperature f
as not only a function of x and y, but also of time, t. Hence, use the
chain rule to determine the value of df
dt t=0 . What are the units on your
answer? What is the practical meaning of this result?
f (x + u1 h, y + u2 h) − f (x, y)
Du f (x, y) = lim (10.6.2)
h→0 h
for those values of x and y for which the limit exists. ♦
The quantity Du f (x, y) is called a directional derivative. When we evaluate
the directional derivative Du f (x, y) at a point (x0 , y0 ), the result Du f (x0 , y0 )
tells us the instantaneous rate at which f changes at (x0 , y0 ) per unit increase
in the direction of the vector u. In addition, the quantity Du f (x0 , y0 ) tells us
the slope of the line tangent to the surface in the direction of u at the point
(x0 , y0 , f (x0 , y0 )).
x = x0 + u1 t and y = y0 + u2 t.
dy
Observe that dx dt = u1 and dt = u2 for all values of t. Since u is a unit
vector, it follows that a point moving along this line moves one unit of distance
per one unit of time; that is, each single unit of time corresponds to movement
of a single unit of distance in that direction. This observation allows us to
use the Chain Rule to calculate the directional derivative, which measures the
instantaneous rate of change of f with respect to change in the direction u.
In particular, by the Chain Rule, it follows that
dx dy
Du f (x0 , y0 ) = fx (x0 , y0 ) +fy (x0 , y0 )
dt (x0 ,y0 ) dt (x0 ,y0 )
1 The symbol ∇ is called nabla, which comes from a Greek word for a certain type of harp
• (x0 , y0 ) = (2, 0)
• (x0 , y0 ) = (0, 2)
• (x0 , y0 ) = (2, 2)
• (x0 , y0 ) = (2, 1)
• (x0 , y0 ) = (−3, 2)
• (x0 , y0 ) = (−2, −4)
• (x0 , y0 ) = (0, 0)
c. What do you notice about the relationship between the gradient at
(x0 , y0 ) and the contour passing through that point?
y y y
∇f (x0 , y0 ) ∇f (x0 , y0 ) ∇f (x0 , y0 )
u
θ θ θ
u
(x0 , y0 ) (x0 , y0 ) u (x0 , y0 )
x x x
b. At the point (x0 , y0 ), the only quantity in Equation (10.6.6) that can
change is θ (which determines the direction u of travel). Explain why
θ = 0 makes the quantity
|hfx (x0 , y0 ), fy (x0 , y0 )i| cos(θ)
as large as possible.
c. When θ = 0, in what direction does the unit vector u point relative to
∇f (x0 , y0 )? Why? What does this tell us about the direction of greatest
increase of f at the point (x0 , y0 )?
d. In what direction, relative to ∇f (x0 , y0 ), does f decrease most rapidly
at the point (x0 , y0 )?
e. State the unit vectors u and v (in terms of ∇f (x0 , y0 )) that provide the
directions of greatest increase and decrease for the function f at the point
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 165
4 y
x
1 2 3 4
h. Suppose you are standing at the point (3, 3). In which direction should
you move to cause f to increase as rapidly as possible? At what rate
does f increase in this direction?
10.6.6 Applications
The gradient finds many natural applications. For example, situations often
arise — for instance, constructing a road through the mountains or planning
the flow of water across a landscape — where we are interested in knowing the
direction in which a function is increasing or decreasing most rapidly.
For example, consider a two-dimensional version of how a heat-seeking mis-
sile might work.(This application is borrowed from United States Air Force
Academy Department of Mathematical Sciences.) Suppose that the tempera-
ture surrounding a fighter jet can be modeled by the function T defined by
100
T (x, y) = ,
1 + (x − 5)2 + 4(y − 2.5)2
where (x, y) is a point in the plane of the fighter jet and T (x, y) is measured
in degrees Celsius. Some contours and gradients ∇T are shown on the left in
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 167
Figure 10.6.6.
5 y 5 y
4 4
3 3
2 2
1 1
x x
1 2 3 4 5 6 7 8 9 10 1 2 3 4 5 6 7 8 9 10
Figure 10.6.6 Contours and gradient for T (x, y) and the missile’s path.
A heat-seeking missile will always travel in the direction in which the tem-
perature increases most rapidly; that is, it will always travel in the direction
of the gradient ∇T . If a missile is fired from the point (2, 4), then its path will
be that shown on the right in Figure 10.6.6.
In the final activity of this section, we consider several questions related to
this context of a heat-seeking missile, and foreshadow some upcoming work in
Section 10.7.
Activity 10.6.6
a. The temperature T (x, y) has its maximum value at the fighter jet’s loca-
tion. State the fighter jet’s location and explain how Figure 10.6.6 tells
you this.
b. Determine ∇T at the fighter jet’s location and give a justification for
your response.
c. Suppose that a different function f has a local maximum value at (x0 , y0 ).
Sketch the behavior of some possible contours near this point. What is
∇f (x0 , y0 )? (Hint: What is the direction of greatest increase in f at the
local maximum?)
d. Suppose that a function g has a local minimum value at (x0 , y0 ). Sketch
the behavior of some possible contours near this point. What is ∇g(x0 , y0 )?
10.6.7 Summary
• The directional derivative of f at the point (x, y) in the direction of the
unit vector u = hu1 , u2 i is
f (x + u1 h, y + u2 h) − f (x, y)
Du f (x, y) = lim
h→0 h
for those values of x and y for which the limit exists. In addition,
Du f (x, y) measures the slope of the graph of f when we move in the
direction u. Alternatively, Du f (x0 , y0 ) measures the instantaneous rate
of change of f in the direction u at (x0 , y0 ).
• The gradient of a function f = f (x, y) at a point (x0 , y0 ) is the vector
∇f (x0 , y0 ) = hfx (x0 , y0 ), fy (x0 , y0 )i .
10.6.8 Exercises
1. Consider the function f (x, y, z) = xy + yz 2 + xz 3 .
Find the gradient of f :
h , ,
i
Find the gradient of f at the point (2, 3, 5).
h , ,
i
Find the rate √
of change√of the function
√ f at the point (2, 3,5) in the
direction u = h4/ 66, −5/ 66, −5/ 66i.
Answer 1. y + z 3
Answer 2. x + z 2
Answer 3. 2yz + 3xz 2
Answer 4. 128
Answer 5. 27
Answer 6. 180
Answer 7. −64.376849782191
2. If f (x, y) = 3x2 − 1y 2 , find the value of the directional derivative at the
point (−1, −4) in the direction given by the angle θ = 2π 6 .
Answer. 3.92820323027551
3. Find the directional derivative of f (x, y, z) = 4xy + z 2 at the point
(−5, 3, −2) in the direction of the maximum rate of change of f .
f<B><I>u</B></I> (−5, 3, −2) = D<B><I>u</B></I> f (−5, 3, −2) =
Answer. 23.6643
100
4. The temperature at any point in the plane is given by T (x, y) = .
x2 + y2 + 3
(a) What shape are the level curves of T ?
hyperbolas
circles
lines
ellipses
parabolas
What is the value of this maximum rate of change, that is, the maxi-
mum value of the directional derivative at (−3, −3)?
(d) Find the direction of the greatest decrease in temperature at the
point (−3, −3).
What is the value of this most negative rate of change, that is, the
minimum value of the directional derivative at (−3, −3)?
Answer 1. B
Answer 2. (0, 0)
Answer 3. 33.3333
Answer 4. h1.36054, 1.36054i
Answer 5. 1.9241
Answer 6. h−1.36054, −1.36054i
2 2 2
5. The temperature at a point (x,y,z) is given by T (x, y, z) = 200e−x −y /4−z /9 ,
where T is measured in degrees Celsius and x,y, and z in meters. There
are lots of places to make silly errors in this problem; just try to keep
track of what needs to be a unit vector.
Find the rate of change of the temperature at the point (-1, 1, 2) in
the direction toward the point (-1, 3, -3).
In which direction (unit vector) does the temperature increase the
fastest at (-1, 1, 2)?
h , ,
i
What is the maximum rate of increase of T at (-1, 1, 2)?
Answer 1. 8.33860305193074
Answer 2. 0.948354065592898
Answer 3. −0.237088516398225
Answer 4. −0.210745347909533
Answer 5. 77.4821403113886
6. If f (x, y, z) = 2zy 2 , then the gradient at the point (2, 2, 4) is
∇f (2, 2, 4) =
Answer. h0, 32, 8i
7. The concentration of salt in a fluid at (x, y, z) is given by F (x, y, z) =
2x2 + 3y 4 + 2x2 z 2 mg/cm3 . You are at the point (−1, 1, −1).
(a) In which direction should you move if you want the concentration
to increase the fastest?
direction:
(Give your answer as a vector.)
(b) You start to move in the direction you found in part (a) at a speed
of 5 cm/sec. How fast is the concentration changing?
rate of change =
Answer 1. −8 ĩ + 12 j̃ − 4 k̃
p
Answer 2. 5 −8 · (−8) + 12 · 12 + −4 · (−4)
8. At a certain point on a heated metal plate, the greatest rate of tempera-
ture increase, 4 degrees Celsius per meter, is toward the northeast. If an
object at this point moves directly north, at what rate is the temperature
increasing?
9. Suppose that you are climbing a hill whose shape is given by z = 902 −
0.07x2 − 0.1y 2 , and that you are at the point (40, 70, 300).
In which direction (unit vector) should you proceed initially in order
to reach the top of the hill fastest?
h ,
i
If you climb in that direction, at what angle above the horizontal will
you be climbing initially (radian measure)?
Answer 1. −0.371390676354104
Answer 2. −0.928476690885259
Answer 3. 1.5045735390203
10. Are the following statements true or false?
(a) The gradient vector ∇f (a, b) is tangent to the contour of f at (a, b).
(b) f~u (a, b) = ||∇f (a, b)||.
(c) f~u (a, b) is parallel to ~u.
(g) If f (x, y) has fx (a, b) = 0 and fy (a, b) = 0 at the point (a, b), then
f is constant everywhere.
(h) ∇f (a, b) is a vector in 3-dimensional space.
11. Let E(x, y) = 1+(x−5)2100
+4(y−2.5)2 represent the elevation on a land mass
at location (x, y). Suppose that E, x, and y are all measured in meters.
a. Find Ex (x, y) and Ey (x, y).
b. Let u be a unit vector in the direction of h−4, 3i. Determine Du E(3, 4).
What is the practical meaning of Du E(3, 4) and what are its units?
c. Find the direction of greatest increase in E at the point (3, 4).
d. Find the instantaneous rate of change of E in the direction of great-
est decrease at the point (3, 4). Include units on your answer.
b. Find the direction of greatest increase in f at the point (1, 2). Ex-
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 171
-2
-4 -2
0 x
-4 -2 0 2 2
4
y -3
-4
d. Find a direction in which the temperature does not change at (1, 2, −1).
16. Figure 10.6.8 shows a plot of the gradient ∇f at several points for some
function f = f (x, y).
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 172
10.7 Optimization
Motivating Questions
• How can we find the points at which f (x, y) has a local maximum or
minimum?
• How can we determine whether critical points of f (x, y) are local maxima
or minima?
• How can we find the absolute maximum and minimum of f (x, y) on a
closed and bounded domain?
We learn in single-variable calculus that the derivative is a useful tool for
finding the local maxima and minima of functions, and that these ideas may
often be employed in applied settings. In particular, if a function f , such
as the one shown in Figure 10.7.1 is everywhere differentiable, we know that
the tangent line is horizontal at any point where f has a local maximum or
minimum. This, of course, means that the derivative f 0 is zero at any such
point. Hence, one way that we seek extreme values of a given function is to
first find where the derivative of the function is zero.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 173
y = f (x)
b. In the same way, the trace given by holding x = x0 constant has a local
maximum at y = y0 . What does this say about the value of the partial
derivative fy (x0 , y0 )?
c. What may we now conclude about the gradient ∇f (x0 , y0 ) at the local
maximum? How is this consistent with the statement “f increases most
rapidly in the direction ∇f (x0 , y0 )?”
d. How will the tangent plane to the surface z = f (x, y) appear at the point
(x0 , y0 , f (x0 , y0 ))?
e. By first computing the partial derivatives, find any points at which
f (x, y) = 2x − x2 − (y + 2)2 may have a local maximum.
• An absolute minimum point is a point such that f (x, y) ≥ f (x0 , y0 ) for all
points (x, y) in the domain of f . The value of f at an absolute minimum
point is the maximum value of f .
♦
We use the term extremum point to refer to any point (x0 , y0 ) at which f
has a local maximum or minimum. In addition, the function value f (x0 , y0 ) at
an extremum is called an extremal value. Figure 10.7.3 illustrates the graphs
of two functions that have an absolute maximum and minimum, respectively,
at the origin (x0 , y0 ) = (0, 0).
doesn’t exist or for which both partial derivatives are simultaneously zero. For
the latter, note that we have to solve the system of equations
fx (x, y) = 0
fy (x, y) = 0.
Activity 10.7.2 Find the critical points of each of the following functions.
Then, using appropriate technology, plot the graphs of the surfaces near each
critical value and compare the graph to your work.
a. f (x, y) = 2 + x2 + y 2
b. f (x, y) = 2 + x2 − y 2
c. f (x, y) = 2x − x2 − 14 y 2
x y
Activity 10.7.3 Recall that the Second Derivative Test for single-variable
functions states that if x0 is a critical point of a function f so that f 0 (x0 ) = 0
and if f 00 (x0 ) exists, then
• if f 00 (x0 ) < 0, x0 is a local maximum,
• if f 00 (x0 ) > 0, x0 is a local minimum, and
You can verify that each function has a critical point at the origin (0, 0). You
should check this.
x y
a. The graphs of these three functions are shown in Figure 10.7.6, with
z = 4 − x2 − y 2 at left, z = x2 + y 2 in the middle, and z = x2 − y 2 at
right. Use the graphs to decide if a function has a local maximum, local
minimum, saddle point, or none of the above at the origin.
b. There is no single second derivative of a function of two variables, so we
consider a quantity that combines the second order partial derivatives.
2
Let D = fxx fyy − fxy . Calculate D at the origin for each of the functions
f1 , f2 , and f3 . What difference do you notice between the values of D
when a function has a maximum or minimum value at the origin versus
when a function has a saddle point at the origin?
c. Now consider the cases where D > 0. It is in these cases that a function
has a local maximum or minimum at a point. What is necessary in these
cases is to find a condition that will distinguish between a maximum and
a minimum. In the cases where D > 0 at the origin, evaluate fxx (0, 0).
What value does fxx (0, 0) have when f has a local maximum value at
the origin? When f has a local minimum value at the origin? Explain
why. (Hint: This should look very similar to the Second Derivative Test
for functions of a single variable.) What would happen if we considered
the values of fyy (0, 0) instead?
Activity 10.7.3 provides the basic ideas for the Second Derivative Test for
functions of two variables.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 177
c. f (x, y) = x3 + y 3 − 3xy.
As we learned in single-variable calculus, finding extremal values of func-
tions can be particularly useful in applied settings. For instance, we can often
use calculus to determine the least expensive way to construct something or to
find the most efficient route between two locations. The same possibility holds
in settings with two or more variables.
Activity 10.7.5 While the quantity of a product demanded by consumers is
often a function of the price of the product, the demand for a product may also
depend on the price of other products. For instance, the demand for blue jeans
at Old Navy may be affected not only by the price of the jeans themselves, but
also by the price of khakis.
Suppose we have two goods whose respective prices are p1 and p2 . The
demand for these goods, q1 and q2 , depend on the prices as
The seller would like to set the prices p1 and p2 in order to maximize
revenue. We will assume that the seller meets the full demand for each product.
Thus, if we let R be the revenue obtained by selling q1 items of the first good
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 178
at price p1 per item and q2 items of the second good at price p2 per item, we
have
R = p1 q1 + p2 q2 .
We can then write the revenue as a function of just the two variables p1
and p2 by using Equations (10.7.1) and (10.7.2), giving us
4000
2000
-2000
-4000
-6000
-8000
-10000
10 20 30 40 20 10 0
p2 50 60 60 50 40 30 p1
a. Find all critical points of the revenue function, R. (Hint: You should
obtain a system of two equations in two unknowns which can be solved
by elimination or substitution.)
b. Apply the Second Derivative Test to determine the type of any critical
point(s).
c. Where should the seller set the prices p1 and p2 to maximize the revenue?
The absolute extremes must occur at either a critical point in the interior
of R or at a boundary point of R. We therefore must test both possibilities,
as we demonstrate in the following example.
Example 10.7.8 Suppose the temperature T at each point on the circular
plate x2 + y 2 ≤ 1 is given by
T (x, y) = 2x2 + y 2 − y.
The domain R = {(x, y) : x2 + y 2 ≤ 1} is a closed and bounded region, as
shown on the left of Figure 10.7.9, so the Extreme Value Theorem assures us
that T has an absolute maximum and minimum on the plate. The graph of
T over its domain R is shown in Figure 10.7.9. We will find the hottest and
coldest points on the plate.
1.5 y
z
2
1.0
0.5
1
x
-1.5 -1.0 -0.5 0.5 1.0 1.5
-0.5 -1
0 0 x
y 1
-1.0
-1
-1.5
Figure 10.7.9 Domain of the temperature T (x, y) = 2x2 + y 2 − y and its
graph.
If the absolute maximum or minimum occurs inside the disk, it will be at
a critical point so we begin by looking for critical points inside the disk. To do
this, notice that critical points are given by the conditions Tx = 4x = 0 and
Ty = 2y − 1 = 0. This means that there is one critical point of the function at
the point (x0 , y0 ) = (0, 1/2), which lies inside the disk.
We now find the hottest and coldest points on the boundary of the disk,
which is the circle of radius 1. As we have seen, the points on the unit circle
can be parametrized as
x(t) = cos(t), y(t) = sin(t),
where 0 ≤ t ≤ 2π. The temperature at a point on the circle is then described
by
T (x(t), y(t)) = 2 cos2 (t) + sin2 (t) − sin(t).
To find the hottest and coldest points on the boundary, we look for the
critical points of this single-variable function on the interval 0 ≤ t ≤ 2π. We
have
dT
= −4 cos(t) sin(t) + 2 cos(t) sin(t) − cos(t)
dt
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 180
This shows that we have critical points when cos(t) = 0 or sin(t) = −1/2.
This occurs when t = π/2, 3π/2, 7π/6, and 11π/6. Since we have x(t) = cos(t)
and y(t) = sin(t), the corresponding points are
π 3π
• (x, y) = t = 2,
(0, 1) when • (x, y) =
(0, −1) when
t= 2 ,
√ √
• (x, y) = 23 , − 12 when t = 11π
6 . • (x, y) = − 23 , − 21 when t = 7π
6 .
These are the critical points of T on the boundary and so this collection of
points includes the hottest and coldest points on the boundary.
We now have a list of candidates for the hottest and coldest points: the
critical point in the interior of the disk and the critical points on the boundary.
We find the hottest and coldest points by evaluating the temperature at each
of these points, and find that
• T 0, 12 = −14 ,
• T (0,
1) = 0, • T (0, −1) = 2,
√ √
3 1 9 3 1 9
• T − 2 , −2 = 4, • T − 2 , −2 = 4.
2 2 9
So the maximum √ value of T on the disk x + y ≤ 1 is 4 , which occurs at
the two points ± 23 , − 12 on the boundary, and the minimum value of T on
the disk is − 14 which occurs at the critical point 0, 21 in the interior of R.
From this example, we see that we use the following procedure for deter-
mining the absolute maximum and absolute minimum of a function on a closed
and bounded domain.
• Find all critical points of the function in the interior of the domain.
• Find all the critical points of the function on the boundary of the domain.
Working on the boundary of the domain reduces this part of the problem
to one or more single variable optimization problems. Note that there
may be endpoints on portions of the boundary that need to be considered.
z
y 5
4
1 2 2 1 x 0
y 3 3
3 -5 4
-10
2
-15
1
-20
x
1 2 3 4
d. Parameterize the hypotenuse of the triangular domain, and find the crit-
ical points of f on the hypotenuse. (Hint: You may need to consider
endpoints.)
e. Find the absolute maximum and absolute minimum values of f on R.
10.7.4 Summary
• To find the extrema of a function f = f (x, y), we first find the critical
points, which are points where one of the partials of f fails to exist, or
where fx = 0 and fy = 0.
10.7.5 Exercises
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 182
1. The function 2
k(x, y) = e−y cos(5x)
has a critical point at (0, 0).
What is the value of D at this critical point? D =
What type of critical point is it? ( maximum minimum sad-
dle point point with unknown behavior)
Answer 1. 2 · 25
Answer 2. maximum
2. Consider the function
Find and classify all critical points of the function. If there are more
blanks than critical points, leave the remaining entries blank.
fx =
fy =
fxx =
fxy =
fyy =
There are several critical points to be listed. List them lexicograhically,
that is in ascending order by x-coordinates, and for equal x-coordinates
in ascending order by y-coordinates (e.g., (1,1), (2, -1), (2, 3) is a correct
order)
The critical point with the smallest x-coordinate is
( , ) Classification:
(local minimum, local maximum, saddle point, cannot be determined)
The critical point with the next smallest x-coordinate is
( , ) Classification:
(local minimum, local maximum, saddle point, cannot be determined)
The critical point with the next smallest x-coordinate is
( , ) Classification:
(local minimum, local maximum, saddle point, cannot be determined)
The critical point with the next smallest x-coordinate is
( , ) Classification:
(local minimum, local maximum, saddle point, cannot be determined)
The critical point with the next smallest x-coordinate is
( , ) Classification:
(local minimum, local maximum, saddle point, cannot be determined)
Answer 1. (18 − 2x) 18y − y 2
Answer 2. 18x − x2 (18 − 2y)
Answer 3. −2 18y − y 2
Answer 4. (18 − 2x)(18 − 2y)
Answer 5. −2 18x − x2
Answer 6. 0
Answer 7. 0
Answer 8. saddle point
Answer 9. 0
Answer 10. 18
Answer 11. saddle point
Answer 12. 9
Answer 13. 9
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 183
f (x, y) = 9 − 6x + 8y
on the closed triangular region with vertices (0, 0), (8, 0) and (8, 12).
List the minimum/maximum values as well as the point(s) at which
they occur. If a min or max occurs at multiple points separate the points
with commas.
Minimum value:
Occurs at
Maximum value:
Occurs at
Answer 1. −39
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 184
Answer 2. (8, 0)
Answer 3. 57
Answer 4. (8, 12)
7. Find the maximum and minimum values of f (x, y) = xy on the ellipse
5x2 + y 2 = 3.
maximum value =
minimum value =
Answer 1. 0.670820393249937
Answer 2. −0.670820393249937
8. Find A and B so that f (x, y) = x2 + Ax + y 2 + B has a local minimum
at the point (7, 0), with z-coordinate 5.
A=
B=
Answer 1. −14
Answer 2. 5 + 7 · 7
9. The contours of a function f are shown in the figure below.
For each of the points shown, indicate whether you think it is a local
maximum, local minimum, saddle point, or none of these.
(a) Point P is ( a local maximum a local minimum a
saddle point none of these)
(b) Point Q is ( a local maximum a local minimum a
saddle point none of these)
(c) Point R is ( a local maximum a local minimum a
saddle point none of these)
(d) Point S is ( a local maximum a local minimum a saddle
point none of these)
Answer 1. a local maximum
Answer 2. a local maximum
Answer 3. none of these
Answer 4. a saddle point
10. Consider the three points (5, 1), (7, 3), and (9, 4).
(a) Supposed that at (5, 1), we know that fx = fy = 0 and fxx = 0,
fyy = 0, and fxy > 0. What can we conclude about the behavior of this
function near the point (5, 1)? ( (5,1) is a local maximum (5,1) is
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 185
x2 y2 z2
+ + =1
16 81 4
Hint: By symmetry, you can restrict your attention to the first octant
(where x, y, z ≥ 0), and assume your volume has the form V = 8xyz.
Then arguing by symmetry, you need only look for points which achieve
the maximum which lie in the first octant.
Maximum volume:
Answer. 110.851251684408
16. Design a rectangular milk carton box of width w, length l, and height h
which holds 500 cm3 of milk. The sides of the box cost 1 cent/cm2 and
the top and bottom cost 2 cent/cm2 . Find the dimensions of the box that
minimize the total cost of materials used.
dimensions =
(Enter your answer as a comma separated list of lengths.)
Answer. 12.5992, 6.2996, 6.29962
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 186
17. Respond to each of the following prompts to solve the given optimization
problem.
a. Let f (x, y) = sin(x) + cos(y). Determine the absolute maximum
and minimum values of f . At what points do these extreme values
occur?
b. For a certain differentiable function F of two variables x and y, its
partial derivatives are
22. The airlines place restrictions on luggage that can be carried onto planes.
• A carry-on bag can weigh no more than 40 lbs.
• The length plus width plus height of a bag cannot exceed 45 inches.
• The bag must fit in an overhead bin.
Let x, y, and z be the length, width, and height (in inches) of a carry on
bag. In this problem we find the dimensions of the bag of largest volume,
V = xyz, that satisfies the second restriction. Assume that we use all
45 inches to get a maximum volume. (Note that this bag of maximum
volume might not satisfy the third restriction.)
a. Write the volume V = V (x, y) as a function of just the two variables
x and y.
a. Calculate Sm and Sb .
b. Solve the system Sm (m, b) = 0 and Sb (m, b) = 0 to show that the
critical point satisfies
Pn Pn Pn
n ( i=1 xi yi ) − ( i=1 xi ) ( i=1 yi )
m= Pn Pn 2
n ( i=1 x2i ) − ( i=1 xi )
Pn Pn 2
P n Pn
( i=1 yi ) i=1 xi − ( i=1 xi ) ( i=1 xi yi )
b= Pn Pn 2 .
2
n ( i=1 xi ) − ( i=1 xi )
• How can we exploit this geometric condition to find the extreme values
of a function subject to a constraint?
We previously considered how to find the extreme values of functions on
both unrestricted domains and on closed, bounded domains. Other types of
optimization problems involve maximizing or minimizing a quantity subject
to an external constraint. In these cases the extreme values frequently won’t
occur at the points where the gradient is zero, but rather at other points that
satisfy an important geometric condition. These problems are often called con-
strained optimization problems and can be solved with the method of Lagrange
Multipliers, which we study in this section.
Preview Activity 10.8.1 According to U.S. postal regulations, the girth
plus the length of a parcel sent by mail may not exceed 108 inches, where
by “girth” we mean the perimeter of the smallest end. Our goal is to find
the largest possible volume of a rectangular parcel with a square end that
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 190
can be sent by mail. (We solved this applied optimization problem in single
variable Active Calculus, so it may look familiar. We take a different approach
in this section, and this approach allows us to view most applied optimization
problems from single variable calculus as constrained optimization problems,
as well as provide us tools to solve a greater variety of optimization problems.)
If we let x be the length of the side of one square end of the package and y
the length of the package, then we want to maximize the volume f (x, y) = x2 y
of the box subject to the constraint that the girth (4x) plus the length (y) is
as large as possible, or 4x + y = 108. The equation 4x + y = 108 is thus an
external constraint on the variables.
a. The constraint equation involves the function g that is given by
g(x, y) = 4x + y.
100
80
A
60
C
40
5000
20000
20 D
1000 10000
B x
10 20 30 40
b. Figure 10.8.1 shows the graph of the constraint equation g(x, y) = 108
along with a few contours of the volume function f . Since our goal is to
find the maximum value of f subject to the constraint g(x, y) = 108, we
want to find the point on our constraint curve that intersects the contours
of f at which f has its largest value.
i. Points A and B in Figure 10.8.1 lie on a contour of f and on the
constraint equation g(x, y) = 108. Explain why neither A nor B
provides a maximum value of f that satisfies the constraint.
ii. Points C and D in Figure 10.8.1 lie on a contour of f and on the
constraint equation g(x, y) = 108. Explain why neither C nor D
provides a maximum value of f that satisfies the constraint.
iii. Based on your responses to parts i. and ii., draw the contour of f
on which you believe f will achieve a maximum value subject to the
constraint g(x, y) = 108. Explain why you drew the contour you
did.
c. Recall that g(x, y) = 108 is a contour of the function g, and that the
gradient of a function is always orthogonal to its contours. With this in
mind, how should ∇f and ∇g be related at the optimal point? Explain.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 191
100
80
60
40
20
x
10 20 30 40
Figure 10.8.2 Contours of f and the constraint contour.
To find this point where the graph of the constraint is tangent to a contour
of f , recall that ∇f is perpendicular to the contours of f and ∇g is perpen-
dicular to the contour of g. At such a point, the vectors ∇g and ∇f are
parallel, and thus we need to determine the points where this occurs. Recall
that two vectors are parallel if one is a nonzero scalar multiple of the other, so
we therefore look for values of a parameter λ that make
∇f = λ∇g. (10.8.1)
∇f = 2xyi + x2 j and ∇g = 4i + j,
So the point (0, 108) is a point we need to consider. Next, provided that λ 6= 0
(from which it follows that x 6= 0 by Equation (10.8.3)), we may divide both
sides of Equation (10.8.2) by the corresponding sides of (10.8.3) to eliminate
λ, and thus find that
2y
= 4, so
x
y = 2x.
4x + 2x = 108
or
x = 18.
Thus we have y = 2x = 36 and λ = x2 = 324 as another point to consider.
So the points at which the gradients of f and g are parallel, and thus at which
f may have a maximum or minimum subject to the constraint, are (0, 108) and
(18, 36). By evaluating the function f at these points, we see that we maximize
the volume when the length of the square end of the box is 18 inches and the
length is 36 inches, for a maximum volume of f (18, 36) = 11664 cubic inches.
Since f (0, 108) = 0, we obtain a minimum value at this point.
We summarize the process of Lagrange multipliers as follows.
Activity 10.8.2 A cylindrical soda can holds about 355 cc of liquid. In this
activity, we want to find the dimensions of such a can that will minimize the
surface area. For the sake of simplicity, assume the can is a perfect cylinder.
a. What are the variables in this problem? Based on the context, what
restriction(s), if any, are there on these variables?
b. What quantity do we want to optimize in this problem? What equation
describes the constraint? (You need to decide which of these functions
plays the role of f and which plays the role of g in our discussion of
Lagrange multipliers.)
c. Find λ and the values of your variables that satisfy Equation (10.8.1) in
the context of this problem.
d. Determine the dimensions of the pop can that give the desired solution
to this constrained optimization problem.
The method of Lagrange multipliers also works for functions of more than
two variables.
Activity 10.8.3 Use the method of Lagrange multipliers to find the dimensions
of the least expensive packing crate with a volume of 240 cubic feet when the
material for the top costs $2 per square foot, the bottom is $3 per square foot
and the sides are $1.50 per square foot.
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 193
The method of Lagrange multipliers also works for functions of three vari-
ables. That is, if we have a function f = f (x, y, z) that we want to optimize
subject to a constraint g(x, y, z) = k, the optimal point (x, y, z) lies on the
level surface S defined by the constraint g(x, y, z) = k. As we did in Preview
Activity 10.8.1, we can argue that the optimal value occurs at the level surface
f (x, y, z) = c that is tangent to S. Thus, the gradients of f and g are parallel
at this optimal point. So, just as in the two variable case, we can optimize
f = f (x, y, z) subject to the constraint g(x, y, z) = k by finding all points
(x, y, z) that satisfy ∇f = λ∇g and g(x, y, z) = k.
10.8.2 Summary
• The extrema of a function f = f (x, y) subject to a constraint g(x, y) = c
occur at points for which the contour of f is tangent to the curve that
represents the constraint equation. This occurs when
∇f = λ∇g.
10.8.3 Exercises
1. Use Lagrange multipliers to find the maximum and minimum values of
f (x, y) = 3x − 4y subject to the constraint x2 + 3y 2 = 129, if such values
exist.
maximum =
minimum =
(For either value, enter DNE if there is no such value.)
Answer 1. 3 · 3 · 3 + 4 · 4
Answer 2. −3 · 3 · 3 − 4 · 4
2. Use Lagrange multipliers to find the maximum and minimum values of
f (x, y) = x2 y+3y 2 −y, subject to the constraint x2 +y 2 ≤ 38.3333333333333
maximum =
minimum =
(For either value, enter DNE if there is no such value.)
Answer 1. 137.926
Answer 2. −59.2593
3. Find the absolute maximum and minimum of the function f (x, y) = x2 +
y 2 subject to the constraint x4 + y 4 = 6561.
As usual, ignore unneeded answer blanks, and list points in lexico-
graphic order.
Absolute minimum value:
attained at ( , ), ( ,
),
( , ), ( , ).
Absolute maximum value:
attained at ( , ), ( ,
),
( , ), ( , ).
Answer 1. 81
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 194
Answer 2. −9
Answer 3. 0
Answer 4. 0
Answer 5. −9
Answer 6. 0
Answer 7. 9
Answer 8. 9
Answer 9. 0
Answer 10. 114.551298552221
Answer 11. −7.56806773728343
Answer 12. −7.56806773728343
Answer 13. −7.56806773728343
Answer 14. 7.56806773728343
Answer 15. 7.56806773728343
Answer 16. −7.56806773728343
Answer 17. 7.56806773728343
Answer 18. 7.56806773728343
4. Find the absolute maximum and minimum of the function f (x, y) = x2 −
y 2 subject to the constraint x2 + y 2 = 361.
As usual, ignore unneeded answer blanks, and list points in lexico-
graphic order.
Absolute minimum value:
attained at ( , ) and ( ,
).
Absolute maximum value:
attained at ( , ) and ( ,
).
Answer 1. −361
Answer 2. 0
Answer 3. −19
Answer 4. 0
Answer 5. 19
Answer 6. 361
Answer 7. −19
Answer 8. 0
Answer 9. 19
Answer 10. 0
5. Find the minimum distance from the point (1, 1, 11) to the paraboloid
given by the equation z = x2 + y 2 .
Minimum distance =
Note: If you need to find roots of a polynomial of degree ≥ 3, you
may want to use a calculator of computer to do so numerically. Also be
sure that you can give a geometric justification for your answer.
√
Answer. 3.53771
6. For each value of λ the function h(x, y) = x2 + y 2 − λ(2x + 8y − 20) has
a minimum value m(λ).
(a) Find m(λ)
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 195
m(λ) =
(Use the letter L for λ in your expression.)
(b) For which value of λ is m(λ) the largest, and what is that maximum
value?
λ=
maximum m(λ) =
(c) Find the minimum value of f (x, y) = x2 + y 2 subject to the con-
straint 2x+8y = 20 using the method of Lagrange multipliers and evaluate
λ.
minimum f =
λ=
(How are these results related to your result in part (b)?)
(2·2+8·8)L2
Answer 1. 20L − 4
2·20
Answer 2. 2·2+8·8
20·20
Answer 3. 2·2+8·8
20·20
Answer 4. 2·2+8·8
2·20
Answer 5. 2·2+8·8
Answer 2. 1
Answer 3. −20
Answer 4. 1
10. Find the maximum and minimum volumes of a rectangular box whose
surface area equals 7000 square cm and whose edge length (sum of lengths
of all edges) is 440 cm.
Hint: It can be deduced that the box is not a cube, so if x, y, and z
are the lengths of the sides, you may want to let x represent a side with
x 6= y and x 6= z.
Maximum value is ,
occuring at ( ,
, ).
Minimum value is ,
occuring at ( ,
, ).
Answer 1. 34481.4814814815
Answer 2. 63.3333333333333
Answer 3. 23.3333333333333
Answer 4. 23.3333333333333
Answer 5. 25000
Answer 6. 10
Answer 7. 50
Answer 8. 50
P3 P3
11. (a) If i=1 xi = 4, find the values of x1 , x2 , x3 making i=1 xi 2 mini-
mum.
x1 , x2 , x3 =
(Give your values as a comma separated list.)
Pn(b) Generalize
2
the
Pnresult of part (a) to find the minimum value of
i=1 x i subject to i=1 xi = 4.
minimum value =
Answer 1. 1.33333, 1.33333, 1.33333
4·4
Answer 2. n
12. The Cobb-Douglas production function is used in economics to model
production levels based on labor and equipment. Suppose we have a
specific Cobb-Douglas function of the form
where x is the dollar amount spent on labor and y the dollar amount spent
on equipment. Use the method of Lagrange multipliers to determine how
much should be spent on labor and how much on equipment to maximize
productivity if we have a total of 1.5 million dollars to invest in labor and
equipment.
13. Use the method of Lagrange multipliers to find the point on the line
x − 2y = 5 that is closest to the point (1, 3). To do so, respond to the
following prompts.
a. Write the function f = f (x, y) that measures the square of the
distance from (x, y) to (1, 3). (The extrema of this function are the
same as the extrema of the distance function, but f (x, y) is simpler
to work with.)
CHAPTER 10. DERIVATIVES OF MULTIVARIABLE FUNCTIONS 197
∇f (x0 , y0 , z0 ) · r0 (t0 ) = 0,
∇g(x0 , y0 , z0 ) · r0 (t0 ) = 0, and
∇h(x0 , y0 , z0 ) · r0 (t0 ) = 0.
g(x, y, z) = c, and
h(x, y, z) = k.
for x, y, z, λ, and µ.
Use this idea to find the maximum and minium values of f (x, y, z) =
x + 2y subject to the constraints y 2 + z 2 = 8 and x + y + z = 10.
16. There is a useful interpretation of the Lagrange multiplier λ. Assume
that we want to optimize a function f with constraint g(x, y) = c. Recall
that an optimal solution occurs at a point (x0 , y0 ) where ∇f = λ∇g. As
the constraint changes, so does the point at which the optimal solution
occurs. So we can think of the optimal point as a function of the parameter
c, that is x0 = x0 (c) and y0 = y0 (c). The optimal value of f subject
to the constraint can then be considered as a function of c defined by
f (x0 (c), y0 (c)). The Chain Rule shows that
df ∂f dx0 ∂f dy0
= + .
dc ∂x0 dc ∂y0 dc
df dg
=λ .
dc dc
df
= λ.
dc
Conclude that λ tells us the rate of change of the function f as the
parameter c increases (or by approximately how much the optimal
value of the function f will change if we increase the value of c by 1
unit).
Multiple Integrals
199
CHAPTER 11. MULTIPLE INTEGRALS 200
2
0 6
y
d
c
x
a b
f (x∗11 , y11
∗
)∆A.
k. Write a double sum using summation notation that expresses the arbi-
trary sum from part (j).
CHAPTER 11. MULTIPLE INTEGRALS 202
♦
If f (x, y) ≥ 0 on the rectangle R, we may ask to find the volume of the
solid bounded above by f over R, as illustrated on the left of Figure 11.1.4.
This volume is approximated by a Riemann sum, which sums the volumes of
the rectangular boxes shown on the right of Figure 11.1.4.
z z
z = f (x, y)
y y
x x
z z
y y
x x
R ♦
Some textbooks use the notation R f (x, y) dA for a double integral. You
will see this in some of the WeBWorK problems.
z = f (x, y)
y
R
RR
RR realize the double integral R f (x, y) dA as a difference in
We can then
volumes: R f (x, y) dA tells us the volume of the solids the graph of f
bounds above the xy-plane over the rectangle R minus the volume of the
solids the graph of f bounds below the xy-plane under the rectangle R.
This is shown on the right of Figure 11.1.7.
• The average of the finitely many mn values f (x∗ij , yij
∗
) that we take in a
double Riemann sum is given by
n m
1 XX
Avgmn = f (x∗ij , yij
∗
).
mn j=1 i=1
d. Give two interpretations for the meaning of the sum you just calculated.
p
Activity 11.1.4 Let f (x, y) = 4 − y 2 on the rectangular domain R =
[1, 7] × [−2, 2]. Partition [1, 7] into 3 equal length subintervals and [−2, 2]
into 2 equal length subintervals. A table of values of f at some points in R is
given in Table 11.1.8, and a graph of f with the indicated partitions is shown
in Figure 11.1.9.
Table 11.1.8
p Table of values of
f (x, y) = 4 − y 2 .
−2 −1 0 1 2
√ √
1 0 3 2 3 0
√ √ 1
2 0 3 2 3 0
√ √
3 0 3 2 3 0 3
√ √
4 0 3 2 3 0 x
√ √
5 0
√
3 2
√
3 0 5
6 0 3 2 3 0
√ √
7 0 3 2 3 0 -2 -1 0 1
y
Figure 11.1.9 Graph of f (x, y) =
p
4 − y 2 on R.
a. Sketch the region R in the plane using the values in Table 11.1.8 as the
partitions.
b. Calculate the double Riemann sum using the given partition of R and
the values of f in the upper right corner of each subrectangle.
RR
c. Use geometry to calculate the exact value of R f (x, y) dA and compare
it to your approximation. Describe one way we could obtain a better
approximation using the given data.
We conclude this section with a list of properties of double integrals. Since
similar properties are satisfied by single-variable integrals and the arguments
for double integrals are essentially the same, we omit their justification.
11.1.4 Summary
• Let f be a continuous function on a rectangle R = {(x, y) : a ≤ x ≤ b, c ≤
y ≤ d}. The double Riemann sum for f over R is created as follows.
b−a
m . Let x0 , x1 , . . ., xm be the endpoints of these subintervals,
where a = x0 < x1 < x2 < · · · < xm = b.
◦ Partition the interval [c, d] into n subintervals of equal length ∆y =
d−c
n . Let y0 , y1 , . . ., yn be the endpoints of these subintervals, where
c = y0 < y1 < y2 < · · · < yn = d.
◦ These two partitions create a partition of the rectangle R into mn
subrectangles Rij with opposite vertices (xi−1 , yj−1 ) and (xi , yj ) for
i between 1 and m and j between 1 and n. These rectangles all have
equal area ∆A = ∆x · ∆y.
∗
◦ Choose a point (x∗ij , yij ) in each rectangle Rij . Then a double Rie-
mann sum for f over R is given by
n X
X m
f (x∗ij , yij
∗
) · ∆A.
j=1 i=1
• With terms defined as in the Double Riemann Sum, the double integral
of f over R is
ZZ n X
X m
f (x, y) dA = lim f (x∗ij , yij
∗
) · ∆A.
R m,n→∞
j=1 i=1
RR
• Two interpretations of the double integral R
f (x, y) dA are:
◦ The volume of the solids the graph of f bounds above the xy-plane
over the rectangle R minus the volume of the solids the graph of f
bounds below the xy-plane under the rectangle R;
◦ Dividing the double integral of f over R by the area of R gives us
the average value of the function f on R. If f (x, y) ≥ 0 on R, we
can interpret this average value of f on R as the height of the box
with base R that has the same volume as the volume of the surface
defined by f over R.
11.1.5 Exercises
1. Suppose f (x, y) = 25 − x2 − y 2 and R is theZ Zrectangle with vertices (0,0),
(6,0), (6,4), (0,4). In each part, estimate f (x, y) dA using Riemann
R
sums. For underestimates or overestimates, consistently use either the
lower left-hand corner or the upper right-hand corner of each rectangle in
a subdivision, as appropriate.
(a) Without subdividing R,
Underestimate =
Overestimate =
(b) By partitioning R into four equal-sized rectangles.
Underestimate =
Overestimate =
Answer 1. −648
Answer 2. 600
Answer 3. −180
Answer 4. 444
CHAPTER 11. MULTIPLE INTEGRALS 207
2. Consider the solid that lies above the square (in the xy-plane) R = [0, 1]×
[0, 1], and below the elliptic paraboloid z = 100 − x2 + 6xy − 3y 2 .
Estimate the volume by dividing R into 9 equal squares and choosing
the sample points to lie in the midpoints of each square.
Answer. 100.203703703704
3. Let R be the
√ rectangle with vertices (0, 0), (2, 0), (2, 2), and (0, 2) and let
f (x, y) = 0.333333xy. R
(a) Find reasonable upper and lower bounds for R f dA without sub-
dividing R.
upper bound =
lower bound =R
(b) Estimate R f dA three ways: by partitioning R into four subrect-
angles and evaluating f at its maximum and minimum values on each
subrectangle, and then by considering the average of these (over and un-
der) estimates. R
overestimate: RR f dA ≈
underestimate:
R R
f dA ≈
average: R f dA ≈
q
Answer 1. 4 13 · 2 · 2
Answer 2. 0
q1 q1 q
3 ·4 3 ·4
Answer 3. 44 4 + 2 2 + 1
3 · 4
q1
·4
Answer 4. 44 34
q q1 q
1
4 3 ·4 3 ·4 1
Answer 5. 8 2 4 +2 2 + 3 ·4
4. Using Riemann sums with four subdivisions in each direction, find upper
and lower bounds for the volume under the graph of f (x, y) = 6 + 3xy
above the rectangle R with 0 ≤ x ≤ 1, 0 ≤ y ≤ 6.
upper bound =
lower bound =
Answer 1. 0.25 · 1.5 · 208.5
Answer 2. 0.25 · 1.5 · 136.5
5. Consider the solid that lies above the square (in the xy-plane) R = [0, 2] ×
[0, 2],
and below the elliptic paraboloid z = 36 − x2 − 2y 2 .
(A) Estimate the volume by dividing R into 4 equal squares and choos-
ing the sample points to lie in the lower left hand corners.
(B) Estimate the volume by dividing R into 4 equal squares and choos-
ing the sample points to lie in the upper right hand corners..
(C) What is the average of the two answers from (A) and (B)?
Answer 1. 138
Answer 2. 114
Answer 3. 126
6. The figure below shows contours of g(x, y) on the region R, with 5 ≤ x ≤
11 and 2 ≤ y ≤ 8.
CHAPTER 11. MULTIPLE INTEGRALS 208
then use an integral to add the volumes of the slices. In what follows, we are
able to use single integrals to generalize this approach to handle even more
general geometric shapes.
y = 4 − x2
What is the geometric meaning of the value of A(x) relative to the surface
defined by f . (Hint: Think about the trace determined by the fixed
value of x, and consider how A(x) is related to the image at left in
Figure 11.2.2.)
z 25 z 25
20 20
15 15
10 10
5 5
-3 -3
-1 -1
1 x 1 x
-4 -2 0 y 2 -4 -2 0 y 2
Figure 11.2.2 Left: A cross section with fixed x. Right: A cross section
with fixed x and ∆x.
b. For a fixed value of x, say x∗i , what is the geometric meaning of A(x∗i ) ∆x?
(Hint: Consider how A(x∗i )∆x is related to the image at right in Fig-
ure 11.2.2.)
c. Since f is continuous on R, we can define the function A = A(x) at every
value of x in [−3, 3]. Now think about subdividing the x-interval [−3, 3]
into m subintervals, and choosing a value
Pm x∗i in each of those subintervals.
What will be the meaning of the sum i=1 A(x∗i ) ∆x?
CHAPTER 11. MULTIPLE INTEGRALS 212
R3
d. Explain why −3 A(x) dx will determine the exact value of the volume
under the surface z = f (x, y) over the rectangle R.
The function A = A(x) determines the value of the cross sectional area (by
area we mean “signed” area) in the y direction for the fixed value of x of the
solid bounded between the surface defined by f and the xy-plane.
z 25 z 25 z 25
20 20 20
15 15 15
10 10 10
5 5 5
-3 -3 -3
-1 -1 -1
1 x 1 x 1 x
-4 2 0 y 2 -4 2 0 y 2 -4 2 0 y 2
where x∗i
is any number in the subinterval [xi−1 , xi ]. Each term A(x∗i )∆x in the
sum represents an approximation of a fixed cross sectional slice of the surface
in the y direction with a fixed width of ∆x as illustrated in Figure 11.2.3. We
add the signed volumes of these slices as shown in the frames in Figure 11.2.3
to obtain an approximation of the total signed volume.
Pm in the x direction approach infinity,
As we let the number of subintervals
we can see that the Riemann sum i=1 A(x∗i )∆x approaches a limit and that
limit is the sum of signed volumes bounded by the function f on R. Therefore,
since A(x) is itself determined by an integral, we have
ZZ m Z b Z b Z d !
X
f (x, y) dA = lim A(x∗i )∆x = A(x) dx = f (x, y) dy dx.
R m→∞ a a c
i=1
is called an iterated integral, and we see that each double integral may be
represented by two single integrals.
CHAPTER 11. MULTIPLE INTEGRALS 213
The fact that integrating in either order results in the same value is known
as Fubini’s Theorem.
Fubini’s Theorem.
If f = f (x, y) is a continuous function on a rectangle R = [a, b] × [c, d],
then
ZZ Z dZ b Z bZ d
f (x, y) dA = f (x, y) dx dy = f (x, y) dy dx.
R c a a c
Note that you will be integrating with respect to y, and holding x con-
stant. Your result should be a function of x only.
b. Next, use your result from (a) along with the Fundamental Theorem of
R3
Calculus to determine the value of −3 A(x) dx.
RR
c. What is the value of R f (x, y) dA? What are two different ways we may
interpret the meaning of this value?
Activity 11.2.3 Let f (x, y) = x + y 2 on the rectangle R = [0, 2] × [0, 3].
RR
a. Evaluate R f (x, y) dA using an iterated integral. Choose an order for
integration by deciding whether you want to integrate first with respect
to x or y.
RR
b. Evaluate R f (x, y) dA using the iterated integral whose order of inte-
gration is the opposite of the order you chose in (a).
11.2.2 Summary
RR
• We can evaluate the double integral R f (x, y) dA over a rectangle R =
[a, b] × [c, d] as an iterated integral in one of two ways:
R b R d
◦ a c f (x, y) dy dx, or
R d R b
◦ c a
f (x, y) dx dy.
CHAPTER 11. MULTIPLE INTEGRALS 214
11.2.3 Exercises
R4R4
1. Evaluate the iterated integral 0 0
12x2 y 3 dxdy
Answer. 16384
R5R3
2. Evaluate the iterated integral 4 2
(3x + y)−2 dydx
Answer. 0.00394481921566777
R 1 R 13
3. Find 0 9 (x + ln y) dydx
Answer. 11.569320450974
R5R4
4. Find 0 2 xyex+y dydx
Answer. 93006.8798911913
RR
5. Calculate the double integral R
(4x + 4y + 16) dA where R is the region:
0 ≤ x ≤ 2, 0 ≤ y ≤ 2.
Answer. 96
RR
6. Calculate the double integral R
x cos(x + y) dA where R is the region:
0 ≤ x ≤ π6 , 0 ≤ y ≤ π4
Answer. 0.0767515510438744
7. Consider the solid that lies above the square (in the xy-plane) R = [0, 2] ×
[0, 2],
and below the elliptic paraboloid z = 49 − x2 − 3y 2 .
(A) Estimate the volume by dividing R into 4 equal squares and choos-
ing the sample points to lie in the lower left hand corners.
(B) Estimate the volume by dividing R into 4 equal squares and choos-
ing the sample points to lie in the upper right hand corners..
(C) What is the average of the two answers from (A) and (B)?
(D) Using iterated integrals, compute the exact value of the volume.
Answer 1. 188
Answer 2. 156
Answer 3. 172
Answer 4. 174.666666666667
Z 3 Z −4 ZZ
8. If f (x)dx = −2 and g(x)dx = 4, what is the value of f (x)g(y)dA
1 −5 D
where D is the rectangle: 1 ≤ x ≤ 3, −5 ≤ y ≤ −4?
Answer. −8
9. Find the average value of f (x, y) = 4x6 y 3 over the rectangle R with ver-
tices (−3, 0), (−3, 6), (3, 0), (3, 6).
Average value =
Answer. 22494.8571428571
√
10. Find the average value of f (x, y) = 7ey x + ey over the rectangle R =
[0, 8] × [0, 5].
Average value =
Answer. 1745.04372328632
CHAPTER 11. MULTIPLE INTEGRALS 215
b y
0.5
a
x
x
0.5 1.0
Figure 11.3.1 Left: The tetrahedron T . Right: Projecting T onto the xy-
plane.
z = 1 − (x + y).
The issue that is new here is how we find the limits on the integrals;
note that the outer integral’s limits are in x, while the inner ones are in
y, since we have chosen dA = dy dx. To see the domain over which we
need to integrate, think of standing way above the tetrahedron looking
straight down on it, which means we are projecting the entire tetrahedron
onto the xy-plane. The resulting domain is the triangular region shown
CHAPTER 11. MULTIPLE INTEGRALS 217
(Hint: Consider the cross sectional slice shown at right in Figure 11.3.1.)
c. Explain why it makes sense to now write the volume integral in the form
Z x=? Z y=? Z x=1 Z y=1−x
1 − (x + y) dy dx = 1 − (x + y) dy dx.
x=? y=? x=0 y=0
and compare to your result from part (a). (As with iterated integrals
over rectangular regions, start with the inner integral.)
2 y 2 y
1 R 1
D D
x x
-2 -1 1 2 -2 -1 1 2
-1 -1
-2 -2
Figure 11.3.2 Left: A non-rectangular domain. Right: Enclosing this domain
in a rectangle.
We can enclose D in a rectangular domain R as shown at right in Fig-
ure 11.3.2 and extend the function f to be defined over R in order to be able
to use the definition of the double integral over a rectangle. We extend f in
such a way that its values at the points in R that are not in D contribute 0 to
the value of the integral. In other words, define a function F = F (x, y) on R
as (
f (x, y), if (x, y) ∈ D,
F (x, y) = .
0, if (x, y) 6∈ D
We then say that the double integral of f over D is the same as the double
integral of F over R, and thus
ZZ ZZ
f (x, y) dA = F (x, y) dA.
D R
2 D 2 D 2 D
1 1 1
x x x
1 2 3 1 2 3 1 2 3
Since we are integrating with respect to y first, the iterated integral has
the form ZZ Z x=b
x2 y dA = A(x) dx,
D x=a
x=2
9 x5
=
8 5 x=0
9 32
=
8 5
36
= .
5
Since we are integrating with respect to x first, the iterated integral has
the form ZZ Z d
x2 y dA = A(y) dy,
D c
where A(y) is a cross sectional area of the solid in the x direction. Several
slices of the domain — perpendicular to the y-axis — are shown at right
in Figure 11.3.4. On a slice with fixed y value, the x values are bounded
below by the x coordinate on the hypotenuse of the right triangle and
above by 2. So h2 (y) = 2; to find h1 (y), we need to write the hypotenuse
as a function of y. Solving the earlier equation we have for the hypotenuse
(y = 23 x) for x gives us x = 23 y. This makes h1 (y) = 23 y. The lowest
horizontal cross section is at y = 0 and the uppermost one is at y = 3,
so we have c = 0 and d = 3. Therefore,
ZZ Z y=3 Z x=2
2
x y dA = x2 y dx dy.
D y=0 x=(2/3)y
c. Evaluate one of the two iterated integrals above. Explain what the value
you obtained tells you.
d. Set up and evaluate a single definite integral to determine the exact area
of D, A(D).
e. Determine the exact average value of f (x, y) = 4x + 10y over D.
R x=4 R y=2 2
Activity 11.3.4 Consider the iterated integral x=0 y=x/2 ey dy dx.
2
a. Explain why we cannot find a simpleR antiderivative for ey with respect to
x=4 R y=2 y2
y, and thus are unable to evaluate x=0 y=x/2 e dy dx in the indicated
order using the Fundamental Theorem of Calculus.
2 R x=4 R y=2 2
b. Given that D ey dA = x=0 y=x/2 ey dy dx, sketch the region of inte-
RR
gration, D.
c. Rewrite the given iterated integral in the opposite order, using dA =
CHAPTER 11. MULTIPLE INTEGRALS 222
11.3.2 Summary
RR
• For a double integral D f (x, y) dA over a non-rectangular region D, we
enclose D in a rectangle R and then extend integrand f to a function
F so that F (x, y) = 0 at all points in R RRoutside of D and F (x, y) =
fRR(x, y) for all points in D. We then define D f (x, y) dA to be equal to
R
F (x, y) dA.
• In an iterated double integral, the limits on the outer integral must be
constants while the limits on the inner integral must be constants or in
terms of only the remaining variable. In other words, an iterated double
integral has one of the following forms (which result in the same value):
Z x=b Z y=g2 (x)
f (x, y) dy dx,
x=a y=g1 (x)
where g1 = g1 (x) and g2 = g2 (x) are functions of x only and the region
D is described by the inequalities g1 (x) ≤ y ≤ g2 (x) and a ≤ x ≤ b or
Z y=d Z x=h2 (y)
f (x, y) dx dy,
y=c x=h1 (y)
where h1 = h1 (y) and h2 = h2 (y) are functions of y only and the region
D is described by the inequalities h1 (y) ≤ x ≤ h2 (y) and c ≤ y ≤ d.
11.3.3 Exercises
ZZ
1. Evaluate the double integral I = xy dA where D is the triangular
D
region with vertices (0, 0), (6, 0), (0, 2).
Answer. 6
ZZ
2. Evaluate the double integral I = xy dA where D is the triangular
D
region with vertices (0, 0), (1, 0), (0, 6).
Answer. 1.5
3. Evaluate
R 1R 3 the integral by reversing the order of integration.
x2
0 3y
e dxdy =
Answer. 1350.34732126256
4. Decide, without calculation, if each of the integrals below are positive,
negative, or zero. Let D be the region inside the unit circle centered at
the origin. Let T, B, R, and L denote the regions enclosed by the top half,
the bottom half, the right half, and the left half of unit circle, respectively.
ZZ
(a) (y 3 + y 5 ) dA
T
CHAPTER 11. MULTIPLE INTEGRALS 223
ZZ
(b) (y 3 + y 5 ) dA
L
ZZ
(c) (y 3 + y 5 ) dA
D
ZZ
(d) (y 3 + y 5 ) dA
R
ZZ
(e) (y 3 + y 5 ) dA
B
2 2
5. The region W lies below the surface f (x, y) = 8e−(x−3) −y and above the
disk x2 + y 2 ≤ 16 in the xy-plane.
(a) Think about what the contours of f look like. You may want to
using f (x, y) = 1 as an example. Sketch a rough contour diagram on a
separate sheet of paper.
(b) Write an integral giving the area of the cross-section of W in the
plane x = 3.
Rb
Area = a d
,
where a = and b =
(c) Use your work from (b) to write an iterated double integral giving
the volume of W , using the work from (b) to inform the construction of
the inside integral.
RbRd
Volume = a c
d d ,
where a = ,b= c= and d =
2
Answer 1. 8e−y
Answer 2. y
Answer 3. −2.64575
Answer 4. 2.64575
2
−y 2
Answer 5. 8e−(x−3)
Answer 6. y
Answer 7. x
Answer 8. −4
Answer 9. 4
√
Answer 10. − 16 − x2
√
Answer 11. 16 − x2
6. Set up a double integral in rectangular coordinates for calculating the
volume of the solid under the graph of the function f (x, y) = 22 − x2 − y 2
and above the plane z = 6.
Instructions: Please enter the integrand in the first answer box. De-
pending on the order of integration you choose, enter dx and dy in either
order into the second and third answer boxes with only one dx or dy in
each box. Then, enter the limits of integration.
Z BZ D
A C
CHAPTER 11. MULTIPLE INTEGRALS 224
A=
B=
C=
D=
p p
Answer. 16 − x2 − y 2 ; dx; dy; −4; 4; − 16 − y 2 ; 16 − y 2
7. Find the volume of the solid bounded by the planes x = 0, y = 0, z = 0,
and x + y + z = 6.
Answer. 36
Z 6 Z √36−y
8. Consider the integral f (x, y)dxdy. If we change the order of
0 0
integration we obtain the sum of two integrals:
Z b Z g2 (x) Z d Z g4 (x)
f (x, y)dydx + f (x, y)dydx
a g1 (x) c g3 (x)
a= b=
g1 (x) = g2 (x) =
c= d=
g3 (x) = g4 (x) =
Answer 1. 0
Answer 2. 5.47722557505166
Answer 3. 0
Answer 4. 6
Answer 5. 5.47722557505166
Answer 6. 6
Answer 7. 0
Answer 8. 36 − x2
9. A pile of earth standing on flat ground has height 36 meters. The ground
is the xy-plane. The origin is directly below the top of the pile and the
z-axis is upward. The cross-section at height z is given by x2 + y 2 = 36 − z
for 0 ≤ z ≤ 36, with x, y, and z in meters.
(a) What equation gives the edge of the base of the pile?
x2 + y 2 = 36
x + y = 36
x2 + y 2 = 6
x+y =6
None of the above
x2 + y 2 = 25
√
x2 + y 2 = 31
x2 + y 2 = 31
x2 + y 2 = 5
None of the above
CHAPTER 11. MULTIPLE INTEGRALS 225
12. The temperature at any point on a metal plate in the xy-plane is given
by T (x, y) = 100 − 4x2 − y 2 , where x and y are measured in inches and
T in degrees Celsius. Consider the portion of the plate that lies on the
region D that is the finite region that lies between the parabolas x = y 2
and x = 3 − 2y 2 .
a. Construct a labeled sketch of the region D.
RR
b. Set up an iterated integral whose value is D T (x, y) dA, using dA =
dxdy. (Hint: It is possible that more than one integral is needed.)
RR
c. Set up an integrated integral whose value is D T (x, y) dA, using
dA = dydx. (Hint: It is possible that more than one integral is
needed.)
d. Use the Fundamental Theorem of Calculus to evaluate the integrals
you determined in (b) and (c).
e. Determine the exact average temperature, TAVG(D) , over the region
D.
13. Consider the solid that is given by the following description: the base is
the given region D, while the top is given by the surface z = p(x, y). In
each setting below, set up, but do not evaluate, an iterated integral whose
value is the exact volume of the solid. Include a labeled sketch of D in
each case.
a. D is the interior of the quarter circle of radius 2, centered at the
origin, that lies in the second quadrant of the plane; p(x, y) = 16 −
x2 − y 2 .
b. D is the finite region between the line y = x + 1 and the parabola
y = x2 ; p(x, y) = 10 − x − 2y.
c. D is the triangular region with vertices (1, 1), (2, 2), and (2, 3);
p(x, y) = e−xy .
√
d. D is the region bounded by the y-axis, y = 4 and x = y; p(x, y) =
p
1 + x2 + y 2 .
Z x=4 Z y=2
14. Consider the iterated integral I = √
cos(y 3 ) dy dx.
x=0 y= x
a. Sketch the region of integration.
d. Determine the exact average value of cos(y 3 ) over the region D that
is determined by the iterated integral I.
CHAPTER 11. MULTIPLE INTEGRALS 227
• Given a mass density function on a lamina, how can we find the lamina’s
center of mass?
• What is a joint probability density function? How do we determine the
probability of an event if we know a probability density function?
So far, we have interpreted the double RR integral of a function f over a do-
main D in two different ways. First, D f (x, y) dA tells us a difference of
volumes — the volume the surface defined by f bounds above the xy-plane on
D minus
RR the volume the surface bounds below the xy-plane on D. In addition,
1
A(D) D
f (x, y) dA determines the average value of f on D. In this section, we
investigate several other applications of double integrals, using the integration
process as seen in Preview Activity 11.4.1: we partition into small regions,
approximate the desired quantity on each small region, then use the integral
to sum these values exactly in the limit.
The following preview activity explores how a double integral can be used to
determine the density of a thin plate with a mass density distribution. Recall
that in single-variable calculus, we considered a similar problem and computed
the mass of a one-dimensional rod with a mass-density distribution. There, as
here, the key idea is that if density is constant, mass is the product of density
and volume.
Preview Activity 11.4.1 Suppose that we have a flat, thin object (called a
lamina) whose density varies across the object. We can think of the density
on a lamina as a measure of mass per unit area. As an example, consider a
circular plate D of radius 1 cm centered at the origin whose density δ varies
depending on the distance from its center so that the density in grams per
square centimeter at point (x, y) is
δ(x, y) = 10 − 2(x2 + y 2 ).
11.4.1 Mass
Density is a measure of some quantity per unit area or volume. For example,
we can measure the human population density of some region as the number
of humans in that region divided by the area of that region. In physics, the
mass density of an object is the mass of the object per unit area or volume.
As suggested by Preview Activity 11.4.1, the following holds in general.
3 y
2
D
1
x
-1 1 2 3 4
-1
-2
-3
Figure 11.4.1 A half disk lamina.
11.4.2 Area
If we consider the situation where the mass-density distribution is constant,
we can also see how a double integral may be used to determine the area of
a region. Assuming that δ(x, y) = 1 over a closed bounded RR region D, where
the units of δ are “mass per unit of area,” it follows that D 1 dA is the mass
of the lamina. But since the density is constant, the numerical value of the
integral is simply the area.
As the following activity demonstrates, we can also see this fact by consid-
ering a three-dimensional solid whose height is always 1.
Activity 11.4.3 Suppose we want to find the area of the bounded region D
between the curves
y = 1 − x2 and y = x − 1.
2 y
1
x
-1 D 1 2 3
-1
-2
-3
b. Evaluate the iterated integral from (a). What does the result tell you?
We now formally state the conclusion from our earlier discussion and Ac-
tivity 11.4.3.
y
(x2 , y2 )
(x3 , y3 )
(x, y)
(x4 , y4 )
(x1 , y1 )
x
δ(x∗ij , yij
∗
)∆A
2 y
x
1
Figure 11.4.4 The lamina bounded by the x-axis and the lines x = 1 and
y = 2x in the first quadrant.
11.4.4 Probability
Calculating probabilities is a very important application of integration in the
physical, social, and life sciences. To understand the basics, consider the game
of darts in which a player throws a dart at a board and tries to hit a particular
target. Let us suppose that a dart board is in the form of a disk D with radius
10 inches. If we assume that a player throws a dart at random, and is not
aiming at any particular point, then it is equally probable that the dart will
strike any single point on the board. For instance, the probability that the dart
1
will strike a particular 1 square inch region is 100π , or the ratio of the area
of the desired target to the total area of D (assuming that the dart thrower
always hits the board itself at some point). Similarly, the probability that the
dart strikes a point in the disk D3 of radius 3 inches is given by the area of D3
divided by the area of D. In other words, the probability that the dart strikes
the disk D3 is ZZ
9π 1
= dA.
100π D3 100π
1
The integrand, 100π , may be thought of as a distribution function, describ-
ing how the dart strikes are distributed across the board. In this case the
distribution function is constant since we are assuming a uniform distribution,
but we can easily envision situations where the distribution function varies. For
example, if the player is fairly good and is aiming for the bulls eye (the center
of D), then the distribution function f could be skewed toward the center, say
2
+y 2 )
f (x, y) = Ke−(x
for some constant positive K. If we assume that the player is consistent enough
so that the dart always strikes the board, then the probability that the dart
strikes the board somewhere is 1, and the distribution function f will have to
CHAPTER 11. MULTIPLE INTEGRALS 233
satisfy1 ZZ
f (x, y) dA = 1.
D
For such a function f , the probability that the dart strikes in the disk D1
of radius 1 would be ZZ
f (x, y) dA.
D1
Indeed, the probability that the dart strikes in any region R that lies within
D is given by ZZ
f (x, y) dA.
R
The preceding discussion highlights the general idea behind calculating
probabilities. We assume we have a joint probability density function f , a
function of two independent variables x and y defined on a domain D that
satisfies the conditions
• f (x, y) ≥ 0 for all x and y in D,
• the probability that x is between some values a and b while y is between
some values c and d is given by
Z b Z d
f (x, y) dy dx,
a c
Note that it is possible that D could be an infinite region and the limits
on the integral in Equation (11.4.1) could be infinite. When we have such a
probability density function f = f (x, y), the probability that the point (x, y)
is in some region R contained in the domain D (the notation we use here is
“P ((x, y) ∈ R)”) is determined by
ZZ
P ((x, y) ∈ R) = f (x, y) dA.
R
or that Z ∞ Z ∞
f (x, y) dy dx = 1.
0 0
Use your knowledge of improper integrals to verify that f is indeed a
probability density function.
b. Assume that the smoke detector fails only if both of the supplied compo-
nents fail. To determine the probability that a randomly selected detector
will fail within one year, we will need to determine the probability that
the life span of each component is between 0 and 1 years. Set up an
appropriate iterated integral, and evaluate the integral to determine the
probability.
c. What is the probability that a randomly chosen smoke detector will fail
between years 3 and 7?
d. Suppose that the manufacturer determines that one of the components is
more likely to fail than the other, and hence conjectures that the proba-
bility density function is instead f (x, y) = Ke−x e−2y . What is the value
of K?
11.4.5 Summary
RR mass of a lamina D with a mass density function δ = δ(x, y) is
• The
D
δ(x, y) dA.
• The area of a region D in the plane has the same numerical value as the
volume ofRRa solid of uniform height 1 and base D, so the area of D is
given by D 1 dA.
• The center of mass, (x, y), of a continuous lamina with a variable density
δ(x, y) is given by
RR RR
D
xδ(x, y) dA yδ(x, y) dA
x = RR and y = RRD .
D
δ(x, y) dA D
δ(x, y) dA
11.4.6 Exercises
1. The masses mi are located at the points Pi . Find the center of mass of
the system.
m1 = 1, m2 = 8, m3 = 6.
P1 = (1, −4), P2 = (−6, 7), P3 = (9, −7).
x̄=
ȳ=
Answer 1. 0.466666666666667
Answer 2. 0.666666666666667
2. Find the centroid (x̄, ȳ) of the triangle with vertices at (0, 0), (5, 0), and
(0, 2).
x̄=
CHAPTER 11. MULTIPLE INTEGRALS 235
ȳ=
Answer 1. 1.66666666666667
Answer 2. 0.666666666666667
3. Find the mass of the rectangular region 0 ≤ x ≤ 2, 0 ≤ y ≤ 1 with density
function ρ (x, y) = 1 − y.
Answer. 1
4. Find the mass of the triangular region with vertices (0, 0), (3, 0), and (0,
5), with density function ρ (x, y) = x2 + y 2 .
Answer. 42.5
5. A lamina occupies the region inside the circle x2 + y 2 = 8y but outside the
circle x2 + y 2 = 16. The density at each point is inversely proportional to
its distance from the orgin.
Where is the center of mass?
( , )
Answer 1. 0
Answer 2. 5.05816623182515
6. A sprinkler distributes water in a circular pattern, supplying water to a
depth of e−r feet per hour at a distance of r feet from the sprinkler.
A. What is the total amount of water supplied per hour inside of a
circle of radius 6?
f t3 per hour
B. What is the total amount of water that goes through the sprinkler
per hour?
f t3 per hour
Answer 1. 6.1741640923836
Answer 2. 6.28318530717959
1
7. Let p be the joint density function such that p(x, y) = 144 xy in R, the
rectangle 0 ≤ x ≤ 6, 0 ≤ y ≤ 4, and p(x, y) = 0 outside R. Find the
fraction of the population satisfying the constraint x + y ≤ 10
fraction =
Answer. 1
8. A lamp has two bulbs, each of a type with an average lifetime of 10
hours. The probability density function for the lifetime of a bulb is f (t) =
1 −t/10
10 e , t ≥ 0.
What is the probability that both of the bulbs will fail within 3 hours?
Answer. 0.0671751947305907
9. For the following two functions p(x, y), check whether p is a joint density
function. Assume p(x, y) = 0 outside the region R.
(a) p(x, y) = 3, where R is −1 ≤ x ≤ 1, 0 ≤ y ≤ 0.5.
p(x, y) ( is a joint density function is not a joint density func-
tion)
(b) p(x, y) = 1, where R is 1 ≤ x ≤ 1.5, 0 ≤ y ≤ 2.
p(x, y) ( is a joint density function is not a joint density func-
tion)
Then, for the region R given by −1 ≤ x ≤ 3, 0 ≤ y ≤ 3, what constant
function p(x, y) is a joint density function?
p(x, y) =
Answer 1. is not a joint density function
CHAPTER 11. MULTIPLE INTEGRALS 236
c. Determine the exact center of mass of the plate. Draw and label the
point you find on your sketch from (a).
d. What is the average density of the plate? Include units on your
answer.
12. Let D be a half-disk lamina of radius 3 in quadrants IV and I, centered
at the origin as in Activity 11.4.2. Assume the density at point (x, y) is
equal to x.
a. Before doing any calculations, what do you expect the y-coordinate
of the center of mass to be? Why?
b. Set up iterated integral expressions which, if evaluated, will deter-
mine the exact center of mass of the lamina.
c. Use appropriate technology to evaluate the integrals to find the cen-
ter of mass numerically.
13. Let x denote the time (in minutes) that a person spends waiting in a
checkout line at a grocery store and y the time (in minutes) that it takes
to check out. Suppose the joint probability density for x and y is
1 −x/4−y/2
f (x, y) = e .
8
b. Part (a) indicates that the two pieces of information completely deter-
mine the location of a point: either the traditional (x, y) coordinates, or
alternately, the distance r from the point to the origin along with the
angle θ that the line through the origin and the point makes with the
positive x-axis. We write “(r, θ)” to denote the point’s location in its
polar coordinate representation. Find polar coordinates for the points
with the given rectangular coordinates.
i. (0, −1) ii. (−2, 0) iii. (−1, 1)
c. For each of the following points whose coordinates are given in polar
form, determine the rectangular coordinates of the point.
√ 5π
i. (5, π4 ) ii. (2, 5π
6 ) iii. ( 3, 3 )
π/2
P 3π/4 π/4
π 0
r y 1 2 3 4
5π/4 7π/4
θ
3π/2
x
Figure 11.5.1 The polar coordinates of a point and the polar coordinate grid.
Trigonometry and the Pythagorean Theorem allow for straightforward con-
version from rectangular to polar, and vice versa.
d. What does the region defined by 1 ≤ r ≤ 3 (where θ can have any value)
look like? (Hint: Compare to your response from part (a).)
e. What does the region defined by 1 ≤ r ≤ 3 and π/4 ≤ θ ≤ π/2 look like?
f. Consider the curve r = sin(θ). For some values of θ we will have r < 0.
In these situations, we plot the point (r, θ) as (|r|, θ + π) (in other words,
when r < 0, we reflect the point through the origin). With that in mind,
what do you think the graph of r = sin(θ) looks like? Plot this curve
using technology and compare to your intuition.
where D is the unit disk. While we cannot directly evaluate this integral in
rectangular coordinates, a change to polar coordinates will convert it to one
we can easily evaluate. ZZ
We have seen how to evaluate a double integral f (x, y) dA as an iter-
D
ated integral of the form
Z b Z g2 (x)
f (x, y) dy dx
a g1 (x)
CHAPTER 11. MULTIPLE INTEGRALS 240
ri+1
ri+1
ri
ri
θj+1
θj
y to being in terms of r and θ, but we also have to change the area element
to dA = r dr dθ in polar coordinates. As we saw in Activity 11.5.3, the reason
the additional factor of r in the polar area element is due to the fact that
in polar coordinates, the cross sectional area element increases as r increases,
while the cross sectional area
RR element in rectangular coordinates is constant.
So, given a double integral D f (x, y) dA in rectangular coordinates, to write a
corresponding iterated integral in polar coordinates, we replace x with r cos(θ),
y with r sin(θ) and dA with r dr dθ. Of course, we need to describe the region
D in polar coordinates as well. To summarize:
2 2
Example 11.5.3RRLet f (x, y) = ex +y on the disk D = {(x, y) : x2 + y 2 ≤ 1}.
We will evaluate D f (x, y) dA. RR
In rectangular coordinates the double integral D f (x, y) dA can be written
as the iterated integral
√
ZZ Z x=1 Z y= 1−x2
2
+y 2
f (x, y) dA = √ ex dy dx.
D x=−1 y=− 1−x2
2 2
We cannot evaluate this iterated integral, because ex +y does not have
an elementary antiderivative with respect to either x or y. However, since
r2 = x2 + y 2 and the region D is circular, it is natural to wonder whether
converting to polar coordinates will allow us to evaluate the new integral. To
do so, we replace x with r cos(θ), y with r sin(θ), and dy dx with r dr dθ to
obtain ZZ ZZ
2
f (x, y) dA = er r dr dθ.
D D
The disc D is described in polar coordinates by the constraints 0 ≤ r ≤ 1
and 0 ≤ θ ≤ 2π. Therefore, it follows that
ZZ Z θ=2π Z r=1
r2 2
e r dr dθ = er r dr dθ.
D θ=0 r=0
1 θ=2π
Z
= (e − 1) dθ
2 θ=0
Z θ=2π
1
= (e − 1) dθ
2 θ=0
θ=2π
1
= (e − 1) [θ]
2 θ=0
= π(e − 1).
While there is no firm rule for when polar coordinates can or should be
used, they are a natural alternative anytime the domain of integration may be
CHAPTER 11. MULTIPLE INTEGRALS 242
expressed
psimply in polar form, and/or when the integrand involves expressions
such as x2 + y 2 .
Activity 11.5.4 Let f (x, y) = x + y and D = {(x, y) : x2 + y 2 ≤ 4}.
a. Sketch the region D and then write the double integral of f over D as
an iterated integral in rectangular coordinates.
b. Write the double integral of f over D as an iterated integral in polar
coordinates.
c. Evaluate one of the iterated integrals. Why is the final value you found
not surprising?
Activity 11.5.5 Consider the circle given by x2 + (y − 1)2 = 1 as shown in
Figure 11.5.4.
2 y
x
-1 1
a. Determine a polar curve in the form r = f (θ) that traces out the circle
x2 + (y − 1)2 = 1. (Hint: Recall that a circle centered at the origin of
radius r can be described by the equations x = r cos(θ) and y = r sin(θ).)
p
b. Find the exact average value of g(x, y) = x2 + y 2 over the interior of
the circle x2 + (y − 1)2 = 1.
c. Find the volume under the surface h(x, y) = x over the region D, where
D is the region bounded above by the line y = x and below by the circle
(this is the shaded region in Figure 11.5.4).
d. Explain why in both (b) and (c) it is advantageous to use polar coordi-
nates.
11.5.3 Summary
• The polar representation of a point P is the ordered pair (r, θ) where r
is the distance from the origin to P and θ is the angle the ray through
the origin and P makes with the positive x-axis.
• The polar coordinates r and θ of a point (x, y) in rectangular coordinates
satisfy p y
r = x2 + y 2 and tan(θ) = ;
x
CHAPTER 11. MULTIPLE INTEGRALS 243
dA = r dr dθ.
RR
• To convert the double integral D f (x, y) dA to an iterated integral in
polar coordinates, we substitute r cos(θ) for x, r sin(θ) for y, and r dr dθ
for dA to obtain the iterated integral
ZZ
f (r cos(θ), r sin(θ)) r dr dθ.
D
11.5.4 Exercises
1. For each set of Polar coordinates, match the equivalent Cartesian coordi-
nates.
√
2. (a) The Cartesian coordinates of a point are (−1, − 3).
(i) Find polar coordinates (r, θ) of the point, where r > 0 and 0 ≤ θ <
2π.
r=
θ=
(ii) Find polar coordinates (r, θ) of the point, where r < 0 and 0 ≤
θ < 2π.
r=
θ=
(b) The Cartesian coordinates of a point are (−2, 3).
(i) Find polar coordinates (r, θ) of the point, where r > 0 and 0 ≤ θ <
2π.
r=
θ=
(ii) Find polar coordinates (r, θ) of the point, where r < 0 and 0 ≤
θ < 2π.
r=
θ=
Answer 1. 2
Answer 2. 4.18879020478639
Answer 3. −2
Answer 4. 1.0471975511966
Answer 5. 3.60555127546399
Answer 6. 2.15879893034246
Answer 7. −3.60555127546399
Answer 8. 5.30039158393226
3. (a) You are given the point (1, π/2) in polar coordinates.
(i) Find another pair of polar coordinates for this point such that r > 0
and 2π ≤ θ < 4π.
r=
θ=
(ii) Find another pair of polar coordinates for this point such that
r < 0 and 0 ≤ θ < 2π.
CHAPTER 11. MULTIPLE INTEGRALS 244
r=
θ=
(b) You are given the point (−2, π/4) in polar coordinates.
(i) Find another pair of polar coordinates for this point such that r > 0
and 2π ≤ θ < 4π.
r=
θ=
(ii) Find another pair of polar coordinates for this point such that
r < 0 and −2π ≤ θ < 0.
r=
θ=
(c) You are given the point (3, 2) in polar coordinates.
(i) Find another pair of polar coordinates for this point such that r > 0
and 2π ≤ θ < 4π.
r=
θ=
(ii) Find another pair of polar coordinates for this point such that
r < 0 and 0 ≤ θ < 2π.
r=
θ=
Answer 1. 1
Answer 2. 7.85398163397448
Answer 3. −1
Answer 4. 4.71238898038469
Answer 5. 2
Answer 6. 10.2101761241668
Answer 7. −2
Answer 8. −5.49778714378214
Answer 9. 3
Answer 10. 8.28318530717959
Answer 11. −3
Answer 12. 5.14159265358979
4. Decide if the points given in polar coordinates are the same. If they are
the same, enter T . If they are different, enter F .
a.) (5, π3 ), (−5, −π 3 )
b.) (2, 35π
4 ), (2, − 35π
4 )
4π
c.) (0, 5π), (0, 4 )
π
d.) (1, 141π
4 ), (−1, 4 )
92π −π
e.) (2, 3 ), (−2, 3 )
f.)(5, 15π), (−5, 15π)
Answer 1. F
Answer 2. F
Answer 3. T
Answer 4. T
Answer 5. T
Answer 6. F
5. A curve with polar equation
17
r=
7 sin θ + 57 cos θ
CHAPTER 11. MULTIPLE INTEGRALS 245
With a = ,b= ,
c= , and d = ,
RbRd
integral = a c d d
(b) The region
With a = ,b= ,
c= , and d = ,
RbRd
integral = a c d d
Answer 1. 0
Answer 2. 4
Answer 3. 4
Answer 4. 7
Answer 5. f (x, y)
Answer 6. y
Answer 7. x
Answer 8. 1.0472
Answer 9. 1.5708
Answer 10. 0
Answer 11. 2
Answer 12. rf (r cos(t) , r sin(t))
Answer 13. r
Answer 14. t
10. A Cartesian equation for the polar equation r = 3 can be written as:
x2 + y 2 =
Answer. 9
11. Using polar coordinates, evaluate the integral which gives the area which
lies in the first quadrant between the circles x2 +y 2 = 36 and x2 −6x+y 2 =
CHAPTER 11. MULTIPLE INTEGRALS 247
0.
Answer. 14.1371669411541
12. (a) Graph r = 1/(4 cos θ) for −π/2 ≤ θ ≤ π/2 and r = 1. Then write
an iterated integral in polar coordinates representing the area inside the
curve r = 1 and to the right of r = 1/(4 cos θ). (Use t for θ in your work.)
With a = ,b= ,
c= , and d = ,
RbRd
area = a c d d
(b) Evaluate your integral to find the area.
area =
Answer 1. − cos−1 14
Answer 2. cos−1 14
1
Answer 3. 4 cos(t)
Answer 4. 1
Answer 5. r
Answer 6. r
Answer 7. t
tan(cos−1( 41 ))
Answer 8. cos−1 1
4 − 42
ZZ
13. Using polar coordinates, evaluate the integral sin(x2 + y 2 )dA where
R
R is the region 4 ≤ x2 + y 2 ≤ 25.
Answer. −5.16743746933879
14. Sketch the region of integration for the following integral.
Z π/4 Z 5/ cos(θ)
f (r, θ) r dr dθ
0 0
The region of integration is bounded by
y = 0, y = x, and y = 5
p
y = 0, x = 25 − y 2 , and y = 5
√
y = 0, y = 25 − x2 , and x = 5
y = 0, y = x, and x = 5
A C
A=
B=
C=
D=
CHAPTER 11. MULTIPLE INTEGRALS 248
(b) Evaluate the integral and find the volume. Your answer will be in
terms of a.
Volume V =
(c) What does the volume approach as a → ∞?
lim V =
a→∞
2
Answer 1. e−r r; dr; dtheta; 0; 2π; 0; a
2
Answer 2. π 1 − e−a
Answer 3. 3.14159
R 0 R 0√ y
17. Consider the iterated integral I = −3 − 9−y 2 x2 +y 2 +1
dx dy.
a. Sketch (and label) the region of integration.
b. Convert the given iterated integral to one in polar coordinates.
c. Evaluate the iterated integral in (b).
d. Find the exact volume of the solid that lies under the surface z =
8 − x2 − y 2 and over the unit disk, D.
19. For each of the following iterated integrals,
• sketch and label the region of integration,
• convert the integral to the other coordinate system (if given in polar,
to rectangular; if given in rectangular, to polar), and
• choose one of the two iterated integrals to evaluate exactly.
R 3π/2 R 3
a. π 0
r3 dr dθ
R 2 R √1−(x−1)2 p
b. 0
√ 2
x2 + y 2 dy dx
− 1−(x−1)
R π/2 R sin(θ) √
c. 0 0
r 1 − r2 dr dθ.
R √2/2 R √1−y2
d. 0 y
cos(x2 + y 2 ) dx dy.
CHAPTER 11. MULTIPLE INTEGRALS 249
form a surface in space. The equations x = x(s, t), y = y(s, t), and z = z(s, t)
are the parametric equations for the surface, or a parametrization of the surface.
In Preview Activity 11.6.1 we investigate how to parameterize a cylinder and
a cone.
Preview Activity 11.6.1 Recall the standard parameterization of the unit
circle that is given by
where 0 ≤ t ≤ 2π.
a. Determine a parameterization of the circle of radius 1 in R3 that has its
center at (0, 0, 1) and lies in the plane z = 1.
b. Determine a parameterization of the circle of radius 1 in 3-space that has
its center at (0, 0, −1) and lies in the plane z = −1.
c. Determine a parameterization of the circle of radius 1 in 3-space that has
its center at (0, 0, 5) and lies in the plane z = 5.
d. Taking into account your responses in (a), (b), and (c), describe the graph
that results from the set of parametric equations
z z
r
a
x
b
y
Now using our earlier parametric equations for x(s) and z(s) for the original
smaller circle, we have an overall parameterization of the torus given by
To trace out the entire torus, we require that the parameters vary through
the values 0 ≤ s ≤ 2π and 0 ≤ t ≤ 2π.
Activity 11.6.2 In this activity, we seek a parametrization of the sphere of
radius R centered at the origin, as shown on the left in Figure 11.6.5. Notice
that this sphere may be obtained by revolving a half-circle contained in the
xz-plane about the z-axis, as shown on the right.
CHAPTER 11. MULTIPLE INTEGRALS 252
z z
R
x
y
In the usual way, we slice the domain into small rectangles. In particular,
we partition the interval [a, b] into m subintervals of length ∆s = b−a
n and let
s0 , s1 , . . ., sm be the endpoints of these subintervals, where a = s0 < s1 <
CHAPTER 11. MULTIPLE INTEGRALS 253
s2 < · · · < sm = b. Also partition the interval [c, d] into n subintervals of equal
length ∆t = d−c n and let t0 , t1 , . . ., tn be the endpoints of these subintervals,
where c = t0 < t1 < t2 < · · · < tn = d. These two partitions create a partition
of the rectangle R = [a, b] × [c, d] in st-coordinates into mn sub-rectangles Rij
with opposite vertices (si−1 , tj−1 ) and (si , tj ) for i between 1 and m and j
between 1 and n. These rectangles all have equal area ∆A = ∆s · ∆t.
Now we want to think about the small piece of area on the surface itself
that lies above one of these small rectangles in the domain. Observe that if we
increase s by a small amount ∆s from the point (si−1 , tj−1 ) in the domain, then
r changes by approximately rs (si−1 , tj−1 )∆s. Similarly, if we increase t by a
small amount ∆t from the point (si−1 , tj−1 ), then r changes by approximately
rt (si−1 , tj−1 )∆t. So we can approximate the surface defined by r on the st-
rectangle [si−1 , si ]×[tj−1 , tj ] with the parallelogram determined by the vectors
rs (si−1 , tj−1 )∆s and rt (si−1 , tj−1 )∆t, as seen in Figure 11.6.6.
rt ∆t
rs ∆s
where the latter equality holds from standard properties of the cross product
and length.
We sum the surface area approximations from Equation (11.6.1) over all
sub-rectangles to obtain an estimate for the total surface area, S, given by
m X
X n
S≈ |rs (si−1 , tj−1 ) × rt (si−1 , tj−1 )|∆s∆t.
i=1 j=1
Taking the limit as m, n → ∞ shows that the surface area of the surface
defined by r over the domain D is given as follows.
Surface area.
Let r(s, t) = hx(s, t), y(s, t), z(s, t)i be a parameterization of a smooth
surface over a domain D. The area of the surface defined by r on D is
CHAPTER 11. MULTIPLE INTEGRALS 254
given by ZZ
S= |rs × rt | dA. (11.6.2)
D
Activity 11.6.3 Consider the cylinder with radius a and height h defined
parametrically by
r(s, t) = a cos(s)i + a sin(s)j + tk
for 0 ≤ s ≤ 2π and 0 ≤ t ≤ h, as shown in Figure 11.6.7.
z
11.6.3 Summary
• A parameterization of a curve describes the coordinates of a point on
the curve in terms of a single parameter t, while a parameterization of a
surface describes the coordinates of points on the surface in terms of two
independent parameters.
• If r(s, t) = hx(s, t), y(s, t), z(s, t)i describes a smooth surface in 3-space
on a domain D, then the area, S, of that surface is given by
ZZ
S= |rs × rt | dA.
D
11.6.4 Exercises
1. Consider the cone shown below.
If the height of the cone is 8 and the base radius is 9, write a parame-
terization of the cone in terms of r = s and θ = t.
x(s, t) = ,
y(s, t) = ,
and
z(s, t) = ,
with
≤s≤ and
≤t≤ .
Answer 1. s cos(t)
Answer 2. s sin(t)
Answer 3. 8 − 89 s
Answer 4. 0
Answer 5. 9
Answer 6. 0
Answer 7. 2π
2. Parameterize the plane through the point (−5, −1, 3) with the normal
vector h−5, 4, 2i
~r(s, t) =
(Use s and t for the parameters in your parameterization, and enter
your vector as a single vector, with angle brackets: e.g., as \lt 1
+ s + t, s - t, 3 - t \gt.)
Answer. h− (5 + 4s + 2t) , − (1 + 5s) , 3 − 5ti
√
3. Parameterize a vase formed by rotating the curve z = 2 x − 1, 1 ≤ x ≤ 2,
around the z-axis. Use s and t for your parameters.
CHAPTER 11. MULTIPLE INTEGRALS 256
x(s, t) = ,
y(s, t) = ,
and
z(s, t) = ,
with
≤s≤ and
≤t≤
Answer 1. s cos(t)
Answer 2. s sin(t)
√
Answer 3. 2 s − 1
Answer 4. 1
Answer 5. 2
Answer 6. 0
Answer 7. 2π
4. Find parametric equations for the sphere centered at the origin and with
radius 4. Use the parameters s and t in your answer.
x(s, t) = ,
y(s, t) = ,
and
z(s, t) = ,
where
≤s≤ and
≤t≤ .
Answer 1. 4 cos(t) sin(s)
Answer 2. 4 sin(t) sin(s)
Answer 3. 4 cos(s)
Answer 4. 0
Answer 5. 2π
Answer 6. 0
Answer 7. π
5. Find the surface area of that part of the plane 5x + 10y + z = 3 that lies
2 2
inside the elliptic cylinder x16 + y81 = 1
Surface Area =
Answer. 1269.51444272221
6. Find the surface area of the part of the circular paraboloid z = x2 + y 2
that lies inside the cylinder x2 + y 2 = 9.
Answer. 117.318700709818
7. Find the surface area of the part of the plane 3x + 4y + z = 5 that lies
inside the cylinder x2 + y 2 = 16.
Answer. 256.304675910625
8. Write down the iterated integral which expresses the surface area of z =
y 6 cos4 x over the triangle with vertices (-1,1), (1,1), (0,2):
Z b Z g(y) p
h(x, y) dxdy
a f (y)
a=
b=
f (y) =
g(y) =
CHAPTER 11. MULTIPLE INTEGRALS 257
h(x, y) =
Answer 1. 1
Answer 2. 2
Answer 3. y − 2
Answer 4. 2 − y
Answer 5. 1 + y 12 · 16 cos6 (x) sin2 (x) + 36y 10 cos8 (x)
9. A decorative oak post is 60 inches long and is turned on a lathe so that
its profile is sinusoidal as shown in the figure below.
x2 y2 z2
+ + = 1.
16 25 9
In Activity 11.6.2, we found that a parameterization of the sphere S
of radius R centered at the origin is
for − π2 ≤ s ≤ π
2 and 0 ≤ t ≤ 2π.
a. Let (x, y, z) be a point on the ellipsoid and let X = x4 , Y = y5 , and
Z = z3 . Show that (X, Y, Z) lies on the sphere S. Hence, find a
parameterization of S in terms of X, Y , and Z as functions of s and
t.
b. Use the result of part (a) to find a parameterization of the ellipse in
terms of x, y, and z as functions of s and t. Check your parametriza-
tion by substituting x, y, and z into the equation of the ellipsoid.
Then check your work by plotting the surface defined by your pa-
rameterization.
11. In this exercise, we explore how to use a parametrization and iterated
integral to determine the surface area of a sphere.
a. Set up an iterated integral whose value is the portion of the surface
area of a sphere of radius R that lies in the first octant (see the
parameterization you developed in Activity 11.6.2).
b. Then, evaluate the integral to calculate the surface area of this por-
tion of the sphere.
d. Observe that the vector u = h2, 0, −4i points from (0, 0, 24) to
(2, 0, 20) along one side of the surface generated by the plane f
over D. Find the vector v such that u and v together span the
parallelogram that represents the surface defined by f over D, and
hence compute |u × v|. What do you observe about the value you
find?
13. A cone with basepradius a and height h can be realized as the surface
defined by z = ha x2 + y 2 , where a and h are positive.
p
a. Find a parameterization of the cone described by z = ha x2 + y 2 .
CHAPTER 11. MULTIPLE INTEGRALS 259
c. Evaluate the iterated integral to find a formula for the lateral surface
area of a cone of height h and base a.
y
x
Figure 11.7.1 A partitioned three-dimensional domain.
c. Let Bijk denote the sub-box [xi−1 , xi ] × [yj−1 , yj ] × [zk−1 , zk ]. Say that
∗
we choose a point (x∗ijk , yijk ∗
, zijk ) in the i, j, kth sub-box for each possible
combination of i, j, k. What is the meaning of δ(x∗ijk , yijk ∗ ∗
, zijk )? What
∗ ∗ ∗
physical quantity will δ(xijk , yijk , zijk )∆V approximate?
d. What final step(s) would it take to determine the exact mass of the piece
of granite?
• Let Bijk be the sub-box of B with opposite vertices (xi−1 , yj−1 , zk−1 )
and (xi , yj , zk ) for i between 1 and m, j between 1 and n, and k between
1 and `. The volume of each Bijk is ∆V = ∆x · ∆y · ∆z.
CHAPTER 11. MULTIPLE INTEGRALS 261
X n X
m X `
f (x∗ijk , yijk
∗ ∗
, zijk ) · ∆V.
i=1 j=1 k=1
♦
If f (x, y, z) represents the mass density of the box B, then, as we saw in
Preview Activity 11.7.1, the triple Riemann sum approximates the total mass
of the box B. In order to find the exact mass of the box, we need to let the
number of sub-boxes increase without bound (in other words, let m, n, and `
go to infinity); in this case, the finite sum of the mass approximations becomes
the actual mass of the solid B. More generally, we have the following definition
of the triple integral.
Definition 11.7.3 With following notation defined as in a triple Riemann
sum, the triple integral of f over B is
ZZZ m X
X n X
`
f (x, y, z) dV = lim f (x∗ijk , yijk
∗ ∗
, zijk ) · ∆V.
B m,n,`→∞
i=1 j=1 k=1
♦
As we noted earlier, if f (x, y, z) represents the density of the solid B at
each point (x, y, z), then
ZZZ
M= f (x, y, z) dV
B
is the mass of B. Even more importantly, for any continuous function f over
the solid B, we can use a triple integral to determine the average value of
f over B, fAVG(B) . We note this generalization of our work with functions
of two variables along with several others in the following important boxed
information. Note that each of these quantities may actually be considered
over a general domain S in R3 , not simply a box, B.
• The triple integral ZZZ
V (S) = 1 dV
S
represents the volume of the solid S.
• The average value of the function f = f (x, y, x) over a solid domain S is
given by ZZZ
1
fAVG(S) = f (x, y, z) dV,
V (S) S
4 y
x
-4 -2 2 4
-3
-2
0
-3 0 3
-4
Figure 11.7.4 Left: The cone. Right: Its projection.
iii. Finally, now that we have thought about slicing up the two-dimensional
domain that is the projection of the cone, what are the limits on z
in the innermost integral? Note that over any point (x, y) in the
plane, a vertical slice in the z direction will involve a range of values
from the cone itself to its flat top. In particular, observe that at
least one of these limits is not constant but depends on x and y.
iv. In conclusion, write an iterated integral of the form (11.7.1) that
represents the mass of the cone S.
Note well: When setting up iterated integrals, the limits on a given variable
can be only in terms of the remaining variables. In addition, there are multiple
different ways we can choose to set up such an integral. For example,
RRR two pos-
sibilities for iterated integrals that represent a triple integral S
f (x, y, z) dV
over a solid S are
R b R g (x) R h (x,y)
• a g12(x) h12(x,y) f (x, y, z) dz dy dx
R s R p2 (z) R q2 (x,z)
• r p1 (z) q1 (x,z)
f (x, y, z) dy dx dz
z 6 y
2
5
4
3
3
y 2
1
x
6
x 2 4 6
Figure 11.7.6 Left: The tetrahedron. Right: Its projection.
We find the mass, M , of the tetrahedron by the triple integral
ZZZ
M= δ(x, y, z) dV,
S
x + 2y + 3z = 6;
CHAPTER 11. MULTIPLE INTEGRALS 264
x + 2y = 6.
The right image in Figure 11.7.6 shows the projection of the tetrahedron
onto the xy-plane.
If we choose to integrate with respect to y for the middle integral in the
iterated integral, then the lower limit on y is the x-axis and the upper limit
is the hypotenuse of the triangle. Note that the hypotenuse joins the points
(6, 0) and (0, 3) and so has equation y = 3 − 12 x. Thus, the bounds on y are
0 ≤ y ≤ 3 − 12 x. Finally, the x values run from 0 to 6, so the iterated integral
that gives the mass of the tetrahedron is
Z 6 Z 3−(1/2)x Z (1/3)(6−x−2y)
M= x + y + z dz dy dx. (11.7.2)
0 0 0
-1
0 x
-1 0 y 1
a. First, set up an iterated double integral to find the volume of the solid S
as a double integral of a solid under a surface. Then set up an iterated
triple integral that gives the volume of the solid S. You do not need to
evaluate either integral. Compare the two approaches.
b. Set up (but do not evaluate) iterated integral expressions that will tell us
the center of mass of S, if the density at point (x, y, z) is δ(x, y, z) = x2 +1.
c. Set up (but do not evaluate) an iterated integral to find the average
density on S using the density function from part (b).
d. Use technology appropriately to evaluate the iterated integrals you de-
termined in (a), (b), and (c); does the location you determined for the
center of mass make sense?
11.7.2 Summary
• Let f = f (x, y, z) be a continuous function on a box B = [a, b] × [c, d] ×
[r, s]. The triple integral of f over B is defined as
ZZZ m X
X n X
l
f (x, y, z) dV = lim f (x∗ijk , yijk
∗ ∗
, zijk ) · ∆V,
B ∆V →0
i=1 j=1 k=1
CHAPTER 11. MULTIPLE INTEGRALS 266
where the triple Riemann sum is defined in the usual way. The definition
of the triple integral naturally extends to non-rectangular solid regions
S.
RRR
• The triple integral S
f (x, y, z) dV can tell us
◦ the volume of the solid S if f (x, y, z) = 1,
◦ the mass of the solid S if f represents the density of S at the point
(x, y, z).
Moreover, ZZZ
1
fAVG(S) = f (x, y, z) dV,
V (S) S
is the average value of f over S.
11.7.3 Exercises
1. Find the triple integral of the function f (x, y, z) = x2 cos(y + z) over the
cube 6 ≤ x ≤ 8, 0 ≤ y ≤ π, 0 ≤ z ≤ π.
Answer. −394.667
2. Evaluate the triple integral
ZZZ
xyz dV
E
Use these definitions to find the moment of inertia about the z-axis of
the rectangular solid of mass 27 given by 0 ≤ x ≤ 3, 0 ≤ y ≤ 3, 0 ≤ z ≤ 3.
Ix =
Iy =
Iz =
CHAPTER 11. MULTIPLE INTEGRALS 267
1·3·3·3·(3·3+3·3)
Answer 1. 3
1·3·3·3·(3·3+3·3)
Answer 2. 3
1·3·3·3·(3·3+3·3)
Answer 3. 3
ZZZ
8. Express the integral f (x, y, z)dV as an iterated integral in six dif-
E
ferent ways, where E is the solid bounded by z = 0, x = 0, z = y − 6x and
y = 12.
Z b Z g2 (x) Z h2 (x,y)
1. f (x, y, z)dzdydx
a g1 (x) h1 (x,y)
a= b=
g1 (x) = g2 (x) =
h1 (x, y) = h2 (x, y) =
Z b Z g2 (y) Z h2 (x,y)
2. f (x, y, z)dzdxdy
a g1 (y) h1 (x,y)
a= b=
g1 (y) = g2 (y) =
h1 (x, y) = h2 (x, y) =
Z b Z g2 (z) Z h2 (y,z)
3. f (x, y, z)dxdydz
a g1 (z) h1 (y,z)
a= b=
g1 (z) = g2 (z) =
h1 (y, z) = h2 (y, z) =
Z b Z g2 (y) Z h2 (y,z)
4. f (x, y, z)dxdzdy
a g1 (y) h1 (y,z)
a= b=
g1 (y) = g2 (y) =
h1 (y, z) = h2 (y, z) =
Z b Z g2 (x) Z h2 (x,z)
5. f (x, y, z)dydzdx
a g1 (x) h1 (x,z)
a= b=
g1 (x) = g2 (x) =
h1 (x, z) = h2 (x, z) =
Z b Z g2 (z) Z h2 (x,z)
6. f (x, y, z)dydxdz
a g1 (z) h1 (x,z)
a= b=
g1 (z) = g2 (z) =
h1 (x, z) = h2 (x, z) =
Answer 1. 0
Answer 2. 2
Answer 3. 6x
Answer 4. 12
Answer 5. 0
Answer 6. y − 6x
Answer 7. 0
Answer 8. 12
Answer 9. 0
y
Answer 10. 6
Answer 11. 0
CHAPTER 11. MULTIPLE INTEGRALS 268
Answer 12. y − 6x
Answer 13. 0
Answer 14. 12
Answer 15. z
Answer 16. 12
Answer 17. 0
y−z
Answer 18. 6
Answer 19. 0
Answer 20. 12
Answer 21. 0
Answer 22. y
Answer 23. 0
y−z
Answer 24. 6
Answer 25. 0
Answer 26. 2
Answer 27. 0
Answer 28. 12 − 6x
Answer 29. 6x + z
Answer 30. 12
Answer 31. 0
Answer 32. 12
Answer 33. 0
z
Answer 34. 2 − 6
Answer 35. 6x + z
Answer 36. 12
9. Calculate the volume under the elliptic paraboloid z = 4x2 + 8y 2 and over
the rectangle R = [−2, 2] × [−3, 3].
Answer. 704
10. The motion of a solid object can be analyzed by thinking of the mass
as concentrated at a single point, the center of mass. If the object has
density ρ(x, y, z) at the point (x, y, z) and occupies a region W , then the
coordinates (x, y, z) of the center of mass are given by
Z Z Z
1 1 1
x= xρ dV y = yρ dV z = zρ dV,
m W m W m W
z = 2.Z Z Z
p
(a) z − x2 + y 2 dV
W
ZZZ
(b) (z − 2) dV
W
ZZZ
(c) e−xyz dV
W
12. Set up a triple integral to find the mass of the solid tetrahedron bounded
by the xy-plane, the yz-plane, the xz-plane, and the plane x/3+y/2+z/6 =
1, if the density function is given by δ(x, y, z) = x + y. Write an iterated
integral in the form below to find the mass of the solid.
ZZZ Z BZ DZ F
f (x, y, z) dV =
A C E
R
dz dy dx
with limits of integration
A=
B=
C=
D=
E=
F=
Answer 1. x + y
Answer 2. 0
Answer 3. 3
Answer 4. 0
Answer 5. 2 − 0.666667x
Answer 6. 0
Answer 7. 6 − 2x − 3y
13. Consider the solid S that is bounded by the parabolic cylinder y = x2 and
the planes z = 0 and z = 1 − y as shown in Figure 11.7.8.
z
1
y
1
x 1
y z z
1.0 1.0 1.0
0.5 0.5
0.5
x x
y
-1 1 -1 1
0.5 1.0
z z
(ρ, θ, φ)
(r, θ, z)
z φ ρ
y y
θ r θ
x x
b. Let ∆A be the area of the projection of the box, B, onto the xy-plane,
which is shaded blue in Figure 11.8.2. Recall that we previously deter-
mined the area ∆A in polar coordinates in terms of r, ∆r, and ∆θ. In
light of the fact that we know ∆A and that z is the standard z coordi-
nate from Cartesian coordinates, what is the volume ∆V in cylindrical
coordinates?
Activity 11.8.3 demonstrates that the volume element dV in cylindrical
coordinates is given by dV = r dz dr dθ, and hence the following rule holds in
general.
1 y
1.0
z 0.5
x
-1.0 -0.5 0.5 1.0
-0.5
-1
-1.0
0
-1 0 1
Figure 11.8.3 The cylindrical cone r = 1 − z and its projection onto the
xy-plane.
-2
0
-2 0 2
p
p 11.8.4 A solid bounded by the cones z =
Figure x2 + y 2 and z =
4 − x2 + y 2 .
z
ρ sin φ
(ρ, θ, φ)
φ ρ ρ cos φ
y
θ
ρ sin φ
π
with θ = 6, 0 ≤ ρ ≤ 1, and 0 ≤ φ ≤ π?
R
R
Q
Q P
S
P
S
y
y
x
x
Figure 11.8.6 Left: A spherical box. Right: A spherical volume element.
d. Use your work in (a), (b), and (c) to determine an approximation for ∆V
in spherical coordinates.
CHAPTER 11. MULTIPLE INTEGRALS 279
Letting ∆ρ, ∆θ, and ∆φ go to 0, it follows from the final result in Ac-
tivity 11.8.6 that dV = ρ2 sin(φ) dρ dθ dφ in spherical coordinates, and thus
allows us to state the following general rule.
in spherical coordinates.
The latter expression is an iterated integral in spherical coordinates.
π
Figure 11.8.7 The solid cut from the sphere ρ = 2 by the cone φ = 4.
11.8.5 Summary
• The cylindrical coordinates of a point P are (r, θ, z) where r is the dis-
tance from the origin to the projection of P onto the xy-plane, θ is the
angle that the projection of P onto the xy-plane makes with the positive
x-axis, and z is the vertical distance from P to the projection of P onto
the xy-plane. When P has rectangular coordinates (x, y, z), it follows
that its cylindrical coordinates are given by
y
r 2 = x2 + y 2 , tan(θ) = , z = z.
x
CHAPTER 11. MULTIPLE INTEGRALS 280
When P has given cylindrical coordinates (r, θ, z), its rectangular coor-
dinates are
x = r cos(θ), y = r sin(θ), z = z.
Given the point P in spherical coordinates (ρ, φ, θ), its rectangular coor-
dinates are
11.8.6 Exercises
1. What are the rectangular coordinates of the point whose cylindrical co-
ordinates are
(r = 4, θ = 2π
7 , z = −1) ?
x=
y=
z=
Answer 1. 2.49395920743493
Answer 2. 3.12732592987212
Answer 3. −1
2. What are the rectangular coordinates of the point whose spherical coor-
dinates are
2, − 23 π, 0π ?
x=
y=
z=
Answer 1. 2 cos −2
3 π sin(0π)
Answer 2. 2 sin −2
3 π sin(0π)
Answer 3. 2 cos(0π)
3. What are the cylindrical coordinates of the point whose spherical coordi-
nates are
CHAPTER 11. MULTIPLE INTEGRALS 281
5π
(3, 3, 6 ) ?
r=
θ=
z=
Answer 1. 1.5
Answer 2. 3
Answer 3. −2.59807621135332
4. Find an equation for the paraboloid z = x2 + y 2 in spherical coordinates.
(Enter rho, phi and theta for ρ, φ and θ, respectively.)
equation:
cos(phi)
Answer. rho = sin2(phi)
5. Match the given equation with the verbal description of the surface:
A. Plane
B. Half plane
C. Circular Cylinder
(a) r = 4
(b) z = r2
(c) ρ = 2 cos(φ)
π
(d) φ = 3
(e) r = 2 cos(θ)
(f) r2 + z 2 = 16
(g) ρ cos(φ) = 4
π
(h) θ = 3
(i) ρ = 4
6. Match the integrals with the type of coordinates which make them the
easiest to do. Put the letter of the coordinate system to the left of the
number of the integral.
ZZZ
(a) z dV where E is: 1 ≤ x ≤ 2, 3 ≤ y ≤ 4, 5 ≤ z ≤ 6
E
Z 1Z y2
1
(b) dx dy
0 0 x
ZZZ
(c) z 2 dV where E is: −2 ≤ z ≤ 2, 1 ≤ x2 + y 2 ≤ 2
E
ZZ
1
(d) dA where D is: x2 + y 2 ≤ 4
D x2 + y 2
ZZZ
(e) dV where E is: x2 + y 2 + z 2 ≤ 4, x ≥ 0, y ≥ 0, z ≥ 0
E
CHAPTER 11. MULTIPLE INTEGRALS 282
A. polar coordinates
B. cylindrical coordinates
C. spherical coordinates
D. cartesian coordinates
7. Evaluate √
the integral.
Z 7 Z 49−x2 Z √49−x2 −z2
1
√ √ dy dz dx =
0 − 49−x2 − 49−x2 −z 2 (x2 + y2+ z 2 )1/2
Answer. 49π
ZZZ p
8. Use cylindrical coordinates to evaluate the triple integral x2 + y 2 dV ,
E
where E is the solid bounded by the circular paraboloid z = 1−1 x2 + y 2
and the xy -plane.
Answer. 0.837758040957278
ZZZ
9. Use spherical coordinates to evaluate the triple integral x2 + y 2 +
E
z 2 dV , where E is the ball: x2 + y 2 + z 2 ≤ 4.
Answer. 80.4247719318987
10. Find the volume of the solid enclosed by the paraboloids z = 1 x2 + y 2
and z = 8 − 1 x2 + y 2 .
Answer. 50.2654824574367
11. FInd the volume of the ellipsoid x2 + y 2 + 6z 2 = 100.
Answer. 1710.06644021582
12. The density, δ, of the cylinder x2 + y 2 ≤ 9, 0 ≤ z ≤ 4 varies with the
distance, r, from the z-axis:
δ = 4 + r g/cm3 .
π
Answer 4. 2
Answer 5. π
Answer 6. 0
Answer 7. 4
Answer 8. π · 42
14. In each of the following questions, set up an iterated integral expression
whose value determines the desired result. Then, evaluate the integral
first by hand, and then using appropriate technology.
a. Find the volume of the “cap” cut from the solid sphere x2 +y 2 +z 2 =
4 by the plane z = 1, as well as the z-coordinate of its centroid.
c. Find the volume of the solid bounded below by the xy-plane, on the
sides by the sphere ρ = 2, and above by the cone φ = π/3.
d. Find the z coordinate of the center qpof mass of the region that is
bounded above by the surface z = x2 + y 2 , on the sides by the
2 2
cylinder x + y = 4, and below by the xy-plane. Assume that the
density of the solid is uniform and constant.
e. Find the volume of the solid that lies outside the sphere x2 +y 2 +z 2 =
1 and inside the sphere x2 + y 2 + z 2 = 2z.
15. For each of the following questions,
• sketch the region of integration,
• change the coordinate system in which the iterated integral is writ-
ten to one of the remaining two,
R π/2 R π R 1
b. 0 0 0
ρ2 sin(φ) dρ dφ dθ
R 2π R 1 R 1
c. 0 0 r
r2 cos(θ) dz dr dθ
16. Consider the solid region S bounded above by the paraboloid z = 16 −
x2 − y 2 and below by the paraboloid z = 3x2 + 3y 2 .
a. Describe parametrically the curve in R3 in which these two surfaces
intersect.
b. In terms of x and y, write an equation to describe the projection of
the curve onto the xy-plane.
c. What coordinate system do you think is most natural for an iterated
integral that gives the volume of the solid?
d. Set up, but do not evaluate, an iterated integral expression whose
value is average z-value of points in the solid region S.
e. Use technology to plot the two surfaces and evaluate the integral in
CHAPTER 11. MULTIPLE INTEGRALS 284
The vertices of the polar rectangle P are transformed into the vertices of
a closed and bounded region P 0 in rectangular coordinates. If we view the
standard coordinate system as having the horizontal axis represent r and the
vertical axis represent θ, then the polar rectangle P appears to us at left in
Figure 11.9.1. The image P 0 of the polar rectangle P under the transformation
given by (11.9.2) is shown at right in Figure 11.9.1. We thus see that there is a
correspondence between a simple region (a traditional, right-angled rectangle)
and a more complicated region (a fraction of an annulus) under the function
T given by T (r, θ) = (r cos(θ), r sin(θ)).
θ y
2 2
P′
1 P 1
r x
1 2 1 2
Let’s see what happens to the rectangle T = [0, 1] × [1, 4] in the st-plane
under this change of variable.
a. Draw a labeled picture of T in the st-plane.
b. Find the image of the st-vertex (0, 1) in the xy-plane. Likewise, find the
respective images of the other three vertices of the rectangle T : (0, 4),
(1, 1), and (1, 4).
are given by the vectors v and w, where v is the vector from (x(a, b), y(a, b))
to (x(a + ∆s, b), y(a + ∆s, b)), and w is the vector from (x(a, b), y(a, b)) to
(x(a, b + ∆t), y(a, b + ∆t)).
An example of an image T 0 in the xy-plane that results from a transforma-
tion of a rectangle T in the st-plane is shown in Figure 11.9.2.
t y
2 2
T′
1 T 1
w v
s x
1 2 1 2
For small ∆s and ∆t, the definition of the partial derivative tells us that
∂x ∂y ∂x ∂y
v≈ (a, b), (a, b), 0 ∆s and w≈ (a, b), (a, b), 0 ∆t.
∂s ∂s ∂t ∂t
Recall that the area of the parallelogram with sides v and w is the length
of the cross product of the two vectors, |v × w|. From this, we observe that
∂x ∂y ∂x ∂y
v×w≈ (a, b), (a, b), 0 ∆s × (a, b), (a, b), 0 ∆t
∂s ∂s ∂t ∂t
∂x ∂y ∂x ∂y
= 0, 0, (a, b) (a, b) − (a, b) (a, b) ∆s ∆t.
∂s ∂t ∂t ∂s
Finally, by computing the magnitude of the cross product, we see that
∂x ∂y ∂x ∂y
|v × w| ≈ 0, 0, (a, b) (a, b) − (a, b) (a, b) ∆s ∆t
∂s ∂t ∂t ∂s
∂x ∂y ∂x ∂y
= (a, b) (a, b) − (a, b) (a, b) ∆s ∆t.
∂s ∂t ∂t ∂s
CHAPTER 11. MULTIPLE INTEGRALS 288
Activity 11.9.3 Find the Jacobian when changing from rectangular to polar
coordinates. That is, for the transformation given by x = r cos(θ), y = r sin(θ),
determine a simplified expression for the quantity
∂x ∂y ∂x ∂y
− .
∂r ∂θ ∂θ ∂r
What do you observe about your result? How is this connected to our
earlier work with double integrals in polar coordinates?
Activity 11.9.4 Let D0 be the region in the xy-plane bounded by the lines
y = 0, x = 0, and x + y = 1. We will evaluate the double integral
√
ZZ
x + y(x − y)2 dA (11.9.4)
D0
CHAPTER 11. MULTIPLE INTEGRALS 289
∂x ∂x ∂x
∂s ∂t ∂u
∂y ∂y ∂y .
∂s ∂t ∂u
∂z ∂z ∂z
∂s ∂t ∂u
(Recall that this determinant was introduced in Section 9.4.) That is, ∂(x,y,z)
∂(s,t,u)
is given by
∂x ∂y ∂z ∂y ∂z ∂x ∂y ∂z ∂y ∂z ∂x ∂y ∂z ∂y ∂z
− − − + − .
∂s ∂t ∂u ∂u ∂t ∂t ∂s ∂u ∂u ∂s ∂u ∂s ∂t ∂t ∂s
To summarize,
conditions
RRR (that are studied in advanced calculus), the triple integral
S0
f (x, y, z) dV is equal to
ZZZ
∂(x, y, z)
f (x(s, t, u), y(s, t, u), z(s, t, u)) ds dt du.
S ∂(s, t, u)
Activity 11.9.5 Find the Jacobian when changing from Cartesian to cylin-
drical coordinates. That is, for the transformation given by x = r cos(θ),
y = r sin(θ), and z = z, determine a simplified expression for the quantity
∂(x, y, z)
.
∂(r, θ, z)
What do you observe about your result? How is this connected to our
earlier work with triple integrals in cylindrial coordinates?
Activity 11.9.6 Consider the solid S 0 defined by the inequalities 0 ≤ x ≤ 2,
x x x
2 ≤ y ≤ 2 + 1, and 0 ≤ z ≤ 6. Consider the transformation defined by s = 2 ,
x−2y z
t = 2 , and u = 3 . Let f (x, y, x) = x − 2y + z.
a. The transformation turns the solid S 0 in xyz-coordinates into a box S in
stu-coordinates. Apply the transformation to the boundries of the solid
S 0 to find stu-coordinatte descriptions of the box S.
∂(x,y,z)
b. Find the Jacobian ∂(s,t,u) .
c. RRR
Use the transformation to perform a change of variables and evaluate
S0
f (x, y, z) dV by evaluating
ZZZ
∂(x, y, z)
f (x(s, t, u), y(s, t, u), z(s, t, u)) ds dt du.
S ∂(s, t, u)
11.9.4 Summary
• If an integral is described in terms of one set of variables, we may write
that set of variables in terms of another set of the same number of vari-
ables. If the new variables are chosen appropriately, the transformed
integral may be easier to evaluate.
• The Jacobian is a scalar function that relates the area or volume element
in one coordinate system to the corresponding element in a new system
determined by a change of variables.
11.9.5 Exercises
∂(x,y)
1. Find the absolute value of the Jacobian, ∂(s,t) , for the change of variables
given by x = 6s + 8t, y = 4s + 8t
∂(x,y)
∂(s,t) =
Answer. |6 · 8 − 8 · 4|
∂(x,y,z)
2. Find the Jacobian. ∂(s,t,u) , where x = 3s + t + 4u, y = 2s − t + u, z =
3s + t − 5u.
∂(x,y,z)
∂(s,t,u) =
Answer. 3 · (−1 · (−5) − 1 · 1) − 1 · (2 · (−5) − 1 · 3) + 4 · (2 · 1 − −1 · 3)
CHAPTER 11. MULTIPLE INTEGRALS 291
Answer 1. 1
Answer 2. 31
Answer 3. 49
Answer 4. −49
Answer 5. 31
Answer 6. −31
Answer 7. −49
Answer 8. 49
Answer 9. −31
Answer 10. 7535362.66666667
Use the change of variables s = y, t = y − x2 to evaluate
RR
4. R
x dx dy over
the region R in the first quadrant bounded by y = 0, y = 36, y = x2 , and
y =RxR2 − 1.
R
x dx dy =
1·36
Answer. 2
5. Use the change of variables s = x + y, t = y to find the area of the ellipse
x2 + 2xy + 2y 2 ≤ 1.
area =
Answer. π
Use the change of variables s = xy, t = xy 2 to compute R xy 2 dA, where
R
6.
R isR the region bounded by xy = 1, xy = 2, xy 2 = 1, xy 2 = 2.
R
xy 2 dA =
Answer. 1 · 1
7. Find positive numbers a andR Rb so that the change of variables s = ax, t =
by transforms the integral R
dx dy into
Z Z
∂(x, y)
ds dt
T ∂(s, t)
10. Let D0 be the region in the xy-plane that is the parallelogram with vertices
(3, 3), (4, 5), (5, 4), and (6, 6).
a. Sketch and label the region D0 in the xy-plane.
RR
b. Consider the integral D0 (x + y) dA. Explain why this integral
would be difficult to set up as an iterated integral.
c. Let a change of variables be given by x = 2u + v, y = u + 2v. Using
substitution or elimination, solve this system of equations for u and
v in terms of x and y.
d. Use your work in (c) to find the pre-image, D, which lies in the uv-
plane, of the originally given region D0 , which lies in the xy-plane.
For instance, what uv point corresponds to (3, 3) in the xy-plane?
e. Use the change of variables in (c) and your other work to write a
RR integral in u and v that is equivalent to the original xy
new iterated
integral D0 (x + y) dA.
f. Finally, evaluate the uv integral, and write a sentence to explain
why the change of variables made the integration easier.
CHAPTER 11. MULTIPLE INTEGRALS 293
e. Explain why this new integral is better, but is still difficult to eval-
uate. What additional change of variables would make the resulting
integral easier to evaluate?
f. Convert the integral from (d) to a new integral in spherical coordi-
nates.
g. Finally, evaluate the iterated integral in (f) and hence determine the
volume of the ellipsoid.
Index
294
INDEX 295