Multi Variable Calculus Notes
Multi Variable Calculus Notes
Lecture Notes 1
Y. L. Wong
Department of Mathematics, National University of Singapore, Singapore 117543
...................................................................................................................................
................................. .................
................. ............
............ ..........
................ ........
...
....... .. ... ... .
. . ..
. ... . .. . .
. . .
.... ... .. . ... . .. ..
. ... .
.. ... ... . .. ... ... ... .....
. ... ... ... . .. ... ... ... ... ... ... ...
.......
......
....... ... ... ... ... ... . ... ... ... ... . .....
...
. ... ... ... .. ...
..
... ... ...
.... ....... . .
........................................................................................................................................................................... ..... ....
...... ..... ................ . ... ......
........... ..................... ..............
....... ... ..
..... ..........
....... .................
.... ..................... ..... .
.... .............. .......... ...... ......
......... .......... . .........
...........
............ ........ ............ .........
................. ....... ................ .... ...
.... ...... ................................
.................................... ............................ ..... ....
... ... ... . .................................... ... ...
... ... ... . ................ ......... ... .
..... .. ... ...
...
. . ..... ...... ...
. ... ... .
.. ... ...
...... ... ... ... ... .
.. ... ... ... ... ... ... ............ ... ... ... ... .. ......... ... ... ... ... ... ........ ... ... ... ... ... ... ... ... ..... .....
......
.......
........
..........
• .
...
.
.....
....
.
..
.... ....
..... ..................
............. ..... ..... ...... ...........
..................
......................................................................... ......... ..............
....................................................................
....
............
.
...
1
This notes is written exclusively for students taking the modules MA1104, Multi-variable Calculus, at the
National University of Singapore. The contents follow closely the reference [1].
Contents
1. Notations 4
2. Vectors in R3 5
3. Cylinders and Quadric Surfaces 17
4. Cylindrical and Spherical Coordinates 20
5. Vector Functions 23
6. Functions of several variables 27
7. Limits and Continuity 30
8. Partial Derivatives 34
9. Maximum and Minimum Values 45
10. Lagrange Multipliers 48
11. Multiple Integrals 52
12. Surface Area 60
13. Triple Integrals 61
14. Change of Variables in Multiple Integrals 66
15. Vector Fields 70
16. Line Integrals 72
17. Line Integrals of Vector Fields 75
18. The Fundamental Theorem for Line Integrals 76
19. Independence of Path 77
20. Green’s Theorem 81
21. The Curl and Divergence of a Vector Field 86
22. Parametric Surfaces and their Areas 90
23. Oriented Surfaces 95
24. Surface Integrals of Vector Fields 97
25. Stokes’ Theorem 99
26. The Divergence Theorem 102
Bibliography 105
3
Semester II(2005/06)
1. Notations
The collection of all real numbers is denoted by R. Thus R includes the integers
. . . , −2, −1, 0, 1, 2, 3 . . . ,
√ rational numbers, p/q, where p and q are integers (q 6= 0), and the irrational numbers, like
the
2, π, e etc. members of R may be visualized as points on the real-number line as shown in
Figure 1.
.............................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................................
√ e 3π
−3 − 52 −2 −1 0 1 1 2 2 4
2
Figure 1 The Number Line
if a ≥ 0
a
|a| =
−a if a < 0.
A function f : A −→ B is a rule that assigns to each a ∈ A one specific member f (a) of B. The
fact that the function f sends a to f (a) is denoted symbolically by a 7→ f (a). For example,
f (x) = x2 /(1 − x) assigns the number x2 /(1 − x) to each x 6= 1 in R. We can specify a function
f by giving the rule for f (x). The set A is called the domain of f and B is the codomain of f .
The range of f is the subset of B consisting of all the values of f . That is, the range of f =
{f (x) ∈ B | x ∈ A}.
Given f : A −→ R. It means that f assigns a value f (x) in R to each x ∈ A. Such a function
is called a real-valued function. For a real-valued function f : A −→ R defined on a subset A
of R, the graph of f consists of all the points (x, f (x)) in the xy-plane.
c Department of Mathematics 4
Semester II(2005/06)
y
....
..........
•
..
..
.. ....
.... • ..
.
..
.. .. ..
... .. ..
... ...
... ..
... ...
graph of f .
.
..
..
..... ... ...
...... ..... ..
.
.
... .. ..
... ..
......... ... ...
... .... .....
... ... ........
... .. ...
... ...
... ....
... ...
(x, f (x)) ·.
..
..
..
..
.
...
....
A = domain of f
2. Vectors in R3
2.1. The Euclidean 3-space. The Euclidean 3-space denoted by R3 is the set
{(x, y, z) | x, y, z ∈ R}.
To specify the location of a point in R3 geometrically, we use a right-handed rectangular
coordinate system, in which three mutually perpendicular coordinate axes meet at the origin.
It is common to use the x and y axes to represent the horizontal coordinate plane and the
z-axis for the vertical height.
z..
..
..........
...
....
..
z• ....
......
... ..
... ...
... ..
... ..
... ..
... ..
... ..
... ..
...
...
...
•(x, y, z) ..
.
..
..
... ..
... ..
... ..
... ..
... ..
... ..
... ..
... ..
.
................... ... .
.....
.
...
.
..
... ..
.
..
................
y
.. .......................
.
..
.....
.....
..
..
.
..
..
..
• ................
...
. ....
... ................................
................ .
.
.
...
.
.
.
...
.
.
..... .
..
..
.
.
.
. .. ... ..
. ... ...................
.. y
x• .
.....
.
..
.
.
..
.
.......
.
. . ... ... ... .
.. ... ... ..
.
.
.. .. .....
.. . ...
.
. . .
.
....
.
.
.
...
.
.
..
(x, y)
.....
..........
.....
x
Example 2.1. (a) z = 3 is the equation of a horizontal plane at level 3 above the xy-plane.
(b) y = 2 is the equation of a vertical plane parallel to the xz-coordinate plane. Every point of
this plane has y coordinate equal to 2. (c) Similarly x = 2 is the equation of a vertical plane
parallel to the yz-coordinate plane.
z.. z.. z..
... ... ...
....... ........ ........
. .... z=3 .... ....
.
. . .
...
. . ................... ... .. ...
.....
.. .... . ........
..... .. ....................... ... ..... .. ...
..
. ... ................ ... ..... .... ...
..... ...
................
.......... ... .
..... ... ...
.
........ ... .... ... ..
.....
.. ... ...
... ..... .....
... .......
.
....................
•3 .
.
.
.
.
.. ...
.
.
.......
.
...
...
. .
...
.
......
.
...
...
...
...
...
. ...
.. ... ... ... ..................
................ ....
.
................
................ ...
. ..... ... ..
..... ... ... ...................
.................. .... ... ...
. ... ................
. ..... .
....
.. ... .. ................
.... ................................. ... .... .
.. ... ...
... ... ...........
.... ... ... ... ... ... ...
.. ... ... ...
... ... .
.
. ...
... ... ... ... ... .
.
.
...
. .
. . ...
.
................................. .... .............................
. . . ... ...
..
. .
... ........................... .
.. ................ .. . ..................... ... .. ................ ....
.. .
.....
.... ................
................ . .. .
...
.
....
. .
.
.. .
. .•
................ ....
...
. ...
... .
.
.
........
. ..........................
....................
........
.
.....
. ....................
.. y .
.. .
...
.
.. .
.... .
....
.
. 2 .....
. .
.
.
.
.
. ...........................
y •2
...
... ....... ...
.
.
.... ...
..
....
........... y
..
.
...
.
......
.
....
.
.
.
.... y=2 ...
.
.
.... ....
.
.......
...........
... ....... .
.
....
.....
... ..... ..... ....... x=2
..... ..... ... ..... ..... ..
....
..... ..... ..... ..... ... ...
.......... .
............ ...
. ..
......
. .......... .... ...
. ....
....... ....... . ....... ...
x x ....
... ........
..... x ................
................
................ .....
... ....... ................
................ ..
....... ................
z.. ...........................................
..
.......... .......... ...... ... ..........
....... ....
... ...... ..... .. ......
.... ..
...... .
. .. .....
....
.. ....
. .
.. .. ...
... .... ..
.
.... ... . .. ... ... .. . . . .
..
... ... ....
...
..
... ..
. .. . . ..
. . .
.. .. .... .....
. .
... .. ...... ... ... .
r •P
. .. ... .... ...
... ........ ..
. .......... ......
... .... ... ........ .....
... . . ........ ..
...
...
.....
......
.......
.....
...
...
.........
O• .
.
.
.
.
....
.....
.
..
...
... . ..
. . . ..
..
... ....... . .. ...
... ... ........... ..... ...... ...
... ............. .. .......... ...
...
... ... .................................................. .
...
...
... ... .. ....
.
..... ... .. .
... ..... ... .. .....
... ...... .. . .........
....... ..... . ..
... ......... ..... ... .......
... ..............................................
...
...
...
...
...
.. .......
..
.... ................................
. ..
...
. ................
..
. ................
..
....
. ................
................
..... ................
.....
.
...
.
.
.
....
.
.. ..................
....... y
..
..
..
....
.
....
.....
.....
.....
..
......
.
..
.....
....
.....
..........
......
x Figure 5 A sphere
Solution. Using the method of completing square, the given equation can be written as (x +
2)2 + (y√− 3)2 + (z + 1)2 = 8. Hence, it is the equation of a sphere centred at (−2, 3, −1) with
radius 8.
c Department of Mathematics 6
Semester II(2005/06)
....
.....
.....
(x + 2)2 + (y − 3)2 = 4, z = 1
.....
............................................
.......... ......
.....
.....
...
is a circle
........ ........ ... .......... ....
........ . .. ...... ........
... . . .... . .
.
..... .. ... ....... ... ... ... ... ... ... ......... ... ... ..................
.... ... ... . .... ................
....... ..
..
. ......
..... .. ... ... ... ... ... ... ... ... ... ... ......... ....
..... ..
.. .. . ... .... ... ....
... ......... ..... ... ... . ... ............ .....
.... ................................... ............... ... ..
... .... ...
... ... ... ..................................................... ... ..
.... .. . .....
..... ... .. ...
.....
....
......
.....
...
..
O• .
.
..
.
.
.
...
..
.......
.
......... ... . ...
.. ..
.
... ....... .. . .
... .......... .... .. ..... .
...... ....
... ............
.................. ................. .
...
... ... .............................................. ...
..
...
.....
....
..
.
.
. ...
.
... .
....
.....
..... ... .. .....
...... . .. .....
....... ....... .. .....
......... ..... ... .......
.........................................................
..
......
.
.
.....
.
. ..................
. y
..
....
.
..
....
.
..
.....
.....
....
..........
.....
x
Figure 7 Vector
Instead of using an arrow on top of P Q or a, we shall suppress the arrow but write PQ or a
→ −→
in bold to denote the vector a =P Q. If P is the origin O, a is called the position vector of the
point Q.
The position vectors of (1, 0, 0), (0, 1, 0 and (0, 0, 1) are denoted by i, j and k respectively. In
other word, i = h1, 0, 0i, j = h0, 1, 0i and k = h0, 0, 1i. i, j and k are called the standard basis
vectors. Therefore, if Q = (x, y, z), then q = hx, y, zi = xi + yj + zk.
Suppose p P = (x1 , y1 , z1 ) and Q = (x2 , y2 , z2 ). The magnitude of a vector PQ is defined to be
|PQ| = (x2 − x1 )2 + (y2 − y1 )2 + (z2 − z1 )2 . A vector PQ also has a direction determined
by the orientation that arrow is pointing.
The zero vector h0, 0, 0i is denoted by 0. Clearly |0| = 0. We say that two vectors are equal
if and only if they have the same direction and the same magnitude. This condition may be
expressed algebraically by saying that if v1 = hx1 , y1 , z1 i and v2 = hx2 , y2 , z2 i, then v1 = v2 if
and only if x1 = x2 , y1 = y2 and z1 = z2 .
If v1 = hx1 , y1 , z1 i and v2 = hx2 , y2 , z2 i, then define the sum v1 + v2 to be the vector hx1 +
x2 , y1 + y2 , z1 + z2 i.
If λ is any real number and v = hx, y, zi, then define the scalar multiple λv to be the vector
hλx, λy, λzi.
c Department of Mathematics 7
Semester II(2005/06)
... ...
... ..............
... ... ... .......
v....2.......................... ............................
... ..
. ... ... ......
....... ......
........ ......
.
..... .. ........
.
.. .
... ......
... ......
... ......
......
v1 ..
.
. ...
.
......
....
....... v1 + v2
... ......
... ......
... ......
..... ..
..
.......
..
... ...... .............
...... ............. ..
... ...... .............
... ....... .............
.
.... ........... ..... ..
.................
. ...
.............................
.................
..... .
v2
O
.
....................
............. ........ A
a-b .
...........
...
...
..
............. ...........
.............
. .
.......
..
....
B .........
.............
.............
......
......
...
......
............ ..
.......
.
.
.. ..
... ......
... ......
... ......
......
.... ..........
b ...
...
...
......
.
......
......
.
a
... ......
..... . .. . . ......
.
..
... ......
......
...
... θ ......
......
.
......................
. ..
... .......
.........
.......
O
Figure 9 Angle between two vectors
Applying cosine rule to 4OAB, we have
|a − b|2 = |a|2 + |b2 | − 2|a||b| cos θ.
As |a − b|2 = (a − b) · (a − b) = |a|2 − 2a · b + |b2 |, it follows that a · b = |a||b| cos θ or
a·b
cos θ = .
|a||b|
Two vectors a and b are said to be orthogonal or perpendicular if the angle between them is
90◦ . In other words,
a and b are orthogonal ⇐⇒ a · b = 0.
Example 2.12. 2i + 2j − k is orthogonal to 5i − 4j + 2k because (2i + 2j − k) · (5i − 4j + 2k) =
(2)(5) + (2)(−4) + (−1)(2) = 0.
Let a=ha1 , a2 , a3 i =
6 0. The angles α, β, γ in [0, π] that a makes with the x, y, z axes respectively
are called the direction angles of a.
z..
..
........
....
...
... a=ha1 , a2 , a3 i
..
..........
... ...
... ...
... .
...
... ...
... ...
... ..
.
... ..
..
γ ...
... ....
...
..............
...............
..... .. ....
α β
.. ..... ....
..... ...........................
... ...... ................
................
... ................
....
..... ................ ...
..
.
.
...
.
.
.. ...............
y
..
.....
..
....
.
...
.....
.....
.....
.............
x .
..
The cosines of these angles, cos α, cos β, cos γ are called the direction cosines of a. We may
express a vector a = ha1 , a2 , a3 i in terms of its magnitude and the direction cosines.
Note that the scalar projection is negative if θ > 90◦ . Moreover, in figure 11. SR = PR−PS =
a·b
b − |a| 2 a. Thus the distance from R to the line P Q is given by
a · b
|RS| = b − a .
|a|2
Example 2.14. Find the scalar and vector projection of b = h1, 1, 2i onto a = h−2, 3, 1i.
√
Solution. |a| = (−2)2 + 32 + 12 = 14. Thus the scalar projection of b onto a is a·b
p
√ |a| =
√1 ((−2)(1) + (3)(1) + (1)(2)) = 3/ 14.
14
The vector projection of b onto a is √3 a = 3
= h− 37 , 14
9 3
, 14 i.
14 |a| 14 a
Exercise 2.15. Find the angle between two long diagonals of a unit cube.
[70.5◦ ].
Exercise 2.16. Prove the Cauchy-Schwarz inequality: |a·b| ≤ |a||b|. Determine when equality
holds.
c Department of Mathematics 10
Semester II(2005/06)
Proof Exercise.
Corollary 2.20. b × c is perpendicular to both b and c.
Proof
b1 b2 b3
b · (b × c) = b1 b2 b3 = 0.
c1 c2 c3
c1 c2 c3
c · (b × c) = b1 b2 b3 = 0.
c1 c2 c3
Theorem 2.21. If θ is the angle between a and b, 0 ≤ θ ≤ π, then |a × b| = |a||b| sin θ.
Proof First we need the following identity
(a2 b3 −a3 b2 )2 +(a3 b1 −a1 b3 )2 +(a1 b2 −a2 b1 )2 = (a21 +a22 +a23 )(b21 +b22 +b23 )−(a1 b1 +a2 b2 +a3 b3 )2
which can be easily verified by direct simplification of both sides.
Using this identity, we have
|a × b|2 = (a2 b3 − a3 b2 )2 + (a3 b1 − a1 b3 )2 + (a1 b2 − a2 b1 )2
= (a21 + a22 + a23 )(b21 + b22 + b23 ) − (a1 b1 + a2 b2 + a3 b3 )2
= |a|2 |b|2 − (a · b)2
= |a|2 |b|2 − |a|2 |b|2 cos2 θ
= |a|2 |b|2 sin2 θ.
c Department of Mathematics 11
Semester II(2005/06)
Proof Exercise.
The relation a · (b × c) = (a × b) · c can be proved by direct expansion in component form. Al-
ternatively, it can be deduced by the property of the determinant: If two rows of a determinant
are switched, the determinant changes sign. Therefore,
c Department of Mathematics 12
Semester II(2005/06)
a1 a2 a3 a1 a2 a3 c1 c2 c3
a · (b × c) = b1 b2 b3 = − c1 c2 c3 = a1 a2 a3
c
1 c2 c3
b
1 b 2 b3
b
1 b2 b3
= c · (a × b) = (a × b) · c.
a1 a2 a3
In fact, a · (b × c) = b1 b2 b3 is the algebraic or sign volume of the parallelepiped
c c2 c3
1
determined by a, b, c.
b×c
.. .
.........
.... .. ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..........
.... .
... .
.. ..... .. ..
... ... .. ... .
... .... ..
. ... ..
.............
...a ...
. .... ...
.
......... ... ..... ... ... ... ... ... ... ... ... ... ... ... ... ... .....
...
.
.... ...
..
.. .......... .. .. ..
... . .. .. ..
... .. .
h = |a| cos θ ... ....
... ...
. .
..
. .
..
.
..... ... ... ... ... ... ... ... ... ... ... ... ... ........... ... .......
.
..
θ ... ...
c
........... ..........
........... ..
.. .........
..
... ... .......
.............
.. ... ............
.. ..... ..
Area of the base
............ .....
................................................................................................... parallelogram = |b × c|
b
Figure 14 Volume = |a||b × c|| cos θ|= |a · (b × c)|.
Corollary 2.24. The vectors a,b,c are coplanar (i.e. they all lie on a plane) if and only if
a · (b × c) = 0.
.....................
..................... ...
....................
..................... ...
.
..
....
.
....
.
............
. b ...
...
...
...
.
a .. .
..
.
....
.
.. .
... ...
...
...
. .. .......... ...
.... ...
... ... ..
....
. ...
..
. ..
. ..
..... ...
... ... .
...... ...
... ... ..
....
. ...
..
. ... ..
..... ...
.. ......... ...
..
.
.
..
..
.
.
.. .......
......... ..
. .
.... .......... .
.... . . . . ..
. .
..
.. . .
.
................ c ...
...
. . . .
... ..
. ...
... ......................................... .
..
.
... . . . . . . . . . . . ....
.
.
.. 0
..................................... ..
...
... ..
........ ..
... ..
.......... ....................... ........ . ....
.. . . . .
Example 2.25. Show that the vectors a = h1, 4, −7i, b = h2, −1, 4i, c = h0, −9, 18i are
coplanar.
1 4 −7
Solution. As a · (b × c) = 2 −1 4 = 0, it follows from 2.22 that a, b and c are coplanar.
0 −9 18
Example 2.26. Suppose a rigid body rotates with angular velocity w about an axis ` through
a point O. The angular velocity w is represented by a vector along `. If r is the position vector
from O of a point inside the rigid body, then the velocity at this point is given by w × r.
Exercise 2.29. Suppose the vectors a, b,c,d are coplanar. Show that (a × b) × (c × d) = 0.
2.5. Lines and Planes. Let L be a line passing through a point P0 (x0 , y0 , z0 ) i the
direction of the vector v = ha, b, ci. Then any point P on L has position vector r = r0 + tv for
some t ∈ R.
c Department of Mathematics 13
Semester II(2005/06)
P0 v P
..................................................................................................................................................................................................
......... .......
........
L
. ........
... ........
... ...............
...
... ........
hx0 , y0 , z0 i = r0 ...
...
r ........
........
........
..
. . ..... .........
... ........
........
... .........
.........
0
Figure 16 Vector equation of a line
The numbers a, b, c are called the direction numbers, or direction cosines of the straight line.
If a, b or c is zero, we may still write the symmetric equation of the line. For example, if a = 0,
we shall write the symmetric equations as
y − y0 z − z0
x = x0 , = ,
a a
which is a line lying on the plane x = x0 .
L1 : x = 1 + t, y = −2 + 3t, z = 4 − t,
L2 : x = 2s, y = 3 + s, z = −3 + 4s,
are skew, i.e. they do not intersect. Hence they do not lie in the same plane.
has a (unique) solution in s and t. The first two equations give t = 11/5, s = 8/5. But these
values of t and s do not satisfy the last equation. Thus, L1 and L2 do not intersect.
.....
..
..... L2
.....
..
.....
.
....
.....
.....
.....
.
..
......
.
.....
....
........................
.......................
..
...................... L1
. .
.. ......
...
...
...
. .
...
.... .
. ......... .
..............
.. ....
.....
.....
.....
..
.......
....
.....
.....
.....
..
......
.
..
.....
Figure 18 n · (r − r0 ) = 0
n · (r − r0 ) = 0
Example 2.31. Find an equation of the plane passing through the points P (1, 3, 2), Q(3, −1, 6)
and R(5, 2, 0).
−→ −→
Solution. P Q= h3 − 1, −1 − 3, 6 − 2i = h2, −4, 4i. P R= h4, −1, −2i. Thus a normal vector n
to the plane is given by
i j k
−→ −→
P Q × P R= 2 −4 4 = h12, 20, 14i.
4 −1 −2
Exercise 2.32. (a) Find the angle θ, (0 ≤ θ ≤ 90◦ ) between the planes x + y + z = 1 and
x − 2y + 3z = 1.
(b) Find the symmetric equations for the line of intersection of the planes in (a).
P1 .
.............
......
..... ....
.. ....
b ...
...
...
...
... ...
.
. .
...
........... ... ... ... ... ...... ..........................................................................
..... ....
. . ..... .....
.
... .... .. ....
.
.. ... .....
..
.
.
..... ...
..... .....
. n ...
.
. ..
..
.....
.....
..... ... ........
. .....
..... ..... ......... .....
..... ............................................... .....
..... .....
.... ....
.
....
.....
.
..
P0 N .
.....
....
.
..
....................................................................................................................
Example 2.34. Find the distance between the parallel planes 10x+2y−2z = 5 and 5x+y−z =
1.
Solution. The planes are parallel because their normal vectors h10, 2, −2i and h5, 1, −1i are
parallel. Pick any point on the plane 10x + 2y − 2z = 5. For example, (1/2, 0, 0) is a point on
10x + 2y − 2z = 5. Then the distance between the two planes is
√
|5(1/2) + 0(1) + 0(−1) − 1| 3
p
2 2 2
= .
5 + 1 + (−1) 6
Example 2.35. Find the distance between the skew lines:
L1 : x = 1 + t, y = −2 + 3t, z = 4 − t
L2 : x = 2s, y = 3 + s, z = −3 + 4s
Solution. As L1 and L2 are skew, they are contained in two parallel planes respectively. A
normal to these two parallel planes is given by
i j k
1 3
−1 = h13, −6, −5i.
2 1 4
Let s = 0 in L2 . We get the point (0, 3, −3) on L2 . Therefore, an equation of the plane
containing L2 is hx − 0, y − 3, z − (−3)i · h13, −6, −5i = 0. That is 13x − 6y − 5z + 3 = 0. Let
= 0 in L1 . We get the point (1, −2, 4) on L1 . Thus, the distance between L1 and L2 is give by
|13(1) − 6(−2) − 5(4) + 3| 8
p
2 2 2
=√ .
13 + (−6) + (−5) 230
Exercise 2.36. Find the equation of the straight line passing through the point P0 (1, 5, −1)
and perpendicular to the lines L1 : x = 5 + t, y = −1 − t, z = 2t and L2 : x = 11t, y = 7t, z =
−2t.
y−5
[ x−1 z+1
2 = −4 = −3 .]
c Department of Mathematics 16
Semester II(2005/06)
........
.. ................................
... .. ........................
................. .
z
....
.
..
...
...........
.
..
.. ................
... ... .........
. .
... ...
.. ..................... . .
.. .. .. .................... .. ..
.. .. ... . ... .
.
.. ..
.
.
.. .. ... . .
.. .. .. ... .... .... ...
... .. . .. .. .... ...
... .. .. .. .. .. ...
... .....
. . .
.. ..
. .
...
... .. .. .. ..
................ ... ... .....
.. ....................... . ..
.. .... . ... ..
... ... ................................... ........ ..
................ ... .... . .......... ..
.
... ........ .
...... ........................................ ... ....
. .. . ...
........................ .. . ...
................. .................... .... ......................... ...
................ ....................
....................... .................... ....
. ... ....
.... ...........................................
..... .. ...............
..... ................
..... ....................
.
..
.........
.
.. ......... y
.......
x
z.. z..
.. ..
.......... ..........
... ...
. ....
.......
............... .... ....................... ..
.. . ...... ....
......
..... ...
....
. ... ... ...
..... ..
. ..
....
..
.
......... .. ......... ..
... ............ . ...
... ... ..
... ... ................................................. .. . ............................
... ... ... .. ... .... ..
... ... ... .. ... .
..
..... .. .
..........
... ... ... .. ... .
..
.
.... .. ......... ........
..... .. ....... ...
... ... . ... .................. ... ... ... . .. ........ ...... ......
... . .... .. .. .. .. ... ... ..
..
. . . ..... ...
....... ... ... .. ... ... ..... . ....... ..... ....
... ...
.
... ....................................................................... ...... ... ....... ... ....
. ...
... ...
.... ..
... .......
........ ...
... ... y ...........
.... . ..
.. . .....
................................................ .................................................................................
.. ...... .... ..... .
...
y
........
......... .......... ..... ...... ...
.. .... ...
.
........ .
. ..
...
...
...
. ....... ...
... ..
.. ... .. .
... ............ .
.
.
..
. .
.. ................. .
. . .. ........ . ... ...
..... ... ... ....
..... ... ..
...... .....
.... ..... ........
...
......
.......... ... ..... ... ... ...... ... ..............
...... ....
x ...
.....
.
..
.....
.
. . .
... .......
..
...
..... ... .....
...
... ......... .. .....
.... ..... ...........
.
. ...
......... .....
.......... ............. ......
....... .............................
x
Figure 21 x2 + y 2 = 1 y2 + z2 = 1
Using translation and rotation, the equation can be expressed in one of the following two
standard forms:
y2 z2
Example 3.3. The graph of the equation x2 + 9 + 4 = 1 is an ellipsoid.
c Department of Mathematics 17
Semester II(2005/06)
z..
..
..........
...
...
...
...
(0, 0, 2) ...
...
.
...
.... ...
.. .. .............
.............. ...... .... ...... ..................
.
.. .
.. .... .. ... ..
. ......
.. . ..
. ..... ... ... .......... ... ..... ... ... ...... ... ... .......
. .. ... ... . ...
... ...... ..
.
... .... . .......(0, 3, 0)
....... .... ... ..
...
..........
...
. .......
..
....................................................................................................
. ...
y
... ............ ..... ........ ... .......
...........
... ... . ..... .
..... .............................................................. ........
...... ..... .. . .....
........... ... .. .......
(1, 0, 0) .....................................................................
.
....
.....
..........
.......
x
y2 z2
Figure 22 x2 + 9 + 4 =1
2 2
The vertical traces (or sections) in x = k are ellipses: y9 + z4 = 1 − k 2 , where −1 < k < 1.
2 2
The vertical traces in y = k are ellipses: x2 + z4 = 1 − k9 , where −3 < k < 3.
2 2
The horizontal traces in z = k are also ellipses: x2 + y9 = 1 − k4 , where −2 < k < 2.
Example 3.4. The graph of the equation z = 4x2 + y 2 is an elliptical paraboloid.
z..
..
..........
.
................. .... ..........................
. .. ......... ... . .......
. ... ... . .
......
.
....
. .. ... ...
...
..... . .. ... ..
.......... .
.
.
. .
.
. ..............
.
.
... ............ . . ........ ..
... ............................................... ... ....
.. .. .. .. ..
..
..
..
.. ... ... .. ..
.. .. . . .. .. ...
..
..
... ... ....... ... ...
. . . .
. .
..
.. .....
. . .. ... .. ...
.. .. ... .. .... ........
... .. . .. .. ....
.
... .. . . ... ... . . . .
.
... ... .. ... .....
..... .. . .. .....
......... ... .......
.................................................................................
.....
.... .. y
..
......
.. .....
....
.....
.....
..........
......
x
Figure 23 z = 4x2 + y 2
The horizontal traces in z = k are ellipses: 4x2 + y 2 = k, where k > 0.
The vertical traces in x = k are parabolas: z = y 2 + 4k 2 .
Similarly, the vertical traces in y = k are parabolas: z = 4x2 + k 2 .
x2 z2
Example 3.5. Sketch the surface 4 + y2 − 4 = 1.
z..
..
..........
..
...
.........
. ........... ... ......................
. . . . . ... . ........
... . .... ... .. ......
... ... .. . ....
.
..... . .. .. ...
..... .
.
.
. .
. .. ........
.
...... .
. .
. ......
.
.
.......... .
. .
...
... .................. .. ....................................
... .................. . ....
... .. .... .. .... ....
... .. ... .. ......
..
... .. ... .. ..
...
... .. .... ... ...
. .. .
. .
.
... ... ... ..... . .
... ...... ... .. ...... ..... ..... ....
..... . .. ...
... ................................................................................
.... .. . ...
........... ....... ... ... ...........
y
.... ................................................. .....
.. ..... .. .. ..
... ..... .. . ...
....... .. . ..
.. .
. ......
.
.. ... ........
.
. ..
............ .. . .. ...
...... .... ................................................ .... .....
x . ................. ..
.
......... .....
. .......
.
.... ... .. ........
....
..... ... .. .....
..
. ..
. ..
...
... .. ...
.....
....... ..
. .. .........
.
.
.......... .
.
. ...
..........
.
....................................
x2 z2
Figure 24 4 + y2 − 4 =1
x2 k2
The traces in z = k are ellipses: 4 + y2 =1+ 4 .
c Department of Mathematics 18
Semester II(2005/06)
z2 2
The traces in x = k are hyperbolas: y 2 − 4 = 1 − k4 .
2 z2
The traces in y = k are hyperbolas: x4 − 4 = 1 − k2 .
....
...
.........
z .
z
...
......... ...
z
...
......... .
.... ... .... ... ... ... ...
.... ... .... .... ... ... ...
.... ... .... ... ... ...
...
... ... .....
. .....
..
...
.... . ......
... .... .... .
.... ....
... ... ..
..... ... ..... ....
....
... ... ...
..... ... ..... ..... .. .....
...
... ... ...
.....
..... ... .......
. ......
......... .... ..............
... ... ..... ... .... ....................
... . ..... ..... ..
... ... ... ..... ... ........ ...
... ... ... ..... .. ..... ...
.. .. ............
...............................................................................................................................................
.. ..
...
... x ...........................................................................................................................................
.... . ....
..... ... .....
. .
x .......................................................................................................................................
..
...
x
... ...
..... ... .........
... ... ...
... ..... ...
... ... ... ...
... .... .....
..... .. .... .. .
........ .... .................
. ......
.... .
. ... ..... ... ..... ....
.
.. .
. ... ..... ... ..... ..... ....
.....
.
.
...
.
..
. ... ..
....... .
.
.....
.... ..
.
.
....
.
.
.....
....
..
. .
.
.
... ... .... ....
.
.... ....
.... . .... ... ...
.
. .
. .... ... .
. .
. ...
..
... .
.
. .... .. ...
. .
.
. ...
.. .... . ...
.... .... .. .... ... ....
.. ... ..
.. . ..
Figure 26 4x2 − y 2 + 2z 2 + 4 = 0
·
..
. . .
..
. ...
..... (3, 1, 0) ..... .....
.
......
. .
. . ..
.
.
.
.
...
. ...
..
...
. ...
.
. .
... .... ... .
..
. ..
...... .
.
.... ...... ...
...
..... ...... ...
..... ...... . ..
..
..... ....... ....
...
..........
.
.
.
............ ..
. . ... . ..
.
. ..................
......
x ..
Figure 27 x2 + 2z 2 − 6x − y + 10 = 0
By the method of completing squares, the equation can be written as y − 1 = (x − 3)2 + 2z 2 .
If we make a change of coordinates: x0 = x − 3, y 0 = y − 1, z 0 = z so that the new origin is at
(3, 1, 0), then the equation becomes y 0 = x02 + 2z 02 . Therefore it is an elliptic paraboloid with
vertex at (3, 1, 0).
Exercise 3.8. Describe the traces of the surface z = xy.
c Department of Mathematics 19
Semester II(2005/06)
Exercise 3.10. Find the equation of the surface obtained by rotating the curve y = x2 about
the y-axis.
The graphs of the quadric surfaces are summarized in the following table.
z.. z..
..
..
........ .........
.. 2
.... ... z2 x2
...
. .
..
........................ ... ............................. c2
= a2
+ yb2 . Horizontal
.
.... .
....
............................................ x 2 y2 z 2 ..
...... .
.....................................................................
......... ..... ... .... ............ + + = 1. The ..... . traces are ellipses. Verti-
.... . . . ...
. ...... a2 b2 c2 ..... .. .....
.
... .. .
. .. ... ..... .. .....
.
... ... ... ............ ... ...... ... ... ........ ... ... ... ...... ..... .... .........
....... ... .. ... .. ...... traces are ellipses. When ..... .. .....
..... .. ..... cal traces in x = k and
...... .... ......................................................................................
.....
...........
...
... .........
... .. .......
.. ... y ........
.....................
......... ........ ....................
... ............... .. ...... ..
.............................................................. ......
.......... a = b = c 6= 0, it is a ................ ..... .............. y = k are hyperbolas for
....
.....
...... ........ .... . . .
.
.
......
.
..
...... .....
..
.. .
...........
.....
.....
..
........... y
.......... .... .. ....... sphere. .
..................... ... ... ... ... ... ... ... ............ k 6= 0 but pairs of lines
.... ........................................................ ..
.
......... ...... .
.....
.......... ..
.
. ................ ..
x ...... x ... ..........................................
for k = 0.
...................
z ............... .... ....................
........ . .... .....
...... z
.......... ...
. −x −
2 y2
+ z2
= 1. Hor-
.. .... .....
.......... ........ .... ........
...................... .. .............................
.. a2 b2 c2
.. .... ....... ..
.. ... .... ....... ... ........... ...
. ..
.. .... .
. .
. .. .
. ... 2 2 .... . ..
. izontal traces in z = k
................. .... ........ .. .............. .... z x y .... . ..
....
...
.
..
.. .
. .
. .
... . ..
.
. c
= a2
− b2
. Hori- .....
...... ... .........
.
... .
... .. .. ...... ..
. .... .... ............ ............
........ are ellipses if k > c or
... ......... ..... ........ .... ........... ..
..... ...
... . ................... ...... ........................ ..... ....
.. .. ..........................................................................................
zontal traces are hyper- ...
.......
.... ... .............................. ..................
.. ...................... .......
...
.......
... ... ....
... ...
.
y bolas. Vertical traces are ...
....... ...............................................
..... .... ........ y k < −c. Vertical traces
........ . .............. ... ... ... ... ... .................................... .
...... ...
... .... ... ..... ............................... are hyperbolas. The two
.. . ...................................................... ... ................ ... ..........
. . ...
........
............................
............. ..
........ parabolas. Here c < 0. ...
.. ..
..... ....
.
..
....
....
.......... .... .. minus signs indicate two
x ..... ........ ... ... ... ......... ... ... ... ........
x .....
.... .
..
......
.
...... ... sheets.
............ ......
Hyperbolic Paraboloid ......................................
Table 1
y.
...
..........
....
(x, y)
...
...
• .....
...
.....
... .....
... ..
......
.
... .....
...
...
r .....
....
.....
... ..
......
.
... ....
.....
... ........
... ..............
θ .
.
..
..
................................................................................................................................................................
. .
x
O ...
...
...
...
...
The relations between Cartesian and polar coordinates are given by the following formulas:
y
r2 = x2 + y 2 , tan θ = x
x = r cos θ, y = r sin θ
y. P (1, π4 )
....
.........
• ....
.....
.....
..
.... ..
......
.
....
.. .....
...
... 1 .....
.....
.....
... .
...
.
..
... .....
... ......... π
.. ....
. .............................
. 4
.......................................................................................................................................
.... .. x
..... ..
..... ...
.....
.....
O ...
..
.. ...
..... ...
..
....... ...
..
..... ...
..
..... ...
....
(−1, π4 ) = (1, 3π ...
..
.
...
4 )•
..
.. ...
..
..... ....
....
Q
Figure 29 Polar coordinates of P = (1, π4 )
Exercise 4.1. Find the equation in polar coordinates of the curve x2 + y 2 = 2x.
[Answer : r = 2 cos θ.]
.....
. .
.
........
... .
.
. (r, θ, 0)
x
Figure 30 Cylindrical coordinates
The relations between Cartesian and cylindrical coordinates are given by the following formulas:
r2 = x2 + y 2 , tan θ = xy , z = z
x = r cos θ, y = r sin θ, z = z
Example 4.2. The surface whose equation in cylindrical coordinates is z = r is a double cone
with the origin as the vertex.
z..
..
.........
...
....
..
...
...
.
.. ... ..
. ... .. .
. ........................ .... ........................................
........... ....
.........
..
......
........... ... ......
..... ........................ . .............
..... ....................................................................
..... . ..
..... .... .....
.....
..... ...• .... .....
.. . .....
.....
.....
.....
.....
...
.
. .
... .. ...... .
P (r, θ, z)
. ...........
.....
..... ... ........... ...
..... ... ....... ..
.............. .
..............
O
.....
.... ..... ........ ..
........... ......................................................
...............
........... ..... ......
...
.. ..
............
.
•Q ......
.....
...............
...............
...............
...... ..... ..... ..................
...
.
...................
.... .......
.....
......
.......
........ y
. ... .... ...... ..
...
..... ....... .. ... ... ... ... ... ... ... ... ... ... ... . ........
..... ..... ... ... . .. ... ... .......
..... ............ . ....
...... .... .....
...... ..........
...........
.....
...
....... ......................
.... .. . ......................... ..
. .. .
.. .. ..... ..
. .. ... .. .
. .
..
.....
.
.
......
........
.........
x
Figure 31 A double cone
Let P be a point on this surface with cylindrical coordinates (r, θ, z). Since z = r, the triangle
OP Q in which ∠OQP is a right angle is isosceles with OQ = r = z = P Q. Thus the cone
opens up an angle of 45◦ with the z-axis. To convert the equation to Cartesian form, we can
square both sides of z = r, thus z 2 = x2 + y 2 is the Cartesian equation of the doublepcone. If
we take positive square root on both sides, the graph of the resulting equation z = x2 + y 2
is the inverted cone on the upper half space z ≥ 0.
Exercise 4.3. Find the equation of the ellipsoid 4x2 + 4y 2 + z 2 = 1 in cylindrical coordinates.
[Answer : z 2 = 1 − 4r2 .]
4.3. Spherical Coordinates. Another coordinate system in 3-dimensional space is the
spherical coordinate system. Given a point P (x, y, z) with P 0 the foot of the perpendicular
from P onto the xy-plane, let ρ ≥ 0 be its distance from the origin, θ the angle that OP 0 makes
with the x-axis and φ the angle that OP makes with the z-axis. Here θ is measured in the
counterclockwise sense from the x-axis with 0 ≤ θ ≤ 2π, and φ is measured from the z-axis
with 0 ≤ φ ≤ π. Note that OP 0 = ρ sin φ and P P 0 = ρ cos φ.
c Department of Mathematics 22
Semester II(2005/06)
z..
..
........
....
...
...
...
...
...
...
• P (ρ, θ, φ)
..
.......
..... ..
... ..... ...
... ...
...... ....
..
......................... ...
φ ... ..........
... .......
ρ ...
...
.........................
.
......
. ...
..
.
... .......
. .....................
.
. . ............................
..
.. .
.
....
. θ
..... ....................... .........
.
.......
.. ... ... ..................
.................. ...
.
.
.
....
.
......
. .....
.....
.
...
...
.
.............. y
.
... . .
..... ...... .
..
.. .. .
..... ...
....
.....
.....
.....
• 0 .....
.....
. .
.
........
... .
.
. P (r, θ, 0)
x
Figure 32 Spherical coordinates ρ ≥ 0, 0 ≤ θ ≤ π
The relations between Cartesian and spherical coordinates are given by the following formulas:
p
ρ = x2 + y 2 + z 2 ,
cos φ = ρz , 0 ≤ φ ≤ π
x
cos θ = ρ sin φ
5. Vector Functions
Definition 5.1. A vector function r(t) is a function whose domain is a set of real numbers
and whose range is a set of vectors.
In other word,
r(t) = hf (t), g(t), h(t)i = f (t)i + g(t)j + h(t)k.
f, g, h are called the component functions of r.
√
Example 5.2. Consider the vector function r(t) = ht3 , ln(3 − t), ti.
For each of the component functions to be defined, we must have 3 − t ≥ 0 and t ≥ 0. Thus
the domain of r is [0, 3). The image of r traces out a curve in R3 .
In general if r(t) = hf (t), g(t), h(t)i is a vector function, then x = f (t), y = g(t), z = h(t) give
the parametric equations of a curve in R3 .
c Department of Mathematics 23
Semester II(2005/06)
z..
..
........
....
...
... ...
... ...
... ....
............ .
... ....
....... ... ...................
R3
r ...
.. ....
... ......
.. ........
•
.......................................................................... .................................................... ... ........ ..
..... ...............................
.....
...
R t .
..
.
.
.....
. ...
...
... ..................
....... ........
.
• ..................
........ y
.....
..
......
.
...
....
.....
.
....
.
. r(t)
.........
.......
x
Figure 33 A vector function
Example 5.3. The vector function r(t) = h1 + t, 2 + 5t, −1 + 6ti defines a curve which is a
straight line in R3 .
Example 5.4. Sketch the curve whose vector equation is r(t) = hcos t, sin t, ti.
Solution. The parametric equations of the curve are x = cos t, y = sin t, z = t. Consider a point
P (x, y, z) on this curve. Since the x, y and z coordinates of P satisfy the relation x2 + y 2 = 1,
it lies on the cylinder x2 + y 2 = 1.
z
....
..........
......................... ... ..................................
......... ....
....... ... .
................ .. ........
... .............................................................. ....
....
... ... . . .
... ..... ..
............................................. ...
................... ............. ...
................. ...................................... .... ..
.....
... ....... ..
.. ..
... .
... .. ...
... .. ..
... ... .
. .......
. ... ..
... . . . . .
............ .. .......................
.
. .. .... ..
................. ..... ........ ...
.
............................................................... .... ..
... . ....
.. ..
... ..
... .. ..
... ................................................... ......
............ . ......... ..
....
......... ..... .
. ...
...
...
.. . .. . . . . . . . ...................................................................
.
.
...... ....
................................................................................................
.. y
.
.....
........
.
.....
.....
.....
.........
.......
x
Figure 34 A helix
Moreover, P lies directly above the point (x, y, 0), which moves counterclockwise around the
circle x2 + y 2 = 1. Since z = t, the curve spirals upward around the cylinder as t increases.
The curve is a Helix.
Example 5.5. Find the vector function that represents the curve of intersection C of the
cylinder x2 + y 2 = 1 and the plane y + z = 2.
Solution. Since C lies on the cylinder which projects onto the circle x2 +y 2 = 1 on the xy-plane,
we can write x = cos t, y = sin t with 0 ≤ t ≤ 2π. Since C also lies on the plane, its x, y, z
coordinates should satisfy the equation of the plane. Thus, z = 2 − y = 2 − sin t. Consequently,
the vector equation of C is r(t) = hcos t, sin t, 2 − sin ti.
c Department of Mathematics 24
Semester II(2005/06)
z ..
......
.................. ..............................
....................................... .........
.. ..
. ..... .
. ..... .
.. ....................... ............ ...... ..
.. ....
.......... ...
..... ................ .. .....
..... .
.
.
......................... .. ......... ........ ....
..... . .........................................
..... .
. ... . . ....
..... .............................
. .. . . ...
.......
.....
..... .... ...... ......... ...
..... ..
. ... .
...... ... . ...
..... ....
.
. .. .. ...
..... . ..
y+z =2 ..... .... ..... .. .......... ...... ...
..... .. ... . ..... ...
. ...
..... ... .. . ..... ...
....... ... . ..... ... ..
..... ... .
. ..... ... ....
............ ......
... ....... .......
..... .....
• ..
..
.... ... ..
... .....
...
... ... .....
.
... ..... ..... .. ... ... .....
... ..... ..... ... ... .. ......
... ..... ..... .
..... ....... .... ... .....
... ..... ....... ...... .....
.....
... ..... . ...... .... .....
... ..... .. ....... .
. .....
...... ....... . .....
... ... ......... .. . .....
... ........... ................... .....
... ... ... ... ... ... ... ... ..................... ... ... ... . .... ....
......... .. ..... .. ... .. .......
....... ... ... ... ......... ... ... ... ... .............................................
x +y =1 2 2 ....
.......
.........
.
. .. ..
. ... .. . .....
... ... .... ...
. ..
... ...
. ... y
..............
....................................................................... .......
.. ..... ............
.. ....... ............
.
.
.....
.....
.....
..........
.......
x
Figure 35 An ellipse
The curve C is an ellipse with centre (0, 0, 2) and it inclines at an angle 45◦ to the horizontal
plane.
Let r(t) = hf (t), g(t), h(t)i. The limit of r(t) as t tends to a is defined by:
lim r(t) = hlim f (t), lim g(t), lim h(t)i.
t→a t→a t→a t→a
Example 5.9. Let r(t) = h1 + t3 , 2t, 1i. Find the unit tangent vector to the curve defined by
r(t) at the point where t = 0.
2j
Solution. First we have r0 (t) = h3t2 , 2, 0i. Thus, r0 (0) = h0, 2, 0i = 2j. Therefore, T(0) = 2 =
j.
Example 5.10. Find parametric equations for the tangent line ` to the helix with parametric
equations x = 2 cos t, y = sin t, z = t at t = π2 .
Solution.
z ..
.........
................ .... ..............................
.... ................ ... .......
. ..
.......... .. .....
... ........................................................................... .....
... .. ..
.......
... . ....
..
... .......................................................... ....
......... .. ............. ........ ...
.................................................... .......
...
... .. ..
... .. ..
...
...
..
....
..
..
....
`
...
..
. .
..... ....
.
.....
.. .....
... .............................................................. .... .....
.
........ . . . . . . . . .. ....... .. .
. ..
.
...... ...
.................................................. ..... .....
... .. ... ........
... .. .... .........
... . π
.
... ...................................................................... • (0, 1, .
... . .
.
2
)
......
......
............ .
.....
.... .. ...
.....
........................................................................................................ y
..........................................................................................
..
.. .. .
..... ....
..... .....
.... .....
..... ....
....... .....
. .
........
....
x
Figure 37 The tangent to the helix
The vector equation of the helix is r(t) = h2 cos t, sin t, ti. Thus, r0 (t) = h−2 sin t, cos t, 1i
and r0 ( π2 ) = h−2, 0, 1i is a tangent vector to the helix at t = π2 . Therefore, the parametric
equations of the tangent line ` are given by: x = 0 + (−2)t, y = 1 + (0)t, z = π2 + (1)t. That is
x = −2t, y = 1, z = π2 + t.
Given a vector function r(t), we may compute successively r0 (t), r00 (t), r000 (t) etc, provided they
exist.
Theorem 5.11. Let u and v be differentiable vector functions of t, c a scalar and f a real-
valued function. Then we have the followings:
d
1. dt (u(t) + v(t)) = u0 (t) + v0 (t).
d
2. dt (cu(t)) = cu0 (t).
3. dt (f (t)u(t)) = f 0 (t)u(t) + f (t)u0 (t).
d
d
4. dt (u(t) · v(t)) = u0 (t) · v(t) + u(t) · v0 (t).
d
5. dt (u(t) × v(t)) = u0 (t) × v(t) + u(t) × v0 (t).
6. (Chain Rule) dt d
(u(f (t))) = f 0 (t)u0 (f (t)).
Exercise 5.12. Suppose |r(t)| = c, where c is a positive constant. Show that r(t) is orthogonal
to r0 (t) for all t.
Let r(t) = f (t)i + g(t)j + h(t)k be a continuous vector function. The definite integral of r(t)
from t = 1 to t = b is defined as:
! ! !
Z b Z b Z b Z b
r(t)dt = f (t)dt i + g(t)dt j + h(t)dt k.
a a a a
Z π
2
Example 5.13. Let r(t) = 2 cos ti + sin tj + 2tk. Find r(t)dt.
0
c Department of Mathematics 26
Semester II(2005/06)
π
π π h iπ π2
Z
2
Solution. r(t)dt = [2 sin t]02 i − [cos t]02 j + t2 2
k = 2i + j + k.
0 0 4
6. Functions of several variables
6.1. Functions of 2 variables.
Definition 6.1. A function f of 2 variables is a rule that assigns to each ordered pair of
real numbers (x, y) in a set D a unique real number denoted by f (x, y). Here D is called the
domain of f . The set of values that f takes on is called the range of f . That is Range of f =
{f (x, y) | (x, y) ∈ D}.
We usually write z = f (x, y) to indicate that z is a function of x and y. Moreover, x, y are
called the independent variables and z is called the dependent variable.
y.
....
........
...
....
..
...
.. .
.....
. ...
............ ........................... ...
. .........
.... ....... ..
...
...
....... ...
.......
.. .....
f
... ................................
... ........ ................ ...........
...
... ... ................. ........
.. ... .............. ........
.. .
. .
....... ......... .......
...... .
·
.. . ....... .......
.. .... ...... .......
... . .....
·
... .
. .....
... (x, y).
.
. ......
. .
.
............................................................................................................................................................................... ..
x ... ...................................................................................
...
...
..
...
....
....
...
. R
..
..
....
.....
.
.
....
..
...
...
f (x, y) 0
...
D ......
.......
........
...
. ....
.
...
............... .....
....
. .......
................................................................
....
..
...
..
Figure 38 f : D −→ R
Example 6.2. Find the domain of f (x, y) = x ln(y 2 − x).
Solution. The expression x ln(y 2 − x) is defined only when y 2 − x > 0. That is y 2 > x. The
curve y 2 = x separates the plane into two regions, one satisfying the inequality y 2 > x, the
other satisfying y 2 < x. To find out which region is determined by the inequality y 2 > x. Pick
any point in one of the regions and test whether it satisfies the inequality. If it does, then by
‘connectivity’, that whole region is the one satisfying y 2 > x, otherwise, it must be the other
region. For example, pick the point (3, 2). Since 22 > 3, the region satisfying y 2 > x is the one
containing (3, 2). Thus, domain of f is {(x, y) ∈ R2 | y 2 > x}.
y.
. . . . . . . ......... . . . . . . . . .
. . . . . . .. . . . . . . . . . .
. . . . . . . .... . . . . . . . . . .
. . . . . . . . . . . . . . . . .
. . . . . . . ... . . . . . . . . . . ...
. . . . . . . . . . . . . . . .......
. . . . . . . .... . . . . . . ....... ...... ......
. . . . . . . . . . . . ..
. . . . . . . ... . . . . ....... ...
. . . . . . . . . ....
. . . . . . . .... . ..........
(3, 2)
·
. . . . . . . .
. . . . . . . ... .......
. . . . . . ...
. . . . . . . .... 0
.............................................................................................................................................................. x
. . . . . . . ....
. . . . . . ..
. . . . . . . .... .......
. . . . . . . ..
. . . . . . . ... . ..........
. . . . . . . . . ...
. . . . . . . .... . . . .... ... ... 2
. . . . . . . . . . . . ....
. . . . . . . ... . . . . . . .... ...... ...... ..
y =x
. . . . . . . . . . . . . . . ....
. . . . . . . .... . . . . . . . . . . ......
. . . . . . .. . . . . . . . . . .
. . . . . . . .. . . . . . . . . . .
. . . . . . . . . . . . . . . . .
. . . . . . . ... . . . . . . . . . .
.
Definition 6.4. Let f be a function of 2 variables with domain D. The graph of f is the set
of all points (x, y, z) ∈ R3 such that z = f (x, y).
z..
..
.........
...
....
..
(x, y, f (x, y))
... .........
......................... ..........
f (x, y) ........
...
.
....
...
...
...
.....
...... · .....
....... ....................................
.. .......................
..
..
..
.. ........................... ......
..
..
.....
....
..... .....
..... .
..
..
... .. .. ..
... .. .. ..
... .. .. ..
... .. .. ..
... .. .. ..
... .. .. ..
....................................................................................................................................
.
.. .
...
.
. . .
..
.
.
..
..
..
... y
.
.....
·
..... ... ... ..
..
.. . . .......
........ ......................
. .
.
.. ..... .......... .
..... ... . ..... .
..... ....... ....................................................
.... .....
. .....
.........
.
. D
x .
Figure 41 z = 6 − 3x − 2y
Example 6.6. The graph of h(x, y) = 4x2 + y 2 is an elliptic paraboloid.
z..
..
..........
.
................. .... .........................
. .. . ........ ... ......
.. ... ...
....
. ...
.... ... .
..... .
. ..
.......... .
. ........
.
... ............ .
.
. ....... .... ....
. .
... ...................................
.. . ...
.. ....... ... ......... ... ... ... ... ....
.. .... .. ... ...
.. .
....
.....
.......
.. ........
.
...
.... ......................
.....
.......
... 4x2 + y 2 = k, k > 0
.. ...........................................
... . ..
... . .
... ... ...
..
.....
...... ... ...........
..............................................................................
.
.... .. y
.....
.....
..
.......
....
.....
.....
..........
......
x
Figure 42 z = 4x2 + y 2
The domain of h is R2 .Since 4x2
+ y2
≥ 0, the range of h is [0,∞). Each horizontal trace is
an ellipse with equation given by 4x + y 2 = k, where k > 0.
2
Thus, a large K gives rise to an ellipse with longer major and minor axes.
Exercise 6.10. Sketch the level curves of p(x, y) = −xy for K = −4, −1, 0, 1, 4. The graph of
p is shown in figure 46.
.......... ...
........ ...
.. ..... .......
.. ... ....... ...
.. ......... ........ ...
............ ... ... ...
..
... . . .. ..............
. .
.
............... .. .
. ... ................ . .
... ... . ..
... .. .. ... ... ...
... .. .. .. .. ..
.
.. ........ .. ..
. ... .... .... ......
... . . . ....
..... ... ... .. ... ... ....... .. ..
... .... ... ... ... .. .. ... ... . ..... ....
... .. ... ... ...... ... ... ... ... ... ...
... .
. ... ........
. ..
......... ... ...
.
. ... ...
. .
....................
... .. . .
........... ...................
... ...
... ...
.... .. ...... .
..... .
. ..
.......... .. ...
... .......................... ..........
.......... ...............
........... .... ... .
.
... .. . ...... .. ... ... ... ... ... ... ... ... ... ... ... ... ... ... ..................
.
.
............... ...
. ................
.
. .
....... ... .............
x ...
...
............................
......................
...............................................................
............. ....
. .
..
.
.........
y
.
... .
.. .......
Figure 46 p(x, y) = xy
....
.
....
. .. .
.
.
.
.....
....
...
f
.. .
. ...
.... . . .
.. ... ... ... .
. .. ... ..
. .
.. . .. ... .........
................. ...............................
.. ... ... ... . .. ....... ........... ... ... ...
...
... ........ ..........
.
.. ........
.. . ..
...... .. ........... ... ...
.... ....
.. ...... .........
. .
....
..
. ..... . ..... .........
.... ...... ... .
.
K=1 .. . ... .. .. ......... .
...... .... . ...... ... .... ... ... ... ..... ... .. .............
..
..... ... ..... .
. ..... ..
..... .
. ...
. .....
... ....... ... ....
...
.....
...
....... .
... .. .
.. ... ..
. .
.. . .. .. .
.
.
. .
.
...
.. ... . .
. ... ... . .. .. . . .. ..
.
.
..
.
.
y
...............................
.. . ...
.................
. ......................................................
R
...... ............ ..... ....... ..
........
... .......
... ................................................ .....
... ...
.....
.... ...
..... ..
2
..
... .........
. ..
......
.
..
.. ...
. ..... ... ......... ...
... .......... . ........ . . .
... ............... .......................... ...........
....... ...
... .... .............................................................. ...
... ..... ...
... ..... ..
..... ........ ...
. ..
........ ..... 2 2 2
x +y +z =2 ....
..... .......... .....
..... ...... ......
.
......
. ........ .. ..........
. ... .......... .. .. .....
.
.. . .. ..
. ........................... ..
....
..
...........
......
x
Exercise 6.12. Sketch the level surfaces of q(x, y, z) = x2 + y 2 − z 2 for K = −4, −1, 0, 1, 4.
c Department of Mathematics 30
Semester II(2005/06)
z..
..
........
.... f (x, y)
...
... .......................................................................................
...........
. ............ .................
L+ ........ ..
..... ... ... ... ......
....
·◦ ............................................ ......
..
..
.
......... ... ... ... ... ... ... ............. ... ... ... ... ... .................................................................................. ......
L ... .
... ....................................... .
.........
... .......................... .
. .....
L− ........... ... .. ... .....
... .. ... ....
.... ... .....
... ... .........................................................................................................
... ......... .. ...
... .. ...
... .. ...
... .. ...
... .. ...
.... .. ...
.. ...
...
...
...
................................ .. ...
...
. . .
.
....
.
. . .. .
.. .. ..
. ........ .....................................
.
...
.
. . ...............................................
.. ...
..
.
.
...
.
.
..
.. ...
....................................
.......... y
··
.
....
. . ..
.
.
...
. ....... ........ .....
.
.. .
...
.
.
(x, y) .
. .. .
. .. ......................... . radius of the disc
....
.
..... ..
... ... δ=
... ..
.
..... ... ... ... .... centred at (a, b)
..........
......
x
Note that f is not required to be defined at (a, b). The idea is that as (x, y) approaches (a, b),
f (x, y) approaches L. In other words, f (x, y) can be made as close to the number L as we
wish by requiring (x, y) sufficiently close to (a, b). This is the meaning of the above definition.
y.
...
..........
...
....
..
...
... ................... ....
.
..
....
. .
...
..
....... ......... ..
..... .........
...
...
.... ............
...
...
......
f
..... ...............................................
...
... ... ..... ................... .............
.
. ..... ............
... .
. ...... ......... ...........
..
..
.
..
. . ... ... ............ ........
.
............
............. .
·
.. .
. . .. ... ....... ..............
... . . .. ..
....... ......
.
.
. .
. ..
... ....... .. ....................... .....
...
... b .. .. ..
..
...... .
..................................................................................................
...
...
.
....
...
δ ...
... ... ... ......
. . ..
.....
....
....
( . ) R
...
..
.. ... ... L− L L+
.. .
. .
.....................................................................................................................................................................................
..... ..
...... .. a
..
.
....
. x
....... ...
...
D ............. .....
... ....................................................................................
....
..
...
..
The implication in definition 7.1 says that all points (x, y) which are inside the disc centred at
(a, b) with radius δ are mapped by f into the interval (L − , L + ). See Figure 49.
z..
..
.........
....
... ....................................................
... ................... ..........................
... .........
..... .. ............. ......
..
.. ..... ........................ ........... ......
........... .... ......
L .... .
.
.
....
.
. ◦ .
. ...
........
.
... ... .. .....
... .. .....
.. ... ....
... ... ...........................................................................
... ...........................
. .
.
... ..
... ..
... ..
... ..
.... ..
.
.................................... ..
...
. ................................
.
...
.
.
. ..... . . .. ..................................
.
.
. .. ................................
.
.....
.
.
.
...
. ..
..........................
. y
·
.....
. ..
..
. .
. .. .
.............
..
....
. ........................ .. ............
...
.....
..... (a, b)
..........
... .
x
Proposition 7.2. Let x = α(t), y = β(t) be the parametric equations of a path in R2 such that
(α(t), β(t)) lies in the domain of f (x, y) for all t in a certain open interval containing to and
lim α(t) = a and lim β(t) = b. Suppose lim f (x, y) = L. Then lim f (α(t), β(t)) = L.
t→to t→to (x,y)→(a,b) t→to
Proof Let any positive number. Since lim f (x, y) = L, there exists a positive number δ
(x,y)→(a,b)
p
such that |f (x, y) − L| < whenever 0 < (x − a)2 + (y − b)2 < δ.
Now because lim α(t) = a and lim β(t) = b, there exists a positive η such that |α(t) − a| < δ/2
t→to t→to
and |β(t) − b| < δ/2 for all t satisfying 0 < |t − to | < η.
Thus for all t satisfying 0 < |t − to | < η, we have
q q
0< (α(t) − a)2 + (β(t) − b)2 < δ 2 /4 + δ 2 /4 < δ
so that
|f (α(t), β(t)) − L| < .
This shows that lim f (α(t), β(t)) = L.
t→to
Exercise 7.3. Prove that if lim f (x, y) exists, then there its value is unique. That is
(x,y)→(a,b)
there is only one number L satisfying the definition 7.1.
x2 − y 2
Example 7.4. Show that lim does not exist.
(x,y)→(0,0) x2 + y 2
x2 − y 2
Solution. Let f (x, y) = . First let’s approach (0, 0) along the x-axis.
x2 + y 2
x2 − 02
lim f (x, y) = lim f (x, 0) = lim 2 = lim 1 = 1.
(x, y) → (0, 0) x→0 x→0 x + 02 x→0
along y = 0
Next let’s approach (0, 0) along the y-axis.
02 − y 2
lim f (x, y) = lim f (0, y) = lim = lim −1 = −1.
(x, y) → (0, 0) y→0 y→0 02 + y 2 y→0
along x = 0
Since f has two different limits along 2 different paths, the given limit does not exist.
xy
Example 7.5. Show that lim does not exist.
(x,y)→(0,0) x2 + y 2
xy
Solution. Let f (x, y) = 2 . First let’s approach (0, 0) along the x-axis.
x + y2
x·0
lim f (x, y) = lim f (x, 0) = lim 2 = 0.
(x, y) → (0, 0) x→0 x→0 x + 02
along y = 0
Next let’s approach (0, 0) along the y-axis.
0·y
lim f (x, y) = lim f (0, y) = lim = 0.
(x, y) → (0, 0) y→0 y→0 02 + y 2
along x = 0
At this point, we cannot conclude anything as the limit may exist or may not exist. Now let’s
approach (0, 0) along the path y = x.
x2 1
lim f (x, y) = lim f (x, x) = lim 2 2
= .
(x, y) → (0, 0) x→0 x→0 x + y 2
along x = y
c Department of Mathematics 32
Semester II(2005/06)
Since f has two different limits along 2 different paths, the given limit does not exist.
xy 2
Example 7.6. Let f (x, y) = . Show that lim f (x, y) does not exist.
x2 + y 4 (x,y)→(0,0)
Solution. Let’s approach (0, 0) along the line y = mx, where m is any real number.
x · (mx)2 x3 m2 xm2
lim f (x, y) = lim 2 = lim = lim = 0.
(x, y) → (0, 0) x→0 x + (mx)2 x→0 x2 (1 + m2 x2 ) x→0 1 + m2 x2
along y = mx
Thus, the limit as (x, y) approaches to the origin along any straight line is zero. However, we
still cannot conclude anything as the limit may exist or may not exist. Now let’s approach
(0, 0) along the curve y = x2 .
y2 · y2
lim f (x, y) = lim 4 = 12 .
(x, y) → (0, 0) x→0 y + y 4
along y = x2
Since f has two different limits along 2 different paths, the given limit does not exist.
y.
....
........
.....
.....
..
...
... .
..... ............
y2 = x
..... ... ..................
..... ..............
..... ... ......
..... .
. .. .... ............... .
..... .... ....... ..
..... ..... ......
....... .......
......... ... ....................... ...................
..... . . . ..
..... ... ..... ......... ..............
..... ..... ..... ........
.... ... ..........
..............................................................................................................................................................
. ...
........ .. ...
............... ........ ........
x
..
....... ..... ... ...... ........
................ .......... . .
.... ......... ..............
........ ...........
....... ... .............
....... ....... ... ...............
....... .
.
...........
......
.. .. ... .
.
. ...............
....... .
.
.
..... .............
..... ...........
.... .....
..
...
Figure 51 Approach (0, 0) along y 2 = x
3x2 y
Example 7.7. Prove that lim = 0.
(x,y)→(0,0) x2 + y 2
Solution. Let a positive number. We wish to find a positive number δ such that
q 3x2 y
0< x2 + y2 < δ =⇒ 2 − 0 < .
x + y2
In order to obtain
the δ that
enables the above implication to hold. We begin by estimating
3x2 y
the expression 2 − 0. Note that
x + y2
3x2 y x2 q
− 0 = 3 |y| ≤ 3|y| ≤ 3 x2 + y 2 .
2
x + y2 x2 + y 2
Thus, if we choose δ = /3, then
q 3x2 y q
0< x2 + y 2 < δ =⇒ 2 − 0 ≤ 3 x2 + y 2 < 3δ = .
x + y2
3x2 y
By the definition of limit, we have = 0. lim
(x,y)→(0,0) x2 + y 2
Remark 7.0.1. We remark that the usual limit theorems hold for limits of functions of two
variables. For example
lim (f (x, y) + g(x, y)) = lim f (x, y) + lim g(x, y).
(x,y)→(a,b) (x,y)→(a,b) (x,y)→(a,b)
c Department of Mathematics 33
Semester II(2005/06)
Remark 7.0.3. For functions of three or more variables, there are similar definition of limits
and continuity. See section 15.1 of [1]. More precisely, for functions of three variables, these
are stated as follows:
Definition 7.13. lim f (x, y, z) = L if for any > 0, there is a corresponding δ > 0
(x,y,z)→(a,b,c)
such that
q
(x, y, z) ∈ D and 0 < (x − a)2 + (y − b)2 + (z − c)2 < δ =⇒ |f (x, y, z) − L| < .
Definition 7.14. A function f is called continuous at (a, b, c) if
lim f (x, y, z) = f (a, b, c).
(x,y,z)→(a,b,c)
8. Partial Derivatives
Definition 8.1. Let f be a function of two variables. The partial derivative of f with respect
to x at (a, b) is
f (a + h, b) − f (a, b)
fx (a, b) = lim .
h→0 h
The partial derivative of f with respect to y at (a, b) is
f (a, b + h) − f (a, b)
fy (a, b) = lim .
h→0 h
c Department of Mathematics 34
Semester II(2005/06)
There are different notations for the partial derivative of a function. If z = f (x, y), we write
∂f ∂ ∂z
fx (x, y) = fx = = f (x, y) = ,
∂x ∂x ∂x
∂f ∂ ∂z
fy (x, y) = fy = = f (x, y) = .
∂y ∂y ∂y
In other words, in order to find fx , we may simply regard y as constant and differentiate f (x, y)
with respect to x. Similarly, to find fy , one can simply regard x as constant and differentiate
d d
f (x, y) with respect to y. That is fx (a, b) = f (x, b)|x=a and fy (a, b) = f (a, y)|y=b .
dx dy
Example 8.2. Let f (x, y) = x3 + x2 y 3 − 2y 2 . Then fx = 3x2 + 2xy 3 and fy = 3x2 y 2 − 4y.
Thus for example, fx (1, 1) = 5 and fy (1, 1) = −1.
Geometrically, fx (a, b) measures the rate of change of f in the direction of i at the point (a, b).
If we consider the line y = b on the xy-plane parallel to the x-axis and passing through the
point (a, b), the image of this line under f is a curve C1 on the surface z = f (x, y). Then
fx (a, b) is just the gradient of the tangent line to C1 at (a, b). Similarly, fy (a, b) is just the
derivative at (a, b) of the curve C2 traced out as the image of the line x = a under f .
· ·
...
.. ........ . ....
.. .
.......
.
........ .........
... ..... ........ ...... ... ..... .... ......
........ ..... ...... ...... ........ ..... ........... ......
f (a, b) ... ..
..... .
......
.
.....
...... f (a, b) ... ...
..... ................... ........
..... ...
......
.....
..
... ...
..
... C .
..
1 ........ ..
...... ..
.....
.
...
.. ..
... ...
..
. C 2 ..
..
.....
...
.
... ... .................................................................... ... ... ....................................................................
.................. ... .................. ...
...
... ..
z = f (x, y) ...
... ..
... .. ... ..
... .. ... ..
. .
.
.......................................................... ..
.
.......................................................... ..
.
...
. ................................. .
...
. .................................
.
. .. ................................ .
. .. ................................
..
.
.
...
. .........................
. y ..
.. .
... .........................
. y
· ·
.. ... ..
..... .. ....... .....
..
......
. . ... ..
......
. ..
.
..
....
. (a, b) ..
.. ... ..
....
. ............................
.... .....
. ....
.......... ..... the line y = b .......... (a, b)
...... ......
x x
∂f ∂f ∂f ∂w ∂w ∂w
fx , fy , fz , , , , or , , .
∂x ∂y ∂z ∂x ∂y ∂z
c Department of Mathematics 35
Semester II(2005/06)
Exercise 8.4. Let f (x, y, z) = exy ln z. Find fx , fy and fz . Show that xfx = yfy .
As in the case of function of one variable, we may also define higher order partial derivatives
of a function of several variables. Let f be a function of x and y. Then fx and fy are also
functions of x and y. Thus we may consider (fx )x , (fx )y , (fy )x and (fy )y . For convenience,
we shall simply denote them by fxx , fxy , fyx and fyy respectively. These are the second order
partial derivatives of f .
There are other notations for the higher order partial derivatives. Suppose z = f (x, y). Then
we also write:
∂ ∂f ∂2f ∂2z
fxx = = =
∂x ∂x ∂x2 ∂x2
∂ ∂f ∂2f ∂2z
fxy = = =
∂y ∂x ∂y∂x ∂y∂x
∂ ∂f ∂2f ∂2z
fyx = = =
∂x ∂y ∂x∂y ∂x∂y
∂ ∂f ∂2f ∂2z
fyy = = 2
= 2
∂y ∂y ∂y ∂y
Example 8.5. Let f (x, y) = x3 + x2 y 3 − 2y 2 . Find fxx , fxy , fyx and fyy .
Solution. First fx = 3x2 + 2xy 3 and fy = 3x2 y 2 − 4y. Thus, fxx = 6x + 2y 3 , fyy = 6x2 y − 4,
fxy = (fx )y = 6xy 2 and fyx = (fy )x = 6xy 2 .
Theorem 8.6. (Clairaut’s Theorem) If f is defined on an open disk containing the point (a, b).
If fxy and fyx are both continuous at (a, b), then fxy (a, b) = fyx (a, b).
As a result of the above theorem, if the function f (x, y) has continuous partial derivatives on
some open disk, then the order of partial differentiation can be changed without affecting the
result. Similar result holds for functions of 3 or more variables.
Proof. We assume fxy and fyx are defined in a small open disk D centered at (a, b). Let (x, y)
be a point in D with (x, y) 6= (a, b). Fix x and consider the function [f (x, y) − f (a, y)] −
[f (x, b) − f (a, b)] in y. Apply Mean Value Theorem with respect to y. We have [f (x, y) −
f (a, y)] − [f (x, b) − f (a, b)] = [fy (x, ζ1 ) − fy (a, ζ1 )](y − b) for some ζ1 between y and b. Next,
we apply Mean Value Theorem to fy (x, ζ1 ) with respect to x. We get
[f (x, y) − f (a, y)] − [f (x, b) − f (a, b)] = fyx (ζ2 , ζ1 )(x − a)(y − b),
for some ζ2 between x and a. Now we can rewrite the expression [f (x, y) − f (a, y)] − [f (x, b) −
f (a, b)] as [f (x, y) − f (x, b)] − [f (a, y) − f (a, b)]. Applying Mean Value Theorem first with
respect to x and then with respect to y, we obtain
[f (x, y) − f (a, y)] − [f (x, b) − f (a, b)]= [f (x, y) − f (x, b)] − [f (a, y) − f (a, b)]
= fxy (ζ3 , ζ4 )(y − b)(x − a),
where ζ3 is between x and a and ζ4 is between y and b. Thus we have fxy (ζ2 , ζ1 ) = fyx (ζ3 , ζ4 ).
Since fxy and fyx are continuous at (a, b), so by taking limit as (x, y) tends to (a, b), we obtain
fxy (a, b) = fyx (a, b).
c Department of Mathematics 36
Semester II(2005/06)
Exercise 8.8. Let f (x, y, z) = sin(3x + yz). Find fxxyz and fxzyx . Show that fxxyy = fxyyx .
8.1. Tangent Plane. Let f be a function of two variables. The graph of f is a surface in
R3 with equation z = f (x, y). Let P (x0 , y0 , z0 ) be a point on this surface. Thus, z0 = f (x0 , y0 ).
Assuming a tangent plane to the surface exists, we shall find its equation.
z..
..
.......... ..................................................
... ..................... .....
............. ......
.... .......... .....
.........
..
... ..
............ ............
.
... .....................
. .
......
.
...
..
..... ........ ....... ....... ...
..
·
.... .... .
... ... .... ........................ ...
.......... ... .... P . .... . ...
f (x , y ) . ....................................................... .
0 0 ... . .. ...
. .
.... ... ..... . ..... C ..
... .... . ........... .......... 2 ...
.... ... .... . ...
... ... .... ...................... ..
... ..
. ................... . .
...
... ... C 1.
...................................................................
...
... ... ...............
... .. ........... .
... .. .................
.
.
..
. z = f (x, y)
... ....... .
..
... ..
.
.
.........................................................
.
..
.
...
. ................................ ..
..
. ... . . . .......................................
.
.... ................................
.
.
.....
.....
.
. ..
..
.. ........
...........................
. y
·
..
..
....
. .
. ..
..
..
.... ... ...
...
.
.. . . .. .. .
.
............................
..... ....
..... (x , y ) 0 0 .........
.....
.......... ....
......
x
Recall that the equation of a plane passing through P (x0 , y0 , z0 ) is of the form A(x − x0 ) +
B(y − y0 ) + C(z − z0 ) = 0. Assuming the plane is not vertical, we have C is not zero. Thus
we may write the equation of the plane as
z − z0 = a(x − x0 ) + b(y − y0 ).
The tangent line to C1 at P is obtained by taking y = y0 in the above equation. That is
z − z0 = a(x − x0 ). Since fx (x0 , y0 ) is the gradient of the tangent line C1 at P , we have
a = fx (x0 , y0 ). Similarly, b = fy (x0 , y0 ). Consequently, the equation of the tangent plane to
the surface z = f (x, y) at P is
z = z0 + fx (x0 , y0 )(x − x0 ) + fy (x0 , y0 )(y − y0 )
Example 8.9. Find the equation of the tangent plane to the elliptic paraboloid z = 2x2 + y 2
at the point (1, 1, 3).
Solution. Let f (x, y) = 2x2 + y 2 . Then fx (x, y) = 4x and fy (x, y) = 2y so that fx (1, 1) = 4
and fx (1, 1) = 2. Hence, the equation of the tangent plane at (1, 1, 3) is given by z = 3 + 4(x −
1) + 2(y − 1). That is z = 4x + 2y − 3.
c Department of Mathematics 37
Semester II(2005/06)
8.2. Linear Approximation. Since the tangent plane to the surface z = f (x, y) at P
is very close to the surface at least when it is near P , we may use the function defining the
tangent plane as a linear approximation to f . Recall that the equation of the tangent plane to
the graph of f (x, y) at P (a, b, f (a, b)) is
z = f (a, b) + fx (a, b)(x − a) + fy (a, b)(y − b).
Definition 8.10. The linear function L whose graph is this tangent plane is given by
L(x, y) = f (a, b) + fx (a, b)(x − a) + fy (a, b)(y − b).
L is called the linearization of f at (a, b). The approximation
f (x, y) ≈ L(x, y) = f (a, b) + fx (a, b)(x − a) + fy (a, b)(y − b)
is called the linear approximation or tangent plane approximation of f at (a, b).
Example 8.11. Let f (x, y) = xexy . Find the linearization of f at (1, 0). Use it to approximate
f (1.1, −0.1).
Solution First we have fx (x, y) = exy + xyexy and fy (x, y) = x2 exy . Thus fx (1, 0) = 1 and
fy (1, 0) = 1. Then, L(x, y) = f (1, 0)+fx (1, 0)(x−1)+fy (1, 0)(y−0) = x+y. The corresponding
linear approximation is xexy ≈ x + y. Therefore, f (1.1, −0.1) ≈ 1.1 + (−0.1) = 1. The actual
value of f (1.1, −0.1) is 0.98542 round up to 5 decimal places.
8.3. The differential. Let z = f (x, y). As in the case of functions of one variable, we
take the differentials dx and dy to be independent variables.
z..
..
..........
... z = f (x, y)
.... .
............ ...
.. ......... ... . .........
... .......... ..... ...
... .... ................. ..
..... ......... .
..
... ..................... ... ........ ............
... ....
. . .. ... ... ... ... ....... ........ ..... 4z
... .. . ... ... ..... ......
... .... ... .. ..... .... dz
... ... ... ..... ..... . .
... ... ... ....... .... ........ .........
... ........ ............. ......... .. ....
... .
. .
...... .......... ... .
. ...... .... .
... ........
. .. ...
... ..... ....... ................. .... .... ..
...... ........ ................... .... ..
·
...... ............................... .... .
.. ..
..... .......................... ... ..
f (a, b) .................. ..
.
. ..
.... ... ..
.
...
.
....
.
.
.
...
...
...
...
...
... .
.
.
... . .
. . . .. ..
.. . .. .... .
. ..... .
............ ............... ..... ......... ... ... ...
. ............
. ........ ...... . ............................................................
.. .. y
·
..
.. .
. ..
..
....
. .
. ..
..
..... .
. ..... ... ... ... . .. ... ... .. . . .. ... .
.. ... ... . .. ..
. . (a + dx, b + dy)
.
.... .. ... ...
..
..... ...
..... .. ... ...
..... .. ... ...
.
..
...... .
. .. . ...dx ..
. ...
.
..... ...
·
.. ..... ...
..... .. .... ...
..... ....... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
.....
.............
.
x ... (a, b) dy
Note that from the definition of differentiability, if f (x, y) is differentiable at (a, b), then fx (a, b)
and fy (a, b) exist. However the converse is not necessarily true.
xy
(
x2 +y 2
if (x, y) 6= (0, 0)
Exercise 8.16. Let f (x, y) = . Show that fx (0, 0) and fy (0, 0)
0 if (x, y) = (0, 0)
exist but f is not differentiable at (0, 0).
Theorem 8.17. Suppose fx (x, y) and fy (x, y) exist in an open disk containing (a, b) and are
continuous at (a, b). Then f is differentiable at (a, b).
Proof. Write
4z = f (a + 4x, b + 4y) − f (a, b)
= [f (a + 4x, b + 4y) − f (a + 4x, b)] + [f (a + 4x, b) − f (a, b)].
By assumption, fx and fy exist near (a, b). Thus when 4x and 4y are sufficiently small, we
may apply the Mean Value Theorem to each of the above differences. Then,
4z = fy (a + 4x, b + c1 4y)4y + fx (a + c2 4x, b)4x,
where c1 , c2 ∈ (0, 1). Moreover, by the assumption that both fx and fy are continuous at (a, b),
we have
fy (a + 4x, b + c1 4y) = fy (a, b) + 2 ,
fx (a + c2 4x, b) = fx (a, b) + 1 ,
where lim 2 = 0 and lim 1 = 0. Thus, 4z = fx (a, b)4x + fy (x, y)4y +
(4x,4y)→(0,0) (4x,4y)→(0,0)
1 4x + 2 4y, with both the limit of 1 and 2 tend to 0 as 4x and 4y tend to 0. Therefore
f is differentiable at (a, b).
That fx and fy are continuous in an open disk containing the point is only a sufficient condition
for the differentiability of f at the point. It is by no means necessary. See the following exercise.
c Department of Mathematics 39
Semester II(2005/06)
(
1
(x2 + y 2 ) sin x2 +y 2 if (x, y) 6= (0, 0)
Exercise 8.18. Let f (x, y) = . Show that f is differ-
0 if (x, y) = (0, 0)
entiable at (0, 0) but fx and fy are not continuous at (0, 0).
Theorem 8.19. (The chain rule, case 1) Suppose z = f (x, y) is a differentiable function of x
and y, where x = g(t), y = h(t) are both differentiable functions of t. Then z is a differentiable
functions of t and
dz ∂f dx ∂f dy
= + .
dt ∂x dt ∂y dt
Proof. Exercise. See also [1] on page 968.
dz
Example 8.20. Let z = x2 y + 3xy 4 , where x = sin 2t, y = cos t. Find .
dt
Solution.
dz ∂z dx ∂z dy
= + = (2xy + 3y 4 )(2 cos 2t) + (x2 + 12xy 3 )(− sin t).
dt ∂x dt ∂y dt
Theorem 8.21. (The chain rule, case 2) Suppose z = f (x, y) is a differentiable function of x
and y, where x = g(s, t), y = h(s, t) are both differentiable functions of s and t. Then z is a
differentiable functions of s and t and
∂z ∂f ∂x ∂f ∂y
= + ,
∂s ∂x ∂s ∂y ∂s
∂z ∂f ∂x ∂f ∂y
= + .
∂t ∂x ∂t ∂y ∂t
∂z ∂z
Example 8.22. Let z = ex sin y, where x = st2 , y = s2 t. Find and .
∂s ∂t
Solution.
∂z ∂z ∂x ∂z ∂y
= + = (ex sin y)t2 + (ex cos y)(2st).
∂s ∂x ∂s ∂y ∂s
∂z ∂z ∂x ∂z ∂y
= + = (ex sin y)(2st) + (ex cos y)(s2 ).
∂t ∂x ∂t ∂y ∂t
Exercise 8.23. Suppose z = f (x, y) has continuous 2nd order partial derivatives and x =
∂z ∂2z
r2 + s2 , y = 2rs. Find and .
∂r ∂r2
8.4. Implicit Differentiation. Suppose F (x, y) = 0 defines y implicitly as a function of
x. That is y = f (x). Then F (x, f (x)) = 0. Now we use the chain rule(case 1) to differentiate
F with respect to x. Thus
dx dy
Fx + Fy = 0.
dx dx
Therefore,
dy Fx
=− .
dx Fy
dy
Example 8.24. Find if x3 + y 3 = 6xy.
dx
c Department of Mathematics 40
Semester II(2005/06)
Solution. Let F (x, y) = x3 + y 3 − 6xy. The given equation is simply F (x, y) = 0. Therefore,
dy Fx 3x2 − 6y
=− =− 2 .
dx Fy 3y − 6x
∂x ∂y ∂z
Fx + Fy + Fz = 0.
∂x ∂x ∂x
∂x ∂y ∂z
Note that = 1 and = 0. Thus Fx + Fz = 0. Hence,
∂x ∂x ∂x
∂z Fx
=− .
∂x Fz
Similarly,
∂z Fy
=− .
∂y Fz
Exercise 8.25. Three ants A, B and C crawl along the positive x, y and z axes respectively.
A and B are crawling at a constant speed of 1 cm/s, C is crawling at a constant speed of 3 cm/s
and they are all travelling away from the origin. Find the rate of change of the area of triangle
ABC when A is 2 cm away from the origin while B and C are p 1 cm away from the origin. [The
1
area of the triangle A(x, 0, 0)B(0, y, 0)C(0, 0, z) is given by 2 x2 y 2 + y 2 z 2 + z 2 x2 .]
[Answer : 4 cm2 /s.]
z..
...
.......
..
...
... .. ... ... .
.. .. ... ... z = f (x, y)
... .. ... ... ... .................................
. .................
... ............ .... ... ... ... ... ...........
... .. .......... ... ... ... ..
... ................ ... ... ... ................
...... ... .. .......
... ..
..
.... ..
.......
... .. ......
·
..
.. ......
... ..... ............................. P (x , y , z ) 0 0 ....
·
......
... ..
....... ..
. . .. . .. ............................ 0
. . . . . . . .. ......
.... .
................................................................................................................................... . .....
... .
.. .
...................................................................................................................................... .................. . . . . .....
... ... ..
·
.
...
..
..................................................................................................................................................... .. . .......
... .. .......
... ..... .....................................................................................................................................................................
. . .
.........
... ... ......................................................................................................................................................................... .......
.
. ... . . . . . . ... ...
... ..
... ................. ................................................................................................................................................................................................................ ......
........... ....................................................................................................................................................................... .....
... ... ............ ............................................................................................................................... .........
... .......... ........................................................................................................................................................................ ......
... .
............................................................................................................................................................................ .
... ........................................................................................................................................................................
... ......................................................................................................................................................................
... .....................................................................................................................................................................
... .....................................................................................................................................................................
... .....................................................................................................................................................................
... .....................................................................................................................................................................
.... .....................................................................................................................................................................
... ...
...
...
...
...
...
................................
..
....
. ... .
.... ..... .
...
. .
..................................................................................
..........................................................................................................................................................................................................................
..
....
. ............................................................................................................................................
..
....
. .........................................................................................................................................
.
..... ..............................................................................................................................................................
....................................
.
.
...
. ............................................................................................................................................................................................................. .......... y
·
.
.... ...........................................................................................................
..
.. .. ... .........................................................
.
....
. . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
..
.. .
... 0
.
. .. . . . .. h
.......................................................................
. ... .. . . . .
.... ........... . . . . . . . . . . .
..... P (x , y ) .
0 0 .........................................................................................................................
.....
..... .. ..
.
.....
....
.
û = ha, bi
.
....
..
..
..
.........
.......
x
Note that Di f (x0 , y0 ) = fx (x0 , y0 ) and Dj f (x0 , y0 ) = fy (x0 , y0 ), where i and j are the standard
basis vectors in R2 .
Theorem 8.27. Let f be a differentiable function of x and y. Then f has a directional
derivative in the direction of any unit vector u = ha, bi and
Du f (x, y) = fx (x, y)a + fy (x, y)b = hfx (x, y), fy (x, y)i · u.
Proof. Consider g(h) = f (x0 + ha, y0 + hb). Clearly g 0 (0) = Du f (x0 , y0 ). By the chain
rule, g 0 (h) = ∂f dx ∂f dy 0
∂x dh + ∂y dh . At h = 0, we have g (0) = fx (x0 , y0 )a + fy (x0 , y0 )b. Hence,
Du f (x0 , y0 ) = fx (x0 , y0 )a + fy (x0 , y0 )b.
Example 8.28. Let f (x, y) = x3 − 3xy + 4y 2 . Find Du f (1, 2), where u is the unit vector
making an angle of π6 with the positive x-axis.
Solution.
y.
..
......... ...
...........
... ....... .
.... .......
....... π π
·
..
..
... ......
.
.......
. û = hcos 6
, sin 6
i
... ................. ... ... ... ... ..
...
(1, 2)
...
...
...
...
...
................................................................................................. x
First, fx = 3x2 − 3y, fy = −3x + 8y. Thus √ fx (1, 2) = −3 and fy (1, 2) = 13. Therefore,
Du f (1, 2) = h−3, 13i · hcos π6 , sin π6 i = (13 − 3 3)/2.
Definition 8.29. Let f be a differentiable function of x and y. The gradient of f is the vector
function
∂f ∂f
∇f (x, y) = hfx (x, y), fy (x, y)i = i+ j.
∂x ∂y
c Department of Mathematics 42
Semester II(2005/06)
Thus we have the following formula for the directional derivative in terms of the gradient of f .
Example 8.30. Let f (x, y) = x2 y 3 − 4y. Find the directional derivative of f at (2, −1) in the
direction 3i + 4j.
Solution. The unit vector along 3i + 4j is u = 35 i + 45 j. The gradient of f is ∇f = h2xy 3 , 3x2 y 2 −
4i. Thus ∇f (2, −1) = h−4, 8i. Consequently, Du f (2, −1) = ∇f (2, −1)·u = h−4, 8i·h 53 , 45 i = 4.
Definition 8.31. Let f be a function of x, y and z. The directional derivative of f at (x0 , y0 , z0 )
in the direction of a unit vector u = ha, b, ci in R3 is
f (x0 + ha, y0 + hb, z0 + hc) − f (x0 , y0 , z0 )
Du f (x0 , y0 , z0 ) = lim
h→0 h
if this limit exists.
Similarly, the gradient of a differentiable function f is defined to be
∂f ∂f ∂f
∇f = hfx , fy , fz i = i+ j+ k.
∂x ∂y ∂z
The formula Du f = ∇f · u is also valid for any function f of more than 2 variables.
Exercise 8.32. Let f (x, y, z) = xyz 2 . Let u be the unit vector √1 i + √2 i − √1 i. Find
6 6 6
Du f (1, 1, 1).
Theorem 8.33. Let f be a differentiable function of 2 or 3 variables. Let P be a point in the
domain of f . The maximum value of Du f (P ) is |∇f (P )| and it occurs where u has the same
direction as the gradient vector ∇f (P ).
Proof. First Du f (P ) = ∇f (P ) · u = |∇f (P )||u| cos θ = |∇f (P )| cos θ, where θ is the angle
between ∇f (P ) and u. Therefore, Du f (P ) attains its maximum value |∇f (P )| when θ = 0,
i.e. when u has the same direction as the gradient vector ∇f (P ).
Exercise 8.34. Let f (x, y, z) = xey , P = (2, 0) and Q = ( 12 , 2).
−→
(b) Find the rate of change of f at P in the direction P Q. In other words, find Du f (P ).
(b) In which direction does f have the maximum rate of change? and what is this maximum
rate of change?
y. .. y
.... ... Contour of f (x, y) = xe
....... ....
.... . ....
... .... ....
... . ... ....
... .... .... ....
......... ... ......
... ...... .... .......
... ... ... .... ........
....... .... ..... ......
.......
... ....... ... ...... .......
... ...... .... ........ ........
.........
... ... .... .... ......
...... ..........
... ... ... ...... ....... ...........
............
... ... .... ....... ........
........ ..............
... .... .... ......... ∇f ................... ................
..................
... ... .... ......
.......
...........
. .............. ......
32
... .... ....... ....... .
. .............
... ... ...... ........ ... ...............
... ... ...... ........... ......
...
... ......
....... .. ..................... 16
...
.....
..... ....... ..... .............
... .....
......
..........
. ... 8
·
... ...... .. ..............
.
4
..........................................................................................................................................................................
2 ... x
P (2, 0)
8.6. Tangent Planes to Level surfaces. Let S be a surface with equation F (x, y, z) =
k, where k is a constant. That is S is a level surface of F . Let P (x0 , y0 , z0 ) be a point in S.
Let’s find the equation of the tangent plane to S at P .
............................................
.................. ............
.......... ..........
........ ∇F (P ) ........
....... .
......... ........
.
..
....... .. . .......
.... .................................................................................. .......
·
.
... ...
. ......
...
. .... .
. .... ......
.... .
.....
.... .
.....
.
. P . .. .
.
............................
.
. ..
...
.
..
.....
.. ......
.....
..
..... ..
. .... ... ... ... .. ... ... .....
.
.. .
. .....
.... ...... ...... ... .
......... .....
.
..
.
...
. .
.................................................................................... .
.....
........
.
... .. . ... ...
.......
...
. ..... . ... ...
.. . ...... .
........... ........
.... ................. .....
... ............... ...........
.....
S .. ............ ......... .....
... ......... ........
.....
......... .......
...... ....
.
...... ...
.
..... ......
..... ...
.....
F (x, y, z) = k
Take any curve r(t) = hx(t), y(t), z(t)i on the surface S such that r(0) = (x0 , y0 , z0 ). It’s
tangent vector r0 (0) shall lie on the tangent plane to S at P . Now if we use the chain rule to
differentiate F (x(t), y(t), z(t)) = k with respect to t, we have
dx dy dz
Fx+ Fy + Fz = 0.
dt dt dt
In other words, ∇F · r0 (t) = 0. At t = 0, we have ∇F (P ) · r0 (0) = 0. Therefore, ∇F (P ) is
perpendicular to the tangent plane.
Example 8.35. Find the equation of the tangent plane and the normal line at the point
(−2, 1, −3) to the ellipsoid
x2 z2
+ y2 + = 3.
4 9
2 2
Solution. Let F (x, y, z) = x4 + y 2 + z9 . Then ∇F = hx/2, 2y, 2z/9i. Thus ∇F (−2, 1, −3) =
h−1, 2, −2/3i. Therefore, the equation of the tangent plane is: hx+2, y−1, z+3i·h−1, 2, −2/3i =
0. That is 3x − 6y + 2z + 18 = 0. Also the equation of the normal line in symmetric form is
given by
x+2 y−1 z+3
= = .
−1 2 − 23
In the special case in which the surface S is the graph of a function z = f (x, y), S can be
regarded as the level surface
F (x, y, z) ≡ f (x, y) − z = 0.
In other words, the graph of z = f (x, y) is simply the level surface of F (x, y, z) at level 0.
In this case, ∇F = hfx , fy , −1i is a normal vector to the tangent plane of S at (x, y, f (x, y)).
Therefore, if we consider a point P (x0 , y0 , f (x0 , y0 )) on the graph of z = f (x, y). The equation
of the tangent plane is given by
hx − x0 , y − y0 , z − f (x0 , y0 )i · hfx , fy , −1i = 0.
That is
z = f (x0 , y0 ) + fx (x0 , y0 )(x − x0 ) + fy (x0 , y0 )(y − y0 ).
This is the same formula obtained in 8.1.
c Department of Mathematics 44
Semester II(2005/06)
.
.
..
....................
............
....... ... ·
.........................................
..................
...
......................................................
..
.
.
..
...
..
................................................................................. .....
...........
......
..
....
........
. ....
.... ..
...
.... ... ... . ..
...... .
.
...
... ... ...
.... ... .. ...
.
..
... .. ... ... .... ... .. ....
... . . .. .. ... ... .... . . ...
.
..
..................... .. .. . . . . .... .... ... . . .. .. . . .. . . . .. .. ....
. . . . ... ............... .....
. . .
...
...
........ .... ..... ... ... ... ... ... ... ... .
.....
.
.... .
... .... ..... ... .... .
..
. ... ...... .
.
.... ...... ...
................................................................
·
.... ... ..
........ ... .... ....... .... . ........
.............. ........ ...............
.....................................
............................................................................
Theorem 9.2. If f has a local maximum or a local minimum at (a, b) and fx (a, b) and fy (a, b)
exist, then fx (a, b) = 0 and fy (a, b) = 0. That is ∇f (a, b) = 0.
Definition 9.3. A point (a, b) is called a critical point of f if fx (a, b) = 0 and fy (a, b) = 0,
or if one of these partial derivatives does not exist.
.
..
· ....
.........
... .. ...
... .. ....
... .... .........
.... ........ ... ... .....................
.................... .... ..... .............
........... ..... . ......
....... .......... ..
.
....... ...........
..... ...
........................... ....
.
.......... ..
.........
... .
..... ..
.....
....... ..................... .......
.....
........... ..........
........................
........................................
Note that if f has a local minimum or a local maximum at (a, b), then (a, b) is a critical point
of f . However not all critical points of a function give rise to local maximum or local minimum.
In other words, at a critical point, a function could have a local maximum, or a local minimum
or neither.
Example 9.4. Let f (x, y) = x2 + y 2 − 2x − 6y + 14. Find the local maxima and local minima
of f .
Example 9.5. Find the local extrema (i.e. local maximum or local minimum) of f (x, y) =
y 2 − x2 .
Solution. First fx = −2x and fy = 2y. Therefore, the only critical point is (0, 0). However, f
has neither a maximum nor a minimum at (0, 0). To see this, consider the function f along
y = 0, f (x, 0) = −x2 < 0 for x 6= 0. So f has a local maximum along y = 0. On the other
hand, if we consider f along x = 0, we have f (0, y) = y 2 > 0 for all y 6= 0. Thus f has a local
minimum along x = 0. Therefore f has neither a maximum nor a minimum at (0, 0). Such a
point is called a saddle point.
c Department of Mathematics 45
Semester II(2005/06)
..........
z . ..
... ........ ...
.. ..... ........ .......
.. ... .. .......
...... ...
.. ......... ... ... ...
............ ... ... ... ..
............. ..
.............. . .
...
... . . . . .
....
...
.
.
....
.
...
. .
............... .
. ..
z = y 2 − x2
.
.
... .. .. ... .. ... ...
... .. .. .. .. ..
... .. ... .. .... ....
... .. ... .
. ...
. ... ...
... .. ...
... .... ... ... ...
... .. .... .... ...
·
... ...
...
... ..
. ..... .. ....... ... ...
.......... ........
. ...........................................................................................................
....
... .
..
.. ...
.. .
. ..
.
..... ... ...
... .....
. y
.. ........ ..... . .. ... ...
... ... ... ... ......... .. .. .....
.... ........... ... .... ... ... ........ ... ... ... ... ... ... ... ... .........
.. . ...
.. ..... ..
... ..... .. ...
..... .................................................................... ..
... ....................... ............ ....
......................................... ...........
. . .
x
Figure 58 A saddle point
Definition 9.6. f is said to have a saddle point at (a, b) if there is a disk centered at (a, b)
such that f assumes its maximum value on one diameter of the disk only at (a, b), and assume
its minimum value on another diameter of the disk only at (a, b).
In other word, f has a saddle point at (a, b) if there are some directions along which f has a
local maximum at (a, b) and some directions along which f has a local minimum at (a, b).
Exercise 9.7. Let f (x, y) = x3 − 3xy 2 . Draw the contours of f near (0, 0). The point (0, 0)
is called a monkey saddle.
Theorem 9.8. (The Second Derivative Test) Suppose fxx , fxy , fyx and fyy are continuous on
a disk with centre (a, b) and suppose fx (a, b) = 0, fy (a, b) = 0. Let
D = fxx (a, b)fyy (a, b) − fxy (a, b)2 .
(a) If D > 0 and fxx (a, b) > 0, then f has a local minimum at (a, b).
(b) If D > 0 and fxx (a, b) < 0, then f has a local maximum at (a, b).
(c) If D < 0, then f has a saddle point at (a, b).
Note that if D = 0, then no conclusion can be drawn from it. The point can be a local
maximum, a local minimum, a saddle point or neither of these. Interested readers are invited
to come up with examples for all these cases.
Example 9.9. Let f (h, k) = Ah2 + 2Bhk + Ck 2 . Suppose A 6= 0. We may write f in the form
h i
f (h, k) = A (h + Bk/A)2 + (AC − B 2 )k 2 /A2 .
Assume AC − B 2 > 0 and A > 0. Then f (h, k) ≥ 0 and f (h, k) = 0 only if (h, k) = (0, 0).
Thus f has a minimum at (0, 0). Similarly, if AC − B 2 > 0 and A < 0, then f has a maximum
at (0, 0). Next suppose AC − B 2 < 0 and A > 0, then f has a minimum along the line k = 0
at (0, 0) but a maximum along the line Ah + Bk = 0, thus giving a saddle point at (0, 0).
continuous functions, so there is a disk B with center (a, b) and radius δ > 0 such that
fxx (x, y) > 0 and D(x, y) > 0 whenever (x, y) is in B. Therefore, by looking at the above
equation, we see that Du2 f (x, y) > 0 whenever (x, y) is in B. This means that if C is the curve
obtained by intersecting the graph of f with the vertical plane through the point P (a, b, f (a, b))
in the direction u, then C is concave upward on an interval of length 2δ. This is true in the
direction of every vector u, so if we restrict (x, y) to lie in B, then the graph of f lies above its
horizontal tangent plane at P . Thus f (x, y) ≥ f (a, b) for all (x, y) in B. This shows that f has
a local minimum at (a, b). The proof of (b) is similar. The proof of (c) is left as an exercise.
Example 9.10. Find the local maximum, local minimum and saddle points (if any) of f (x, y) =
x4 + y 4 − 4xy + 1.
Solution. First, fx = 4x3 − 4y and fy = 4y 3 − 4x. Now we proceed to solve 4x3 − 4y = 0 and
4y 3 − 4x = 0 for the critical points. The two equations are equivalent to y = x3 and x = y 3 .
Substituting one into the other, we obtain x9 −x = 0. That is x(x+1)(x−1)(x2 +1)(x4 +1) = 0.
Thus the real solutions are x = 0, −1, 1. Therefore, the critical points are (0, 0), (−1, −1) and
(1, 1). To apply the second derivative test, we compute the second order partial derivatives.
fxx = 12x2 , fyy = 12y 2 , fxy = −4. Thus D(x, y) = fxx fyy − fxy
2 = 144x2 y 2 − 16.
At (0, 0), D(0, 0) = −16 < 0. Hence, f has a saddle point at (0, 0). At (−1, −1), D(−1, −1) =
128 > 0 and fxx (−1, −1) = 12 > 0. Hence f has a local minimum at (−1, −1). At (1, 1),
D(1, 1) = 128 > 0 and fxx (1, 1) = 12 > 0. Hence f has a local minimum at (1, 1).
Definition 9.11. A bounded set in R2 is one that is contained in some disk. A closed set in
R2 is one that contains all its boundary points.
...................... .
........
.....................
...................... .... ...
.... . . . . . . ..... ... . ..
................................. ................
.
.. . . . .
.......................................................... ......................
... . . . . . . . . ... ................
.......................
...................... ..................
....................... ... . . . . . ...
....................... .....................
....... . . . ...... ..............................
......................................................... .......................... . .
.........................................................
................ ................................................
........ . . ....... ... . . . . . . . ...
..... ................ ....................
.................... ..................
.......................... .... . . . . ....
................................................. .................................... .... ..............
.............
... . . . . . . . ... ... ................... .... . ....
.................... ... . . . . . . . ... ... ......
..... . . . . .... ......................
...... . . . ...... ...... . . . . ....... .... ...
..................... .......................... .......
.................. ..
Figure 59 Sets in R2
Theorem 9.12. (Extreme Value Theorem) If f is continuous on a closed, bounded set D in R2 ,
then f attains an absolute maximum value f (x1 , y1 ) and an absolute minimum value f (x2 , y2 )
at some points (x1 , y1 ) and (x2 , y2 ) in D.
The following is a procedure to find the absolute maximum and the absolute minimum value
of a function defined on a closed and bounded set.
c Department of Mathematics 47
Semester II(2005/06)
L ·
.....................................
.....................................
4 ....................................................... 2
(1, 1)
.......................................................
.....................................
.....................................
L
.
...................................................................................................................
(0, 0) L1 x
(3, 0)
Figure 60
Along L1 : y = 0. That is f (x, 0) = x2 for x ∈ (0, 3) which is increasing, thus giving no critical
point along L1 .
Along L2 : x = 3. That is f (3, y) = 9 − 6y + 2y = 9 − 4y for y ∈ (0, 2) which is decreasing,
thus giving no critical point along L2 .
Along L3 : y = 2. That is f (x, 2) = (x − 2)2 for x ∈ (0, 3). It has a critical point (a local
minimum) at x = 2. That is at the point (2, 2).
Along L4 : x = 0. That is f (0, y) = 2y for y ∈ (0, 2) which is increasing, thus giving no critical
point along L4 .
Now let’s compute the values of f at all the critical points (including the four vertices of the
rectangle).
f (1, 1) = 1, f (2, 2) = 0, f (0, 0) = 0, f (3, 0) = 9, f (3, 2) = 1, f (0, 2) = 4. Thus the absolute
maximum value of f is 9 and the absolute minimum value is 0.
Exercise 9.14. Find the maximum and minimum values of f (x, y) = x2 + y 2 − x − y + 1 on
the triangular region R with vertices (0, 0), (2, 0), (0, 2). See Figure 61.
[Answer : Maximum value = 3, Minimum value = 1/2.]
y.
....
........
....
...
(0, 2) ...
..........
............
..............
................
L
.................. 2
.. . . . . ...
L 3 ................................
.........................
..........................
R
............................
. ..
................................................................................................................................
(0, 0) L1
x
(2, 0)
Figure 61
z..
.
.........
...
....
..
...
... ... ..
. .
. .. ...
. ... ..........................
... .. . .......... ... ........
. ...... ...
......
... z = f (x, y)
.... ... .
....... . ...
.. .........
.. ..............
.
..
. .. .. .. .
............
............................................. . .
... ..
.
...
... ... ...
...... ... ...
.... ... ... ........
... .... .
. .. .
....... .....
.
... .... . ... .
...
.. .... .... ..........
.... .... .. ...... ... ...
........ .... ... . .
.. ........................................ ..................... .. ...
... .
... .... .. .
... ... ......
... ....
......... ... ............ ... ...........
.. . . .
....... ....... .. ....... ...........
..... ...... ... ............. ...
... .
.................................................................................................................................................................................................
...
...... .
..
...
.
.....
.. y
........ ......... ..
....
.... ..........
.................................
..... g(x, y) = 0
.....
.. .......
....
.....
..........
....
x
Figure 62
If we confine the point (x, y) to lie on the curve g(x, y) = 0 on the xy-plane, its image under
f gives a curve on the graph of z = f (x, y). We are looking for the highest and lowest points
of this curve. Suppose the extreme value of f (x, y) subject to the constraint g(x, y) = 0 is k
and is attained at the point (x0 , y0 ). By examining at the contour of f , we see that at the
extreme point, the curve g(x, y) = 0 must touch the level curve f (x, y) = k, because if the curve
g(x, y) = 0 cuts across the level curve f (x, y) = k, one can still move the point along g(x, y) = 0
so as to increase or decrease the value of f . In other words, the gradients of f and g must be
parallel at the extreme point (x0 , y0 ). Consequently, we must have ∇f (x0 , y0 ) = λ∇g(x0 , y0 )
if ∇g(x0 , y0 ) 6= 0.
y. ∇f
.
.. ..
.......... ...........
... ..
.... ..
.
.. ... ..... ... ...
... .... ... ... ... ...
... ... ... ... ..
... ... ... ... ... .....
... .... ...
...
...
...
..... .........
..... ..... ∇g . .
... ...
... ... ...
...
...
.....
..........
............
.. ...........................
...... f (x, y) = 12
... ... . ...
. ....
...
... ...
.
......
... ......
.
. .... ...... ....... .... .
y 0 ....... ... .....
..... • ..
............................................................................ f (x, y) = 11
...
...
..... ...... ............. .......
.
... .
.....
..... . . ................. ......
.....
... ....... ..... ........ ........
.......... ....
.......... . ..
..... ...
... .......
..... .
.
..............
..... . . . ... ........................................
.....
f (x, y) = 10
g(x, y) = 0 ..... .. .....
.... . ....... .......
. .......... ......
....... ......
.. ............. ......................... ...........
.................. ..
... ....... ......
....... f (x, y) = 9 . . .............
... ........
... .........
...........
... .............. f (x, y) = 8
.....
..
...........................................................................................................................................................
...
x x
... 0
Example 10.1. Find the extreme values of f (x, y) = x2 + 2y 2 on the circle x2 + y 2 = 1. See
figure 62.
Solution Since the circle is a closed and bounded set and f is a continuous function, there
is always an absolute maximum and an absolute minimum of f over it. They are among
the extreme values of f defined over the circle. To find them, first ∇f (x, y) = h2x, 4yi and
∇g(x, y) = h2x, 2yi. Thus ∇f (x, y) = λ∇g(x, y) is equivalent to
2x = λ2x,
4y = λ2y.
Together with the constraint x2 + y 2 = 1, we need to solve the following system of equations
in x, y and λ:
2x(λ − 1) = 0,
2y(λ − 2) = 0,
2
x + y2 = 1.
The first equation gives either x = 0 or λ = 1.
If x = 0, then the constraint equation gives y = ±1. Thus we have the solutions (0, −1), (0, 1).
If λ = 1, then the second equation gives y = 0. Thus, by the constraint equation, we have the
solutions (−1, 0), (1, 0).
Consequently, we have four solutions (0, 1), (0, −1), (−1, 0), (1, 0). Now f (0, −1) = 2, f (0, −1) =
2, f (−1, 0) = 1, f (1, 0) = 1. Therefore, the absolute maximum value is 2 and the absolute min-
imum value is 1.
Exercise 10.2. Find the rectangular box with the largest volume that can be inscribed in the
ellipsoid
x2 y 2 z 2
+ 2 + 2 = 1.
a2 b c
8abc
[Answer: 3 ]
32
z..
..
..........
..
....
..
...
...
.. ..
.... .
.. ..
...... ..................
................. ... ........... ...................... .......................................
.
·
. .
.... ..
. . . ...
...
. .
.............. ...
.... .. .. .. ...... ........
....... ............................................................................................. .... ..........
. . .
.. . . . . . .....
.... ... . ....... ... ... ...... ... ... ......... ... ... .................. .. ..
... ... ... ... .............. ... ... ... .
. .. .
. . . ... ... .......
... .... .. . .. .. .. ..
........ ..... ................................................................................................................
.
............... .. . .
. .
... .......................................... ........
. .. ..
. .
.. .
... ... .....................
. y
..... . .... ... ... ... .................................................................................................................
. . .
..
..... ....... .. .. .....
...... .................................................................................................. .
......
........ .. . . .
..
.......... ....... ...
..... .. .........
...
.............. .. .............
..... ...................................................
.....
.
............
.
x
Figure 64
Exercise 10.3. Find the point on the sphere x2 + y 2 + z 2 = 4 that are closest to and farthest
from the point (3, 1, −1). (Consider the line passing through (3, 1, −1) and the centre of the
sphere, it intersects the sphere diametrically at two points.)
√ √
[Answer: Min = 11 − 2 at ( √611 , √211 , √−2
11
), Max = 11 + 2 at ( √−6
11
, √−2
11
, √211 ) ]
multipliers as in example 10.1, we obtain 4 critical points (−1, 0), (1, 0), (0, −1), (0, 1). Now,
f (0, 0) = 0, f (0, −1) = 2, f (0, −1) = 2, f (−1, 0) = 1, f (1, 0) = 1. Therefore, the absolute
maximum value is 2 and the absolute minimum value is 0.
Remark 10.5. If we define a new function L(x, y; λ) = f (x, y)−λg(x, y). Then ∂L∂x = fx −λgx ,
∂L ∂L
∂y = fy − λgy and ∂λ = g. Therefore, the critical points of L correspond to the extreme points
of the original problem. The L is a called a Lagrangian.
The method of Lagrange Multipliers can be applied to the case of more than one constraints.
Consider the problem of maximizing or minimizing f (x, y, z) subject to the constraints g(x, y, z)
= 0 and h(x, y, z) = 0. If f attains an extreme value at (x0 , y0 , z0 ) , then
∇f (x0 , y0 , z0 ) = λ∇g(x0 , y0 , z0 ) + µ∇h(x0 , y0 , z0 ).
(For this linear combination to be valid, we need to assume ∇g(x0 , y0 , z0 ) 6= 0 and ∇h(x0 , y0 , z0 )
6= 0 and that they are not parallel.)
Solving this vector equation and the constraint equations give all the possible extreme points.
Equating components, these equations are equivalent to the following system:
fx (x, y, z) = λgx (x, y, z) + µhx (x, y, z)
y (x, y, z) = λgy (x, y, z) + µhy (x, y, z)
f
f (x, y, z) = λgz (x, y, z) + µhz (x, y, z)
z
g(x, y, z) = 0
h(x, y, z) = 0
In this case, we have two multipliers.
Example 10.6. Find the maximum value of f (x, y, z) = x+2y +3z on the curve of intersection
of the plane x − y + z = 1 and the cylinder x2 + y 2 = 1.
Solution. We wish to maximize f (x, y, z) = x + 2y + 3z subject to the constraints g(x, y, z) =
x − y + z − 1 and h(x, y, z) = x2 + y 2 − 1.
...............................................
..................................... ........
...................... ........... ......
......................... ........... ..
...
..
..... . ...................... ..........
.....
..... .
.. ....... .................................................. ...
..... .
...... .......... ... ..
..... .
. ....... .. ...
..... .... ...... ..... ...
..... . ...... ...
. . ...... ... ...
.....
..... .......
. ....... ... . ...
x−y+z =1 ..... ... ... .
..
..... ... .. ...
..... .. ... ..... ... .
..... .. ... ..... ... ....
....... .. ..... ... ..
...... ... .
... . . .... ......
.... ......... ........ ...
... .... ....
... ....... ....... ... .. .....
... ..... ..... ... ... ......
... ..... ..... ... .. .....
..... ...... .....
... ..... .....
..... ......
.....
.... .....
... ..... ....... .. .....
... ..... ....... .....
... ..... .......... ........ .....
.....
..... ........ ..
... ... ... ... ... ... ........... ... ... ... ... .. . ...
...... ... ... ... ......
2 2 ....
.....
.....
......
.
.
...........
x +y =1 ........
..................
.....
.. .... .... . .
............
..................................................... .......
..... ............
......
Figure 65
First we have ∇f = h1, 2, 3i, ∇g = h1, −1, 1i and ∇h = h2x, 2y, 0i. Thus we need to solve the
system of equations: ∇f = λ∇g + µ∇h, x − y + z = 1, x2 + y 2 = 1. That is
1 = λ + 2xµ (1)
−λ + 2yµ
2 = (2)
3 = λ+0 (3)
x−y+z = 1 (4)
x2 + y 2
= 1 (5)
c Department of Mathematics 51
Semester II(2005/06)
From (3), λ = 3. Substituting this into (1) and (2), we get x = − µ1 and y = 5
2µ . Note that
µ 6= 0 by (2) and (3).
From (4), we have z = 1 − x + y = 1 + µ1 + 2µ
5 7
= 1 + 2µ . (6)
√
Using (5), we have (− µ1 )2 + ( 2µ ) = 1. From this, we can solve for µ, giving µ = ± 229 .
5 2
Thus x = − √229 or x = √229 . The corresponding values of y are √529 , − √529 . Using (6), the
corresponding values of z are 1 + √7 , 1 − √7 .
29 29
Therefore, the two possible extreme values are at P1 = (− √229 , √529 , 1 + √729 ) and P2 =
√ √
( √229 , − √529 , 1 − √729 ). As f (P1 ) = 3 + 29 and f (P2 ) = 3 − 29, the maximum value is
√ √
3 + 29 and the minimum value is 3 − 29.
z.. z = f (x, y)
...
.......... ..........
....... .......
... .......
.... . ...
......... .
....
.
........ .. .
.. .
.. ..........
... ............. .. ..
... ..... .. ..
.. ..
...
...
...
...
.......
.
.. .
.. ..
.. ..
.. ..
· .........
...................................
........................... ..............
..................................... ..
.
.
... ...
.
..
..
... .. .. . .. .
......................................................
. . . . . . .
. .
... ..... .
..
. .
. ..
... .
. ............... ...
.
. .
. ..
...
.. .
. .....................
.
. .
. ..
...
.. .
. .....................
.
. .
. ..
...
.. .
. .....................
.
. .
. ..
.. . ..................... . ..
... . c . .
. ..... . . d
........................................................................................................................................................................................................................
.
...
.
.
...
.
... .. .
.
.. . .
. .............. ..
.
.
.
.
.
.
.
.
. .
.
.
.
.
.
.
.
.
.
.. ..... ... y
..... .. .. ........... .. . ..
............. .. ...........................................................................................
a . . . . . . .
... .........
·
.
...
.....
... .
....................................... .. .........
. ..... .. .....
..... .. ...............
b ........... .........................................................................
..
.....
....
.....
.....
..... R = [a, b] × [c, d]
. ..
......
.
......
...... 2
x = {(x, y) ∈ R | a ≤ x ≤ b, c ≤ y ≤ d}
Figure 66
To do so, we first subdivide R into mn small rectangles Rij , each having area 4A, where
i = 1, · · · , m and j = 1, · · · , n. For each pair (i, j), pick an arbitrary point (x∗ij , yij
∗ ) inside
∗ ∗
Rij . We then use the value f (xij , yij ) as the height of a rectangular solid erected over Rij .
Thus its volume is f (x∗ij , yij
∗ )4A. The sum of the volume of all these small rectangular solids
approximates the volume of the solid under the graph of z = f (x, y) over R. This sum
m X
n
f (x∗ij , yij
∗
X
)4A
i=1 j=1
is called a Riemann sum of f . We define the double integral of f over R as the limit of the
Riemann sum as m and n tend to infinity. In other words,
Definition 11.1. The double integral of f over R is
ZZ m X
n
f (x∗ij , yij
∗
X
f (x, y) dA = lim )4A
R m,n→∞
i=1 j=1
If f (x, y) ≥ 0, then the volume V of the solid lies above the rectangle R and below the surface
z = f (x, y) is ZZ
V = f (x, y) dA.
R
11.2. Iterated Integrals. Let f (x, y) be a function defined on R = [a, b] × [c, d]. We
Z d
write f (x, y)dy to mean that x is regarded as a constant and f (x, y) is integrated with
c Z d
respect to y from y = c to y = d. Therefore, f (x, y)dy is a function of x and we can
c
integrate it with respect to x from x = a to x = b. The resulting integral
Z bZ d
f (x, y)dydx
a c
Z dZ b
is called an iterated integral. Similarly one can define the iterated integral f (x, y)dxdy.
c a
Z 3Z 2 Z 2Z 3
2
Example 11.3. Evaluate the iterated integrals (a) x y dydx, (b) x2 y dxdy.
0 1 1 0
" #y=2 " #x=3
3Z 2 3 x2 y 2 3 3x2 x3
Z Z Z
Solution. (a) x2 y dydx = dx = dx = = 27/2.
0 1 0 2 y=1 0 2 2 x=0
" #x=3 " #y=2
Z 2Z 3 Z 3
2 x y Z 2
9y 2
2
(b) x y dxdy = dy = 9y dy = = 27/2.
1 0 1 3 x=0 1 2 y=1
Consider a positive function f (x, y) defined on a rectangle R = [a, b] × [c, d]. Let V be the
volume of the solid under the graph of f over R. We may compute V by means of either one
of the iterated integrals:
Z bZ d Z dZ b
f (x, y)dydx, or f (x, y)dxdy.
a c c a
z.. z..
.. ..
.......... z = f (x, y) ..........
.... ...
...... ... ......
... ............. .............
.. ........ .... ... ... ........ .... ...
... ........ ... ........
... . .. . .. ....... ..... ... ... .
.. . .
. .
........ ..... ..
... ....... . .
... .. ... ....... .... . ..
... ..
........ ........ ..........
... .......... ... .. ... .......... . ... ..
... ...... ........ ... ...... ... .... ...
... .
. ...... .. ... ..................... ..
.. ... .
....... .
............ ...
.
..
..
... .
... .. ....
....................... ... .. ... ... ..
. .
........... .....
. . ..
... .. ... . .................................. ... ... .
... .
. .
............. ......
.
. . .
.............................................................. ... .. ... . .
. ................................ . ..
... .... .. ... ... . ..
... ...... .................................................... ... .. ... ...... ..................................... ..
. ..
... ..
....
............................... .. ..................... ... .. ... ..
..
..
................................................. .. ..
.... . .. . . .. . . . . . . . . . . .. ... . . .................. ..
... . . . . . . . . . . . . . . . . . . .. ... .. . ..
... .. .......... ..................... ... ... .. . ................... .. y
. .
.
. c
.......... ....................... ... d ..
. .
.
. c
. ................... .. d ..
..
.................................................................................................................................................................................................................
. .
y ..
............................................................................................................................................................................................................
. . .
y
..
....
. .
. ........ ..................... ..
.
. . . .. ..
....
. .
. . .................. ..
. ...
..... .. ................................................................................................................. ..... .. .......
.. ... . . . . . ... ..
a ........... .. ...................................... ... ....... a ........... .. ........................................................................................
.
..
......
. .
. .
. . ...
. . . . . . . . . . . . . . .
. ... ..
......
. .
. ..
..... ...................... ..
. ...
..
.. .... . . . . . . . . . . .. ... ... .................. .. ...
..... .. ......................................... ............ .
...
.. .
. .
...
.
..
............... .
.....
.
x ............ .. ................................................................................ ..... .. ..... .. .....
..... .. ........ ..... .. ........ . . .... .. ........
.... . .....
. ....
. ..... .... ................
..
.. . . ..... .. ..... ..... .. .....
b .......... ................................................................................ b ............. ....................................................................................
..... ..
..... .....
.......... ..........
....... .......
x x
RbRd RdRb
Figure 67 a c
f (x, y)dydx Figure 68 c a
f (x, y)dxdy
Theorem 11.4. (Fubini’s Theorem) If f (x, y) is continuous on R = [a, b] × [c, d], then
Z bZ d Z dZ b
f (x, y)dydx = f (x, y)dxdy.
a c c a
c Department of Mathematics 53
Semester II(2005/06)
Example 11.5. Find the volume of the solid S that is bounded by the elliptic paraboloid
x2 + 2y 2 + z = 16, the planes x = 2, y = 2, and the 3 coordinate planes. See figure 69.
ZZ Z 2Z 2
2 2
Solution. Volume = 16 − x − 2y dA = 16 − x2 − 2y 2 dxdy = 48.
R 0 0
z..
..
.........
...
. .............
..
...
..... .. .............
...... ....
. ......
. .....
... .... .....
... ... .....
..
. .
. z = 16 − x2 − 2y 2
.....
.....
... ..
. ....
... .
. ....
. ...
.... ..
. ...
................. .... ...
...
... ........... .. .
....
........... ..........
.
..
......
.... .............. .......... ....
..
... .
. ....... ...
..... .....
.............
.... .... ... ....
... . . 2
...............................................................................................................
...
.. .... .
...
.
..
.
.
. ..
..... y
... .... .... ....
..... R .....
.... ......... ... .........
... ....... .. .....
2 ..................................................................................
.....
.....
.....
..
......
.
....
.....
..........
.......
x
16 − x2 − 2y 2 dA
RR
Figure 69 R
ZZ
Example 11.6. Let R = [0, π2 ] × [0, π2 ]. Evaluate sin x cos y dA.
R
ZZ Z π Z π Z π Z π
2 2 2 2
Solution. sin x cos y dA = sin x cos y dydx = sin x dx cos y dy = 1 × 1 = 1.
R 0 0 0 0
! Z !
ZZ Z b d
g(x)h(y) dA = g(x) dx h(y) dy ,
R a c
11.3. Doubles Integral over General Regions. Let f (x, y) be ZZ a continuous function
defined on a closed and bounded region D in R2 . The double integral f (x, y) dA can be
D
defined similarly as the limit of a Riemann sum.
c Department of Mathematics 54
Semester II(2005/06)
z.. z = f (x, y)
.
.........
...
....
..
... ...............................
......
... ..... ... ... ... ... ............
... ......... ...
....... .
......
... .............................................. ..
... .. ..
... .. ..
... .. ..
.
..
...........................................................................................................
.
y
..
....
. ..
. ...... ... ... .... ..... .......
..........
.
.. .
..... .............. .
.....
.....
...D
.............................................
...
..
..
..........
......
x
Figure 70 Double integral over a more general region
In particular, if D is one of the following two types of regions in R2 . We may set up the
corresponding iterated integral to compute it.
If D is the region bounded by two curves y = g1 (x) and y = g2 (x) from x = a to x = b, where
g2 (x) ≥ g1 (x) for all x ∈ [a, b], we called it a type 1 region. In this case, the double integral of
f over D can be expressed as an iterated integral as given in figure 71.
Similarly, if D is the region bounded by two curves x = h1 (y) and x = h2 (y) from y = c to
y = d, where h2 (y) ≥ h1 (y) for all y ∈ [c, d], we called it a type 2 region. In this case, the
double integral of f over D can be expressed as an iterated integral as given in figure 72.
y.
..
.........
y = g (x) 2
.... ....... .......................
.. .. .................................... .....
... .
. .
. ..
... .... .... ...
... .
. .
. ...
... ..
. ..
. ...
... ..
. ..
. ...
... ..
. D ..
. ...
ZZ Z b Z g2 (x)
... ..
. ..
. ...
...
...
.. ...............
.
.........
.
...... ..
...... .
.
.
.
. ...
...
. f (x, y) dA = f (x, y) dydx
... ...... .
...
..
D a g1 (x)
... ..... .
.......................
..
...
... y = g (x) 1
...
... . . .
..........................................................................................................................................
a x
x
b
y.
.
...
........ x = h (y) 1 x = h (y) 2
..
d .......... .................................................
... ..... ...
.. ..... .
...
... .
.. ...
... .... .
... ...
...
D ..
..
.
. . .
.
y ........
..
...........................................................
..... .
... .... ... d Z h2 (y)
ZZ Z
... ...
... .....
...
...
.. ...
..
...
...
...
...
f (x, y) dA = f (x, y) dxdy
c ......... ..................................................... D c h1 (y)
....
..
.....................................................................................................................................
x
ZZ
Example 11.7. Evaluate (x + 2y) dA, where D is the region bounded by the parabolas
D
y = 2x2 and y = 1 + x2 .
y
.. ...
... ......... ..
... ... ...
.
... ... ... ...
..... .
. ........
......
.... ....
y = 1 + x (1, 2) 2..........
(−1, 2) .... ... ..
....... ... ....
....... .
. ..
.........
.. ..... .. .
.... ...
.. ..... . .
.
.. .
.. ......
.. ...... .... ...... ..
... ........ . ....... ....
... ............. .... .................... . y = 2x 2
... ......................
. ....
...
D .... ...
...
... .
. .....
.
....
.... .... ....
..... . ....
. ..... .... .......... .
...................................................................................................................................................................................................................................
... x
−1 1
Figure 73
ZZ Z 1 Z 1+x2 Z 1 h iy=1+x2
(x + 2y) dA = (x + 2y) dydx = xy + y 2 dx
2x2 y=2x2
D Z−1
1
−1
y.
....
........
...
...
...
..
1 .......................................................................
.... .... ..
..... 1Z 1
.....
... ...
Z
.....
...
...
..
... ....
... .........
.. ......
I= sin(y 2 ) dydx
x ...... ...... 0 x
... .....
y=x
... ..
......
.
... ....
.....
... ........
... .......
............................................................................................................
. x
x 1
y.
...
..........
....
...
...
.
1 .......................................................................
...
.... .....
... .....
... ..
......
.
... ....
..... 1Z y
.....
Z
..
y ...........................................
....
...
y=x .
.....
.
..
...
I= sin(y 2 ) dxdy
... .....
..... 0 0
... .........
..........
...........................................................................................................
x
y 1
Figure 74
1Z y 1 1h 1
1
Z Z ix=y Z
sin(y 2 ) dxdy = y sin(y 2 ) dy = − cos(y 2 ) = 12 (1 − cos 1). x sin(y 2 ) dy =
0 0 0 x=0 0 2 0
Example 11.10. Find the volume of the solid above the xy-plane and is bounded by the
cylinder x2 + y 2 = 1 and the plane z = 0 and z = y. See figure 75.
c Department of Mathematics 56
Semester II(2005/06)
z..
..
........
....
...
...
.
..
. ........ ................ .... ............................
. . ..
........ ....
.........
......
..
.... ... ..
......... ... ....
... ............. .
. . ...........
... ......................... ... ................................. ...
........
... .. ....
.. ... ...
.... ... ...................
... ... ...... .... ......
... ... the plane ... ...
... ... ... ...
... ... z=y ........ .
. ........
... ... .... .......... .......
... ... .
. . .
... ... . . ... .... ....
. .
... ... ... .... ..
... .
. .. ... ....... ....
. ... ... ... ... .
... ... ... ... ... ... .. ... .......................... ...... ....
........ • ... .... ....... ....... ........
.. ..
..... .
........................ ................. ... .........................................................
....
.......
...........
....
... ....
............
...
..... .......
... y
.............................................................................
•
...
.
........
.
.....
.....
.....
. .....
...........
.
x
Figure 75
Solution. Since the plane z = y is the top face of the solid, we may use the function defining this
plane as the height function of this solid. The function whose graph is the plane z = y is simply
f (x, y) = y. Therefore, the volume of the solid can be computed by integrating this f over the
bottom face of the solid which is the semi-circular disk D = {(x, y) | x2 + y 2 ≤ 1, y ≥ 0}.
ZZ Z 1 Z √1−y2 Z 1 √
x= 1−y 2
Volume = f (x, y) dA= √ y dxdy = [xy] √ dy
2 x=− 1−y
D 0 − 1−y 2 0
1 1
2
Z q
3
= 2y 1 − y 2 dy = − (1 − y 2 ) 2 = 2/3.
0 3 0
............................................
.............. ... .........
............. . .....
........ ...
....
.. ........ ..
.
..
...
....
...... .. ...
........ .
.
.
2 D .........
.
.....
.. D 1 ... .
.......
.
.... D
....... .. ......
............ .. ......
........................ .. .......
...............................
ZZ ZZ Figure 76
(5) dA = 1 dA = A(D), the area of D.
D D ZZ
(6) If m ≤ f (x, y) ≤ M for all (x, y) ∈ D, then mA(D) ≤ f (x, y) dA ≤ M A(D).
D
11.4. Double Integrals in Polar Coordinates. Consider a point (r, θ) on the plane
in polar coordinates as in figure 77. An increment dr in r and dθ in θ give rise to an area
dA = rdθdr. This is the area differential in polar coordinates.
c Department of Mathematics 57
Semester II(2005/06)
y. y.
...
.. ..........
........ ....
.... ...
.. ..
... ... .....
... ........... ... .. ..........
dr b ... .....
... ... .......... ... ... .....
... ... ..... ... ....... .....
....
... ......
. .... ... R ... .......... ...
... . ..... ........................................
Area dA = rdθdr ... ... ..... ...
... .. ..........
. ... ... ... ..... ..
... . .... .. ... ... .......
r rdθ . ... ............. ... . ..... ....... ..
..
·
.. . . .
... .. ....... ... .... ...
... . . .
. ... .
.. .........
.
.. . .....
... ... . .. ...
.
.......... .......
... .. ..... (r, θ)
.. .
... ..... ...... ............ a
dθ ... ... ....... β ... ... .............
... ... ........... ... .... .......... .... ....
θ .......... ...
α ............... . ..
..................................................................................................... ..............................................................................................................
x x
Figure 77 Figure 78 A polar rectangle
y.
...
..........
... ......
... ... ............
... ... .....
... ...
. .....
.....
... .... ...
... ............. ........
....
D
...
. ... ...... .....
... .... .......... .
.. .
...
. .
... ...... .......
... r = h (θ)
... .......
... .. ......... ....... 2
... ..... .......... ................
. ... ....
r = h (θ)
... ........ ...... ....... 1
... .... .............................
β ... .... .....................
... ... .................... ...
........................ .... ....
α ...........................................................................................................................
x
Figure 80 A general polar region
c Department of Mathematics 58
Semester II(2005/06)
Example 11.12. Find the volume of the solid that lies under the paraboloid z = x2 + y 2 ,
above the xy-plane, and inside the cylinder x2 + y 2 = 2x.
Solution. The cylinder x2 + y 2 = 2x lies over the circular disk D which can be described in
polar coordinates as
π π
D = {(r, θ) | − ≤ θ ≤ , 0 ≤ r ≤ 2 cos θ}.
2 2
z..
...
........
....
..
..
. ... .. . ...............................................................................
.......................................................... .... ........
.....
. .....
....... ...
..................... ...
........ ...
........................................................... .... ..... .
. .........
.......................................................................... ... .
.........
............ .... ...
.............. ... ...
................... ... ...
......................... ... .
....
.
....
........................... .. ......
............................................................. .... ......
.......
......... . . . . . . . ...... .. ........
x ......................................................................................................................................
......... . . . . . . .........
............................................
....
.....
.....
...
......
....
.....
...........
......
y Figure 81
The height of the solid is the z-value of the paraboloid. Hence the volume V of the solid is
π
ZZ Z
2
Z 2 cos θ
2 2
V = (x + y ) dA = r2 rdrdθ = 3π/2.
D − π2 0
Exercise 11.13. Show that the volume of the solid √ region bounded by the three cylinders
x2 + y 2 = 1, y 2 + z 2 = 1 and x2 + z 2 = 1 is 16 − 8 2.
Solution. The object is symmetric with respect to the 3 coordinate axes. Thus, it suffices to
compute the volume of the portion E with nonnegative x, y, z coordinates. In other words, E
lies in the first√octant of the coordinate
p system. For this object E, it is bounded on top by the
graphs of z = 1 − x2 and z = 1 − y 2 . On its side and bottom, it is bounded by√the cylinder
x2 +py 2 = 1 and the three coordinate planes. Furthermore, the graphs of z = 1 − x2 and
z = 1 − y 2 intersect along a curve√ on the plane x = y. Thus we can simply find the volume
of solid under the graph of z = 1 − x2 over the region D = {(r, θ) | 0 ≤ r ≤ 1, 0 ≤ θ ≤ π4 }.
Then the result volume is 16 times this value. Using polar coordinates in D, we have
c Department of Mathematics 59
Semester II(2005/06)
z..
..
..........
..
....
..
....
..
.....
.. E
......................................
. ........
..
..... .... .........
. .......
..... ... ... ......
. ... .....
... ... ... ..... p
.. . 2
..
. .
.
..
.
.
... ..........................
z= 1−y
√ . .
.
....
.
....
.
...
...
.
...
...
...
z= 1−x 2 ........................... . . .
....
... ...
... ... .. ...
.. ...
... ... .. ...
.... ... . ....... ...
... ... .... ...... ...
... .... .... .... ...... ...
. . .
.. ...... ..... ........
. .
...
... .
.
..
...................
.
. .
.. ....... ......................... ... ... .....
.
.
... .
... .......... ... . . . ......... ..... ... .................................
.
..
. . ......................................
...
... ..
.. .
... ........ ...
.... .....
... ....................
.. ...
.
.
.
.
......
. .
. ... . ..............................
.. y
.. ... .............................. 2 2
..
...............
.. ..
.
. D .. .. ............
. . ....
.. .... .
x +y =1
.... ...............................................................
....
.....
.....
..
......
.
.....
.....
.....
.....
..
......
..
..........
......
x
Intersection of 3 cylinders
π
ZZ p Z
4
Z 1p
Volume = 16 1 − x2 dA = 16 1 − r2 cos2 θ r drdθ
D
π 0 0
r=1
1
Z
4
= 16 − (1 − r2 cos2 θ)3/2 dθ
Z0 π
3 cos θ r=0
3
4 1 − sin θ
= 16 dθ
0 3 cos2 θ
16 π √
= [tan θ − sec θ − cos θ]04 = 16 − 8 2.
3
z.. z = f (x, y)
..
.........
.... ................
...............................................
...................
... ........ ...................
... ...... .....
... ........ ...... ..
... P . ...
R ..... ..
... .. ...
.
.. ....................................................... .......
. ..
. .. ......... .... .
... ... ..... .. .......... .... ..
... ... ........ . .... .. ..
... ... . .
.. .......... ...
. .
....... .........
. .... ..
... ... . .............. .. .
...... ...... ..
. ..
.. .......................................... .... .
... . . . .
. ..
... Q ... dS ... ....................... ..... ..
. ..
... ..
. ........... .. ............ ..
. ..
... ... .. ...
. ..
... ..
. .. ... ..
... .. .. .. ..
. .......................................... .
... ..
. . . . .. . . .. .. ..
. ..
... ..
. .. ... ... . ........ .
. .
. ..
... ...
. .
.. ... ....... .
. .. .
. ..
... .. ...............
. .
. .. .
.
. ..
... ............ .
. .. .
. ..
... . .. .
. .. . ..
... . .
. .. . ..
.
. .. .
.
... .
.
.
. .
. ..
... . c .
. .. . d ..
...............................................................................................................................................................................................................................................................................
.
..
.. .
...
.
.... .
.
.
.
.
.
.
.
.
....
.
..
.
.
.
.
.
..
..
.... .... y
..
..... .. .. .. .. ..
..
..... .
. .
. .
. .
. .
.
a ................ .. .. .
......................................................................................................................
..
..... .. ..... . .
..... .. ..... dy 0... .. .....
..
......
. .
. P ...
...... ... ...
...
...
...
.... .. .. .
. .... ..... .. ..
.
.. .
. ... .....
.
..... .. .....
....
.....
.... ....
..... ..
....
.....
....
..... .. dx ..... dA ..... ..... .. .....
..... .. ..... ..... ..... .. .....
..
......
. .
. ...
...... .... ..... ..
.. ..... .
. ..
......
.
..... .. ........................................ ..
..... D .. ......... .. .......
............... .................................................................................................................
b .
.
.....
..
....
.....
.........
.......
x
Example 12.1. Find the area of the part of the paraboloid z = x2 + y 2 that lies under the
plane z = 9.
Solution. The paraboloid lies above the circular disk D = {(r, θ) | 0 ≤ θ ≤ 2π, 0 ≤ r ≤ 3}.
RR q
Surface area = fx2 + fy2 + 1 dA
RRD p
= 2 2
D R 1√+ 4(x + y ) dA
R 2π 3
= 0 √ 0 1 + 4r2 rdrdθ (change to polar coordinates)
π
= 6 (37 37 − 1).
Denote the volume of Cijk by 4V . Then we may form the Riemann sum:
l X
m X
n
f (x∗ijk , yijk
∗ ∗
X
, zijk )4V.
i=1 j=1 k=1
The limit exists if f is continuous. The triple integral of a continuous function defined on a
more general closed and bounded solid in R3 can be defined in a similar way.
Theorem 13.1. (Fubini’s Theorem for triple integrals) If f (x, y, z) is continuous on B =
[a, b] × [c, d] × [r, s], then
ZZZ Z sZ dZ b Z sZ bZ d
f (x, y, z) dV = f (x, y, z)dxdydz = f (x, y, z)dydxdz = etc.
B r c a r a c
(Note that there are 3!=6 such iterated integrals involved and they are all equal.) Furthermore,
the theorem is valid for a general closed and bounded solid.
ZZZ
Example 13.2. Evaluate xyz 2 dV , where B = [0, 1] × [−1, 2] × [0, 3].
B
Solution. ZZZ Z 3Z 2 Z 1
2
xyz dV = xyz 2 dxdydz = 27/4.
B 0 −1 0
13.1. Triple Integrals over a General Bounded Region. For each of the following
three types of solid regions, we may write down the triple integral as an iterated integral of a
double integral and a simple integral.
z.. z = u2 (x, y)
..
.......... ................................................
... ...... ...
......... ..
... ... ................ ........
... ... .............................. ....
... .
...
...
.
... ....... . .. ... ... ... ... ... ... ... ... ....
... ..... ...
...
...
.........
..
. E .
.... .
............................ .......... .
......
... .. ..
... .. ZZ "Z #
... z = u (x, y)
..
.. 1 ...
ZZZ u2 (x,y)
... .. .
...
...
.. ..
.......................................................................................
f (x, y, z) dV = f (x, y, z) dz dA
..
..
.
.
...... ...
...
...
...
...
...
...
...
...
..
.
..
.
... y E D u1 (x,y)
..... .. ..
..
...... .
. ..
.... .
.................
.................................... ..
..... ......
..... .... D .
..... ........ ....
.
......
. ............................................
.
...
.
.
...
.
.....
..........
.....
x
Figure 84 Type 1 solid region
c Department of Mathematics 62
Semester II(2005/06)
z..
.
........
...
.... .............................
.. ............... ......
... ......
... D ...
... ....... ......
... .... ..................................................
... . ... ...
.
... .. ..
... ...
... ... ...
... ... ...
................................. ... ZZ "Z #
... .
.
..
..
. . ....
. x = u (y, z)
....... ... .
.. 1
ZZZ u2 (y,z)
... ......
..... ....
... ........ ...... f (x, y, z) dV = f (x, y, z) dx dA
... ..... ... ... ... ... ... ... ... ..............
. ......
. .........................................
..... .....
.
.... .... ..
.... ................................................. ........
..... .......
E
........
.....
.. ... y E D u1 (y,z)
..... ........ ...
..
......
. ........
.................................... ..
.. ..
..
.....
..... x = u (y, z)
..... 2
.....
..
.....
.
..
....
.
..
.........
.......
x
Figure 85 Type 2 solid region
z
....
......... y = u (x, z) 2
.
............... ... ... ... .... ... ... ... ... ............................................
... ..... .... ... ... ... .....
.... ... ... .... . ..
. ...
.. ... ... E .. .. ... ...
... ... ... ... .. ... ...
..... D .. ... ..... .. ...
.
...
... .... ... ... .. .... ...
... ... .
. ... . ..
. ...
... .. ..
. ... . ....
. ... ZZ "Z #
... ... ... ... . ... ... u2 (x,z)
. .
... ... . .
ZZZ
........... ... ... ... ...... ... ... ... ... ...............................................
y = u (x, z)
....
..
1
f (x, y, z) dV = f (x, y, z) dy dA
. .................................................................
.....
.....
.. .
...
...
...
...
...
... . . ..
...
...
...
.... .
...
...
...
.. .
y E D u1 (x,z)
..
......
.
....
.....
.....
.....
. . .......
.
.....
.....
.....
..........
......
x
Figure 86 Type 3 solid region
ZZZ
Note that if f (x, y, z) = 1 for all (x, y, z) ∈ E, then 1 dV is just the volume of E.
E
ZZZ p
Example 13.3. Evaluate x2 + z 2 dV , where E is the solid region bounded by the
E
paraboloid y = x2 + z 2 and the plane y = 4. See figure 87.
ZZZ p ZZ Z 4 p
x2 + z 2 dV = x2 + z 2 dy dA
E ZZD h x2p
+z 2 i4
= y x2 + z 2 dA
ZZD p x2 +z 2
= x2 + z 2 (4 − x − z 2 ) dA
2
D
Z 2π Z 2
= r(4 − r2 )r drdθ (change to polar coordinates)
0 0
= 128π/5.
c Department of Mathematics 63
Semester II(2005/06)
z..
..
........
u (x, z) = x + z 2 2
1
.... .......
... .......
... ....... .
..... ... ... ... ... ... ... ... ... ... ... ... ... ... ...................... ... ... ..............................
........ ..... ............... ..... ...
... .. .... ............... ...
... .... ..... ............... .. ...
D .... ... ... ........... .. ...
... ... ... .
..
.............. ... ...
... ... ... ................. ... ..
.
... u (x, z) = 4 ... 2
... ... ...........
...
.
.
. E
.. .
.. ...
...
... ... ............ .... .. ... ...
... ... ...... .. ... .. ..
. .
...
.. ....... .. .. .. ... ...
.. . . ... .. ..
.... ........................................................................................................................................................................................ y
... .............. .. ... . .
. ..
... ..... . .
.... . ... . ... ..
...... .. ...... . . .
... .
. .... .
.
...... ..... .
........ ... .. ... ...
..... ... . . .. ... ...
..... .... .. ............... .... . ..
..
......
. . . .... ........... ... ...
. ...
.
.... ... .......... ...
..... ... ... ......... ... ...
..... ... .. ............... . ..
..... ... ... ...............
............... ..... ..
.
............ ..... ....... ................... .......
.
.
x . .
...... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...... . .
z..
..
.........
.... z = u (r, θ)
......................................... 2
... ........ .... .........
.......... ... ....... ..
... .............. ........ ....
... ...................
. ...
.......................
... .. .......
. ... ...
. ...... . ...
... .... ................................. ..........................
... ... ... ...... .....................
... ..
. .. ..................... .... . ...
... ........................... ...........
.
... ..... .. ..........
..... ... ............ z = u (r, θ) 1
... .........................................
...
...
.
...... ................................................................................................................................................................
.
........... .......... ... ... ... ... ...
..... .......... ... ... ... ...
y
.
......
. ... ... ... ..........
.
..... .. ......................
... ....... ...........
..... ... .......... .
.
....
.
.. ...................................
. .......
. D ........
.
...
. . ..... ............
..
.. .
...
.
.. r = u (θ) 1 ..
.....
.. ....
. ..................... ......... r = u (θ) 2
.... ..............
.....
.........
.......
x
Figure 88 A cylindrical rectangle
Example 13.6. Let E be the solid within the cylinder 2 + y 2 = 1, below the plane z = 4,
ZZZx q
and above the paraboloid z = 1 − x2 − y 2 . Evaluate x2 + y 2 dV .
E
c Department of Mathematics 64
Semester II(2005/06)
[Answer=16π/5]
13.3. Triple Integrals in Spherical Coordinates. Consider the volume element in
spherical coordinates. To do so, take any point P (ρ, θ, φ). Make an increment in each of
the coordinates. See figure 90. Let’s calculate the volume of the solid arising from these
increments. The projection of OP onto the xy-plane has length ρ sin φ. Thus the thickness of
this volume element is ρ sin φdθ. It opens up a sector of width of ρdφ. Thus, the volume is
dV = ρ2 sin φdρdθdφ. .... .. ..
..... .........
... .. ..........
... ........ ...
z ... .... .......... .........
..... ....
... ... dρ ..... ...
........
. .. ....
. ..... ..
..... ...
... ... ... ..... ...
.... .. ... .........
... ....... .. .
... .
.
.
...... .. .... .........
. ρ sin φdθ
... ............ ... . .....
... ........ ....... .
.. ... ............
.
. ....... ... ... ...... .
... .. . ................. ...... ............
... ..... .... ...... ........
...
...
....
....
... ............. ...............
... ..... ...............
... ..... ............
......... .......
.... .........
dφ ·P (ρ, θ, φ)
.. .. ......
.......... .
..
.................................................................................................................................................. y
...
.
..
.
.
...
.
.
. 0
.....
.....
.....
..... dV = ρ sin φdρdθdφ 2
..
......
.....
..
......
.
..
.........
.......
x
Figure 90 The volume element in spherical coordinate
Now consider a spherical rectangle
E = {(ρ, θ, φ) | a ≤ ρ ≤ b, α ≤ θ ≤ β, c ≤ φ ≤ d},
where a ≥ 0, β − α ≤ 2π, d − c ≤ π. The triple integral of f over E can be expressed as follow:
ZZZ Z dZ β Z b
f (x, y, z) dV = f (ρ sin φ cos θ, ρ sin φ sin θ, ρ cos φ) ρ2 sin φdρdθdφ.
E c α a
c Department of Mathematics 65
Semester II(2005/06)
3
ZZZ
2 +y 2 +z 2 ) 2
Example 13.8. Evaluate e(x dV , where B is the unit ball {(x, y, z) | x2 + y 2 +
B
z 2 ≤ 1}.
Solution. Using spherical coordinates, we have
ZZZ
4 3
2 +y 2 +z 2 ) 2
Z π Z 2π Z 1 3
2) 2
e(x
ρ2 sin φdρdθdφ = π(e − 1). dV = e(ρ
B 0 0 0 3
Note that the corresponding triple integral formulated in Cartesian coordinates is very hard
to evaluate.
Exercisep13.9. Use spherical coordinates to find the volume of the solid that lies above the
cone z = x2 + y 2 and below the sphere x2 + y 2 + z 2 = z. See figure 91.
z..
.
..........
..
.
. ... ... ...................................
.
. . ....... ....
.......
......
....... ... .....
... ...
... ... ...
. .. .
. ...
... ..
. ...
. ..
.... .. .. .
.
... ... ... . ... ... ... ... ... ... .. ....
. ... .
..
...... ...... ... . .......
..... ....
...... .
. .......
............... . ....
..... ...........................................................................
..... .
..... ..
.
. .........
..... ..
..... .... π .........
..... . .
..... .........4 .....
..... .. ............
..... .. .....
............
.....
....................................................................
y [Answer: π/8]
.....
........
.
....
.....
.....
..........
......
x
Figure 91
Exercise 13.10. Let B be the solid ball in R3 of radius R centered at the origin. The 4-
dimensional ball H in R4 of radius R centered at the origin is the 4-dimensional solid given by
{(x, y, z, 2 2 2 2 2
ZZZ qw) | x + y + z + w ≤ R }. The volume of H can be expressed as the triple integral
2 R2 − x2 − y 2 − z 2 dV. Using spherical coordinates, show that this volume is π 2 R4 /4.
B
Figure 92
For example if T is the transformation to polar coordinates T (r, θ) = (r cos θ, r sin θ), then T
maps a rectangle [r1 , r2 ] × [θ1 , θ2 ] in the rθ-plane to a polar rectangle in the xy-plane.
c Department of Mathematics 66
Semester II(2005/06)
θ.. y.
.. ....
........ ........
.. ..
... ...
... ...
... ...
... ...
... ... ........
.........
.. S ...
... ..
R
... ..........
θ 2 ...
.......................................................................
... ... ... r 2... .
....
.....
.....
.....
... ... ... ... ....
..
...
...
... ...
...
T ...
... • .
...
... ...
...
... ...
... • ... ....................................................................................................................... ... (x, y) ....
......... .......
.
....
...
... ...
... (r, θ) ...
...
r ...
... 1 ..
.
. .......
.... .. .
........ ..... .
... .
... ... .. ... . ... ..... .....
θ ........ ..................................................................... ... .
...... .....
1 .. ... ... ....... ........
.. ... .. .........
...
.
..........................................................................................................................................
θ θ
.......................................1
... .. ....... ....
.. . . 2
.....................................................................................................
.....
r .....
x
r1 r2
Figure 93 Polar coordinates
Example 14.1. Consider T (u, v) = (u2 − v 2 , 2uv). Find the image of the square S =
{(u, v) | 0 ≤ u ≤ 1, 0 ≤ v ≤ 1}.
Solution. First let’s find out the boundary of the image. Label the edges of the square S by
S1 , S2 , S3 and S4 as shown in figure 94.
v.. y.
... ....... ...
(0, 1) ...............................................S 3
..........................................................
(1, 1)
.
.
.........
.
...................
(0, 2)
... .. .
...
. ..
... ... ..... .... ........
1 2 ....
...
....
...
...
...
...
x = −1 + y 4 .
....
.........
. ...
...
.
.
...
...
... x = 1 − 14 y 2
... ... ... .
.
. ....
0 ..... .........
... .... .. 0
.........
..
.
.
..........
. S3... . .
. .. 2 S
S 4 ...
....
..
. S ... S 2 T ..
.
.
.
.
.
.
.
..
..
..
.
. ..
. ..
. ..
. . ..
.... .... ....................................................................................... ... .... ..
...
... . .
...
...
...
...
... .
.
.
....
.
....
.
.
R ...
...
.. .. ...
... ... . . ...
... ... ..
. ..
. ...
. .
................................................................................................................................
. . . . .
..........................................................................................................................................
. ..
.
..
..
u . .. .
.
.
.
x
0 .. S1 (1, 0) (−1, 0) S40 0 S10 (1, 0)
Figure 94
Similarly, we find out S30 as x = −1 + 41 y 2 for y from 2 to 0 and S40 as y = 0 for x from −1 to 0.
The boundary of the image of S encloses a region R. We are going to show that T maps S
bijectively onto R. We leave it the reader to verify that T is a bijective function for u, v ≥ 0.
As we traverse the boundary of S in the counterclockwise direction, the above calculation
shows that the image of the boundary of S also traverses in the counterclockwise direction.
This means that T preserves orientation. In other words, points on the left hand side of the
boundary of S go under T to points on the left hand side of the boundary of R. Therefore, T
maps S onto R. Another easy way to confirm this is to pick a point P , say (1/2, 1/2) inside S
and check that T (P ) is inside R. Then the region S must be mapped by T into R.
Before we derive the formula for change of variables in a multiple integral, let’s review the
formula for functions of 1 variable. Let the continuous function f (x) be integrated over the
c Department of Mathematics 67
Semester II(2005/06)
interval [a, b]. Suppose we make a substitution x = g(u) so that a = g(c) and b = g(d). Thus
we obtain: Z Z b d
f (x) dx = f (g(u))g 0 (u) du.
a c
Here the formula is valid provided g is differentiable and g 0 (u) 6= 0, except possibly at a finite
number of points. The function g is also required to be bijective so that g −1 exists. Observe
that c may not be less than d. More precisely, if g 0 (u) > 0 for all u between c and d, then
g and hence g −1 is increasing. Thus g −1 preserves orientation or ordering. This means that
c < d and [c, d] is an interval. On the other hand, if g 0 (u) < 0 for all u between c and d, then
g and hence g −1 is decreasing. Thus g −1 reverses orientation or ordering. This means that
c > d and it does not make sense to write [c, d] though we could still integrate from c to d. In
this case, the formula can be rewritten as:
Z b Z c
f (x) dx = f (g(u))(−g 0 (u)) du,
a d
so as to keep the lower limit of integration smaller than the upper limit.
Therefore, if the interval [c, d], (c < d) is mapped onto the interval [a, b] under x = g(u), then
the formula for change of variables can be stated as:
Z Z
f (x) dx = f (g(u))|g 0 (u)| du.
[a,b] [c,d]
How does a change of variables affect a double integral? Let T be a transformation mapping
a point (u0 , v0 ) to a point (x0 , y0 ). Consider a small increment du and dv at the point (u0 , v0 )
along the u and v directions respectively. These increments generate a rectangle of area dudv
whose image under T is a curved parallelogram in the xy-plane. The area of this curved
parallelogram up to the first order approximation is given by the area of the parallelogram
generated by the two tangent vectors adu and bdv at (x0 , y0 ), where a is the derivative of the
curve T (u, v0 ) at (u0 , v0 ), and b is the derivative of curve T (u0 , v) at (u0 , v0 ). That is
dT (u, v0 ) ∂x ∂y
a= = h (u0 , v0 ), (u0 , v0 )i,
du
u=u0 ∂u ∂u
dT (u0 , v) ∂x ∂y
b= = h (u0 , v0 ), (u0 , v0 )i.
dv
v=v0 ∂v ∂v
v.. y.
.. ....
......... ........
.... ....
... ...
... ...
... ...
... ... .............................
... ... ............................
......... .......... .... ....
.....
... ... ...... ..........
... ... ............... ..
.. .....
... .......................................................... ... . ............ .. .....
..
...
... ....
.. ...
...
T (u, v) = (x, y) ...
... ....
.
............. ........
...........
...
...
...
dv ... .........
dudv ... ... b ...... a .........
...... .........................................................
...
... ... .. ....................................................................................... ... .
. ........
... ................................................................ ... ............
• ...
...
.... • ...
... T (u, v ) = (x(u, v0 ), y(u, v0 )) 0
(u , v ) du ...
...
0 0 (x , y ) ...
...
0 0
.. ...
.................................................................................................................................. .................................................................................................................................
....
..
u ....
..
x
. .
Figure 95
Therefore, the area element dA in the xy-plane is dudv times the magnitude of
∂x ∂y ∂x ∂y ∂x ∂y ∂x ∂y
h , , 0i × h , , 0i = − k.
∂u ∂u ∂v ∂v ∂u ∂v ∂v ∂u
c Department of Mathematics 68
Semester II(2005/06)
∂x ∂y ∂x ∂y
That is dA = | ∂u ∂v − ∂v ∂u |dudv.
Definition 14.2. The Jacobian of the transformation T given by x = x(u, v), y = y(u, v) is
∂(x, y) ∂x ∂x ∂x ∂y ∂x ∂y
= ∂u ∂v =
− .
∂y ∂y
∂(u, v) ∂u ∂v
∂u ∂v ∂v ∂u
Therefore,
∂(x, y)
dA = dudv.
∂(u, v)
Theorem 14.3. Let T (u, u) be a bijective C 1 -transformation whose Jacobian is nonzero except
possibly at a finite number of points. Suppose T maps a region S in the uv-plane onto a region
R of the xy-plane. Suppose f is continuous on R. Then
∂(x, y)
ZZ ZZ
f (x, y) dA = f (x(u, v), y(u, v))
dudv.
R S ∂(u, v)
Example 14.4. Find the Jacobian of the transformation from polar coordinates to Cartesian
coordinates.
Solution. x = r cos θ and y = r sin θ. Thus,
∂(x, y) ∂x ∂x
∂r ∂θ cos θ
−r sin θ
= ∂y ∂y = = r.
∂(r, θ) ∂r ∂θ
sin θ r cos θ
Therefore, ZZ ZZ
f (x, y) dA = f (r cos θ, r sin θ)rdrdθ.
R S
Example 14.5. Use the change of variables x = u2 − v 2 , y = 2uv to evaluate the integral
ZZ
y dA,
R
where R is the region bounded by the parabolas y 2 = 4 − 4x and y 2 = 4 + 4x, and the x-axis.
Solution. The transformation is the one discussed in example 10.1. First, let’s compute the
Jacobian of T .
∂(x, y) 2u −2v
= = 4u2 + 4v 2 .
∂(u, v) 2v 2u
Therefore,
1Z 1
∂(x, y)
ZZ ZZ Z
y dA = (2uv)
dudv = (2uv)4(u2 + v 2 )dudv = 2.
R S ∂(u, v) 0 0
ZZ x+y
Example 14.6. Evaluate the double integral e x−y dA, where R is the trapezoidal region
R
with vertices (1, 0), (2, 0), (0, −2), (0, −1).
Solution. Under the change of variables u = x + y and v = x − y, the trapezoid R = ABCD
is mapped bijective onto the trapezoid R0 = A0 B 0 C 0 D0 , where A0 = (−2, 2), B 0 = (2, 2), C 0 =
(−1, 1) and D0 = (1, 1). The inverse transformation is x = 21 (u + v) and y = 12 (u − v) and its
Jacobian is
∂(x, y) 12 1
2
1
= 1 =− .
1
∂(u, v) 2 − 2 2
c Department of Mathematics 69
Semester II(2005/06)
v y
... ...
......... .........
... ...
.... ....
A0 (−2, 2) ..
..............................................................................................................
.............................................
..............................................
B (2, 2)0 x = 12 (u + v) ..
...
.............................................
..........................................
.......................................... .. .. ...
.........................................
........................................
......................................
...................................... .
. ..
...... ...
..................................... ....
....................................
..................................
..................................
................................. .... ..... ...
................................
..............................
..............................
.............................
............................
..........................
..........................
...
.
................................................................
.....
..... y = 1 (u − v) ...
...
..........................2
.........................
0 ... .. 0 ............................................................. ...
D (−1, 1) ... .... C (1, 1)... ...
... . .. ...
... ... ..... ...
... .. ...
... ... ... .. C(1, 0)
............................................................................................................................................................................... u .......................................................................................................................................................................................................
x
... ... ..........................
........................................
...
...
... B(2, 0)
.............................
... ..................................................................
... ... .....................................................
... ......................................................
D(0, −1) ...........................................
....................................
............................
Figure 96 ....................
.........
A(0, −2) ......
Therefore,
Z 2Z v
∂(x, y) u/v 1 1
ZZ ZZ ZZ
e(x+y)/(x−y) dA = eu/v dA = e dA = eu/v dudv
R ∂(u, v) 2 1 −v 2
R0 Z R 0Z
1 1 2h iu=v 2
ueu/v dv = v(e − e−1 )dv = 34 (e − e−1 ). =
2 1 u=−v 2 1
Note that the above transformation reverses orientation.
For the case of triple integrals, we have a completely analogous formula for change of variables.
Suppose T (u, v, w) = (x(u, v, w), y(u, v, w), z(u, v, w)) is a C 1 -transformation from R3 to R3
mapping a solid region S bijective onto the solid region R. First, the Jacobian of T is defined
to be
∂x ∂x ∂x
∂(x, y, z) ∂u ∂v ∂w
∂y ∂y ∂y
= .
∂(u, v, w) ∂u ∂v ∂w
∂z ∂z ∂z
∂u ∂v ∂w
If f (x, y, z) is a continuous function defined on R. Then
∂(x, y, z)
ZZ ZZ
f (x, y, z) dV = f (x(u, v, w), y(u, v, w), z(u, v, w)) dudvdw.
R S ∂(u, v, w)
∂(x,y,z)
Exercise 14.7. Show that the Jacobian ∂(ρ,θ,φ) of the transformation from spherical coordi-
2
nates to Cartesian coordinates is −ρ sin φ.
Figure 97
or simply F = P i + Qj.
Definition 15.2. Let E ⊆ R3 . A vector field on E is a function F that assigns to each point
(x, y, z) in E a three dimensional vector F(x, y, z).
That is F(x, y, z) = P (x, y, z)i + Q(x, y, z)j + R(x, y, z)k = hP (x, y, z), Q(x, y, z), R(x, y, z)i.
Example 15.3. A vector field on R2 is defined by F(x, y) = −yi + xj. Show that F(x, y) is
always perpendicular to the position vector of the point (x, y).
Solution.
y.
...
..........
...
..
................................. • ...
........
.........
.....
...
.
.
.
...
.
.. • .
.
.
.
• ....
....
..
. .
.
....... .
...................• .
....
....
..........
... .. ...
... ........... ..
......................................................................................................................................................
• • • • •
...
...
.........
....
.
.
....
x
........ ... ..
.... .................. .........
• ... ... ......
....
•
..... ...
. • ....
..... .
.......
......... ....
• ...............................
...
...
...
...
Figure 98
Figure 98 shows
p the vector field F. Note that hx, yi · F(x, y) = hx, yi · h−y, xi = 0. Also
|F(x, y)| = x2 + y 2 . The vector assigned by F to the origin is the zero vector.
Solution. ∇f (x, y) = 2xyi + (x2 − 3y 2 )j. The gradient field and the contours of f are drawn
on the diagram in figure 99.
c Department of Mathematics 71
Semester II(2005/06)
y.
....
........
....
............... ... ..
............... ..............
............... ... ...............
................ ... ...............
................ .
. ..........................
.
...............
.... ...............
................ ...............
............... ... ..............
............... ... ...............
.................
................ ... ..........................
................ ... .............
........... ...... ..... .... ...
... ..... ...... ... ..... ..... .....
... .... ..... ... ...... ..... ...
... ..... ...... .
. ........................
... ...... ....... ... .. .
... ..... ............ -1 .... ...... ...... ...
... ............ ............ . ....... .... ...
... .......... ........................................................ ........ ....
... .................. .
. ................. ..
........ ... ... ....... . .. .... ...
.
........ .. 0
..... 0 ....
... .....
..... ...
1 .
. .
•
........ ..
......
..
.....
.....
.....
..
.
.
. ......
.
. .
........
.
....
.....
......
............... 1
... . . . . . .. . . .. .. . . . . .. . ........... ..... .. .....
. . ..........................
........................ . ............. . . ......................
.........................................................................................................................................................................................................................................................................
. .
..................................
.................... ..
.
. .
.
... ... .......
.
.
.
.
.
............................
.................. x
...........
..... . .
.
.... ... ........ .............. .. . ..
-1 .....
... .
...
.....
. ...
.
..... ............. ....
.
.
.
.....
......
..... ...
. .
. .......
-1
... ..... ..... .....
. ..... ... .....
∇f (x, y)• ... ..... .. .. ..
.............................. ............................................................ ................... ....
. .. .... ....
.......... .... ...
... ..... ....... ....... ..... ..
... ..... ....... 1 ... ...... ..... ...
... ...... ....... ...
...... .... ...
....................... .
.
...... ..... ..
...... ..... ....
.................... .
.
. ..... ...........
.
............... .
. ..... .... ....
.. . . . ..............
....................... .
.
. ..............
..............
.... .............. .
.
. ...............
.. . . ................. .
.
. ..............
..................... .
. ..............
..
............ .
. ...............
. . . .
. . .
. ..............
. . . .
................. .
.
.
...............
...............
. . ..................
. .
.
. ...............
............... ....
.............
...
...
...
Figure 99
Notice that the gradient vectors are perpendicular to the level curves as is proved in 8.6 using
the chain rule.
p
Exercise 15.6. Find the gradient vector field ∇f of f (x, y) = x2 + y 2 . Sketch ∇f .
Definition 15.7. A vector field F is called a conservative vector field if it is the gradient of
some scalar function, that is there exists a differentiable function f such that F = ∇f . In this
situation, f is called a potential function for F.
For example, F(x, y) = 2xyi + (x2 − 3y 2 )j is conservative since it has a potential function
f (x, y) = x2 y − y 3 .
Not all vector fields are conservative, but such fields do arise frequently in physics. For instance,
the gravitational field given by
−mM Gx −mM Gy −mM Gz
F= 3 i + 3 j + 3 k
2 2 2
(x + y + z ) 2 2 2 2
(x + y + z ) 2 (x + y 2 + z 2 ) 2
2
Figure 100 C = C1 + C2 + C3
y.
...
..........
....
(0, 2) ...
...
• ..
.. . ..
... .... .......................
. .....
..........
....... ... .........
..
..... .
.
.
........
.......
C ... ....
1........... .
.
......
.....
....
............ ..
. ..
.........................................................................................................................................................................
.. ... .
.
. ...
.. x
.
....
..
. ... 0 .
.
.
.... ......
.....
..... .......
..... ... ........
..... ... .........
..
...... . ..
...
.............
.
....
.
.....
.
...............
.............. ...
.
. ...... C
.
.
2
..... ...............
..... ..................... ....
........................................ ... ....
• ..
.......... ..
...
(−5, −3)
Figure 101
Z Z 1 Z 1
2 2
Thus, y dx + xdy = (5t − 3) 5dt + (5t − 5)5dt = −5/6.
C1 0 0
(b) C2 : x = 4 − t2 , y = t, −3 ≤ t ≤ 2.
Z Z 2 Z 2
Thus y 2 dx + xdy = t2 (−2t)dt + (4 − t2 )dt = 245/6.
C2 −3 −3
•B.
...
−C ..•
... B
C... .............
....................... .
...
.
. ..
. ..
..
..
..
...
...
..
..
....
.....
..... .......
..... .....
.... ....
... ...
... ...
... ...
... ...
... ...
A• A•
Figure 102
For example the upper semicircle C in the xy-plane centered at the origin with radius 1 joins
the point (1, 0) to (−1, 0). It has a vector equation in the form r1 (t) = hcos(πt), sin(πt)i,
t ∈ [0, 1]. Then −C can be parametrized by r2 (t) = hcos(π(1 − t)), sin(π(1 − t))i, t ∈ [0, 1] and
−C joins (−1, 0) to (1, 0).
c Department of Mathematics 74
Semester II(2005/06)
Note that because the sign of x0 (t) and y 0 (t) reverses in −C, we have
Z Z Z Z
f (x, y) dx = − f (x, y) dx and f (x, y) dy = − f (x, y) dy.
−C C −C C
But because the arclength differential is always positive,
Z Z
f (x, y) ds = f (x, y) ds.
C −C
For line integral of a function f (x, y, z) along a space curve C, we have the similar definitions:
s
b b 2 2 2
dx dy dz
Z Z Z
0
f (x, y, z) ds = f (r)|r (t)|dt = f (x(t), y(t), z(t)) + + dt,
C a a dt dt dt
Z Z b
f (x, y, z) dx = f (x(t), y(t), z(t))x0 (t) dt,
C a
Z Z b
f (x, y, z) dy = f (x(t), y(t), z(t))y 0 (t) dt,
C a
Z Z b
f (x, y, z) dz = f (x(t), y(t), z(t))z 0 (t) dt.
C a
Z
Example 16.7. Evaluate y sin z ds, where C is the circular helix r(t) = hcos t, sin t, ti,
C
t ∈ [0, 2π].
Solution. √ Z
Z Z 2π 2 2π
q
2
y sin z ds = (sin t)(sin t) sin t + cos2 t + 1 dt = (1 − cos(2t)) dt
C √0 2 0
2 1
2π
√
= t − sin(2t) = 2π.
2 2 0
= P dx + Qdy + Rdz.
a
Sometimes, it is helpful to think of F · dr as hP, Q, Ri · hdx, dy, dzi = P dx + Qdy + Rdz.
Theorem 18.1. Let C be a smooth curve given by r(t), t ∈ [a, b]. Let f be a function of 2 or
3 variables whose gradient ∇f is continuous. Then
Z
∇f · dr = f (r(b)) − f (r(a)).
C
Proof.
.
...........
.....
. .
............ ............ ∇f
...... ...........................................................
..... . ........ .. ...
...... ..............
. . .
................... ........... .......... .......
.............
..
...... ...... ....• r(b) ..
.
. . ..... ... ..
. . .. . ... .
......... ... ......... ........ ........ ............. ..........
.... . .... .... . .. .
. . ...... ...... . ......... ......
.
......... ..
.. ..
. .... ..
.
... ........... ...........
. ........................ ...........
... ... ...... ..... ...... .
......
... ..... .
.... . . .
...
. ......
........... ........... . .
..... ..........
.. . . ..
.
r(a) • .......
..........
.... .... ............... .....
...............................
..
Figure 103
Z Z
∇f · dr = ∇f (r(t)) · r0 (t) dt
C C
Z b ∂f
dx ∂f dy ∂f dz
= + + dt
Zab ∂x dt ∂y dt ∂z dt
d
= f (r(t)) dt by Chain Rule
a dt
=f (r(b)) − f (r(a)). by Fundamental Theorem of Calculus
c Department of Mathematics 76
Semester II(2005/06)
mM G
Example 18.2. Consider the gravitational (force) field F(r) = − r, where r = hx, y, zi.
|r|3
mM G
Recall that F = ∇f , where f (x, y, z) = p . Find the work done by the gravita-
+ y2 + z2 x2
tional field in moving a particle of mass m from the point (3, 4, 12) to the point (1, 0, 0) along
a piecewise smooth curve C.
Z Z
Solution. W ≡ F · dr = ∇f · dr = f (1, 0, 0) − f (3, 4, 12) = 12mM G/13.
C C
C 1
.............................................
C 1
.............................................
.......... ........ .......... ........
....... ...... ....... ......
....... ...... ....... ......
.....
.
.....
.
..
•B ...
... .....
.
.....
.
..
•B ...
...
...
. ...
. ...
. ...
.
...
. ... ...
. ...
... .. ...
.
...
... ... ...
.......... ...
... .... ... .........
..
. ......
. ..
. ........
.
... ..... ... .....
... ...... ... ......
...... ......
... ....... ... .......
........ ........
A•.......
..........
...............................
...... .
. ..
...
C2 A•.......
..........
...............................
...... .
. ..
...
C2
Z Z Z Z Z
We have F · dr = F · dr + F · dr = F · dr − F · dr = 0, since both C1 and
C C1 C2 C1 −C2
−C2 are paths from A and B and the line integral is by assumption independent of path.
To prove the sufficiency, consider 2 paths C1 and C2 having the same initial point
Z A and
terminal point B. See figure 105. Then C = C1 − C2 is a closed path. Thus, 0 = F · dr =
Z Z Z Z Z C
F · dr = F · dr − F · dr. Hence, F · dr = F · dr.
C1 −C2 C1 C2 C1 C2
By 18.3 and the fundamental theorem for line integrals, we have the following implication and
equivalence: (1) =⇒ (2) ⇐⇒ (3).
In fact, the implication (2) =⇒ (1) is true with some suitable assumptions on the domain D.
Definition 19.4. A subset D in R2 (or R3 ) is said to be open if for any point p in D, there
is a disk (ball) with center at p that lies entirely in D. (This means D does not contain any
boundary points.)
..........
........ .............
..... ....
p ... ...
..
. ...
..
.............................................................................................
............................................... ....................................................... •
.................................................................................................................................................
.
.................................................. .........................................................
................................................................... ...........................................................
........................................................ ........................................................................
..................................................... . .
............................................... ...... p• ..........................................................................
...................................................... ..................................................
............................................................ ..................................................
.............. ... ... ... ............................................. ............................................................
............ ............. ...............................
............ ............................................. ............. ...............................
............ .............................. ............. ...............................
.............. ... ... ... ............................................... ............. .. . . . . . . . . . . . . ..
.. . . . . . . . . . . . . . . . . . . . . . . . ...........................................................................
. . . . . . . .
.. . . . . . . . . .............................
. ..................................................
....................................................................... ..................................................
..................................................
.. . . . . . . . . . . . . . . . . . . . . . . .
.............................................................................................
D D
...........................................................................................................................................
Definition 19.5. A subset D in R2 (or R3 ) is said to be connected if any two points in D can
be joined by a path that lies in D.
....................................................................................................................................................
...............................................
................................................
................................................
.. . . . . . . . ... . . . . . . . . . . . . . . .
........................................................................................................
.............................................................
.........................................................
......................................................
............... ... ... ... .................................. ....................................... .........................................................
............... .
................................................. .................... .......................
•
...............
............ ..............................
.............. ... ... ... ...............................................
• .....................
.....................
. .... . . . . . . . . ..
......................
..................................................................................................................
.. . . . . . . . . . . . . . . . . . . . . . . .
................................................
.................... ..
....................
• .
.
. .. . . . . . . . . .
..........................
•
....................................................................... ..................... .......................
.. . . . . . . . . . . . . . . . . . . . . . . . D
.............................................................................................
.....................
.
... ... ... ... ... ... ... ... ... .
.......................
........................................................
Theorem
Z 19.6. Suppose F is a vector field that is continuous on an open connected region D.
If F · dr is independent of path in D, then F is conservative. That is there exists a function
C
f such that ∇f = F.
Proof. Let’s prove the case in R2 . The case in R3 is similar. Fix a point A(a, b) in D. Let
Z (x,y)
f (x, y) = F · dr,
(a,b)
Z
where (x, y) ∈ D. Since F · dr is independent of path in D, f is well-defined. As the domain
C
D is open, there exists a disk centered at (x, y) that lies entirely in D. Pick a point (x1 , y)
in the disk with x1 < x. Let C consist of any path C1 from (a, b) to (x1 , y) followed by the
horizontal line segment C2 from (x1 , y) to (x, y). See figure 111. Then
Z Z Z (x1 ,y) Z
f (x, y) = F · dr + F · dr = F · dr + F · dr.
C1 C2 (a,b) C2
∂f ∂
Z
Thus, =0+ F · dr because the first line integral along C1 is independent of x.
∂x ∂x C2
c Department of Mathematics 78
Semester II(2005/06)
Z Z
Let’s write F = P i + Qj. Then F · dr = P dx + Qdy. On C2 , y = constant so that
Z C2 C2
Qdy = 0. Hence,
C2
∂f ∂
Z
∂
Z x
= F · dr = P dx = P (x, y).
∂x ∂x C2 ∂x x1
.........................................
........... ........
........ .......
....... ......
.... ...... .....
.....
..
..... .....
....... ...
...
.............................
. ....
... . .
.......................... ...
...
. .
.... ..............................................
. ...
... ......................................... ...
..
.
...
...
.
. •
............................. . . . . . ..
........................................................
...
...
...
..
. ...... .........................
. (x, y) ...
.... C ..
....
. ...........................
. . ...
.... ........................ ...
... .. .
....
. . .
............ .. .... ......
................ . .................
... .......... ...
... ......................... ........
.... ... .....................................
... .... ...
. .
.............................................................
.. ... .. . . . . . . . . . . . . . . . . . ..
...
..
.
.. .......................................................... ...
... . . ... ....................................... ...
...
... ...
.... ... ............................................................ ....
... .
.. .. .... ................................................................. ....
C
...
. ... (x , y) .. . . . . . . .2. . . . . . . . . . . . . . ..
...
...
... • ...
...
...
1
• •
.......................................................................................
.......................................................
..... ... ........................................... ...
...
A(a, b)
... .... . (x, y)
.... ......................................................................
... ....................................... ...
... ... ... ... ..................................... ..
... ..
... ... ...........................................
....
..... ..... ... ........................................
........................................
..... .
....
... .... .....................................
.....
D ..... ..................................
......
....... ...... C ................ . ................
......
........ ....... 1
........... .. ... .. ... ........
..............................
Figure 111
Figure 110
Similarly, by considering a path from A(a, b) to a point (x, y1 ) with y1 < y inside the disk
∂f
followed by the vertical segment from (x, y1 ) to (x, y), we can prove that = Q(x, y). There-
∂y
fore,
∂f ∂f
F = P i + Qj = i+ j = ∇f.
∂x ∂y
The openness of D is to ensure the points (x1 , y) and (x, y1 ) exist corresponding to every (x, y)
in D.
Theorem 19.7. Let F(x, y) = P (x, y)i + Q(x, y)j be a vector field on D ⊂ R2 , where P and Q
have continuous partial derivatives in D. If F is conservative, then
∂P ∂Q
= .
∂y ∂x
Proof. As F is conservative in D, there exists a differentiable function f (x, y) in D such that
∇f = F. That is fx = P and fy = Q. By Clairaut’s Theorem,
∂P ∂2f ∂2f ∂Q
= = = .
∂y ∂y∂x ∂x∂y ∂x
Definition 19.8. A simple curve is a curve which does not intersect itself.
c Department of Mathematics 79
Semester II(2005/06)
............................
• ............. . •
.. ........
.. .......
..
.. . ........
.... ...
. . . ...
..
. .. . .
.. ........ .......... .....
....
...
...
....
....
...
....
.....
... ...
•
...
. .. ... ....
... ...
. ... ...
... ... ...
..... ...
....... ... ... ...
.
..
...
.................................................. ..
. ..
. ..
..
..
...... ..
. .. ..
.... .
. ... ..
... .
.. ...
. ..
... . . .
.... ... ... ...
... ... ... ..
... ... ..
. ...
... .
... ..
.. ... ...
.... C .. C ... ...
.
..
1 ..
..
2 ....... .........
.........
•
..
•
..
Figure 112
Definition 19.9. A simply-connected region in the plane is a connected region such that every
simple closed curve in D encloses only points that are in D.
................................................................................................................................................... ....................................................................................................................................................................................................
........................................................................................... .........................................................................................................................
............................................................... ...................................................................................
............................................................................................ .........................................................................................................................................
........................................................................................................... .................................................................................................................................................
..................................................................................................................... ........................................................................................................................................
............................................................................. ................................................................................................
................................................................................................... .........................................................................................................................................................................
................................................................................................. ..................................... .....................................
..................................................................... ........................... ............................
.............................................................................................. .................................... ...................................
...................................................................................................... ..................................... .........................................
.................................................................. ............................ ...........................
.................................................................................................. ..................................... .........................
.................................................................................................... ..................................................................................................................................................................
............................................................................... ...................................................................................................................................
............................................................................................................ ................................................................................................
................................................................................................................ ...............................................................................................................................................
........................................................................................... ................................................................................................................................................
D ..............................................................
....................1 ................................................. D ...................................................................................
.....................2
...................................................................
................................................................................................................................................................................................ . . . . . ....................................................................
..........................................................................................................................................................................................
Figure 113
Theorem 19.10. Let F(x, y) = P (x, y)i+Q(x, y)j be a vector field on an open simply-connected
∂P ∂Q
region D ⊂ R2 , where P and Q have continuous partial derivatives in D. If = , then
∂y ∂x
F is conservative.
This is a consequence of Green’s Theorem in the next section.
Example 19.11. Determine whether the vector field F(x, y) = (3 + 2xy)i + (x2 − 3y 2 )j is
conservative.
Solution. As
∂(x2 − 3y 2 ) ∂(3 + 2xy)
= 2x = ,
∂x ∂y
and the domain of F is R2 which is open and simply-connected, F is conservative by the
Theorem 19.10.
2 2
Z Let F(x, y) = (3 + 2xy)i + (x − 3y )j. Find a function f such that ∇f = F.
Example 19.12.
Also evaluate F · dr, where C is the curve given by r(t) = et sin ti + cos tj, t ∈ [0, π].
C
c Department of Mathematics 80
Semester II(2005/06)
Exercise 19.13. If F(x, y, z) = y 2 i + (2xy + e3z )j + 3ye3z k, Find a function f such that
∇f = F.
[Answer: f (x, y, z) = xy 2 + ye3z + C]
................................................................................................................................................................................ ................................................................................................................................................................................
.............................................................. ..............................................................
. . . . . . . . . . . . . . . . . . . . . . . . . . . . .. . . . . . . . . . . . . . . . . . . . . . . . . . . . . ..
............................................................................................. ..............................................................................................
............................................................................................. .............................................................................................
........................................................................................... ...........................................................................................
........................................................................................... ............................................................................................
............................................................... ...............................................................
........................................................................................... ...........................................................................................
............................................................................................. ............................................................................................
.............................................................. ...............................................................
........................................................................................... ...........................................................................................
............................................................................................. ............................................................................................
D
..............................................................
............................................................................................
D
.............................................................
.............................................................................................
.................................................................................................
.............................................................. C .
..................................................................................................
..
........................................................ . ...
−C
........................................................................................... ...........................................................................................
............................................................................................. ............................................................................................
.............................................................. ...............................................................
........................................................................................... ...........................................................................................
............................................................................................ ............................................................
............................................................................................................ ............................................................................................................................................
Figure 114
In other words, if D is a simply connected region on the plane, then the boundary C of D
oriented in the counterclockwise sense is the positive orientation. We use the notation ∂D to
denote the boundary of D with the positive orientation.
Theorem 20.1. (Green’s Theorem) Let C be a positively oriented, piecewise-smooth, simple
closed curve in the plane and let D be the region bounded by C. If P (x, y) and Q(x, y) have
continuous partial derivatives on an open simply connected region that contains D, then
∂Q ∂P
Z ZZ
P dx + Qdy = − dA.
C D ∂x ∂y
Z
The line integral P dx + Qdy has other notations as
C
Z I Z
P dx + Qdy, P dx + Qdy, or ..............
........ P dx + Qdy.
C C ∂D
c Department of Mathematics 81
Semester II(2005/06)
They all indicate that the line integral is calculated using the positive orientation of C.
Proof. We shall first verify that Green’s Theorem is true for D being a region which is both
of type I and Type II. (See figure 71 for type I and type II regions.) The general case can be
proved by cutting the region into a finite number of regions of both type I and type II and
applying the result to each of them. Let’s consider a type I region. The proof for type II region
is similar.
y.
...
..........
C = ∂D
...
... ........... • ......................................
...........
y = Y (x) .... ..............
...
...
C .....
.......... 3 ...
...
......
.....
.....
2
... .. .. ... ....
... .... .
. ....
....
... .. ..
...
... ..
..
.
..
.
....
...
• ....
..
...
... .. ... ...
... ..
. .
. ...
...
... • ....
.
.
..
D ..
.
..
.
..
C .....
......
..
2
... . .
... ..
. ..
. ....
... .. .. ..
C ...
... 4 ............
.
. .
•
....
... ..
..
.
... ..
... ... ..
... ..
. .
. ...
• ...
...
......
.....
........
....
...
....
.....
.
....
... ................ ... ......
...
...
C ..........
.............1
.
. .
.. .........
..
...
...
...
• y = Y (x) ......................................
...
1
... ..
... .
..............................................................................................................................................................................................................
a ...
.. x
x
b
Figure 115
In this case, the lower and upper boundaries of D consists of simple smooth curves C1 and C3
respectively, and the left and right boundaries are vertical lines C4 : x = a and C2 : x = b.
See figure 115. Let C1 and C3 be parametrized as the graphs of y = Y1 (x) and y = Y2 (x)
respectively for x ∈ [a, b]. Here we assume Y2 (x) > Y1 (x) for all x ∈ (a, b) so that C3 is higher
than C1 . Thus C1 is given the orientation which goes from left to right, whereas C3 is given
the orientation which goes from right to left so that C = C1 + C2 + C3 + C4 . Then
"Z #
ZZ
∂P
Z b Y2 (x) ∂P
dxdy = dy dx
D ∂y a Y1 (x) ∂y
Z b
y=Y (x)
= [P (x, y)]y=Y21 (x) dx
Zab
= P (x, Y2 ) − P (x, Y1 ) dx
aZ
b Z a
=− P (x, Y1 ) dx − P (x, Y2 ) dx
Za Z b
=− P dx − P dx
ZC1 ZC3 Z Z
=− P dx − P dx − P dx − P dx
ZC1 C2 C3 C4
=− P dx.
C
Z Z
Note that P dx and P dx are in fact zero because C2 and C4 are vertical segments:
C2 C4
x = b and x = a, so that dx = 0.
c Department of Mathematics 82
Semester II(2005/06)
∂Q
ZZ Z
Similarly, using the fact that D is also a type II region, dxdy = Q dy. Therefore,
D ∂x C
∂Q ∂P
I ZZ
P dx + Qdy = − dA.
C D ∂x ∂y
Now to extend this result to the general case, first consider a region D which is a union of two
regions D1 and D2 meeting along a curve L in their common boundaries. See figure 116. Thus
∂D1 = L1 + L, ∂D2 = −L + L2 and ∂D = L1 + L2 . Suppose Green’s Theorem holds for the
regions D1 and D2 . Then
∂D.......................... L1.............. ............L
..
2
........ . ............. ..... . . .............
.......... .... ........ ........... .... .... ........
....... ....... ........ .......
...... ... ...... ...... ... ... ......
..... ... ..... ..... ... .. .....
....
. .
. ... .....
. ... ... ...
.. .. ... ... ...
. . ... . L ... −L ... ...
.... .... ... ........ ... ........ ...
.... ... .. .... . . .... ..
D ...
...
D1 ...
.
.
D2 .........
..
...
...
...
D1 .
........
....
....
..
D2 .........
.
.
... .
.
. ... ... ... ... ...
....
.... ..... ..... ... ... ...
.....
..... ..... ....
...... ... ...... ...... ... ... .....
.......
....... ........ ... ......
.......... ..
...................................................... ........... .. ........
.................... .............................
Figure 116
Then, suppressing the terms involving P , Q etc, the following calculation shows that Green’s
theorem holds for the region D.
Z Z Z Z Z Z Z Z
= = + = + + +
∂D ZL1 +L2 Z L1 L2Z L1 Z L −L L2
= + = +
L1 +L ZZ−L+L2ZZ
ZZ ∂D1 ∂D2
= + = .
D1 D2 D
Now any simple closed curve in the plane bounds a region which can be cut into regions both
of type I and type II. See figure 117. Thus, by the above consideration, Green’s Theorem is
valid for any simple closed curve in the plane.
.........................................
........... ........
........ ... ......
...... ... ......
...
...... ... .....
....... .
.
......
........
..
.
...
. D4 ...............
..
.
.
. ...
...
...
.........
D3 ...
...
D
.... ...
. . ...
.
.........
...
...
... ....... .
.... ...
.... .... .... ...
... ..
. .
... ...
... .... ... ...
..................................... ....................................
... ... .. ..
... ... ... ..
.
... ... .. ..
... ...
. ..... ...
.
..... .
... ...... ..... ...
... ......... ..... ...
...
... ......................... ...
... ...
...
.....
.....
D1 .
.
....
D2 .....
.....
. .
..... ... .....
...... ......
....... ... ......
......... . ........
........................................
Figure 117
Lastly, the proof of Theorem 19.10 follows from Green’s Theorem and Theorem 19.6.
Z
Example 20.2. Evaluate x4 dx + xydy, where C is the triangular curve consisting of the
C
line segments from (0, 0) to (1, 0), from (1, 0) to (0, 1) and from (0, 1) to (0, 0).
Solution. The functions P (x, y) = x4 and Q(x, y) = xy have continuous partial derivatives on
the whole of R2 , which is open and simply connected.
c Department of Mathematics 83
Semester II(2005/06)
y.
..
.........
...
....
..
(0, 1) ...
..
............
............
. ..
....................... C
.................
...................
. ...............................
.
............................... y =1−x
............................
..........................
................................................
D
..................................
.....................................
. . . . . . . . . . . . . . . ..
·
........................................
(0, 0)
....
..
· (1, 0)
..................................................................................................................................................................
x x
Figure 118
ZZ " #
∂(xy) ∂x4
Z
4
By Green’s Theorem, x dx + xydy = − dA
C ZZ D ∂x ∂y
= y dydx
ZD
Z1 1−x
= y dydx
0 0
= 16 .
I q
Example 20.3. Evaluate (3y − esin x )dx + (7x + y 4 + 1)dy, where C is the circle x2 +y 2 =
C
9, oriented in the counterclockwise sense.
Solution. C bounds the circular disk D = {(x, y) | x2 + y 2 ≤ 9} and is given the positive
orientation. By Green’s Theorem,
ZZ " p #
∂(7x + y 4 + 1) ∂(3y − esin x )
I q
sin x 4
(3y − e )dx + (7x + y + 1)dy = − dA
C ZZ D ∂x ∂y
= 4 dA
D
= 4(π32 ) = 36π.
20.1. Application of Green’s Theorem to Find Area. Recall that the area of a region
∂Q ∂P
ZZ
D in R2 is 1 dA. Therefore, if we choose P (x, y) and Q(x, y) such that − = 1,
D ∂x ∂y
then by Green’s Theorem we have
ZZ I
Area of D = 1 dA = P dx + Qdy.
D ∂D
There are various choices of P and Q that satisfy this requirement. For example:
Therefore,
1
I I I
Area of D = x dy = − y dx = xdy − ydx.
∂D ∂D 2 ∂D
x2 y 2
Example 20.4. Find the area of the ellipse + 2 = 1.
a2 b
c Department of Mathematics 84
Semester II(2005/06)
Solution. Let the parametric equations for the ellipse be x = a cos t, y = b sin t, for t ∈ [0, 2π].
Then, Area = 12 ∂D xdy − ydx = 12 02π (a cos t)(b cos t) − (b sin t)(−a sin t) dt = 12 02π ab dt =
H R R
πab.
I
Exercise 20.5. Evaluate by Green’s Theorem e−x sin y dx + e−x cos y dy, where C is the
C
rectangle with vertices at (0, 0), (π, 0), (π, π/2), (0, π/2), oriented in the counterclockwise sense.
20.2. Non Simply-Connected Regions. Green’s Theorem is also valid for non simply-
connected regions, that is regions with holes. Consider the region D in figure 119 in which
∂D = C1 + C2 . We may cut the region D by two line segments L1 and L2 into two simply
connected regions D0 and D00 .
∂D = C = C1 + C2
........................................................................................................................................................
............................... ..............................................................................................................................................................................................
............................................................................................... ..................................................................
...............................................................
...............................................................................................
D0
..................................................................................................
.................................................................................................
................................................................................................................................ L
..............1 ........................................................................................................2 L .......
..................... ...................... ........................................ .......................................
................ .................................... ....................................
.................
...................... C2
.............................................................................................................................
............... ..................................... C2
........................................................................................................................
.....................................
................................................................... .......................................................................
............................................................................................... .................................................................................................
D
.............................................................................................
................................................................. D00
................................................................................................
....................................................................
............................................................................................................................................................. ...........................................................................................................................................................
C1 C1
Figure 119
∂Q ∂P ∂Q ∂P ∂Q ∂P
ZZ ZZ ZZ
Then − dA = − dA + − dA
D ∂x ∂y Z D 0 ∂x ∂y Z D 00 ∂x ∂y
= P dx + Qdy + P dx + Qdy
Z∂D0 Z ∂D00
= P dx + Qdy + P dx + Qdy
ZC1 C2
= P dx + Qdy
ZC1 +C2
= P dx + Qdy.
C
Here the third equality is obtained by re-grouping the line integrals and cancelling the line
integrals along L1 and L2 .
−y x
Z
Example 20.7. Let F(x, y) = 2 2
i+ 2 j. Show that F · dr = 2π for every simple
x +y x + y2 C
closed curve that encloses the origin.
Solution. Note that the vector field F is defined on R2 \{(0, 0)}. Let C be any closed curve that
encloses the origin. Choose a circle C 0 centered at the origin with a small radius a such that
C 0 lies inside C. We can parametrize C 0 by x = a cos t, y = a sin t, t ∈ [0, 2π]. Let D be the
region bounded between C and C 0 . We give both C and C 0 the counterclockwise orientation.
Thus ∂D = C − C 0 is given the positive orientation with respect to the region D.
c Department of Mathematics 85
Semester II(2005/06)
y.
....
.........
.......................... ..
........ .......... ...
.. .............. ..
.... ................
. . . ...............
... ........................................................
...
...
...
...
...
C
..........
.....
...
... ...
C0
..
... ..
. ..
.. ... .. .. .. . ..
...
...
.. ...
..
.. .
....... ... ........ .
..
.
.. .
. . .
.
. . . ..
.. . . .. ...
............................................................................................................................................................................
...
...
...
...
..
...
.......
...
.. x
... ..... .
. ..
...
. .
..
... ......... .. ......... .
... ...... ...
....
.... D ... ....
.... ... ....
.... .
. .
.....
.
..... .
..... .... .....
...... .. .....
....... .....
.......... .... .................
.............................
..
Figure 120
Therefore, F · dr = F · dr
C 0
ZC2π
= F(r(t)) · r0 (t) dt
Z02π
(−a sin t)(−a sin t) + (a cos t)(a cos t)
= dt
0 (a2 cos2 t + a2 sin2 t)
= 2π.
∂R ∂Q ∂P ∂R ∂Q ∂P
curl F = − i+ − j+ − k.
∂y ∂z ∂z ∂x ∂x ∂y
The curl of a vector field F is a vector field which measures the rotational effect of F. The
geometric meaning of curl F can be seen after we learn Stokes’ Theorem. At this point, let’s
introduce the del operator ∇. We let
∂ ∂ ∂
∇=i +j +k .
∂x ∂y ∂z
We regard ∇ as a 3-dimensional vector consisting of the operators of partial differentiations
with respect to x, y, z. We can multiple ∇ by a scalar function (on the right), take the dot
product with a function, or the cross product with a vector field. For example, we may regard
the gradient of a function f as being the scalar multiplication of ∇ and f . That is
∂f ∂f ∂f
grad f = ∇f = i +j +k .
∂x ∂y ∂z
The curl of a vector field F = P i + Qj + Rk can be regarded as the cross product between ∇
and F.
i j k
∂ ∂ ∂
curl F = ∇ × F = ∂x ∂y ∂z
P Q R
∂R ∂Q ∂P ∂R ∂Q ∂P
= − i+ − j+ − k.
∂y ∂z ∂z ∂x ∂x ∂y
c Department of Mathematics 86
Semester II(2005/06)
Theorem 21.2. If f (x, y, z) has continuous 2nd order partial derivatives, then curl (∇f ) = 0.
Proof.
i j k
∂ ∂ ∂
curl ∇f = ∂x
∂y ∂z
∂f ∂f ∂f
∂x ∂y ∂z
! ! !
∂2f ∂2f ∂2f ∂2f ∂2f ∂2f
= − i+ − i+ − i
∂y∂z ∂z∂y ∂z∂x ∂x∂z ∂x∂y ∂y∂x
= 0. by Clairaut’s Theorem
Corollary 21.3. If F is conservative (i.e. F = ∇f ), then curl F = 0.
Remark 21.4. If F = P i + Qj is a vector field on R2 , we may regard F as the vector field
∂Q ∂P
F = P i + Qj + 0k in R3 with zero k component. Then curl F = − k. Thus in this
∂x ∂y
∂Q ∂P
case, if F is conservative, then = which is 19.7.
∂x ∂y
For example F(x, y, z) = xzi + xyzj − y 2 k is not conservative because curl F = −(2y + xy)i +
xj + yzk 6= 0.
Using Stokes’ Theorem in the next section, one can prove the following:
Theorem 21.5. Let F be a vector field on R3 whose component functions have continuous
partial derivatives. If curl F = 0, then F is conservative.
Exercise 21.6. Show that the vector field F = x2 i + y 2 j + z 2 k is conservative. Find a function
f such that ∇f = F.
∂ ∂ ∂
= i +j +k · (P i + Qj + Rk) .
∂x ∂y ∂z
= ∇ · F.
Proof.
∂ ∂R ∂Q ∂ ∂P ∂R ∂ ∂Q ∂P
div curl F= − + − + −
∂x ∂y ∂z ∂y ∂z ∂x ∂z ∂x ∂y
= 0.
because the terms cancel in pairs by Clairaut’s Theorem.
For a velocity vector field F, div F measures the amount of flow radiating at a point. If the
flow is uniform and without compression or expansion, then div F = 0. Thus, if div F = 0, we
say that F is incompressible. Whereas curl F measures the rotational effect of the vector field
F. Therefore, if curl F = 0, then we say that F is irrotational.
............................. ...... ................ ................ ................ ................ ................ ................ ................
........
......................... ...... ......... ................ .............. ................ .............. ................ .............. ..............
... .
..
.. .
... . . . ..
..
..
......... .. ........... .
.
.. . .. ........ .
.......... ..... ... ..
.... .............. ...... ................ ...... ... .. ....... ................ ................ ................ ................ ................ ................ ................
. ........ .... ........ . ...... ..
. .......
.. .
. .. .......... ...... ....
.... .... .... ........
... ... ....... ....
.. ........ ....
. .
· .......... .... .... ....
.
........... ......
................ ................ ................ ................ ................ ................ ................
........ .... . .
..... . ... ..
... ....... .... ................ ....... .............. ........ ............... ................ ............... ................ ............... ................ ................
..... ...... .................. ..... ...........
. ....
. .. ....
..... ............. .
..... .. .
....................... ................ ................ ................ ................ ................ ................ ................
..............
. ............................ ................ ................ ................ ................ ................ ................ ................
curl F 6= 0 curl F = 0
irrotational
div F = 0
div F 6= 0
incompressible
Figure 121
Another differential operator occurs when we compute the divergence of a gradient vector field
∇f . If f is a function of three variables, we have
∂2f ∂2f ∂2f
div (∇f ) = ∇ · (∇f ) = + + .
∂x2 ∂y 2 ∂z 2
We abbreviate this expression as ∇2 f . The operator ∇2 = ∇ · ∇ is called the Laplace operator
because of its relation to Laplace’s equation:
∂2f ∂2f ∂2f
∇2 f = + + = 0.
∂x2 ∂y 2 ∂z 2
We can also apply the Laplace operator ∇2 to a vector field F = P i + Qj + Rk in terms of its
components:
∇2 F = ∇2 P i + ∇2 Qj + ∇2 Rk.
c Department of Mathematics 88
Semester II(2005/06)
p
Exercise 21.9. Let r = x2 + y 2 + z 2 . Find ∇2 (r3 ).
Exercise 21.10. Prove that div(f F) = f divF + F · ∇f.
21.1. Green’s Theorem in Vector Forms. Let F = P i + Qj be a vector field. Green’s
Theorem says that
∂Q ∂P
Z ZZ
F · dr = − dA.
∂D D ∂x ∂y
∂Q ∂P
3
Regard F as a vector field in R . That is F = P i + Qj + 0k. Then curl F = − k,
∂x ∂y
∂Q ∂P
so that − = (curl F) · k. Therefore, we may state Green’s Theorem in the following
∂x ∂y
form: Z ZZ
F · dr = (curl F) · k dA.
∂D D
To get a better meaning of this equation, let ∂D be parametrized by the vector equation
r(t) = hx(t), y(t)i for t ∈ [a, b].
We assume the pararmetrization gives the positive orientation of ∂D. Then the unit tangent
vector is
r0 (t) x0 (t) y 0 (t)
T(t) = 0 =h 0 , i.
|r (t)| |r (t)| |r0 (t)|
Thus
r0 (t) 0
Z Z b Z b
F · dr = F(r(t)) · r0 (t) dt = F(r(t)) · 0 |r (t)| dt
∂D a a |r (t)|
Z b Z
= (F(r(t)) · T(t)) |r0 (t)| dt = F · T ds,
a ∂D
where ds = |r0 (t)| dt is the arc length differential. Then, we can also state Green’s Theorem in
the following form:
Z ZZ
F · T ds = (curl F) · k dA.
∂D D
In this form, the equation expresses the line integral of the tangential component of F along
∂D as the double integral of the vertical component of curl F over the region D. This is a
special case of Stoke’s Theorem in which D is not necessarily a planar region but is a surface
in R3 with a boundary curve ∂D.
We could also derive a formula involving the normal component of F along ∂D. In that way,
Greens’ theorem will be stated in terms of the divergence of the vector field F. Using the above
parametrization of C, one can easily verify (by taking dot product with T) that the outward
unit normal vector to ∂D is given by
y 0 (t) x0 (t)
n(t) = h , − i.
|r0 (t)| |r0 (t)|
(It is the outward pointing normal because C is given the positive orientation.) Now we
consider the line integral of the normal component of F along ∂D.
c Department of Mathematics 89
Semester II(2005/06)
Z Z b
F · n ds = (F(r(t)) · n(t)) |r0 (t)| dt
∂D Zab
P (x(t), y(t))y 0 (t) Q(x(t), y(t))x0 (t)
= − |r0 (t)| dt
Zab
|r0 (t)| |r0 (t)|
= P (x(t), y(t))y 0 (t)dt − Q(x(t), y(t))y 0 (t)dt
Za
= P dy − Qdx
C
∂P ∂Q
ZZ
= + dA by Green’s Theorem
ZZD ∂x ∂y
= divF dA.
D
Therefore, Z ZZ
F · n ds = divF dA.
∂D D
This version says that the line integral of the normal component of F along ∂D is equal to
the double integral of the divergence of F over the region D. This result can be generalized to
the case of closed surface enclosing a solid region in R3 which is the content of the Divergence
Theorem.
In the next two exercises, we assume D satisfy the hypotheses of Green’s Theorem and the
appropriate partial derivatives of f and g exist and are continuous. The first exercise is a
consequence of 20.10.
ZZ Z ZZ
2
Exercise 21.11. Prove that f ∇ g dA = f (∇g) · n ds − ∇f · ∇g dA.
D ∂D D
ZZ Z
Exercise 21.12. Prove that (f ∇2 g − g∇2 f ) dA = (f (∇g) − g(∇f )) · n ds.
D ∂D
Figure 122
c Department of Mathematics 90
Semester II(2005/06)
Example 22.2. Identity the surface with vector equation r(u, v) = h2 cos u, v, 2 sin ui.
Solution. The point (x, y, z) = (2 cos u, v, 2 sin u) lies on this surface if and only if x2 + z 2 =
4 cos2 u + 4 sin2 u = 4. Therefore, the surface is the cylinder x2 + z 2 = 4. The domain of r can
be taken as the infinite strip D = {(u, v) | 0 ≤ u ≤ 2π, −∞ < v < ∞}. The function r simply
identifies the two vertical sides of this strip to form the cylinder. Here we omit the line u = 2π
so that r is injective. We could take the domain of r to be the whole xy-plane. In that case, r
takes the whole xy-plane and wraps it up around the cylinder infinitely many times.
v.. z..
.. ..
......... .........
.. ..
.................. ...
.................. ...
.................. ...
.................. ...
..................
..................
r(u, v) ...
...
.................. ...
.................. ...
.......................................................................... ...
D ..................
.................. ... x +z =4 2 2
.................. ..
.................. ..
............................................................................................
. . . ...
................. ... .. ... ... ...
................................................................................................................. ... .... ... ...
................. u ... ... .. ..
..
...
...
2π
.................. .... ... . ..
..................
.................. ........... .
.
.
... ............................................................................................................................................................
... ..
y
.................. .
...
......... .
. ... .
...
.. ... ...
....... .. ..
.
...
......... ...
..... .
.
. ...
.. ...
.
....... ..........................................................................................
...........
x
Figure 123
Example 22.3. Find a vector function that represents the plane that passes through the point
P0 with position vector r0 and contains two non-parallel vectors a and b.
Solution. Let O be the origin. For any point P on the plane, its position vector r can be
expressed as
−→ −→
r =OP0 + P0 P = r0 + ua + vb,
for some numbers u and v.
..
...
... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... .............
.......... ·P
..... .
....... ..
....... ..
.. ....... ..
vb .
...
..
.
.
.......
.......
.....
.....
..
..
...
. .
...
........ .
..
.. .......... .
.......... ........ ..
....... ..
b ..
.
...
.
. . .. . .
.......
....... ..
..
...
... .......... ..
... ........ ..
... ...... ..
............................................................................................................
P0 a ua
Figure 124
Example 22.5. Find a vector function that represents the elliptic paraboloid z = x2 + y 2 .
c Department of Mathematics 91
Semester II(2005/06)
In general, if the surface S is the graph of a function z = f (x, y), then a natural parametric
representation of S is
r(x, y) = hx, y, f (x, y)i.
p
Example 22.6. Find a parametric representation of the cone z = 2 x2 + y 2 .
p
Solution. Since the cone is the graph of the function z = 2 x2 + y 2 , we can simply take
q
r(x, y) = hx, y, 2 x2 + y 2 i.
p
Alternatively, we can consider cylindrical coordinates. The equation of the cone z = 2 x2 + y 2
in cylindrical coordinates is z = 2r. Thus if we use polar coordinates (r, θ) of the xy-plane, we
can write r(r, θ) = hr cos θ, r sin θ, 2ri.
Figure 125
Similarly, the image of the vertical line u = u0 under r is a curve C2 whose tangent vector at
P0 is given by
∂x ∂y ∂z
rv ≡ h (u0 , v0 ), (u0 , v0 ), (u0 , v0 )i.
∂v ∂v ∂v
Both vectors ru and ru lie in the tangent plane to S at P0 . Therefore, the cross product ru ×rv ,
assuming it is nonzero, provides a normal vector to the tangent plane to S at P0 . Therefore,
(r−r0 )·(ru ×rv ) = 0 is the equation of the tangent plane. At this point, lets’ make a definition.
Definition 22.7. The surface S is said to be smooth if ru × rv 6= 0 for all points (x, y) ∈ D.
c Department of Mathematics 92
Semester II(2005/06)
Thus a smooth surface always has a tangent plane at each of its points. Basically, a smooth
surface is one which has no corners and no breaks.
Example 22.8. Find the equation of the tangent plane to the surface with parametric equations
x = u2 , y = v 2 , z = u + 2v at the point (1, 1, 3).
Solution. The vector equation of the surface is
r(u, v) = hu2 , v 2 , u + 2vi.
Therefore, ru = h2u, 0, 1i and rv = h0, 2v, 2i. Thus, a normal vector to the tangent plane is
ru × rv = h−2v, −4u, 4uvi. At the point (1, 1, 3), we have (u, v) = (1, 1). Then, the normal
vector at (u, v) = (1, 1) is h−2, −4, 4i. Therefore, the equation of the tangent plane to the
surface at (1, 1, 3) is (hx, y, zi − h1, 1, 3i) · h−2, −4, 4i = 0, or equivalently, x + 2y − 2z + 3 = 0.
22.3. Surface Area of the Graph of a Function. Let S be a surface which is the graph
of a function f (x, y) defined on a domain D ⊂ R2 . Then a parametric representation of S is
r(x, y) = hx, y, fq
(x, y)i. Thus, rx = h1, 0, fx i and ry = h0, 1, fy i so that rx × ry = h−fx , −fy , 1i,
and |rx × ry | = 1 + fx2 + fy2 . Therefore, the surface area of S over D is given by
ZZ q
A(S) = 1 + fx2 + fy2 dA.
D
This is the same formula derived in section 12.
c Department of Mathematics 93
Semester II(2005/06)
Exercise 22.10. Suppose S is the surface obtained by rotating the curve y = f (x) for x ∈ [a, b]
through an angle 2π about the x-axis, where f (x) ≥ 0, and f 0 (x) is continuous. Show that S
has a parametric representation given by r(x, θ) = hx, f (x) cos θ, f (x) sin θi, for x ∈ [a, b] and
θ ∈ [0, 2π]. Show that the area of S is given by
Z b q
2π f (x) 1 + [f 0 (x)]2 dx.
a
Exercise 22.11. On the xz-plane, the circle (x − b)2 + z 2 = a2 (b > a) is rotated through an
angle 2π about the z-axis to form a torus T . Find the surface area of T .
[Answer: 4π 2 ab]
22.4. Surface Integrals. Let S be a parametric surface with vector equation r(u, v) =
hx(u, v), y(u, v), z(u, v)i, where (u, v) ∈ D. Let f (x, y, z) be a continuous function defined on
S.
where 0 ≤ φ ≤ π and 0 ≤ θ ≤ 2π. From example 22.9, we have |rφ × rθ | = sin φ. Therefore,
ZZ ZZ
2
x dS = (sin φ cos θ)2 |rφ × rθ | dA
S D
Z 2π Z π
= sin3 φ cos2 θ dφdθ
Z0 π
0 Z 2π
3
= sin φ dφ cos2 θ dθ
0 0
= 4π/3.
ZZ
Example 22.14. Evaluate z dS, where S is the surface whose side face S1 is part of the
S
cylinder x2 + y2 = 1 bounded by the bottom face S2 which is the xy-plane and the top face
S3 which is part of the plane z = x + 1 above the xy-plane. See figure 126.
c Department of Mathematics 94
Semester II(2005/06)
z.. S3 : z = 1 + x
.
........ ........
... ........... .........
.......
... ....... ..
.. ........... ...
............. ......
... .....
.
.
. .
.
.
.... ... .
... ....
.
..... .... . .
... ...
..... ... ..
..... ...
... ....
..... ... ...
....... .
. .. ...
.. ....
.
..... ... ...
... ... ... ...
... .....
.
.... .
.
. . .
.... ...
.. .
. ... ...
.... ..
. ..
..... ... 2 2
.. . .
.
. .
. .........
.... S :x +y =1 ... 1
.... .
. ...
. .. ... ......
... .. ... ... ... ... ... ... . ............ .... ... ... ... ... ... ... .
.
.. ... . ........ ..... ....
...... ...... .... .....
..... ...... ...
......
...... ..................
.......
.....
...
.
.....
..
y
.......................................................................................................................
................. ......
......
..
........... ..
.. .....
. ..
.
...................... ................
.... ...................................
.....
..... S :z=0
......... 2
.
.....
.....
.....
..........
.......
x
Figure 126
Let’s first calculate the surface integral over S1 . The surface S1 is a cylinder. By example 22.2,
it has a parametric representation r(φ, z) = hcos θ, sin θ, zi, where 0 ≤ θ ≤ 2π, 0 ≤ z ≤ 1 + x =
1 + cos θ. Thus, rθ × rz = hcos θ, sin θ, 0i and |rθ × rz | = 1. Therefore,
ZZ Z 2π Z 1+cos θ Z 2π 1 3π
z dS = z dzdθ = (1 + cos θ)2 dθ = .
S1 0 0 0 2 2
ZZ
On S2 , we have z = 0. Thus the integrand of z dS is zero so that the integral has value
S2
zero. Therefore,
ZZ
z dS = 0.
S2
The surface S3 is the graph of the function z = 1 + x. Therefore, using polar coordinates, we
have ZZ ZZ q ZZ √
z dS= 2 2
(1 + x) 1 + zx + zy dA = (1 + x) 2 dA
S3 D D
Z 2π
Z 1 √ √
= (1 + r cos θ) 2 rdrdθ = 2π.
0 0
3π √
ZZ
Consequently, z dS = + 2π.
S3 2
ZZ
z 2 dS, where S is the portion of the cone z =
p
Exercise 22.15. Evaluate x2 + y 2 for
S
which 1 ≤ x2 + y 2 ≤ 4.
√
[Answer: 15π 2/2]
............................................................................................................................................
................ ...........
.......... ........
........ .......
..
........ .. ... ... ... ... ... ... ... ... ......... ... ... ... ... ... ... ... ... ...... ... ... ... ... ... ... ... ... ... .....
....
.
... . .. ... . ... ... .
.. ... ...
.
... ..... ...
..... ...... ................... ..... . . . . ..
.....
. ...................
. .....
. ......
...
...
...
..
..
...
............
....
... ...........
.. .........
..................... ... ....
............ ......... ..... ............. ...
...... ....... ..............
.... ....................
.......... ............. ....... .
........ ...........
......... ............
......
. ........... .
.... ......
....................
............................................... .................. .... ...
... ... ... ....... .
...
. ........................................... ... ... .. ....
... ... ... .. ........ ....... ... ....
..... ... ... ... ... ... ... . .. ... ... ... .....
......
........ ... ... ............... ... ... .. ........... ... ... ... .... ... ... ... ... ... ... ... ...
• . ...
..
.......... .
.... ... ... .........
..............
................................... ...... ............ ...........
....... ..............................................................
...........................
.............
.
If S is orientable, then it is possible to choose a unit normal vector n at every point S so that
n varies continuously over S. In that case, S is called an oriented surface and the choice of n is
called an orientation of S. There are only two orientations of an orientable surface S, namely
one for each side of the surface S which corresponds to the choice where all the unit normal
vectors point away from that side of the surface.
n.
..
.........
...
.... n
n .. .
........
......
..........
...................... • .
....................................................... ...........
............ ...........
•
.............................................
.................. ... ...........
........
.....
..... .......... .
. ..
. .. .
...
.. ......... ........ ... .......
..... ...... ....... ... ......
..... ......
..... ...... ....... . ......
....... .....
.
. .....
.
..... .........
• ..
. ... ... ... ... ...
•
... ... ... ... ... ... ... ... ...
.....
.....
... .........
.....
.
• ..... ... . .. ... ... ... ...
.......
.
.
.
. • ..
...
... ... ... ... .............. ... ...
.....
.....
...
... .. . ... ... ... ... ... . ... .......... ...
n ...
. ... ... ... ..
.
.. .... ... ...
.
.. ... ... ... .... .
.. .... ... ...
... ......
..... .
.
....
. ... ... ... ....
. ...... ... . .
...
. ..... ...
......... ....... ......... .....
...... ..
....
......
.
......
....
...... .... ...... . .
......... .
. ...
... ..
...
... ..
...
...
. n ..
.....
...... ......
. .....
......
n .
..
....
........ . .
.. .....
. ........ .... .
.....
.......... ..... .......... .....
............... .......... ............... ..........
.......................................................................................... ..........................................................................................
S −S
Figure 128 There are two possible orientations of S.
z.. z..
.. ..
......... .........
.. ..
... ...
... ...
... ...
... n ...
... .. ...
...
... ........... ...
...
. .. . ..
.
. ..
.. ......... ... .........................
.... .
.
.
. ..
.. ......... ... ......................
.... .
.
..
.....
. .
.
. .
... ....... ..
.....
. ....
......
.....
.... . ..... ...
..... .... ... .... ... ... .....
...
... • ... ....
.... ...
...
...
... ..
.
...
...
•... ....
.. ..
.
...
...
...
... . . . ... .. . .
.....
. ... ... ... ... .. ... . . . ... .. . . .. ...
.
.. ..
. ... ... ... ..
. ... ... .... .... ... ... .... . ... ... ......... ... ... ....
.......... ... .. .......... ... .. −n
...... ....... ...... ...... .............. ......
....
......... ...
. . ..
...
...
...
...
.. ..
....... . . ..
.... .. ... .. .. .. .......................................................
......
y ....
......... ...
.. ..
. . . ..
...
...
....
...
. ..... ....... .. ... .. .. .. .......................................................
......
y
... ............. ... ...... .. ... ............. ... ...... ..
... ................................................................... .... ... ................................................................... ....
... ..... .. ... ..... ..
... ..... ... ... ..... ...
... ......... ... ... ......... ...
........ ... ........ ...
......... ........ .
. ....... ........
.
. .. .. .
.. ......... .. .
..... .............. ...... ..... .......
..... ................................... .....
..........
..............................
..... ....
.......... ..........
...... .....
x The positive orientation x The negative orientation
rφ ×rθ
Therefore, n = |rφ ×rθ | = hsin φ cos θ, sin φ sin θ, cos φi = r(φ, θ), which is the outward pointing
normal.
For a closed surface (i.e. a surface which is the boundary of a solid region E), the convention
is that the positive orientation is the one for which the normal vectors point outward from E,
and the inward-pointing normals give the negative orientation. For an orientable surface which
is the boundary of a solid region in R3 , it is usually understood that it is given the positive
orientation.
z..
..
.........
...
..
...
...
....
. ...
.
.........
....... .... ............... 1 S : z = 1 − x − y2
......
.....
2
.
...... .... .....
...... .
. ..
..
... ....
...
... ...
... ... ...
... ... ...
...
... .
. ...
... ..
. ...
... ..
. ...
... ...
. ...
. ... .. . ...
.. . .
. .. . ... ... . ..
. .
.. .
.. .
. . ... ... ... ....
. ... ... ... ... ... ..
. ..........
........ .... ......
....
....
...... . ...
.
....
. ..
....................................................................................................................... y
......... .. .... ...
............... ...... ..
.....
. .
.
.........................................................
... 2 2
..... S : x +y ≤1
..... 2
.
......
.
..
....
.
....
.....
...........
......
x
Figure 130
First let’s compute the surface integral over S1 . The surface S1 is the graph of the function
g(x, y) = 1−x2 −y 2 over the disk D = {(x, y) | x2 +y 2 ≤ 1}. We have gx = −2x and gy = −2y.
ZZ ZZ
Therefore, F · dS = (−P gx − Qgy + R) dA
S1 ZZD
= [(−y)(−2x) − x(−2y) + (1 − x2 − y 2 )] dA
ZZD
= (1 + 4xy − x2 − y 2 ) dA
D
Z 2π Z 1 π
= (1 + 4r2 cos θ sin θ − r2 ) rdrdθ = .
0 0 2
ZZ
The disk S2 is oriented downward, so its unit normal vector is −k. Then, F · dS =
ZZ ZZ S2
F · (−k) dS = −z dS = 0, since z = 0 on S2 .
S2 S2
ru × rv
ZZ ZZ
F · dS = F· dS
S ZZS |ru × rv |
ru × rv
= F(r(u, v)) · |ru × rv | dA
ZZD |ru × rv |
= F(r(u, v)) · (ru × rv ) dA.
D
Therefore, ZZ ZZ
F · dS = F(r(u, v)) · (ru × rv ) dA.
S D
ZZ
Example 24.2. Let F(x, y, z) = hz, y, xi. Evaluate F · dS, where S is the unit sphere
S
x2 + y2 + z2 = 1.
Solution. A parametric representation of the unit sphere is given by
r(φ, θ) = hsin φ cos θ, sin φ sin θ, cos φi,
c Department of Mathematics 98
Semester II(2005/06)
rφ × rθ = hsin2 φ cos θ, sin2 φ sin θ, sin φ cos φi, and F(r(φ, θ)) = hcos φ, sin φ sin θ, sin φ cos θi.
n.
...
..........
...
... n
n .... .
...........
......
..........
..... ......
...
...
.. .
.......................... .. . • ............................
........... ..........
...........
..... ........ ... .
.. . . .......
..... ........... ....
..... ......
..
......... • ..... .....
.....
..... .......
...... •
... ... ... ... ...
... ... .
.. ... ... ... ... ... ... ... ... ... ... ... ... ... ... ... .
.. ... ...
.....
... ... ....
...
.
.. ...... ..... ...
.
........
..
S ......
....
...
..... .
..
... ...
.....
...... ...
.....
...
......
• ...........
.
. .... .
....... .........................................
........... .........................................................................................
...........
.........
z.. z = g(x, y)
.. ........
......... n ....... ...
....... ......
.... ... ...
......... ...
... ... ....
........ ...................
. ....
.. .
... .................... .. ..
... ... ....
. .
.
.. ..
... ..... ...
... S .. ..
... ... ... ..
...
...
.. .
.. ..
.. ..
• .
..
.
... ..
.....
. ...
..... ..... ..
... ... .. ... .............. ..
... ... ...................... .
.. ..
... ........................ C .
. ..
... .. .
. .
. ..
... .. .
. .
. ..
... .. .
. .
. ..
... .. .
. .
. ..
... .. .
. .
. ..
. . . ......................................................................
.
.
.
...
.
.
.. .
... .
...
...
...
...
...
...
...
...
...
...
.
.
.
.
..
.. .....
.... .. .
.
.
.
.
...... . .. ..
.. ....
...
...
... . . ...
.... ..
.......... .
.
.
.
.
.
..... .... ... ..
..
..
y
.
...
. .. .
.
. . . . . . . . .
..
.. .
... .
.
. ...
........................ ..............................................
. .
.. .
..... . ........ . .....
..
.
.
...
.
. .
.
.
...
... D .
. .
.
................................................................
.
.
..
..
.....
.
....
.....
..... C
.......... 1
....
x
Figure 132
c Department of Mathematics 99
Semester II(2005/06)
Z Z
dxb dy dz
F · dr = P +Q +R dt
C Zab dt dt dt
dx dy dx dy
= P + Q + R(gx + gy ) , dt
Zab dt dt dt dt
dx dy
= (P + Rgx ) + (Q + Rgy ) dt
Za dt dt
= (P + Rgx )dx + (Q + Rgy )dy
C1
∂ ∂
ZZ
= (Q + Rgy ) − (P + Rgx ) dA by Green’s Theorem on D
D ∂x ∂y
∂z ∂z
ZZ Qx + Qz ∂x + (Rx + Rz ∂x )gy + Rgyx
= dA
D ∂z ∂z
ZZ
−P y − P z ∂y − (R y + R z ∂y )g x − Rg xy
Note that Q(x, y, z) = Q(x, y, g(x, y)) is a function of x and y in D. By chain rule, we have
∂Q ∂x ∂y ∂z ∂Q ∂z
= Qx + Qy + Qz . That is = Qx + Qz = Qx + Qz gx . Similarly, we have
∂x ∂x ∂x ∂x ∂x ∂x
∂R
= R x + R z gx .
∂x
Z
Example 25.2. Evaluate F · dr, where F = −y 2 i+xj+z 2 k and C is the curve of intersection
C
of the plane y + z = 2 and the cylinder x2 + y 2 = 1. (C is oriented in the counterclockwise
sense when viewed from above.)
Solution. Let S be the surface enclosed by C on the plane y + z = 2. S is the graph of
z = g(x, y) = 2 − y over the disk D = {(x, y) | x2 + y 2 ≤ 1}.
.....................
... ...
z..
....... ......
..... ...... ....
... ........
. ......
......... ...........
... .... ... .......
... .... ... .....
... ...... C ... .....
.....
... ......... ... ....
...
... ..... .
.
..... . ...
... ..... . ...
... ..... ... ...
...
... ..... .
..... .. ...
... .....
... ......... y+z =2 ..
..
...
... .. .........
... S .. .......
.......... .........
... .. ........ ..
... .. .
.
... .
... ... ... ......... ... ... ... ... . ...
... ... ... ... ... .. ... ....
........ ... .....
...
......
...
...
. .
...
...
...
...
. .. .
. .... ......... ...... ..
......
................................................
... ... y
......... D
..........................
. ................................
.........
........ 2 2
.....
.
. x +y =1
.....
.....
...........
......
x
Figure 133
c Department of Mathematics 100
Semester II(2005/06)
i j k
∂ ∂ ∂
Also curl F = ∂x = (1 + 2y)k. By Stokes’ Theorem,
∂y ∂z
−y 2 z2
x
Z ZZ ZZ Z 2π Z 1
F · dr = curl F · dS = (1 + 2y) dA = (1 + 2r sin θ)rdrdθ
C S D 0 0
Z 2π 1 2
= ( + sin θ) dθ = π.
0 2 3
ZZ
Example 25.3. Use Stokes’ Theorem to compute curl F · dS, where F(x, y, z) = yzi +
S
xzj + xyk and S is the part of the sphere x2 + y 2 + z 2 = 4 that lies above the cylinder
x2 + y 2 = 1 and above the xy-plane.
p
Solution. The cylinder
√ x2 +y 2 = 1 intersects the upper hemisphere z = 4 − x2 + y 2 in√a curve
C at height z = 3. The curve C has a√vector equation √ given by r(t) = hcos t, sin t, 3i and
0
r (t) = h− sin t, cos t, 0i. Also F(r(t)) = h 3 sin t, 3 cos t, cos t sin ti. By Stokes’ Theorem,
z..
..
..........
...
...
....
..
...
.. ... . S ......
................. ........ .............................
...
.....
..
..
............... .......
· ..
.... ... .............................................................. ..........
.... C ..
.
.
..........
....
.....
...
... .
. . ... ...
... .... .. .. ...
. .. . ... ... ... ........ ... ... ... ... ...
.
...
...
.... .. ... ... ... ...... .. ...........
. .
........ ... .. . .. ...... ....
.
.
. ..
...... ........... ... ... ..................................................................................................................
...
.....
........... ...
............................................ . .
....
.
.. ... y
....... .
... .......
...........
............................ ...... ..
. . .
...
............
.
..
.... .......................
.....
.....
.....
..
......
....
.....
.....
..........
.....
x
Figure 134
ZZ Z Z
0
Z 2π √ √
curl F · dS = F · dr = F(r(t)) · r (t) dt = (− 3 sin2 t + 3 cos2 t) dt = 0.
S C C 0
Exercise 25.5. Let S be the capped cylindrical surface shown in figure 135. S is the union of
the surfaces S1 and S2 where S1 = {(x, y, z) | x2 + y 2 = 1, 0 ≤ z ≤ 1} and S2 = {(x, y, z) | x2 +
2 + (z − 1)2 = 1, z ≥ 1}. Let F(x, y, z) = (zx + z 2 y + x)i + (z 2 yx + y)j + z 4 x2 k. Evaluate
yZZ
curl F · dS.
S
[Answer. 0]
c Department of Mathematics 101
.........
............. ....................
Semester II(2005/06)
....... ......
...... .....
..
...... .....
. ...
...
S 2 ... ...
.
... ... ... ... ... ... ... ... ...
...
...
..
... ... ... ... ....
.... ........ ...
...... ......
.... ...
............ ....
......
... ................... ........ ....
... ................................................... ...
... ...
S 1 ..... ... ... ... ... ... ... ... ... ...
...
..
... ...
.... ......... ... ... ... ....
...... .....
... ...
..... ..
........ .....
.............. ........
...................................................
Figure 135
26. The Divergence Theorem
Theorem 26.1. (The Divergence Theorem or Gauss’ Theorem) Let E be a solid region which
is both of type I, II and III, and let S be the boundary of E, given with the positive (out-
ward) orientation. Let F be a vector field whose component functions have continuous partial
derivatives on an open region containing E. Then
ZZ ZZZ
F · dS = div F dV.
S E
Proof. Let F = P i + Qj + Rk. Then div F = Px + Qy + Rz . Thus,
ZZZ ZZZ ZZZ ZZZ
div F dV = Px dV + Qx dV + Rx dV.
E E E E
Let n be the unit outward normal of S. Then
ZZ ZZ ZZ ZZ ZZ
F · dS = (P i + Qj + Rk) · n dS = P i · n dS + Qj · n dS + Rk · n dS.
S S S S S
We shall show
ZZ ZZZ ZZ ZZZ ZZ ZZZ
P i · n dS = Px dV , Qj · n dS = Qy dV , and Rk · n dS = Rz dV .
S E S E S E
z..
...
........ S : z = u2 (x, y)
................................................. 2
........ ....
.... ....... .
... ................. .......
... .. ......................................... .....
... ... ...
... E
... ... ... ... ... ... ... ... ... . .... 3 S
... .......... .. ... ...
...
... .... .
... .................. .......
... .. ...................................... ...
... .. S : z = u1 (x, y) .. 1
... .. ..
... .. ..
... .. ..
... .. ..
... .. .
................................................................................................................................
.
.
...
.
.
....
.
.. .
.
. ... y
.. . ..
..... .. ..
......
. .
.
...
.
. .
. ............ ....................................... ..
..
. . ... ........
....
. ....
.
.....
.
. D ........
...................... ...........................
..
..... ...
.....
..........
......
x
To prove the third equation, we use the fact that E is a type I solid region:
E = {(x, y, z) | u1 (x, y) ≤ z ≤ u2 (x, y), (x, y) ∈ D},
where D is the projection of E onto the xy-plane. The boundary of E consists of S1 , S2 , and
S3 .
ZZ
On S3 , n is perpendicular to k. Thus Rk · n dS = 0.
S
c Department of Mathematics 102
Semester II(2005/06)
The surface S2 is given as the graph of z = u2 (x, y) with (x, y) ∈ D. The outward normal n
points upward. Thus, ZZ ZZ
Rk · n dS = R(x, y, u2 (x, y)) dA.
S D
The surface S1 is given as the graph of z = u1 (x, y) with (x, y) ∈ D. The outward normal n
points downward. Thus,
ZZ ZZ
Rk · n dS = − R(x, y, u1 (x, y)) dA.
S D
ZZ ZZ
Therefore, Rk · n dS = R(x, y, u2 (x, y)) − R(x, y, u1 (x, y)) dA
S D
ZZ ZZ Z u2 (x,y) ∂R
z=u (x,y)
= [R(x, y, z)]z=u21 (x,y) dA = dz dA
D D u1 (x,y) ∂z
ZZZ
= Rz dV .
E
Similarly, using the fact that E is a type II and III solid region, one can prove the second and
first equation. Combining the three equations, we get the Divergence Theorem.
ZZ
Example 26.2. Let F(x, y, z) = zi + yj + xk. Evaluate F · dS, where S is the unit sphere
S
x2 + y 2 + z 2 = 1.
ZZ ZZZ ZZZ
Solution. By the Divergence Theorem, F · dS = div F dV = 1 dV = volume of
S E E
the unit ball = 4π/3.
ZZ
2
Example 26.3. Evaluate F · dS, where F = xyi + (y 2 + exz )j + sin(xy)k and S is the
S
surface of the region E bounded by the parabolic cylinder z = 1−x2 and the planes z = 0, y = 0
and y + z = 2.
Solution. The solid region E can be described as
E = {(x, y, z) | − 1 ≤ x ≤ 1, 0 ≤ z ≤ 1 − x2 , 0 ≤ y ≤ 2 − z}.
z..
..
.........
....
...
...
...
(0, 0, 1) ...
.
y =2−z
................................................................................... ...
.
.... ... ....... ...... ........
... ... ....... ....
..
. ....
.......
.
.
... .
. .
.
... .... .
. .
...... ............
... .... . ......
... .... ... .. ..... ......
... ... ... ... .....
.....
... ... ... ... ...
.... ... ... .... ... ... ... ... ... ... ... ... ... ... ... ... ... ... ............
... ... ... ... ..
... ... ....
. ........
.
.. ... ..... ....
.. . .....
.... ... ..... ... .........
.. ...... ... ....
. ..... ...
....
...
.....
........................................................................................................................................................
..... ... .
.....
y
...
... .
.
.
......
...
..
...... (0, 2, 0) .. .
... . .
... ...... .
. . ...... ...
........ ......... .
..... ........
.
(1, 0, 0) ............................................................................................
..... .
..
......
.
..
.......... z =1−x 2
.......
x
Figure 137
c Department of Mathematics 103
Semester II(2005/06)
Exercise 26.6. Let S be a closed surface and let F a vector field onZZR3 whose component
functions have continuous second order partial derivatives. Show that curl F · dS = 0.
S
ZZ ZZZ
Exercise 26.7. Prove that Dn (f ) dS = ∇2 f dV , where Dn (f ) is the directional
∂E E
derivative of f along the direction of the unit normal n to the surface ∂E.
(Recall that Dn (f ) = ∇f · n.)
c Department of Mathematics 104
Bibliography
105