Vector Calculus Book
Vector Calculus Book
Vector Calculus Book
Mathesis iuvenes tentare rerum quaelibet ardua semitasque non usitatas pandere docet.
ii
Preface v
2 Differentiation 97
2.1 Some Topology . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
2.2 Multivariable Functions . . . . . . . . . . . . . . . . . . . . . . . . 99
2.3 Limits . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
2.4 Definition of the Derivative . . . . . . . . . . . . . . . . . . . . . . . 104
2.5 The Jacobi Matrix . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
2.6 Gradients and Directional Derivatives . . . . . . . . . . . . . . . . . 116
2.7 Levi-Civitta and Einstein . . . . . . . . . . . . . . . . . . . . . . . . 122
2.8 Extrema . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 124
2.9 Lagrange Multipliers . . . . . . . . . . . . . . . . . . . . . . . . . . 129
3 Integration 133
3.1 Differential Forms . . . . . . . . . . . . . . . . . . . . . . . . . . . 133
3.2 Zero-Manifolds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 136
3.3 One-Manifolds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 137
3.4 Closed and Exact Forms . . . . . . . . . . . . . . . . . . . . . . . . 141
3.5 Two-Manifolds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 146
3.6 Change of Variables . . . . . . . . . . . . . . . . . . . . . . . . . . 155
3.7 Change to Polar Coordinates . . . . . . . . . . . . . . . . . . . . . . 162
3.8 Three-Manifolds . . . . . . . . . . . . . . . . . . . . . . . . . . . . 168
iii
iv
These notes started during the Spring of 2003. They are meant to be a gentle
introduction to multivariable and vector calculus.
David A. SANTOS
[email protected]
v
1 Vectors and Parametric Curves
W e start with a naı̈ve introduction to some linear algebra necessary for the
course. Those interested in more formal treatments can profit by reading
[BlRo] or [Lan].
Notice that infinitely many different choices of departure and arrival points
may give the same vector.
a1 = (1, 2), b1 = (3, −4), a2 = (3, 5), b2 = (5, −1), O = (0, 0) b = (2, −6).
1
Some authors use the terminology “fixed vector” instead of “bi-point.”
1
y axis
b
r′
b r = (x, y)
Points and Vectors on the Plane b
r
b
x axis
b
a2
2
Figure 1.1: A point
a in R .
b
Figure 1.2: A bi-point in R2 .
1
O b
b
b2
b
b1
b
b
Though the bi-points [a1 , b1 ], [a2 , b2 ] and [O, b] are in different locations on the
plane, they represent the same vector, as
" # " # " # " #
3−1 5−3 2−0 2
= = = .
−4 − 2 −1 − 5 −6 − 0 −6
The instructions given by the vector are all the same: start at the point, go two
units right and six units down. See figure 1.3.
3Z = {. . . , −6, −3, 0, 3, 6, . . .}, 3Z+1 = {. . . , −5, −2, 1, 4, 7, . . .}, 3Z+2 = {. . . , −4, −1, 2, 5, 8, . . .}.
The equivalence class 3Z comprises the integers divisible by 3, and for exam-
ple, −18 ∈ 3Z. Analogously,
" #
in example "2, the # bi-point [a1 , b1 ] belongs to the
2 2
equivalence class , that is, [a1 , b1 ] ∈ .
−6 −6
−→
3 Definition The vector Oa that corresponds to the point a ∈ R2 is called the
position vector of the point a.
←→
4 Definition Let a 6= b be points on the plane and let ab be the line passing
through a and b. The direction of the bi-point [a, b] is the direction of the line
←→
L, that is, the angle θ ∈ [0; π[ that the line ab makes with the positive x-axis
(horizontal axis), when measured counterclockwise. The direction of a vector
→
− →
−
v 6= 0 is the direction of any of its bi-point representatives. See figure 1.4.
←→
5 Definition We say that [a, b] has the same direction as [z, w] if ab = ← → We
zw.
say that the bi-points [a, b] and [z, w] have the same sense if they have the same
direction and if when translating one so as to its tail is over the other’s tail, both
their heads lie on the same half-plane made by the line perpendicular to then at
their tails. They have opposite sense if they have the same direction and if when
translating one so as to its tail is over the other’s tail, their heads lie on different
half-planes made by the line perpendicular to them at their tails. . See figures 1.5
and 1.6 . The sense of a vector is the sense of any of its bi-point representatives.
Two bi-points are parallel if the lines containing them are parallel. Two vectors
are parallel, if bi-point representatives of them are parallel.
B B B
D
A A A
b θ C C
6 Definition The Euclidean length or norm of bi-point [a, b] is simply the distance
between a and b and it is denoted by
È
||[a, b]|| = (a1 − b1 )2 + (a2 − b2 )2 .
A bi-point is said to have unit length if it has norm 1. The norm of a vector is the
norm of any of its bi-point representatives.
☞ A vector is completely determined by three things: (i) its norm, (ii) its direc-
tion, and (iii) its sense. It is clear that the norm of a vector satisfies the following
properties:
−
1. →
a ≥ 0.
− →
−
2. →
a = 0 ⇐⇒ → −a = 0.
" #
− È √ √
1
7 Example The vector →
−
v = √ has norm →
v = 12 + ( 2)2 = 3.
2
8 Definition If →
−
u and →
−
v are two vectors in R2 their vector sum →
−
u +→
−
v is defined
by the coordinatewise addition
" # " # " #
→
− u1 v1 u1 + v1
u +→
−
v = + = . (1.1)
u2 v2 u2 + v2
−2→
−
u
B →
− →
−
→
− v u
u
A
→
− C
u +→
−
v
1→−
u
2
9 Definition If α ∈ R and →
−
a ∈ R2 we define scalar multiplication of a vector and a
scalar by the coordinatewise multiplication
" # " #
a1 αa1
α→
−
a =α = . (1.3)
a2 αa2
It is easy to see that vector addition and scalar multiplication satisfies the
following properties.
→
− →
−
➊ α(→ −a + b ) = α→−a +αb
➋ (α + β)→ −a = α→−a + β→−a
➌ 1→
−a =→−
a
➍ (αβ)→
−
a = α(β →
−
a)
−v
→
−v
→
−w
+ →
→
−w
→
−
→
−
v
+ →
→
v +→
−
−
−v −
w
−
v
+ →
w
w−
w→
−
→
u→
→
→
−
−
−
w
v
→
→
−
−
→
−
v
−
→
−
→
v )+ w
(u +
−)
→
(
−
→v +w
−v
→ −
→
u + −
→
w
→
−
10 Definition Let →−
u =6 0 . Put"R→
−
u# = {λ→−
u : λ ∈ R} and let a ∈ R2 , The affine
u1
line with direction vector →
−
u = and passing through a is the set of points on
u2
the plane
( ! )
x
a + R→
−
u = 2
∈ R : x = a1 + tu1 , y = a2 + tu2 , t∈R .
y
u2
that is, the affine line is the Cartesian line with slope . Conversely, if y =
u1
mx + k is the equation of a Cartesian line, then
! " # !
x 1 0
= t+ ,
y m k
that# is, every Cartesian line is also an affine line and one may take the vector
"
1
as its direction vector. It also follows that two vectors →
−
u and →
−v are parallel
m
if and only if the affine lines R→−u and R→ −
v are parallel. Hence, →
−
u k →−
v if there
→
−
exists a scalar λ ∈ R such that u = λ v . →
−
☞ →
−
Because 0 = 0→
−
v for any vector →
− →
−
v , the 0 is parallel to every vector.
u→
−
b
r = a + t→
−
u
b
a
" #
1
11 Example Find a vector of length 3, parallel to →
−
v = √ but in the opposite
2
sense.
− √ →
−
v
Solution: ◮ Since → v = 3, the vector → has unit norm and
−
v
has the same direction and sense as →
−
v , and so the vector sought is
" # " √ #
→
−
v 3 1 − 3
−3 → = − √ √ = √ .
−
v 3 2 − 6
!
1
12 Example Find the parametric equation of the line passing through and
−1
" #
2
in the direction of the vector .
−3
Some plane geometry results can be easily proved by means of vectors. Here
are some examples.
−→ −→
13 Example Given a pentagon ABCDE, determine the vector sum AB + BC +
−→ −→ −→
CD + DE + EA.
.◭
14 Example Consider a △ABC. Demonstrate that the line segment joining the
midpoints of two sides is parallel to the third side and it is in fact, half its length.
as we were to shew. ◭
−−−→ −−−→
Solution: ◮ As AMC = MC B, we have,
−→ −→ −−−→ −−−→ −−−→ −−−→
CA + CB = CMC + MC A + CMC + MC B
−−−→ −−−→ −−−→
= 2CMC − AMC + MC B
−−−→
= 2CMC ,
←−−→
Solution: ◮ Since the triangle is non-degenerate, the lines AMA and
←−→ −→
BMB are not parallel, and hence meet at a point G. Therefore, AG and
−−−→ −→ −−−→
GMA are parallel and hence there is a scalar a such that AG = aGMA .
−→ −−−→
In the same fashion, there is a scalar b such that BG = bGMB . From
example 14,
−−−−→ −→
2MA MB = BA
−→ −→
= BG + GA
−−−→ −−−→
= bGMB − aGMA
−−−→ −−−−→ −−−→
= bGMA + bMA MB − aGMA ,
and thus
−−−−→ −−−→
(2 − b)MA MB = (b − a)GMA .
−−−−→ −−−→
Since △ABC is non-degenerate, MA MB and GMA are not parallel,
whence
2 − b = 0, b − a = 0, =⇒ a = b = 2.
◭
MB b b
MA Chapter 1
b
G
b b b
A B
MC
Figure 1.13: Example 17.
It follows that
−−→′ −→ −−→
GG = GB + BG′
−−−→ −−−→
= −2GMB + 2G′ MB
−−−→ −−−→
= 2(MB G + G′ MB )
−−→
= 2G′ G.
Therefore
−−→′ −−→ −−→ →
− −−→ →
−
GG = −2GG′ =⇒ 3GG′ = 0 =⇒ GG′ = 0 =⇒ G = G′ ,
Homework
Problem 1.1.1 Is there is any truth to the such that
statement “a vector is that which has mag- −
→ −
→ −→ 1−→
nitude and direction”? IA = −3IB, JA = − JB,
3
and then demonstrate that for any arbitrary
Problem 1.1.2 ABCD is a parallelogram. point M on the plane
E is the midpoint of [B, C] and F is the mid- −−→ −−
→ −
→
point of [D, C]. Prove that MA + 3MB = 4MI
−→ −→ −→ and
−−→ −−→ −→
AC + BD = 2BC. 3MA + MB = 4MJ.
Problem 1.1.3 (Varignon’s Theorem) Use Problem 1.1.5 Find the Cartesian equa-
vector algebra in order to prove that in any tion corresponding to the line with paramet-
quadrilateral ABCD, whose sides do not ric equation
intersect, the quadrilateral formed by the
x = −1 + t, y = 2 − t.
midpoints of the sides is a parallelogram.
→
−
a →
−
a →
−
a
→
− →
−
c d
→
− →
−
c →
−
c →
−
b b
→
−
b
Figure 1.14: [A]. Problem Figure 1.15: [B]. Problem Figure 1.16: [C]. Problem
1.1.8. 1.1.8. 1.1.8.
→
−
a →
−
a →
− →
−
b c
→
− →
− →
− →
− →
− →
−
c d c d a d
→
− →
− →
− →
−
e
b b f
Figure 1.17: [D]. Problem Figure 1.18: [E]. Problem Figure 1.19: [F]. Problem
1.1.8. 1.1.8. 1.1.8.
18 Definition Let →
−
x ∈ R2 and →−y ∈ R2 . Their scalar product (dot product, inner
product) is defined and denoted by
→
−
x •→
−
y = x 1 y1 + x 2 y2 .
" # " #
1 →
− 3 →
−
19 Example If →
−
a = and b = , then →
−
a • b = 1 · 3 + 2 · 4 = 11.
2 4
The following properties of the scalar product are easy to deduce from the defini-
tion.
SP1 Bilinearity
(→
− y )•→
x +→
− −
z =→
−
x •→
−
z +→
−
y •→
−
z, →
−
x •(→
−
y +→
−
z)=→
−
x •→
−
y +→
−
x •→
−
z (1.4)
(α→
−
x )•→
−
y =→
−
x •(α→
−
y ) = α(→
−
x •→
−
y ), α ∈ R. (1.5)
SP3 Commutativity
→
−
x •→
−
y =→
−
y •→
−
x (1.6)
SP4
→
−
x •→
−
x ≥ 0 (1.7)
SP5
→
− →
−
x •→
−
x = 0 ⇐⇒ →
−
x = 0 (1.8)
SP6 → p
−
x = →−
x •→
−
x (1.9)
→
−
b −→
−
a
→
−
a
→
−
b
→
−
20 Definition Given vectors → −a and b , we define the (convex) angle between them,
\ →
− →
−
denoted by (→ −
a , b ) ∈ [0; π], as the angle between the affine lines R→
−
a and R b .
→
−
21 Theorem Let →
−
a and b be vectors in R2 . Then
→
− →
− →
− \ →
−
a • b = ||→
−
a |||| b || cos (→
−
a , b ).
Proof: From figure 1.20, using Al-Kashi’s Law of Cosines on the length
of the vectors, and (1.4) through (1.9) we have
→
− →
− →
− \ →
−
|| b − →
−
a ||2 = ||→
−a ||2 + || b ||2 − 2||→−
a |||| b || cos (→
−
a, b)
→
− →
− →
− →
− \ →
−
⇐⇒ ( b − → −a )•( b − →−a ) = ||→ −
a ||2 + || b ||2 − 2||→−
a |||| b || cos (→−a, b)
→
− → − →
− →
− →
− \ →
−
⇐⇒ b • b − 2→ −a•b +→ −a •→
−a = ||→−a ||2 + || b ||2 − 2||→
−a |||| b || cos (→−a, b)
→
− →
− →
− →
− \ →
−
⇐⇒ || b ||2 − 2→ −a • b + ||→−a ||2 = ||→
−a ||2 + || b ||2 − 2||→−a |||| b || cos (→
−a, b)
→
− →
− \ →
−
⇐⇒ →
−
a • b = ||→
−
a |||| b || cos (→
−
a , b ),
as we wanted to shew. ❑
\ →
− π
Putting (→
−
a, b)= in Theorem 21 we obtain the following corollary.
2
22 Corollary Two vectors in R2 are perpendicular if and only if their dot product
is 0.
☞ →
−
It follows that the vector 0 is simultaneously parallel and perpendicular
to any vector!
− →
→
− − →
−
24 Definition If →
−
a ⊥ b and →
a = b = 1 we say that →
−
a and b are or-
thonormal.
→
−
Equality occurs if and only if →
−
a k b.
" # " #
a1 →
− b1
If →
−
a = and b = , the CBS Inequality takes the form
a2 b2
Proof:
→
− →
− →
−
||→
−
a + b ||2 = (→
−a + b )•(→ −a + b)
→
− →
− →−
= →
−a •→
−a + 2→ −a•b + b•b
→
− →
−
≤ ||→
−a ||2 + 2||→−
a |||| b || + || b ||2
→
−
= (||→
−a || + || b ||)2 ,
from where the desired result follows. ❑
Also,
→ → − È È p
− −
a + b + →
c = x2 + y 2 + y 2 + z 2 + z 2 + x2 ,
29 Definition Let A and B be points on the plane and let → −u be a unit vector.
−→ →
−
If AB = λ u , then λ is the directed distance or algebraic measure of the line
segment [AB] with respect to the vector →
−u . We will denote this distance by AB −
u,
→
→
−
or more routinely, if the vector u is patent, by AB. Observe that AB = −BA.
→ ←
← →
30 Theorem (Thales’ Theorem) Let D y D′ be two distinct lines on the plane. Let
←→ ←
→
A, B, C be distinct points of D , and A′ , B′ , C′ be distinct points of D′ , A 6= A′ ,
←→ ←→
B 6= B′ , C 6= C′ , A 6= B, A′ 6= B′ . Let AA′ k BB′ . Then
←−
→ ←→ AC A′ C ′
AA′ k CC′ ⇐⇒ = ′ ′.
AB AB
D D′ D′ D
b
C
A b b
A′
B b b
B′ A b b
A′
C b b
C′ B b b
B′
Proof: Refer to figure 1.2. On the one hand, because they are unit
vectors in the same direction,
−→ −→ −− → −− →
AB AC A′ B′ A′ C′
= ; = ′ ′.
AB AC A′ B ′ AC
On the other hand, by Chasles’ Rule,
−−→′ −→ −−→ −−→ −−→ −→ −−→
BB = BA + AA′ + A′ B′ = (A′ B′ − AB) + AA′ .
←−
→ ←→
Since AA′ k BB′ , there is a scalar λ ∈ R such that
−− → −→ −−→
A′ B′ = AB + λAA′ .
←
→ ←
→
31 Corollary Let D and D′ are distinct lines, intersecting in the unique point C.
←
→ ←
→
Let A, B, be points on line D , and A′ , B ′ , points on line D′ . Then
←−
→ ←→ CB CB ′
AA′ k BB′ ⇐⇒ = .
CA CA′
→
−
and so both vectors are parallel to j and hence →
−
x k→
−
y . We have demonstrated
the following theorem.
→
−
v = a→
−
x + b→
−
y, a ∈ R, b∈R
if and only if →
−
x is not parallel to →
−
y . In this last case we say that →
−
x is linearly
→
−
independent from vector y . If two vectors are not linearly independent, then we
say that they are linearly dependent.
" # " #
1 1
33 Example The vectors and are clearly linearly independent, since one
0 1
" #
a
is not a scalar multiple of the other. Given an arbitrary vector we can express
b
it as a linear combination of these vectors as follows:
" # " # " #
a 1 1
= (a − b) +b .
b 0 1
{a→
−
x + b→
−
y : a ∈ [0; 1], b ∈ [0; 1]}.
(" # " #)
1 1
Figure 1.23 shews the fundamental parallelogram of ,
, coloured in
0 1
brown, and the respective tiling of the plane(by various
! !translations
! !)of it. Observe
0 1 1 2
that the vertices of this parallelogram are , , , . In essence
0 0 1 1
then, linear independence of two vectors on the plane means that we may obtain
every vector on the plane as a linear combination of these two vectors and hence
Homework
1 −1 Problem 1.3.3 Consider the line with
Problem 1.3.1 Prove that and Cartesian equation L : ax + by = c, where
1 1
are linearly independent, and draw their not both of a, b are zero. Let t be a point
fundamental parallelogram. not on L. Find a formula for the distance
from t to L.
Problem
1.3.2 Write an arbitrary vector
Problem 1.3.4 Prove that two non-zero
a
on the plane, as a linear combination perpendicular vectors in R2 must be lin-
b early independent.
1 −1
of the vectors and .
1 1
2 2
36 Definition A function T→
v : R → R is said to be a translation if it is of the
−
form T→ →
− →
−
v (x) = x + v , where v is a fixed vector on the plane.
−
A translation simply shifts an object on the plane rigidly (that is, it does not
distort it shape or re-orient it), to a copy of itself a given amount of units from
where it was. See figure 1.24 for an example.
→
−
v
It is clear that the composition of any two translations commutes, that is, if
2 2
T→
v 1 , T→
− v 2 : R → R are translations, then T→
− v 1 ◦ T→
− − v 2 ◦ T→
v 2 = T→
− v 1 . For let
−
→
−
T1 (a) = a + v 1 and T→ →
−
v 2 (a) = a + v 2 . Then
−
− ◦ T→
(T→ − )(a) = T→ − (T→− (a)) = T→
→
− =a+→
− (a) + v
−v +→ −v ,
v 1 v 2 v 1 v 2 v 2 1 2 1
and
→
− →
− →
−
v 2 ◦ T→
(T→
− − v 2 (T→
v 1 )(a) = T→
− v 1 (a)) = T→
− v 1 (a) + v 2 = a + v 1 + v 2 ,
−
!! !
x 2x
Figure 1.25 shews the scaling S2,0.5 = .
y 0.5y
It is clear that the composition of any two scalings commutes, that is, if
Sa,b , Sa′ ,b′ : R2 → R2 are scalings, then Sa,b ◦ Sa′ ,b′ = Sa′ ,b′ ◦ Sa,b . For
!! !
a′ x a(a′ x)
(Sa,b ◦ Sa′ ,b′ )(r) = Sa,b (Sa′ ,b′ (r)) = Sa,b = ,
b′ y b(b′ y)
and !! !
ax a′ (ax)
(Sa′ ,b′ ◦ Sa,b )(r) = Sa′ ,b′ (Sa,b (r)) = Sa′ ,b′ = ,
by b′ (by)
from where the commutativity claim is deduced.
but
!! !!! !! !
−1 −1 0 0
(S2,1 ◦ T→
−)
i
= S2,1 T→
−
i
= S2,1 = .
0 0 0 0
y axis
x axis
Figure 1.26: Reflexions. The original object (in the first quadrant) is yellow. Its reflexion
about the y-axis is magenta (on the second quadrant). Its reflexion about
the x-axis is cyan (on the fourth quadrant). Its reflexion about the origin is
blue (on the third quadrant).
and
ρ sin(α + θ) = ρ sin α cos θ + ρ sin θ cos α = y cos θ + x sin θ.
! !
x x cos θ − y sin θ
Hence the point is mapped to the point .
y x sin θ + y cos θ
y-axis
ρ cos(α + θ)
ρ sin(α + θ)
b
b
ρ cos α
ρ sin α
θ
α
b b
x-axis
Figure 1.27: Rotation by an angle θ in the levogyrate (counterclockwise) sense from the
x-axis.
It is easy to prove that scalings, reflexions, and rotations are linear transfor-
mations from R2 to R2 , but not so translations.
A(x) = L(x) + →
−
v.
It is easy to see that translations are then affine transformations, where for
the linear transformation L in the definition we may take I : R2 → R2 , the
identity transformation I(x) = x.
We have seen that scalings, reflexions and rotations are linear transforma-
tions. If L : R2 → R2 is a linear transformation, then
→
− →
− →
− →
−
L(r) = L(x i + y j ) = xL( i ) + yL( j ),
and
!! ! !
0 0·1 0
Sa,b = = .
1 b·1 b
44 Example (Reflexion
" Matrices)
# It is easy to verify that the matrix for the
" trans-
#
−1 0 1 0
formation RH is , that the matrix for the transformation RV is
0 1 0 −1
" #
−1 0
, and the matrix for the transformation RO is .
0 −1
45 Example
" (Rotating#Matrices) It is easy to verify that the matrix for a rotation
cos θ − sin θ
Rθ is .
sin θ cos θ
2
47 Example (Zero Matrix) # matrix for the null linear transformation N : R →
" The
0 0
R2 , N(x) = O is 02 = .
0 0
From problem 1.4.7 we know that the composition of two linear transfor-
mations is also linear. We are now interested in how to codify the matrix of a
composition of linear transformations L1 ◦ L1 in terms of their individual matri-
ces.
" #
2 2 a b
48 Theorem Let L : R → R have the matrix representation AL = and
c d
" #
′ 2 2 r s
let L : R → R have the matrix representation AL′ = . Then the compo-
t u
sition L ◦ L′ has matrix representation
" #
ar + bt as + bu
.
cr + dt cs + du
!! !!
′ 1 ′ 0
Proof: We need to find (L ◦ L ) and (L ◦ L ) .
0 1
We have
!! !!! !! ! ! !
′ 1 ′ 1 r →
− →
− a b ar + bt
(L◦L ) =L L =L = rL( i )+tL( j ) = r +t = ,
0 0 t c d cr + dt
and
!! !!! !! ! ! !
0 0 s →
− →
− a b as + bu
(L◦L′ ) = L L′ =L = sL( i )+uL( j ) = s +u = ,
1 1 u c d cs + du
" #
′ ar + bt as + bu
whence we conclude that the matrix of L ◦ L is , as
cr + dt cs + du
we wanted to shew.❑
51 Example Find a 2 × 2 matrix that will transform the square in figure 1.28 into
the parallelogram in figure 1.29. Assume in each case that the vertices of the
figures are lattice points, that is, coordinate points with integer coordinates.
" #
a b
Solution: ◮ Let be the desired matrix. Then since
c d
" #" # " #
a b 0 0
= ,
c d 0 0
!
0
the point is a fortiori, transformed to itself. We now assume, with-
0
out loss of generality, that each vertex of the square is transformed in
the same order, counterclockwise, to each vertex of the rectangle. Then
" #" # " # " # " #
a b 1 2 a 2
= =⇒ = =⇒ a = c = 2.
c d 0 2 c 2
3 3
b
2 2
b
1 1
b
0 0
b
-1 -1
-2 -2
-3 -3
-4 -4
-4 -3 -2 -1 0 1 2 3 4 5 -4 -3 -2 -1 0 1 2 3 4 5
Homework
1 −1 Problem 1.4.7 Let L : R2 → R2 and L′ :
Problem 1.4.1 If A = , B = R2 → R2 be linear transformations. Prove
2 3
that their composition L ◦ L′ is also a linear
a b transformation.
and
1 −2
(A + B)2 = A2 + 2AB + B 2 , 5
4
find a and b. 3
2
Problem
1.4.2 Consider
△ABC
with A = 1
−1 0 2 0
, B = , C = , as in fig-
2 −2 1 -1
Problem 1.4.3 Find the effects of the re- Figure 1.30: Problems 1.4.2, 1.4.8, and 1.4.3.
flexions R π , R π , R− π , and R− π on the
2 4 2 4
triangle in figure 1.30.
Problem 1.4.8 Find the effects of the re-
Problem 1.4.4 Prove that the composition flexions RH , RV , and RO on the triangle
of two reflexions is commutative. Prove that in figure 1.30.
the composition of a reflexion and a scaling
is commutative. Prove that the composition Problem 1.4.9 Find all matrices A ∈
of a reflexion and a translation is not nec- M2×2 (R) such that A2 = 02
essarily commutative.
Problem 1.4.5 Prove that the composition Problem 1.4.10 Find the image of the fig-
of two rotations on the plane Rθ and Rθ ′ ure below (consisting of two
circles
and a
satisfies 1 1
triangle) under the matrix .
−1 3
Rθ ◦ Rθ ′ = Rθ+θ ′ = Rθ ′ ◦ Rθ ,
2
and so the composition of two rotations on
the plane is commutative.
1
with common formulæ of areas of plane figures, in particular, the area of a par-
allelogram should be as we learn in elementary geometry and the area of a unit
square should be 1.
c
d b b b
b
b b
a
d
b b
b b b
a c
" #
a
From figures (1.32) and (1.33), the area of a parallelogram spanned by ,
b
" #
c
and is
d
" # " #!
a c
D , = ad − bc.
b d
1
D(→
−
r1 , →
−
r2 ) + D(→
−
r2 , →
−
r ) + D(→
−
r ,→
−
r3 ) + D(→
−
r3 , →
−
r1 ) .
4
This is equivalent to
1 1
D(→
−
r1 , →
−
r2 ) + D(→
−
r2 , →
−
r3 ) + D(→
−
r3 , →
−
r1 ) − D(→
−
r1 , →
−
r2 ) − D(→
−
r2 , →
−
r ) − D(→
−
r ,→
−
r3 ) + D(→
−
r3 , →
−
r1 ) + 2D(→
−
r2 , →
−
r3 )
2 4
We will prove that
D(→
−
r1 , →
−
r2 ) − D(→
−
r2 , →
−
r ) − D(→
−
r ,→
−
r3 ) + D(→
−
r3 , →
−
r1 ) + 2D(→
−
r2 , →
−
r3 ) = 0.
D(→
−
r1 , →
−
r2 ) − D(→
−
r2 , →
−
r ) − D(→
−
r ,→
−
r3 ) + D(→
−
r3 , →
−
r1 ) + 2D(→
−
r2 , →
−
r3 ) = D(→
−
r1 , →
−
r2 ) − D(→
−
r2 , →
−
r3 + →
−
r2 − →
−
r1 ) + 2
−D( r + r − r , r ) + D(→
→
− →
− →
− →
− −r ,→−
r ) 3 2 1 3 3 1
= D(→
−
r1 , →
−
r2 − →
−
r3 ) + D(→
−
r3 + →
−
r2 − →
−
r1 , →
−
r
→
− →
−
+2D( r , r ) 2 3
= D(→
−r3 + → −
r2 , →
−
r2 − →
−
r3 ) + 2D(→−r2 , →
−
r3 )
= D( r3 , r2 ) − D( r2 , r3 ) + 2D( r2 , →
→
− →
− →
− →
− →
− −
r3 )
→
− →
− →
−
= D( r , r ) − D( r , r )→
−
3 2 2 3
= 0,
as claimed. We have proved then that the area of a triangle, whose vertices
(x1 , y1 ), (x2 , y2 ), (x3 , y3 ) are listed in counterclockwise order, is
" # " # " #!
1 x1 x2 x2 x3 x3 x1
det + det + det . (1.12)
2 y1 y2 y2 y3 y3 y1
r4 (x3 , y3 )
(x4 , y4 )b r3 b
r2
Chapter 1
r
b
(x2 , y2 ) b
(x2 , y2 )
r2
(x, y) b
r1 b r1 b
(x3 , y3 ) (x , y ) (x1 , y1 )
r3 1
b 1
r2
b
(x2 , y2 )
Figure 1.34: Area of a quadrilateral. Figure 1.35: Area of a triangle.
r1 b
(x1 , y1 )
that has a reflex angle. A ray produced from this vertex must hit another
vertex, creating a diagonal, otherwise the polygon would have infinite
area. This diagonal divides the polygon into two sub-polygons. These
two sub-polygons are either both convex or at least one is not convex. In
the latter case, we repeat the argument, finding another diagonal and
creating a new sub-polygon. Eventually, since the number of vertices
is infinite, we end up triangulating the polygon. Moreover, the polygon
can be triangulated in such a way that all triangles inherit the positive
orientation of the original polygon but each neighbouring pair of triangles
have opposite orientations. Applying (1.12) we obtain that the area is
" #
X xi xj
det ,
yi yj
where the sum is over each oriented edge. Since each diagonal occurs
twice, but having opposite orientations, the terms
" # " #
xi xj xj xi
det + det = 0,
yi yj yj yi
" #
1 2
. Let us compute their dot product, find a unit vector in the direction of →
−
a,
3 4
and the angle between the vectors. (There must be either a colon or a semicolon
at the end of each statement. The result will not display if a colon is chosen.)
> with(linalg):
> a:=vector([2,1]);
a := [2, 1]
> b:=vector([1,2]);
b := [1, 2]
> normalize(a);
√ √
2 5 5
[ , ]
5 5
> dotprod(a,b);
4
> angle(a,b);
4
arccos
5
> A:=matrix([[1,2],[3,4]]);
" #
1 2
A :=
3 4
> det(A);
−2
Homework
→
−
Find all vectors →
Problem 1.5.1 −
a ∈ R2 then→
−
a = b.
−3 √
such that →
−
a ⊥ and ||a|| = 13.
2 Problem 1.5.5 (Polarisation Identity) Let
→
−
u ,→
−
v be vectors in R2 . Prove that
Problem 1.5.2 (Pythagorean Theorem) If →
− 1
→
− →
− u •→
−
v = ||→
−
u +→
−
v ||2 − ||→
−
u −→
−
v ||2 .
a ⊥ b , prove that 4
→
− − 2 →
− 2
a + b ||2 = →
||→
− a + b . Problem 1.5.6 Consider two lines on the
plane L1 and L2 with Cartesian equations
L1 : y = m 1 x + b1 and L2 : y = m 2 x + b1 ,
where m 1 6= 0, m 2 6= 0. Using Corollary
22, prove that L1 ⊥ L2 ⇐⇒ m 1 m 2 = −1.
Problem 1.5.3 Let a, b be arbitrary real
.
numbers. Prove that
→−
v •→
−
Problem 1.5.8 Let →
−v,→−
w , be vectors on the →
−
a =→
− w −
v − 2 →w is perpendicular to →
−
w.
→
− →
− →
−
plane, with w 6= 0 . Prove that the vector w
and such that r([a; b]) = Γ. r(a) is the initial point of the curve and r(b) its
terminal point. A curve is closed if its initial point and its final point coincide.
The trace of the curve r is the set of all images of r, that is, Γ. If there exist
t1 6= t2 such that r(t1 ) = r(t2 ) = p, then p is a multiple point of the curve. The
curve is simple if its has no multiple points. A closed curve whose only multiple
points are its endpoints is called a Jordan curve.
1 − t2
t/10 Figure 1.39: x = ,
Figure 1.37: x = sin 2t, y = Figure 1.38: x = 2 cos t, 1 + t2
cos 6t. y = 2t/10 sin t. t − t3
y= .
1 + t2
Figure 1.40: x = (1 +
cos t)/2, y = (sin t)(1 +
cos t)/2.
Graphing parametric equations is a difficult art, and a theory akin to the one
studied for Cartesian equations in a first Calculus course has been developed.
Our interest is not in graphing curves, but in obtaining suitable parametrisa-
tions of simple Cartesian curves. We mention in passing however that Maple
has excellent capabilities for graphing parametric equations. For example, the
commands to graph the various curves in figures 1.37 through 1.40 follow.
> with(plots):
> plot([sin(2*t),cos(6*t),t=0..2*Pi],x=-5..5,y=-5..5);
> plot([2ˆ(t/10)*cos(t),2ˆ(t/10)*sin(t),t=-20..10],x=-5..5,y=-5..5);
> plot([(1-tˆ2)/(1+tˆ2),(t-tˆ3)/(1+tˆ2),t=-2..2],x=-5..5,y=-5..5);
> plot([(1+cos(t))/2,sin(t)*(1+cos(t))/2,t=0..2*Pi],x=-5..5,y=-5..5);
Our main focus of attention will be the following. Given a Cartesian curve
with equation f (x, y) = 0, we wish to find suitable parametrisations for them.
That is, we want to find functions x : t 7→ a(t), y : t 7→ b(t) and an interval I
such that the graphs of f (x, y) = 0 and f (a(t), b(t)) = 0, t ∈ I coincide. These
parametrisations may differ in features, according to the choice of functions and
the choice of intervals.
x = t, y = t2 , t∈R
works for the whole curve. Notice that as t increases, the curve is traversed
from left to right.
√ √
2. If x = t and y = t, then again y = t = ( t)2 = x2 . This works only for
t ≥ 0, and hence the parametrisation
√
x = t, y = t, t ∈ [0; +∞[
gives the half of the curve for which x ≥ 0. As t increases, the curve is
traversed from left to right.
√ √
3. Similarly, if x = − t and y = t, then again y = t = (− t)2 = x2 . This
works only for t ≥ 0, and hence the parametrisation
√
x = − t, y = t, t ∈ [0; +∞[
gives the half of the curve for which x ≤ 0. As t increases, x decreases, and
so the curve is traversed from right to left.
4. If x = cos t and y = cos2 t, then again y = cos2 t = (cos t)2 = x2 . Both x
and y are periodic with period 2π, and so this parametrisation only agrees
with the curve y = x2 when −1 ≤ x ≤ 1. For t ∈ [0; π], the cosine decreases
from 1 to −1 and so the curve is traversed from right to left in this interval.
√ √
Figure 1.41: x = t, y = t2 , Figure 1.42: x = t, y = t, Figure 1.43: x = − t, y =
t ∈ R. t ∈ [0; +∞[. t, t ∈ [0; +∞[.
The identities
(x − 1)2 (y + 2)2
56 Example Give two distinct parametrisations of the ellipse + =
4 9
1.
1. The first parametrisation must satisfy that as t traverses the values in the
interval [0; 2π], one starts at the point (3, −2), traverses the ellipse once
counterclockwise, finishing at (3, −2).
x−1
= cos t =⇒ x = 1 + 2 cos t,
2
and
y+2
= sin t =⇒ y = −2 + 3 sin t.
3
Then
x = 1 + 2 cos t, y = −2 + 3 sin t, t ∈ [0; 2π]
is the desired first parametrisation.
begins and ends at the point (0, 1). To begin at the point (1, 0) we must
π π
make a shift: sin 4πt + , cos 4πt + will start at (1, 0) and
2 2
travel clockwise twice, as t traverses [0; 1]. Hence we may take
π π
x = 1 + 2 sin 4πt + , y = −2 + 3 cos 4πt + , t ∈ [0; 1]
2 2
as our parametrisation. ◭
O′ b b
φ B
b
P
b
θ b
O A
π
Let A be the projection of P on the x-axis. Then ∠OAP = ∠OP O′ = ,
2
′ π
∠OO P = π − φ − θ, ∠P OA = − φ, and OP = (R − ρ) sin(π − φ − θ).
2
Hence
π
x = (OP ) cos ∠P OA = (R − ρ) sin(π − φ − θ) cos( − φ),
2
π
y = (R − ρ) sin(π − φ − θ) sin( − φ).
2
Now
π
x = (R − ρ) sin(π − φ − θ) cos( − φ)
2
= (R − ρ) sin(φ + θ) sin φ
(R − ρ)
= (cos θ − cos(2φ + θ))
2
(R − ρ)
= (R − ρ) cos θ − (cos θ + cos(2φ + θ))
2
= (R − ρ) cos θ − (R − ρ)(cos(θ + φ) cos φ).
ρ ρ
Also, cos(θ + φ) = − cos(π − θ − φ) = − =− and cos φ =
OO′ R−ρ
(R − ρ)θ
cos and so
ρ
(R − ρ)θ
x = (R−ρ) cos θ−(R−ρ)(cos(θ+φ) cos φ = (R−ρ) cos θ+ρ cos ,
ρ
b b d→
−
x
b
b
Parametric
b Curves on the Plane
→
−
x
b
→
−
x + d→
−
x Figure 1.48: Length of a curve.
Given a curve Γ how can we find its length? The idea, as seen in figure 1.48
is to consider the projections dx, dy at each point. The length of the vector
" #
dx
dr =
dy
is È
||dr|| = (dx)2 + (dy)2 .
Hence the length of Γ is given by
Z Z È
||dr|| = (dx)2 + (dy)2 . (1.13)
Γ Γ
and to sum over the closed curve, obtaining a total oriented area of
I " # I
1 x dx 1
det = (xdy − ydx). (1.14)
2 Γ y dy 2 Γ
I
Here denotes integration around the closed curve.
Γ
Furthermore, find an integral expression for the perimeter of this ellipse and find
the area it encloses.
x2 y2 A2 cos2 t B 2 sin2 t
+ = + = cos2 t + sin2 t = 1.
A2 B2 A2 B2
Notice that this parametrisation goes around once the ellipse counter-
clockwise. The perimeter of the ellipse is given by
Z Z È
→ 2π
d−
r = A2 sin2 t + B 2 cos2 t dt.
Γ 0
= πAB.
in figure 1.51. Find the perimeter of the astroid and the area it encloses.
= 6.
We can use Maple (at least version 10) to calculate the above integrals.
For example, if (x, y) = (cos3 t, sin3 t), to compute the arc length we use
the path integral command and to compute " the
# area, we use the line
−y/2
integral command with the vector field .
x/2
> with(Student[VectorCalculus]):
> PathInt( 1, [x,y]=Path( <(cos(t))ˆ3,(sin(t))ˆ3>,0..2*Pi));
> LineInt( VectorField(<-y/2,x/2>), Path( <(cos(t))ˆ3,(sin(t))ˆ3>,0..2*Pi));
◭
We include here for convenience, some Maple commands to compute various arc
lengths and areas.
60 Example To obtain the arc length of the path in figure 1.52, we type
> with(Student[VectorCalculus]):
> PathInt(1, [x,y] = LineSegments( <0,0>, <1,1>, <1,2> ,<2,1>,<3,3>,<4,1>);
To obtain the arc length of the path in figure 1.53, we type
> with(Student[VectorCalculus]):
> PathInt(1, [x,y] = Arc( Circle( <0,0>, 3), Pi/6, Pi/5 ) ;
To obtain the area inside the curve in 1.54
> with(Student[VectorCalculus]):
> LineInt( VectorField(<-y/2,x/2>),
> Path( <(1+cos(t))*(cos(t))+1,(1+cos(t))*(sin(t))+2>,0..2*Pi));
5 5
5 4 4
4 3 b 3
b
3 2 2
b b
2 1 1
b b b
1 b
b
1 2 3 4 5 1 2 3 4 5
1 2 3 4 5
Homework
Problem 1.6.1 A curve is represented Problem 1.6.5 Describe the trace of the
parametrically by x(t) = t3 − 2t, y(t) = parametric curve
t3 + 2t. Find its Cartesian equation.
x sin t
= , t ∈ [0; 4π].
Problem 1.6.2 Give an implicit Cartesian y 2 sin t + 1
equation for the parametric representation
t2 t3
x= , y = . Problem 1.6.6 Consider the plane curve
1 + t5 1 + t5 √ √
defined implicitly by x + y = 1. Give a
suitable parametrisation of this curve, and
Problem 1.6.3 Let a, b, c, d be strictly posi- find its length. The graph of the curve ap-
tive real constants. In each case give an im- pears in figure 1.55.
plicit Cartesian equation for the parametric
representation and describe the trace of the
parametric curve.
1. x = at + b, y = ct + d 1
2. x = cos t, y = 0
3. x = a cosh t, y = b sinh t −1 1
π π −1
4. x = a sec t, y = b tan t, t ∈] − ; [
2 2
Figure 1.55: Problem 1.6.6.
Problem 1.6.4 Parametrise the curve y =
π
log cos x for 0 ≤ x ≤ . Then find its arc
3 Problem 1.6.7 Consider the graph given
length.
parametrically by x(t) = t3 + 1, y(t) =
1 − t2 . Find the area under the graph, over
the x axis, and between the lines x = 1 and ρ. The curve traced out by P as the circle
x = 2. rolls along a straight line, without slipping,
is called a cycloid. Find a parametrisation
of the cycloid.
Problem 1.6.8 Find the arc length of the
curve given parametrically by x(t) = 3t2 ,
y(t) = 2t3 for 0 ≤ t ≤ 1.
t2 (2t + 1)3/2 1 1 P d
x(t) = , y(t) = , t ∈ [− ; + ]. b b
C
2 3 2 2
θ ρ
Find the length of this curve.
b b
R3 = {r = (x, y, z) : x ∈ R, y ∈ R, z ∈ R} .
x y
Having oriented the z axis upwards, we have a choice for the orientation of
the the x and y-axis. We adopt a convention known as a right-handed coordinate
system, as in figure 1.60. Let us explain. In analogy to R2 we put
2 3 2 3 2 3
1 0 0
→
− 6 7 →
− 6 7 →
− 6 7
i = 405 , j = 415 , k = 40 5 ,
0 0 1
→
− →
−
j j
→
− →
−
i i
Figure 1.60:
ℓ3 Right-handed Figure 1.62: Left-handed
system. ℓ system.
2
ℓ1
→
−
62 Definition The dot product of two vectors →
−
a and b in R3 is
→
− →
−
a • b = a1 b1 + a2 b2 + a3 b3 .
→
− − →
−
Just as in R2 , the dot product satisfies →
−
a • b = →
a b cos θ, where θ ∈ [0; π]
is the convex angle between the two vectors.
The Cauchy-Schwarz-Bunyakovsky Inequality takes the form
→
− − →
−
|→
−
a • b | ≤ →
a b =⇒ |a1b1 +a2 b2 +a3 b3 | ≤ (a21 +a22 +a23 )1/2 (b21 +b22 +b23 )1/2 ,
Cauchy’s Inequality,
1/2 √
1 1 1 1 √ 7 7 3
|x+y+z| = x + 2y + 3z ≤ (x +4y +9z )
2 2 2 1/2
1+ + = 27 = .
2 3 4 9 6 2
→
−
64 Definition Let a be a point in R3 and let →−v 6= 0 be a vector in R3 . The
parametric line passing through a in the direction of →
−
v is the set
©
r ∈ R3 : r = a + t→
−
v .
1
65 Example Find the parametric equation of the line passing through 2 and
3
−2
−1 .
0
☞ Given two lines in space, one of the following three situations might arise:
(i) the lines intersect at a point, (ii) the lines are parallel, (iii) the lines are skew
(non-parallel, one over the other, without intersecting, lying on different planes).
See figure 1.63.
→
− →
−
Consider now two non-zero vectors → −a and b in R3 . If →
−
a k b , then the set
→
−
{s→
−a + t b : s ∈ R, t ∈ R} = {λ→−a : λ ∈ R},
→
−
which is a line through the origin. Suppose now that → −
a and b are not parallel.
We saw in the preceding chapter that if the vectors were on the plane, they would
span the whole plane R2 . In the case at hand the vectors are in space, they still
span a plane, passing through the origin. Thus
→
− →
−
{s→
−
a + t b : s ∈ R, t ∈ R, →
−
a 6k b }
is a plane passing through the origin. We will say, abusing language, that two
vectors are coplanar if there exists bi-point representatives of the vector that lie
on the same plane. We will say, again abusing language, that a vector is parallel
to a specific plane or that it lies on a specific plane if there exists a bi-point
representative of the vector that lies on the particular plane. All the above gives
the following result.
66 Theorem Let → −
v ,→
−
w in R3 be non-parallel vectors. Then every vector →
−u of the
form
→
−u = a→−v + b→
−w,
→
−
a, b arbitrary scalars, is coplanar with both v and →
→
− −
w . Conversely, any vector t
→
− →
−
coplanar with both v and w can be uniquely expressed in the form
→
−
t = p→
−
v + q→
−
w.
See figure 1.64.
From the above theorem, if a vector →−a is not a linear combination of two other
→
− →−
vectors b , c , then linear combinations of these three vectors may lie outside the
→
− −
plane containing b , → c . This prompts the following theorem.
→
− −
67 Theorem Three vectors →
−
a , b ,→c in R3 are said to be linearly independent if
→
− →
−
α→−
a + β b + γ→ −c = 0 =⇒ α = β = γ = 0.
Any vector in R3 can be written as a linear combination of three linearly indepen-
dent vectors in R3 .
−
→
arbitrary point on this plane. Since a, b, and c are non-collinear, ab and → −
ac,
−
→
which are coplanar, are non-parallel. Since ax also lies on the plane, we have by
Lemma 66, that there exist real numbers p, q with
→
− −
→
ar = pab + q →
−
ac.
By Chasles’ Rule,
→
− →
−
r =→
−
a + p( b − →
−
a ) + q(→
−
c −→
−
a ),
is the equation of a plane containing the three non-collinear points a, b, and c,
→
−
where →−a , b , and →
−
c are the position vectors of these points. Thus we have the
following theorem. 2 3
a
→
− 6 7
n = 4b 5
c
→
−
v
p→
−
v + q→
−
w →
−
x
ax + by + cz = d
→
−
w
68 Theorem Let → −
u and →−
v be linearly independent vectors. The parametric equa-
tion of a plane containing the point a, and parallel to the vectors →
−
u and →
−v is
given by
→
−r −→−a = p→−u + q→
−
v.
Componentwise this takes the form
x − a1 = pu1 + qv1 ,
y − a2 = pu2 + qv2 ,
z − a3 = pu3 + qv3 .
Adding gives,
(u2 v3 − u3 v2 )(x − a1 ) + (u3 v1 − u1 v3 )(y − a2 ) + (u1 v2 − u2 v1 )(z − a3 ) = 0.
Put
a = u2 v3 − u3 v2 , b = u3 v1 − u1 v3 , c = u1 v2 − u2 v1 ,
and
d = a1 (u2 v3 − u3 v2 ) + a2 (u3 v1 − u1 v3 ) + a3 (u1 v2 − u2 v1 ).
Since v is linearly independent from →
→
− −
u , not all of a, b, c are zero. This gives the
following theorem.
69 Theorem The equation of the plane in space can be written in the form
ax + by + cz = d,
which is the Cartesian equation of the plane. Here a2 + b2 + c2 6= 0, that is, at
2 3
a
→
− 6 7
least one of the coefficients is non-zero. Moreover, the vector n = 4 b 5 is normal
c
to the plane with Cartesian equation ax + by + cz = d.
Proof: We have already proved the first statement. For the second
statement, observe that if →−u and → −v are non-parallel vectors and →−r −
→
− →
− →
−
a = p u + q v is the equation of the plane containing the point a and
parallel to the vectors →
−
u and →−v , then if →
−
n is simultaneously perpendic-
ular to u and v then ( r − a ) n = 0 for →
→
− →
− →
− →
− • →
− −
u •→
−
n =0 = → −
v •→
−n . Now,
since at least one of a, b, c is non-zero, we may assume a 6= 0. The
argument is similar if one of the other letters is non-zero and a = 0. In
this case we can see that
d b c
x = − y − z.
a a a
Put y = s and z = t. Then
2
d b3 c3 2
x− − −
a 6 a7 6 a7
y = s6 7 6 7
4 1 5 + t4 0 5
z 0 1
1
70 Example The equation of the plane passing through the point −1 and
2
2 3
−3
6 7
normal to the vector 4 2 5 is
4
71 Example Find both the parametric equation and the Cartesian equation of
2 3 2 3
1 1
6 7 6 7
the plane parallel to the vectors 415 and 415 and passing through the point
1 0
0
−1 .
2
This gives
s = z − 2, t=y+1−s=y+1−z+2=y−z+3
and
x = s + t = z − 2 + y − z + 3 = y + 1.
Hence the Cartesian equation is x − y = 1. ◭
72 Definition If →
−
n is perpendicular to plane Π1 and →−n ′ is perpendicular to plane
Π2 , the angle between the two planes is the angle between the two vectors →−n and
→
− ′
n .
73 Example
Vectors in Space
y y x y
x x
Figure 1.68: Solid bounded
Figure 1.66: The plane z = Figure 1.67: The plane z =
by the planes z = 1 − x and
1 − x. 1 − y.
z = 1 − y in the first octant.
Solution: ◮
1. This appears in figure 1.66.
2. This appears in figure 1.67.
2 3 2 3
1 0
6 7 6 7
3. The vector 405 is normal to the plane x + z = 1, and the vector 415
1 1
is normal to the plane y + z = 1. If θ is the angle between these
two vectors, then
1·0+0·1+1·1 1 π
cos θ = √ √ =⇒ cos θ = =⇒ θ = .
2 2 2
1 +1 · 1 +1 2 2 3
Homework
→
−
Ǒ
→
−
Problem 1.7.1
→ Vectors
a , b satisfy 0
− →
− →− −
→
a = 13, b = 19, a + b = 24. point 0 .
− −
→ 1
Find → a − b .
74 Definition Let →
−
x ,→
−
y ,→
−
z be vectors in R3 , and let α ∈ R be a scalar. The cross
product × is a closed binary operation satisfying
➊ Anti-commutativity: →
−
x ×→
−
y = −(→
−
y ×→
−
x)
➋ Bilinearity:
(→
−
x +→
−
z )×→
−
y =→
−
x ×→
−
y +→
−
z ×→
−
y and →
−
x × (→
−
z +→
−
y)=→
−
x ×→
−
z +→
−
x ×→
−
y
➎ Right-hand Rule:
→
− →
− →
− →− →
− →
− → − →
− →
−
i × j = k, j × k = i , k × i = j .
It follows that the cross product is an operation that, given two non-parallel
vectors on a plane, allows us to “get out” of that plane.
75 Example Find
2 3 2 3
1 0
6 7 6 7
4 0 5 × 415 .
−3 2
Solution: ◮ We have
→
− →
− →
− →
− →
− →
− →
− →
− →
− →
− →
− →
−
( i −3k)×( j +2k) = i × j +2 i × k −3k × j −6k × k
→
− →
− →
− →
−
= k −2 j +3 i +60
→
− →
− →
−
= 3 i −2 j + k.
Hence 2 3 2 3 2 3
1 0 3
6 7 6 7 6 7
4 0 5 × 415 = 4−25 .
−3 2 1
◭
2 3 2 3
x1 y1
6 7 6 7
76 Theorem Let →
−
x = 4x2 5 and →
−
y = 4y2 5 be vectors in R3 . Then
x3 y3
→
− →
− →
− →
−
x ×→
−
y = (x2 y3 − x3 y2 ) i + (x3 y1 − x1 y3 ) j + (x1 y2 − x2 y1 ) k .
−
×→
y
−
−
y
x →
Cross Product
→
−
θ
−
→
x →
− →
y −
x
→
− →
− →
− →
− →
− →
− →
−
Proof: Since i × i = j × j = k × k = 0 , we only worry about
the mixed products, obtaining,
→
− →
− →
− →
− →
− →
− →
−
x ×→
−
y = (x1 i + x2 j + x3 k ) × (y1 i + y2 j + y3 k )
→
− →
− →
− →
− →
− →
− →
− →
−
= x 1 y2 i × j + x 1 y3 i × k + x 2 y1 j × i + x 2 y3 j × k
→
− →
− →
− →
−
+x3 y1 k × i + x3 y2 k × j
→
− →
− →
− →
− →
− →
−
= (x1 y2 − y1 x2 ) i × j + (x2 y3 − x3 y2 ) j × k + (x3 y1 − x1 y3 ) k × i
→
− →
− →
−
= (x1 y2 − y1 x2 ) k + (x2 y3 − x3 y2 ) i + (x3 y1 − x1 y3 ) j ,
Using the cross product, we may obtain a third vector simultaneously per-
pendicular to two other vectors in space.
77 Theorem → −x ⊥ (→
−
x ×→ −
y ) and →
−
y ⊥ (→−
x ×→ −
y ), that is, the cross product of two
vectors is simultaneously perpendicular to both original vectors.
Proof: We will only check the first assertion, the second verification is
analogous.
→
− →
− →
− →
− →
−
x •(→
−
x ×→
−
y) = (x1 i + x2 j + x3 k )•((x2 y3 − x3 y2 ) i
→
− →
−
+(x3 y1 − x1 y3 ) j + (x1 y2 − x2 y1 ) k )
= x 1 x 2 y3 − x 1 x 3 y2 + x 2 x 3 y1 − x 2 x 1 y3 + x 3 x 1 y2 − x 3 x 2 y1
= 0,
Although the cross product is not associative, we have, however, the following
theorem.
78 Theorem
→
− →
− →
− →
− −
a ×(b ×→
−
c ) = (→
−
a •→
−
c ) b − (→
−
a • b )→
c.
Proof:
→
− →
− →
− →
− →
− →
−
a ×(b ×→
−
c) = (a1 i + a2 j + a3 k ) × ((b2 c3 − b3 c2 ) i +
→
− →
−
+(b3 c1 − b1 c3 ) j + (b1 c2 − b2 c1 ) k )
→
− →
− →
−
= a1 (b3 c1 − b1 c3 ) k − a1 (b1 c2 − b2 c1 ) j − a2 (b2 c3 − b3 c2 ) k
→
− →
− →
−
+a2 (b1 c2 − b2 c1 ) i + a3 (b2 c3 − b3 c2 ) j − a3 (b3 c1 − b1 c3 ) i
→
− →
− →
−
= (a1 c1 + a2 c2 + a3 c3 )(b1 i + b2 j + b3 i )+
→
− →
− →
−
(−a1 b1 − a2 b2 − a3 b3 )(c1 i + c2 j + c3 i )
→
− →
−
= (→
−
a •→−c ) b − (→−a • b )→
−
c,
b b
b
P
D′ b
b
→
−
a × b
A′
→
−
b
b
b
C′
→
−
D
c
b
N
θ
−b
→ b
B′
x A b
b
→
−
a y
C
b
b
M
B
79 Theorem Let (→
−
x ,→
\ −
y ) ∈ [0; π] be the convex angle between two vectors →
−
x and
→
−
y . Then
||→
−
x ×→
−
y || = ||→
−
x ||||→
−
y || sin (→
−
x ,→
\ −
y ).
Proof: We have
||→
−
x ×→
−
y ||2 = (x2 y3 − x3 y2 )2 + (x3 y1 − x1 y3 )2 + (x1 y2 − x2 y1 )2
= y 2 y32 − 2x2 y3 x3 y2 + z 2 y22 + z 2 y12 − 2x3 y1 x1 y3 +
+x2 y32 + x2 y22 − 2x1 y2 x2 y1 + y 2 y12
= (x2 + y 2 + z 2 )(y12 + y22 + y32 ) − (x1 y1 + x2 y2 + x3 y3 )2
= ||→
−
x ||2||→
−
y ||2 − (→
−x •→
−y )2
= ||→
−x ||2||→
−y ||2 − ||→
−
x ||2 ||→
−
y ||2 cos2 (→
−
x ,→
\ −
y)
→
− 2 →− 2 2 →
−\ →
−
= || x || || y || sin ( x , y ),
→
−
80 Corollary Two non-zero vectors →
−
x ,→
−
y satisfy →
−
x ×→
−
y = 0 if and only if they
are parallel.
The following result mixes the dot and the cross product.
→
− −
82 Theorem Let →−a, b, → c , be linearly independent vectors in R3 . The signed
→
−
volume of the parallelepiped spanned by them is (→−
a × b)•→ −c.
Proof: See figure 1.71. The area of the base of the parallelepiped is the
→
−
area of the parallelogram
determined by the vectors →−a and b , which
− →
− −
has area → a × b . The altitude of the parallelepiped is →c cos θ
→
−
where θ is the angle between → −
c and →−
a × b . The volume of the paral-
lelepiped is thus
→− − →
→ − →
− →
− →
−
a × b c cos θ = ( a × b ) • c ,
☞
that
Since we may have used any of the faces of the parallelepiped, it follows
→
− →
− →
−
(→
−
a × b)•→
−
c =(b ×→
−
c)•→
−
a = (→
−
c ×→
−
a)• b.
In particular, it is possible to “exchange” the cross and dot products:
→
− →
− →
−
a •(b ×→ −c ) = (→
−
a × b)•→ −
c
1. Find the Cartesian equation of the plane containing the points A, D ′ , and
M.
4. Suppose that a line through M is drawn cutting the line segment [AC ′ ] in
←−→ ←→
N and the line DD′ in P . Find the parametric equation of MP.
Solution: ◮
0
4. Since P is on the z-axis, P = 0 for some real number z ′ > 0.
z′
←→
The parametric equation of the line MP is thus
2 3
x 1 −1
6 7
y = 3 +s 4−35 =⇒ x = 1−s y = 3−3s, z = sz ′ .
z 0 z′
←→ ←→
Since N is on both MP and AC′ we must have
2 − 2t = 1 − s, 3t = 3 − 3s, t = sz ′ .
1 2
Solving the first two equations gives s = ,t = . Putting this into
3 3
0
′
the third equation we deduce z = 2. Thus P = 0 and the
2
desired equation is
2 3
x 1 −1
6 7
y = 3 +s 4−35 =⇒ x = 1−s, y = 3−3s, z = 2s.
z 0 2
Homework
2 3 2 3
Problem 1.8.1 Prove that 1 1
6 7 6 7
parallel to the vectors 415 and 415 and
→
− →
− →
−
(→
−
a − b ) × (→
−
a + b ) = 2→
−
a × b. 1 0
passing through the point (0, −1, 2), by
→
− finding a normal to the plane.
Problem 1.8.2 Prove that →−
x × → −x = 0
follows from the anti-commutativity of the
cross product.
Problem 1.8.4 Redo example 71, that is, Problem 1.8.6 Let a ∈ R. Find a vector
find the Cartesian equation of the plane of unit length simultaneously perpendicu-
2 3 2 3
0 1 Problem 1.8.15 Consider the plane Π
6 7 6 7
lar to →
−
v = 4−a5 and →
−
w = 4a5. passing through the points A(6, 0, 0),
B(0, 4, 0) and C(0, 0, 3), as shewn in fig-
a 0 ure 1.73 below. The plane Π intersects a
3 × 3 × 3 cube, one of whose vertices is at
Problem 1.8.7 (Jacobi’s Identity) Let →−a , the origin and that has three of its edges on
→
− →
b,− c be vectors in R3 . Prove that the coordinate axes, as in the figure. This
→
− →
− →
− →
− intersection forms a pentagon CP QRS.
→
−
a ×( b ×→−
c )+ b ×(→ −
c ×→ −a )+→
−c ×(→−
a×b)= 0.
−→ −→
Problem 1.8.8 Let →−
x ∈ R3 , ||x|| = 1. Find 1. Find CA × CB.
→
− →
− →
−
||→
−
x × i ||2 + ||→
−
x × j ||2 + ||→
−x × k ||2 . −→ −→
2. Find CA × CB.
→
−
Problem 1.8.9 The vectors → −
a , b are con- 3. Find the parametric equation of the
stant vectors. Solve the equation line LC A joining C and A, with a pa-
→
− →
− →
−
a × (→−x × b ) = b × (→ −
x ×→ −
a ). rameter t ∈ R.
→
− − →
− 4. Find the parametric equation of the
Problem 1.8.10 If →
−
a + b +→
c = 0 , prove line LD E joining D and E, with a pa-
that rameter s ∈ R.
→
− →
− →
−
a × b = b ×→ −
c =→−
c ×→
−a.
5. Find the intersection point between
→
− − the lines LC A and LD E .
Problem 1.8.11 Assume →
−a •( b ×→
c ) 6= 0
and that
→
− →
− 6. Find the equation of the plane Π.
x = α→−
a + β b + γ→−c.
→
−
Find α, β, and γ in terms of →
−
a •(b ×→
−
c ). 7. Find the area of △ABC.
→
− − 8. Find the coordinates of the points P ,
Problem 1.8.12 The vectors → −a , b ,→
c are Q, R, and S.
constant vectors. Solve the system of equa-
tions
→
− 9. Find the area of the pentagon
2→−
x +→−
y ×→ −a = b , 3→ −y +→ −
x ×→ −
a =→−c, CP QRS.
→
− − → −
Problem 1.8.13 Let → −a , b ,→c , d be vec-
3
tors in R . Prove the following vector iden- z
tity, C
→
− →
− →
− → − →
− → − −
(→
−a × b )•(→
−c × d ) = (→
−a •→
−c )( b • d )−(→
−
a • d )( b •→
c ).
P
S
→
− − → −
Problem 1.8.14 Let → −a, b, →c , d , be vec-
3
tors in R . Prove that Q R B y
→
− →
− xA
(b ×→ −c ) • (→
−
a × d)
→
− →
−
+(→
−c ×→ −a)•(b × d)
→
− →
−
+(→
−a × b ) • (→−c × d)
= 0. Figure 1.73: Problem 1.8.15..
for all points a, b in R3 and all scalars λ. Such a linear transformation has a
3 × 3 matrix representation whose columns are the vectors T (i), T (j), and T (k).
Solution: ◮
and also
αu1
L(αu) = L αu2
αu3
α(u1 ) − α(u2 ) − α(u3 )
= α(u1 ) + α(u2 ) + α(u3 )
αu3
u1 − u2 − u3
= α u1 + u2 + u3
u3
u1
= αL u2
u3
= αL(u),
2 3
1 −1 −1
6 7
41 1 1 5.
0 0 1
2 3 2 3
1 2 3 a b c
6 7 6 7
87 Example If A = 44 5 05, and B = 4a b 05, then
6 0 0 a 0 0
2 3 2 3
1+a 2+b 3+c 3 6 9
6 7 6 7
A + B = 44 + a 5 + b 0 5, 3A = 412 15 05 ,
6+a 0 0 18 0 0
2 3 2 3
a + 4b + 6c 2a + 5b 3a 6a 3b c
6 7 6 7
BA = 4 a + 4b 2a + 5b 3a5 , AB = 49a 9b 4c5 .
a 2a 3a 6a 6b 6c
It is an easy exercise to prove that the product of two scaling matrices commutes.
Easy to find counterexamples should convince the reader that the product of two
rotations in space does not necessarily commute.
Homework
2 3
1 0 0 Problem 1.9.3 Consider the n × n matrix
6 7
Problem 1.9.1 Let A = 41 1 05. Find 2 3
1 1 1 1 ... 1 1
1 1 1 60
6 1 1 1 ... 1 17 7
A2 , A3 and A4 . Conjecture and, prove by 6 7
induction, a general formula for An . A=6 6
0 0 1 1 ... 1 1 7.
7
6 . . .. 7
4. . . . . . .. .. . . . . . . .5
0 0 0 0 ... 0 1
Problem 1.9.2 Let A ∈ M3×3 (R) be given
Describe A2 and A3 in terms of n.
by
0 0 1
Proof:
→
− − →
−
1. If D(→−
a , b ,→
c) = → −a •(b ×→ −c ), linearity of the first component
follows by the distributive law for the dot product:
→
− − →
−
D(→
−
a +→
−
a ′, b , →
c) = (→
−a +→−a ′) • ( b × →
−
c)
→
− →
− →
−
= a •(b × c)+ a′•(b ×→
→
− →
− −
c)
→
− →
− →
− →
− ′ →
− →
−
= D( a , b , c ) + D( a , b , c ),
and if λ ∈ R,
→
− − →
− − →
− − →
− −
D(λ→
−
a , b,→
c ) = (λ→
−
a )•( b ×→
c ) = λ((→
−
a )•( b ×→
c )) = λD(→
−
a , b,→
c ).
❑
Observe that
2 3
a1 b1 c1
6 7 →
− →
−
det 4a2 b2 c2 5 = a •(b ×→
−
c) (1.16)
a3 b3 c3
→
− →
− −
→
= →
−
a • (b2 c3 − b3 c2 ) i + (b3 c1 − b1 c3 ) j + (b1 c2 − b2 c1(1.17)
)k
= a1 (b2 c3 − b3 c2 ) + a2 (b3 c1 − b1 c3 ) + a3 (b1 c2 − b2 c1 ) (1.18)
" # " # " #
b2 c2 b1 c1 b1 c1
= a1 det − a2 det + a3 det , (1.19)
b3 c3 b3 c3 b2 c2
which reduces the computation of 3 × 3 determinants to 2 × 2 determinants.
2 3
1 2 3
6 7
92 Example Find det 44 5 65 .
7 8 9
> with(linalg):
> a:=vector([-2,0,1]);
a := [−2, 0, 1]
> b:=vector([-1,3,0]);
b := [−1, 3, 0]
> crossprod(a,b);
[−3, −1, −6]
> dotprod(a,b);
2
> angle(a,b);
√
50
arccos
25
Homework
Problem 1.10.1 Prove that Problem 1.10.2 Prove that
2 3 2 3
1 1 1 a b c
6 7 6 7
det 4 a b c 5 = (b − c)(c − a)(a − b). det 4 b c a5 = 3abc − a3 − b3 − c3 .
a2 b2 c2 c a b
D′
b
A′ C′ N M
b b
b b
′
B b
N
b
b
D
b
b M b
A C D′ b b b b
A
b
Q P
B
−−→ √
Solution: ◮ By the Pythagorean Theorem, AD′ = a 2. Because
−−→ −−→
they are diagonals that belong to parallel faces of the cube, AD′ k BC′ .
Now, M and N are the midpoints of the sides [B ′ B] and [B ′ C ′ ] of
−−→ −−→
△B ′ C ′ B, and hence MN k BC′ by example √ 14. The aforementioned
−−→ 1 −−→′ a 2 −−→
example also gives MN = AD = . In consequence, AD′ k
−−→ 2 2
MN. This means that the four points A, D′ , M, N are all on the same √
′
√ a 2
plane. Hence MND A is a trapezoid with bases of length a 2 and
2
(figure 1.75). From the figure
−−→ √
′ −→ 1 −−→′ −−→ a 2
D Q = AP = AD − MN = .
2 4
Also, by the Pythagorean Theorem,
r Ê √
−−→ −−→2 −−→2 a2 a 5
′ ′ ′
D N = D C + C′ N = a2 + = .
4 2
94 Theorem (Thales’ Theorem) Of two lines are cut by three parallel planes, their
corresponding segments are proportional.
←→ ←→
Proof: See figure 1.76. Given the lines AB and CD, we must prove
that
AE CF
= .
EB FD
←→
Draw line AD cutting plane P 2 in G. The plane containing points A, B,
←→
and D intersects plane P 2 in the line EG. Similarly the plane containing
←→
points A, C, and D intersects plane P 2 in the line GF. Since P 2 and P 3
←→ ← →
are parallel planes, EG k BD, and so by Thales’ Theorem on the plane
(theorem 30),
AE AG
= .
EB GD
←→ ← →
Similarly, since P 1 and P 2 are parallel, AC k GF and
CF AG
= .
FD GD
It follows that
AE CF
= ,
EB FD
as needed to be shewn. ❑
P1 b
C
A D′
b
b
A′
b
C′
b
′
B
G F b b
P2E b b
b
b
D b
b M N b
A C
b
P 3B b b
D B
D′
b
A′ D ′′ C′
b b
b
B ′ Theorem in 3D.
Figure 1.76: Thales’ Figure 1.77: Example 95.
b
′′
A b
D C ′′ b
b
b
b
N b
M B ′′
b b
A C
b
Homework
Proof: We only provide the proof for the second statement, as the proof
for the first is similar. Cut any two such solids by horizontal planes
that produce cross sections of area A(x) and cA(x), where c > 0 is the
Z height x above Za fixed base.
constant of proportionality, at an arbitrary
x2 x2
From elementary calculus, we know that A(x)dx and cA(x)dx
x1 x1
give the volume of the portion of each solid
Z cut by all horizontal
Z planes
x2 x2
as x runs over some interval [x1 ; x2 ]. As A(x)dx = c A(x)dx
x1 x1
the corresponding volumes must also be proportional. ❑
97 Example Use Cavalieri’s Principle in order to deduce that the area enclosed
x2 y2
by the ellipse with equation 2 + = 1, a > 0, b > 0, is πab.
a b
b b
b b b
x b
a
x
b
b b
a b
b
r b
a
x
2πa3 4πa3
It follows that the volume of the sphere is 2 = .◭
3 3
Essentially the same method of proof as Cavalieri’s Principle gives the next
result.
99 Theorem (Pappus-Guldin Rule) The area of the lateral surface of a solid of rev-
olution is equal to the product of the length of the generating curve on the side
of the axis of revolution and the length of the path described by the centre of
gravity of the generating curve under a full revolution. The volume of a solid of
revolution is equal to the product of the area of the generating plane on one side
of the revolution axis and the length of the path described by the centre of gravity
of the area under a full revolution about the axis.
r b
b
R
100 Example Since the centre of gravity of a circle is at its centre, by the Pappus-
Guldin Rule, the surface area of the torus with the generating circle having radius
r, and radius of gyration R (as in figure 1.82) is (2πr)(2πR) = 4π 2 rR. Also, the
volume of the solid torus is (πr 2 )(2πR) = 2π 2 r 2 R.
Homework
Problem 1.12.1 Use the Pappus-Guldin of a right circular cone with base radius r
Rule to find the lateral area and the volume and height h.
102 Definition The rectilinear angle of a dihedral angle is the angle whose sides
are perpendicular to the edge of the dihedral angle at the same point, each on
each of the faces. See figure 1.84.
All the rectilinear angles of a dihedral angle measure the same. Hence the
measure of a dihedral angle is the measure of any one of its rectilinear angles.
103 Definition The opening of three or more planes that meet at a common point
is called a polyhedral angle or solid angle. In the particular case of three planes,
we use the term trihedral angle. The common point is called the vertex of the
polyhedral angle. Each of the intersecting lines of two consecutive planes is
called an edge of the polyhedral angle. The portion of the planes lying between
consecutive edges are called the faces of the polyhedral angle. The angles formed
by adjacent edges are called face angles. A polyhedral angle is said to be convex
if the section made by a plane cutting all its edges forms a convex polygon.
104 Theorem The sum of any two face angles of a trihedral angle is greater than
the third face angle.
A D P
b b b
C A5 b
b
b
b b b
B
X b
A2
b VW Z
A1 b
b
Y
Ak b Ak−1
Ak+1
105 Theorem The sum of the face angles of any convex polyhedral angle is less
than 2π radians.
Proof: Let the polyhedral angle have n faces and vertex V . Let the
faces be cut by a plane, intersecting the edges at the points A1 , A2 ,
. . . An , say. An illustration can be seen in figure 1.86, where for con-
venience, we have depicted only five edges. Observe that the polygon
A1 A2 · · · An is convex and that the sum of its interior angles is π(n − 2).
Thus
X X
V Ak Ak−1 + ∠V Ak Ak+1 > ∠Ak−1 Ak Ak+1 = π(n − 2).
1≤k≤n 1≤k≤n
Also, X
V Ak Ak+1 + ∠V Ak+1 Ak + ∠Ak V Ak+1 = πn,
1≤k≤n
since this is summing the sum of the angles of the n triangles of the
faces. But clearly
X X
V Ak Ak+1 = ∠V Ak+1 Ak ,
1≤k≤n 1≤k≤n
since one sum adds the angles in one direction and the other in the
opposite direction. For the same reason,
X X
V Ak Ak−1 = ∠V Ak Ak+1 .
1≤k≤n 1≤k≤n
Hence
X X
∠Ak V Ak+1 = πn − (V Ak Ak+1 + ∠V Ak+1 Ak )
1≤k≤n 1≤k≤n
X
= πn − (V Ak Ak+1 + ∠V Ak Ak−1 )
1≤k≤n
< πn − π(n − 2)
= 2π,
as we needed to shew. ❑
Since n ≥ 3 and m ≥ 3, the above inequality only holds for five pairs (n, m).
Appealing to Euler’s Formula for polyhedrons, which states that V + F = E + 2,
where V is the number of vertices, F is the number of faces, and E is the number
of edges of a polyhedron, we obtain the values in the following table.
4 3 8 12 6 Hexahedron or Cube.
3 4 6 12 8 Octahedron.
5 3 20 30 12 Dodecahedron.
3 5 12 30 20 Icosahedron.
Thus there are at most five Platonic solids. That there are exactly five can
be seen by explicit construction. Figures 1.88 through 1.92 depict the Platonic
solids.
107 Definition Given a point (x, y, z) in Cartesian coordinates, its spherical coor-
dinates are given by
Here φ is the polar angle, measured from the positive z-axis, and θ is the az-
imuthal angle, measured from the positive x-axis. By convention, 0 ≤ θ ≤ 2π
and 0 ≤ φ ≤ π.
108 Definition If a plane intersects with a sphere, the intersection will be a circle.
If this circle contains the centre of the sphere, we call it a great circle. Otherwise
we talk of a small circle. The axis of any circle on a sphere is the diameter of the
sphere which is normal to the plane containing the circle. The endpoints of such
a diameter are called the poles of the circle.
☞ The radius of a great circle is the radius of the sphere. The poles of a great
circle are equally distant from the plane of the circle, but this is not the case in
a small circle. By the pole of a small circle, we mean the closest pole to the
plane containing the circle. A pole of a circle is equidistant from every point of
the circumference of the circle.
109 Definition Given the centre of the sphere, and any two points of the surface
of the sphere, a plane can be drawn. This plane will be unique if and only if the
points are not diametrically opposite. In the case where the two points are not
diametrically opposite, the great circle formed is split into a larger and a smaller
arc by the two points. We call the smaller arc the geodesic joining the two points.
If the two points are diametrically opposite then every plane containing the line
forms with the sphere a great circle, and the arcs formed are then of equal length.
In this case we take any such arc as a geodesic.
b
B
φ
ρ
ρ cos φ
ρ
θ
si
n
φ
☞ A spherical triangle has then six angles: three vertex angles ∠A, ∠B, ∠C,
and three arc angles, ∠a, ∠b, ∠c. Observe that if O is the centre of the sphere
then
−→ −→ −→ −→ −→ −→
∠a = ∠ OB, OC , ∠b = ∠ OC, OA , ∠c = ∠ OA, OB ,
and
−→ −→ −→ −→ −→ −→ −→ −→ −→ −→ −→ −→
∠A = ∠ OA × OB, OA × OC , ∠B = ∠ OB × OC, OB × OA , ∠C = ∠ OC × OA, OC × OB .
x1 x2 + y1 y2 + z1 z2 = ρ2 cos ∠ (AOB) .
Therefore we obtain
that is,
cos a cos b + sin a sin b cos C = cos c.
❑
112 Theorem Let I be the dihedral angle of two adjacent faces of a regular poly-
hedron. Then
I cos π
n
sin = π .
2 sin m
Proof: See figure 1.94. Let AB be the edge common to the two adjacent
faces, C and D the centres of the faces; bisect AB at E, and join CE
and DE; CE and DE will be perpendicular to AB, and the angle CED
is the angle of inclination of the two adjacent faces; we shall denote it
by I. In the plane containing CE and DE draw CO and DO at right
angles to CE and DE respectively, and meeting at O; about O as centre
describe a sphere meeting OA, OC, OE at a, c, e respectively, so that
113 Theorem Let r and R be, respectively, the radii of the inscribed and circum-
scribed spheres of a regular polyhedron. Then
a π I a I π
r= cot tan , R= tan tan .
2 m 2 2 2 n
Here a is the length of any edge of the polyhedron, and I is the dihedral angle of
any two faces.
mF ra2 π
therefore the volume of the polyhedron is . cot
m12
ma2 π
Furthermore, the area of one face of the polyhedron is cot , and therefore
4 m
mF a2 π
the surface area of the polyhedron is cot .
4 m
Homework
Problem 1.14.1 The four vertices of a reg- Problem 1.14.4 Consider an octahedron
ular tetrahedron are whose edge√measures a. Shew that its vol-
Ǒ
Ǒ √
a3 2
1 −1/2 ume is , its surface area is 2a2 3,
√ 3
V1 = 0 , V2 = 3/2 , and √
that the radius of the inscribed sphere
0 0 a 6
is .
Ǒ
Ǒ 6
−1/2 0
√
V3 = − 3/2 , V4 = 0 . Problem 1.14.5 Consider a dodecahedron
√ whose edge measures a. Shew that its vol-
0 2
a3 √
What is the cosine of the dihedral angle be- ume is 15 + 7 5 , its surface area is
tween any pair of faces of the tetrahedron? È 4 √
3a2 25 + 10 5, and that the radius of the
Ê É
Problem 1.14.2 Consider a tetrahedron a 1
inscribed sphere is 10 + 22 .
4 5
√ measures a. Shew that its vol-
whose edge
a3 2 √
ume is , its surface area is a2 3, and
12 Problem 1.14.6 Consider an icosahedron
that
√ the radius of the inscribed sphere is whose edge measures a. Shew that its vol-
5a3 √
a 6
.
12 ume is 3 + 5 , its surface area is
√ 12
5a2 3, and that
a √ √ of the inscribed
the radius
Problem 1.14.3 Consider a cube whose sphere is 5 3 + 15 .
edge measures a. Shew that its volume is 12
a3 , its surface area is 6a2 , and that the ra-
a
dius of the inscribed sphere is .
2
Let us start with the plane. Recall that if a, b, c are real numbers, not all
2 3
a
6 7
zero, then the Cartesian equation of a plane with normal vector 4 b 5 and passing
c
through the point (x0 , y0 , z0 ) is
x − x0 = pu1 + qv1 ,
y − y0 = pu2 + qv2 ,
z − z0 = pu3 + qv3 .
114 Definition A surface S consisting of all lines parallel to a given line ∆ and
passing through a given curve Γ is called a cylinder. The line ∆ is called the
directrix of the cylinder.
x y y
x
Figure 1.95: Circular cylinder x2 + y2 = 1. Figure 1.96: The parabolic cylinder z = y2 .
115 Example Figure 1.95 shews the cylinder with Cartesian equation x2 +y 2 = 1.
One starts with the circle x2 + y 2 = 1 on the xy-plane and moves it up and down
the z-axis. A parametrisation for this cylinder is the following:
The method of parametrisation utilised above for the cylinder is quite useful
when doing parametrisations in space. We refer to it as the method of cylindrical
coordinates. In general, we first find the polar coordinates of x, y in the xy-plane,
and then lift (x, y, 0) parallel to the z-axis to (x, y, z):
x = r cos θ, y = r sin θ, z = z.
z
(r cos θ, r sin θ, z)
b
r
θ
x y
116 Example Figure 1.96 shews the parabolic cylinder with Cartesian equation
z = y 2 . One starts with the parabola z = y 2 on the yz-plane and moves it up and
down the x-axis. A parametrisation for this parabolic cylinder is the following:
x = u, y = v, z = v2 , u ∈ R, v ∈ R.
117 Example Figure 1.98 shews the hyperbolic cylinder with Cartesian equation
x2 − y 2 = 1. One starts with the hyperbola x2 − y 2 on the xy-plane and moves
it up and down the z-axis. A parametrisation for this parabolic cylinder is the
following:
x = ± cosh v, y = sinh v, z = u, u ∈ R, v ∈ R.
We need a choice of sign for each of the portions. We have used the fact that
cosh2 v − sinh2 v = 1. The Maple commands to graph this surface are:
> with(plots):
> implicitplot3d(xˆ2-yˆ2=1,x=-10..10,y=-10..10,z=-10..10,numpoints=5001);
> plot3d({[-cosh(s),sinh(s),t],[cosh(s),sinh(s),t]},
> s=-2..2,t=-10..10,numpoints=5001,axes=boxed);
118 Definition Given a point Ω ∈ R3 (called the apex) and a curve Γ (called the
generating curve), the surface S obtained by drawing rays from Ω and passing
through Γ is called a cone.
☞ A B
In practice, if the Cartesian equation of a surface can be put into the form
f( , ) = 0, where A, B, C, are planes secant at exactly one point, then the
C C
surface is a cone, and its apex is given by A = 0, B = 0, C = 0.
z 2 = x2 + y 2
x 2 y 2
is a cone, as its equation can be put in the form −1 = 0. Considering
+
z z
the planes x = 0, y = 0, z = 0, the apex is located at (0, 0, 0). The graph is shewn
in figure 1.100.
121 Example Find the equation of the surface of revolution generated by revolv-
ing the hyperbola
x2 − 4z 2 = 1
about the z-axis.
x2 y2 z2
Figure 1.100: Cone + 2 = 2.
a2 b c
È
(0, 0, z), that is x2 + y 2 . We must have
È p
x2 + y 2 = 1 + 4z 2 ,
which is to say
x2 + y 2 − 4z 2 = 1.
This surface is called a hyperboloid of one sheet. See figure 1.104.
Observe that when z = 0, x2 + y 2 = 1 is a circle on the xy plane. When
x = 0, y 2 − 4z 2 = 1 is a hyperbola on the yz plane. When y = 0,
x2 − 4z 2 = 1 is a hyperbola on the xz plane.
122 Example The circle (y − a)2 + z 2 = r 2 , on the yz plane (a, r are positive real
numbers) is revolved around the z-axis, forming a torus T . Find the equation of
this torus.
Rearranging
p
x2 + y 2 + z 2 − a2 − r 2 = 2asgn(y − a) r 2 − z 2 ,
or
(x2 + y 2 + z 2 − (a2 + r 2 ))2 = 4a2 r 2 − 4a2 z 2
z z
x y y y
x x
Figure 1.101: Paraboloid z = Figure 1.102: Hyperbolic Figure 1.103: Two-sheet hy-
x2 y2 x2 y2 z2 x2 y2
+ 2. paraboloid z = 2 − 2 perboloid 2 = 2 + 2 +1.
a2 b a b c a b
x = u, y = v, z = u2 − v 2 , u ∈ R, v ∈ R.
x2 y2
127 Example Let a, b, c be strictly positive real numbers. The surface + +
a2 b2
2 2 2
z x y
= 1 is called an ellipsoid. For z = 0, + 1 is an ellipse on the xy
c2 a2 b2
2 2
x z
plane.When y = 0, + = 1 is an ellipse on the xz plane. When x = 0,
a2 c2
2 2
z y
+ = 1 is an ellipse on the yz plane. See figure 1.105. We may parametrise
c2 b2
the ellipsoid using spherical coordinates:
Homework
Problem 1.15.1 Find the equation of the Problem 1.15.7 Shew that the surface in
surface of revolution S generated by revolv- R3 implicitly defined by
ing the ellipse 4x2 + z 2 = 1 about the z-
axis. x4 + y 4 + z 4 − 4xyz(x + y + z) = 1
Problem 1.15.6 Demonstrate that the sur- Problem 1.15.11 (Putnam Exam 1970)
face in R3 Determine, with proof, the radius of the
2
+y 2 +z 2 largest circle which can lie on the ellipsoid
S : ex − (x + z)e−2xz = 0,
x2 y2 z2
implicitly defined, is a cylinder. 2
+ 2 + 2 = 1, a > b > c > 0.
a b c
and such that r([a; b]) = Γ. r(a) is the initial point of the curve and r(b) its
terminal point. A curve is closed if its initial point and its final point coincide.
The trace of the curve r is the set of all images of r, that is, Γ. The length of the
curve is Z
→
d−r .
Γ
x y
130 Example Find a parametric representation for the curve resulting by the in-
tersection of the plane 3x + y + z = 1 and the cylinder x2 + 2y 2 = 1 in R3 .
131 Example Let a, b, c be strictly positive real numbers. Consider the the region
©
R = (x, y, z) ∈ R3 : |x| ≤ a, |y| ≤ b, z = c .
A point P moves along the ellipse
x2 y2
+= 1, z =c+1
a2 b2
once around, and acts as a source light projecting a shadow of R onto the xy-
plane. Find the area of this shadow.
b
b
b
Chapter 1
b b
b b
b
x y
Resizing to a region
©
R = (x, y, z) ∈ R3 : |x| ≤ a, |y| ≤ b, z = c ,
and an ellipse
x2 y2
+
= 1, z = c + 1
a2 b2
we use instead of c + 1, a(c + 1) (parallel to the x-axis) and b(c + 1)
(parallel to the y-axis), so that the area shadowed is
πab(c + 1)2 + 4ab(c + 1)2 + 4abc(c + 1) = c2 ab(π + 12) + 16abc + 4ab.
◭
Homework
Problem 1.16.1 Let C be the curve in R3 Problem 1.16.4 Give a parametrisation for
defined by the part of the ellipsoid
Problem 1.16.2 Consider the surfaces in Problem 1.16.5 Let P be the point (2, 0, 1)
R3 implicitly defined by and consider the curve C : z = y 2 on the
yz-plane. As a point Q moves along C , let
z−x2 −y 2 −1 = 0, z+x2 +y 2 −3 = 0. R be the point of intersection of ← →
PQ and
Describe, as vividly as possible these sur- the xy-plane. Graph all points R on the
faces and their intersection, if they at all in- xy-plane.
tersect. Find a parametric equation for the
curve on which they intersect, if they at all Problem 1.16.6 Let a be a real number pa-
intersect. rameter, and consider the planes
P1 : ax + y + z = −a,
Problem 1.16.3 Consider the space curve
2 3 P2 : x − ay + az = −1.
t4
6 1 + t2 7 Let l be their intersection line.
→
− 6 t3 7
r : t 7→ 6 7
6 1 + t2 7 . 1. Find a direction vector for l.
4 5
t2 2. As a varies through R, l describes a
1 + t2 surface S in R3 . Let (x, y, z) be the
point of intersection of this surface
Let tk , 1 ≤ k ≤ 4 be non-zero real num- and the plane z = c. Find an equa-
bers. Prove that → −
r (t1 ), →
−
r (t2 ), →
−
r (t3 ), and tion relating x and y.
→
−r (t4 ) are coplanar if and only if
3. Find the volume bounded by the two
1 1 1 1 planes, x = 0, and x = 1, and the
+ + + = 0.
t1 t2 t3 t4 surface S as c varies.
→
− →
−
133 Definition If →
−
a and b are two vectors in Rn their vector sum →
−
a + b is defined
by the coordinatewise addition
2 3
a1 + b1
6 7
→
− 6 a2 + b2 7
→
− 6 7
a + b =6 .. 7. (1.20)
6 . 7
4 5
an + bn
135 Definition The standard ordered basis for Rn is the collection of vectors
{→
−
e 1, →
−
e 2, . . . , →
−
e n, }
with 2 3
0
6.7
6.7
6.7
→
− 6 7
ek=6 7
6 17
6.7
6.7
4.5
0
(a 1 in the k slot and 0’s everywhere else). Observe that
2 3
α1
6 7
X
n 6 α2 7
6 7
αk →
−
e k = 6 . 7.
6 .. 7
k=1 4 5
αn
→
−
136 Definition Given vectors →
−
a , b of Rn , their dot product is
→
− X
n
→
−
a•b = ak bk .
k=1
☞ The above proof works not just for Rn but for any vector space (cf. below)
that has an inner product.
→
−
138 Corollary (Triangle Inequality) Let → −
a and b be any two vectors in Rn . Then
we have − → −
→
− → − →
a + b ≤ a + b .
Proof:
→
− →
− →
−
||→
−
a + b ||2 = (→
−a + b )•(→ −a + b)
→
− →
− →−
= →
−a •→
−a + 2→ −a•b + b•b
→
− →
−
≤ ||→
−a ||2 + 2||→−
a |||| b || + || b ||2
→
−
= (||→
−a || + || b ||)2 ,
from where the desired result follows. ❑
Again, the preceding proof is valid in any vector space that has a norm.
This expression agrees with the geometry in the case of the dot product for R2
and R3 .
X
n
Solution: ◮ Using CBS on (ak bk )ck once we obtain
k=1
!1/2 !1/2
X
n X
n X
n
ak bk ck ≤ a2k b2k c2k .
k=1 k=1 k=1
!1/2
X
n
Using CBS again on a2k b2k we obtain
k=1
!1/2 !1/2
X
n X
n X
n
ak bk ck ≤ a2k b2k c2k
k=1 k=1 k=1
!1/4 !1/4 !1/2
Xn X
n X
n
≤ a4k b4k c2k ,
k=1 k=1 k=1
X
n
141 Lemma Let ak > 0, qk > 0, with qk = 1. Then
k=1
!1/x
X
n X
n
lim log qk a x
k = qk log ak .
x→0
k=1 k=1
√ a1 + a2 + · + an
n
a1 a2 · · · an ≤ .
n
giving
√ a1 + a2 + · + an
n
a1 a2 · · · an ≤ ,
n
as wanted. ❑
1 · 3 · 5 · · · · (2n − 1) < nn .
For, by AMGM,
n n
1 + 3 + 5 + · · · + (2n − 1) n2
1 · 3 · 5 · · · · (2n − 1) < = = nn .
n n
Notice that since the factors are unequal we have strict inequality.
144 Definition Let a1 > 0, a2 > 0, . . . , an > 0. Their harmonic mean is given by
n
.
1 1 1
+ + ··· +
a1 a2 an
145 Corollary (Harmonic Mean-Geometric Mean Inequality) Let b1 > 0, b2 > 0, . . . , bn >
0. Then
n
≤ (b1 b2 · · · bn )1/n .
1 1 1
+ + ··· +
b1 b2 bn
1
Proof: This follows by putting ak = in Theorem 142 . For then
bk
1 1 1
1/n + + ··· +
1 1 1 b b2 bn
··· ≤ 1 .
b1 b2 bn n
146 Corollary (Harmonic Mean-Arithmetic Mean Inequality) Let b1 > 0, b2 > 0, . . . , bn >
0. Then
n b1 + b2 + · · · + bn
≤ .
1 1 1 n
+ + ··· +
b1 b2 bn
X
n
s n2
≥
k=1
s − ak n−1
and
X
n
ak n
≥ .
k=1
s − ak n−1
s
Solution: ◮ Put bk = . Then
s − ak
X
n
1 X
n
s − ak
= =n−1
k=1
bk k=1
s
s ak
Since −1= , we have
s − ak s − ak
X
n n
X
ak s
= −1
k=1
s − ak k=1
s − ak
Xn
s
= −n
k=1
s − ak
2
n
≥ −n
n−1
n
= .
n−1
◭
Homework
Problem 1.17.1 The Arithmetic Mean Geo- Problem 1.17.3 (USAMO 1978) Let
metric Mean Inequality says that if ak ≥ 0 a, b, c, d, e be real numbers such that
then
a1 + a2 + · · · + an a+b+c+d+e = 8, a2 +b2 +c2 +d2 +e2 = 16.
(a1 a2 · · · an )1/n ≤ .
n
Equality occurs if and only if a1 = a2 = Maximise the value of e.
. . . = an . In this exercise you will follow the
steps of a proof by George Pólya. Problem 1.17.4 Find all positive real num-
1. Prove that ∀x ∈ R, x ≤ ex−1 . bers
2. Put a1 ≤ a2 ≤ . . . ≤ an
nak
Ak = , such that
a1 + a2 + · · · + an
and Gn = a1 a2 · · · an . Prove that X
n
X
n
X
n
b b
b2 − a2
ε
b
b
(a1 , a2 )
b b
b1 − a1
97
z z
ε
b (a1 , a2 , a3 ) Some Topology
x y x y
Figure 2.3: Open ball in R3 . Figure 2.4: Open box in R3 .
©
153 Example The ellipsoidal region (x, y) ∈ R2 : x2 + 4y 2 < 4 is open in R2 .
The reader will recognise that open boxes, open ellipsoids and their unions and
finite intersections are open sets in Rn .
Homework
Problem 2.1.1 Determine whether the fol- 4. D = {(x, y) ∈ R2 : x2 ≤ y ≤ x}
lowing subsets of R2 are open, closed, or
5. E = {(x, y) ∈ R2 : xy > 1}
neither, in R2 .
6. F = {(x, y) ∈ R2 : xy ≤ 1}
1. A = {(x, y) ∈ R2 : |x| < 1, |y| < 1}
2. B = {(x, y) ∈ R2 : |x| < 1, |y| ≤ 1} 7. G = {(x, y) ∈ R2 : |y| ≤ 9, x < y 2 }
Problem 2.1.2 (Putnam Exam 1969) Let finite collection of open disks in R2 whose
p(x, y) be a polynomial with real coeffi- union contains a set E ⊆ R2 . Shew that
cients in the real variables x and y, defined there is a pairwise disjoint subcollection
over the entire plane R2 . What are the pos- Dk , k ≥ 1 in F such that
sibilities for the image (range) of p(x, y)?
[
n
E⊆ 3Dj .
Problem 2.1.3 (Putnam 1998) Let F be a j =1
f : A → Rm .
158 Example The level curves of the surface f (x, y) = x2 +y 2 (an elliptic paraboloid)
are the concentric circles
x2 + y 2 = c, c > 0.
Homework
Problem 2.2.1 Sketch the level curves for 6. (x, y) →
7 sin(x2 + y 2 )
the following maps. 7. (x, y) → 7 cos(x2 − y 2 )
1. (x, y) 7→ x + y
2. (x, y) 7→ xy Problem 2.2.2 Sketch the level surfaces
for the following maps.
3. (x, y) 7→ min(|x|, |y|)
1. (x, y, z) 7→ x + y + z
4. (x, y) 7→ x3 − x 2. (x, y, z) 7→ xyz
5. (x, y) 7→ x2 + 4y 2 3. (x, y, z) 7→ min(|x|, |y|, |z|)
5. (x, y, z) 7→ x2 + 4y 2 7. (x, y, z) 7→ x2 + y 2 + z 2
2.3 Limits
We will start with the notion of limit.
The notions of infinite limits, limits at infinity, and continuity at a point, are
analogously defined. Limits in more than one dimension are perhaps trickier to
find, as one must approach the test point from infinitely many directions.
x2 y
160 Example Find lim .
(x,y)→(0,0) x2 + y 2
and hence
x2 y
lim = 0.
(x,y)→(0,0) x2 + y 2
The Maple commands to graph this surface and find this limits ap-
pear below. Notice that Maple is unable to find the limit and so returns
unevaluated.
> with(plots):
> plot3d(xˆ2*y/(xˆ2+yˆ2),x=-10..10,y=-10..10,axes=boxed,color=xˆ2+yˆ2);
> limit(xˆ2*y/(xˆ2+yˆ2),x=0,y=0);
◭
x5 y 3
161 Example Find lim .
(x,y)→(0,0) x6 + y 4
x2 y x5 y3
Figure 2.5: (x, y) 7→ . Figure 2.6: (x, y) 7→ .
x2 + y2 x6 + y4
Solution: ◮ Either |x| ≤ |y| or |x| ≥ |y|. Observe that if |x| ≤ |y|,
then
x5 y 3 y8
≤ = y4 .
x6 + y 4 y4
If |y| ≤ |x|, then
x5 y 3 x8
2
x6 + y 4 ≤ x6 = x .
Thus
x5 y 3
4 2 4 2
x6 + y 4 ≤ max(y , x ) ≤ y + x −→ 0,
1+x+y
162 Example Find lim .
(x,y)→(0,0) x2 − y 2
Solution: ◮ When y = 0,
1+x
→ +∞,
x2
as x → 0. When x = 0,
1+y
→ −∞,
−y 2
as y → 0. The limit does not exist. ◭
xy 6
163 Example Find lim .
(x,y)→(0,0) x6 + y 8
xy 6 1
= → +∞,
x6 + y8 2t2
sin(x4 ) + sin(y 4 )
165 Example Find lim p .
(x,y)→(0,0) x4 + y 4
sin x − y
166 Example Find lim .
(x,y)→(0,0) x − sin y
both exist. These are called the iterated limits of f as (x, y) → (x0 , y0 ). The
following possibilities might occur.
1. If lim f (x, y) exists, then each of the iterated limits lim lim f (x, y)
(x,y)→(x0 ,y0 ) y→y0 x→x0
and lim lim f (x, y) exists.
x→x0 y→y0
2. If the iterated limits exist and lim lim f (x, y) 6= lim lim f (x, y)
y→y0 x→x0 x→x0 y→y0
then lim f (x, y) does not exist.
(x,y)→(x0 ,y0 )
3. It may occur that lim lim f (x, y) = lim lim f (x, y) , but that
y→y0 x→x0 x→x0 y→y0
lim f (x, y) does not exist.
(x,y)→(x0 ,y0 )
4. It may occur that lim f (x, y) exists, but one of the iterated limits
(x,y)→(x0 ,y0 )
does not.
Homework
Problem 2.3.1 Sketch the domain of defi- sin xy
p Problem 2.3.5 Find lim .
nition of (x, y) 7→ 4 − x2 − y 2 . (x,y )→(0,2) x
Problem 2.3.2 Sketch the domain of defi- Problem 2.3.6 For what c will the function
nition of (x, y) 7→ log(x + y). ¨ p
1 − x2 − 4y 2 , if x2 + 4y 2 ≤ 1,
f (x, y) =
c, if x2 + 4y 2 > 1
Problem 2.3.3 Sketch the domain of defi-
1
nition of (x, y) 7→ 2 . be continuous everywhere on the xy-plane?
x + y2
x−y
Does lim exist?. x2 y 2
(x,y )→(0,0) x+y but still lim does
(x,y )→(0,0) x 2 y 2 + (x − y)2
not exist.
g(x) − g(a)
lim = g ′ (a)
x→a x−a
exists. The limit condition above is equivalent to saying that
identically.
Solution: ◮ We have
→
− − →
− →
− →
−
f (→
−
x + h,→
y + k ) − f (→
−
x ,→
−
y) = (→
−x + h )•(→
−
y + k)−→−
x •→
−
y
→
− →
− →
− → − →
→ −
= x•y + x•k + h•y + h•k −→
→
− →
− − −
x •→
−
y
→
− →
− →
− →
− →
− →
−
= x•k + h•y + h•k.
→
− →− →
− → −
As ( h , k ) → ( 0 , 0 ),
→we
have by the Cauchy-Buniakovskii-Schwarz
− •→
→ − − →− →
−
inequality, | h k | ≤ h k = o h , which proves the asser-
tion. ◭
Just like in the one variable case, differentiability at a point, implies continuity
at that point.
and so
f (x) − f (a) → 0,
as x → a, proving the theorem. ❑
Homework
Problem 2.4.1 Let L : R3 → R3 be a linear Problem 2.4.2 → Let
f : Rn → R, n ≥
transformation and →
− −
1, f ( x ) = x be the usual norm in Rn ,
− 2
with → x = → −x •→−
x . Prove that
R3 → R3
F : → . →
−
−
x 7→ →
−
x × L(→
−
x) x •→
−v
Dx (f )(→
−
v ) = →
− ,
x
Shew that F is differentiable and that
→
−
→
− →
− →
− for →
−
x 6= 0 , but that f is not differentiable
Dx (F )( h ) = →
−
x × L( h ) + h × L(→
−
x ). →
−
at 0 .
∂fi
Here fi : Rn → R. The partial derivative (x) is defined as
∂xj
To find partial derivatives with respect to the j-th variable, we simply keep
the other variables fixed and differentiate with respect to the j-th variable.
∂f
(x, y, z) = 1 + 3y 2 z 3 ,
∂x
∂f
(x, y, z) = 2y + 6xyz 3 ,
∂y
and
∂f
(x, y, z) = 3z 2 + 9xy 2 z 2 .
∂z
column of the Jacobi matrix. Let f ′ (x) = (Jij ), and observe that
2 3
J1j
6 7
6 J2j 7
6 7
Dx (f )(→
−
e j) = 6 . 7 .
6 .. 7
4 5
Jnj
and put y = x + ε→
−
e j , ε ∈ R. Notice that
176 Example Let f (ρ, θ, z) = (ρ cos θ, ρ sin θ, z) be the function which changes
from cylindrical coordinates to Cartesian coordinates. We have
2 3
cos θ −ρ sin θ 0
6 7
f ′ (ρ, θ, z) = 4 sin θ ρ cos θ 05 .
0 0 1
177 Example Let f (ρ, φ, θ) = (ρ cos θ sin φ, ρ sin θ sin φ, ρ cos φ) be the function
which changes from spherical coordinates to Cartesian coordinates. We have
2 3
cos θ sin φ ρ cos θ cos φ −ρ sin φ sin θ
6 7
f ′ (ρ, φ, θ) = 4 sin θ sin φ ρ sin θ cos φ ρ cos θ sin φ 5 .
cos φ −ρ sin φ 0
Observe that f is not continuous at (0, 0) (f (0, 0) = 0 but f (x, y) = 1 for values
arbitrarily close to (0, 0)), and hence, it is not differentiable there. We have
∂f ∂f
however, (0, 0) = (0, 0) = 1. Thus even if both partial derivatives exist at
∂x ∂y
(0, 0) is no guarantee that the function will be differentiable at (0, 0). You should
also notice that both partial derivatives are not continuous at (0, 0).
∂2 π
180 Example Let f (u, v, w) = eu v cos w. Determine f (u, v, w) at (1, −1, ).
∂u∂v 4
Solution: ◮ We have
∂2 ∂
(eu v cos w) = (eu cos w) = eu cos w,
∂u∂v ∂u
√
e 2
which is at the desired point. ◭
2
Solution: ◮ We have
∂ ∂ xy log z
z xy = e
∂x ∂x
xy ∂z
= y log z + z xy ,
z ∂x
∂ ∂
(xy)z ez log xy
=
∂x ∂x
∂z z
= log xy + (xy)z ,
∂x x
∂ ∂z
xy 2 z 3 = y 2 z 3 + 3xy 2 z 2 ,
∂x ∂x
Hence, at (1, 1, 1) we have
∂z ∂z ∂z 1
+1+1+3 = 0 =⇒ =− .
∂x ∂x ∂x 2
Similarly,
∂ xy ∂ xy log z
z = e
∂y ∂y
xy ∂z
= x log z + z xy ,
z ∂y
∂ ∂
(xy)z =ez log xy
∂y ∂y
∂z z
= log xy + (xy)z ,
∂y y
∂ ∂z
xy 2 z 3 = 2xyz 3 + 3xy 2 z 2 ,
∂y ∂y
∂z ∂z ∂z 3
+1+2+3 = 0 =⇒ =− .
∂y ∂y ∂y 4
◭
Just like in the one-variable case, we have the following rules of differentia-
tion. Let A ⊆ Rn , B ⊆ Rm be open sets f, g : A → Rm , α ∈ R, be differentiable
on A, h : B → Rl be differentiable on B, and f (A) ⊆ B. Then we have
Solution: ◮ We have
2 3
ev uev
6 7
f ′ (u, v) = 4 1 1 5,
v u
and " #
′ 2x 1 0
h (x, y) = .
0 1 1
Observe also that
2 3
ey+z (x2 + y)ey+z
6 7
f ′ (h(x, y)) = 4 1 1 5.
y+z x2 + y
Hence
" # " #
3 2 u(x, y) xz
Solution: ◮ Put g : R → R , g(x, y) = = . Observe
v(x, y y+z
that h = f ◦ g. Now,
" #
′ z 0 x
g (x, y) = ,
0 1 1
f ′ (u, v) = 2u ev ,
f ′ (h(x, y)) = 2xz ey+z .
Thus
∂h ∂h ∂h
(x, y) (x, y) (x, y) = h′ (x, y)
∂x ∂y ∂z
= (f ′ (g(x, y)))(g ′ (x, y))
2 3
z 0 x7 .
6
= 2xz ey+z 4 5
0 1 1
= 2xz 2 ey+z 2x2 z + ey+z
Under certain conditions, the interval of integration in the above theorem need
not be compact.
Z Z
+∞ sin x π +∞ sin2 x
186 Example Given that dx = , compute dx.
0 x 2 0 x2
Z +∞ sin2 ax
Solution: ◮ Put I(a) = dx, with a ≥ 0. Differentiating
x2 0
both sides with respect to a, and making the substitution u = 2ax,
Z +∞ 2x sin ax cos ax
I ′ (a) = dx
Z0+∞ x2
sin 2ax
= dx
Z0+∞ x
sin u
= du
u
π0
= .
2
Integrating each side gives
π
I(a) = a + C.
2
π
Since I(0) = 0, we gather that C = 0. The desired integral is I(1) = .
2
◭
Homework
Problem 2.5.1 Let f : [0; +∞[×]0; +∞[→ Problem 2.5.5 Suppose g : R → R is con-
2
R, f (r, t) = tn e−r /4t , where n is a con- tinuous and a ∈ R is a constant. Find the
stant. Determine n such that partial derivatives with respect to x and y
of
∂f 1 ∂ ∂f Z x2 y
= 2 r2 .
∂t r ∂r ∂r f : R2 → R, f (x, y) = g(t) dt.
a
Rn → Rm
f :
x 7→ f (x)
Rn → R
f :
x 7→ f (x)
Proof: Let
R → Rn
c:
t 7→ c(t)
be a curve lying on this surface. Choose t0 so that c(t0 ) = x. Then
(f ◦ c)(t0 ) = f (c(t)) = K,
which translates to
(∇f (x))•(c′ (t0 )) = 0.
Since c′ (t0 ) is tangent to the surface and its dot product with ∇f (x) is
0, we conclude that ∇f (x) is normal to the surface. ❑
190 Example Find the direction of the greatest rate of increase of f (x, y, z) =
xyez at the point (2, 1, 2).
f (x, y, z) = x + y 2 − z 2 .
2 3
1
6 7
∇f (x, y, z) = 4 2y 5
−2z
which is
2 3
1
6 7
4 4 5
−6
or
x + 4y − 6z = −9.
2
193 Example If f (x, y, z) = (x2 , y 2 , yez ) then
2
divf (x) = 2x + 2y + 2yzez .
2
195 Example If f (x, y, z) = (x2 , y 2 , yez ) then
2
curlf ((x, y, z)) = ∇ × f (x, y, z) = (ez )i.
f (x + t→
−
v ) − f (x)
v f (x) = lim
D→
− .
t→0 t
198 Theorem Let A ⊆ Rn be open and let f : A → R be a scalar field, and assume
→
−
that f is differentiable in A. Let →
−
v ∈ Rn \ { 0 } be such that → −
x + t→−v ∈ A for
sufficiently small t ∈ R. Then the directional derivative of f in the direction of →
−
v
→
−
at the point x is given by
∇f (x)•→−
v.
Solution: ◮ We have
2 3
3x2
6 7
∇f (x, y, z) = 4 3y 2 5
−2z
and so
∇f (x, y, z)•→
−
v = 3x2 + 6y 2 − 6z.
Homework
Problem 2.6.1 Let f (x, y, z) = xey z . Find Problem 2.6.7 Use a linear approximation
of the function f (x, y) = ex cos 2y at (0, 0)
(∇f )(2, 1, 1). to estimate f (0.1, 0.2).
2 3
xz Problem 2.6.8 Prove that
6 xy 7
Problem 2.6.2 Let f (x, y, z) = 4e 5.
∇ • (u × v) = v • (∇ × u) − u • (∇ × v).
z
Find
(∇ × f )(2, 1, 1). Problem 2.6.9 Find the point on the sur-
face
Problem 2.6.3 Find the tangent plane to
2x2 + xy + y 2 + 4x + 8y − z + 14 = 0
x2
the surface − y 2 − z 2 = 0 at the point
2 for which the tangent plane is 4x + y − z =
(2, −1, 1).
0.
201 Example Using Einstein’s Summation convention, the dot product of two
vectors →
−
x ∈ Rn and →
−
y ∈ Rn can be written as
X
n
→
− • →
−
x y = xi yi = <xt yt = .
i=1
→
− − → −
202 Example Given that ai , bj , ck , dl are the components of vectors in R3 , →
−
a , b ,→
c, d
respectively, what is the meaning of
<ai bi ck dk =?
Solution: ◮ We have
X
3
→
− − X
→
3
→
− − •→ −
<ai bi ck dk = = ai bi <ck dk = = →
−
a • b <ck dk = = →
−
a•b ck dk = (→
−
a • b )(→
c d ).
i=1 k=1
203 Example Using Einstein’s Summation convention, the ij-th entry (AB)ij of
the product of two matrices A ∈ Mm×n (R) and B ∈ Mn×r (R) can be written as
X
n
(AB)ij = Aik Bkj = <Ait Btj = .
k=1
204 Example Using Einstein’s Summation convention, the trace tr (A) of a square
X
n
matrix A ∈ Mn×n (R) is tr (A) = Att = <Att =.
t=1
tr (AB) = tr (BA) .
Solution: ◮ We have
and
from where the assertion follows, since the indices are dummy variables
and can be exchanged. ◭
X
3
207 Example It is easy to see that <δik δkj = = δik δkj = δij .
k=1
209 Definition (Levi-Civitta’s Alternating Tensor) The symbol εjkl is defined as fol-
lows: 8
>
> 0 if {j, k, l} =
6 {1, 2, 3}
> !
>
>
< −1 1 2 3
if is an odd permutation
εjkl = j k l
>
> !
>
> 1 2 3
>
: +1 if is an even permutation
j k l
☞
Also,
In particular, if one subindex is repeated we have εrrs = εrsr = εsrr = 0.
210 Example Using the Levi-Civitta alternating tensor and Einstein’s summation
→
− →
− →
−
convention, the cross product can also be expressed, if i = →−
e1 , j = →−
e2 , k =
→
−
e3 , then
→
−
x ×→ −
y = <εjkl (ak bl )→
−
ej = .
211 Example If A = [aij ] is a 3×3 matrix, then, using the Levi-Civitta alternating
tensor,
det A = <εijk a1i a2j a3k = .
Homework
Problem 2.7.1 Let →
−
x ,→
−
y ,→
−
z be vectors in R3 . Demonstrate that
<xi yi zj = = (→
− y )→
x •→
− −
z.
2.8 Extrema
We now turn to the problem of finding maxima and minima for vector functions.
As in the one-variable case, the derivative will provide us with information about
the extrema, and the “second derivative” will provide us with information about
the nature of these extreme points.
To define an analogue for the second derivative, let us consider the following.
Let A ⊂ Rn and f : A → Rm be differentiable on A. We know that for fixed
x0 ∈ A, Dx0 (f ) is a linear transformation from Rn to Rm . This means that we
have a function
A → L (Rn , Rm )
T : ,
x 7→ Dx (f )
Rn × Rn → L (Rn , Rm )
Bx 0 :
(x1 , x2 ) 7→ Dx20 (f )(x1 , x2 )
is well defined, and linear in each variable x1 and x2 , that is, it is a bilinear
function.
Just as the Jacobi matrix was a handy tool for finding a matrix representa-
tion of Dx (f ) in the natural bases, when f maps into R, we have the following
analogue representation of the second derivative.
f (x, y, z) = xy 2 z 3 .
Then
2 3
0 2yz 3 3y 2 z 2
6 7
H(x,y,z) f = 4 2yz 3 2xz 3 6xyz 2 5
3y 2 z 2 6xyz 2 6xy 2 z
From the preceding example, we notice that the Hessian is symmetric, as the
∂2 ∂2
mixed partial derivatives f = f , etc., are equal. This is no coinci-
∂x∂y ∂y∂x
dence, as guaranteed by the following theorem.
We are now ready to study extrema in several variables. The basic theorems
resemble those of one-variable calculus. First, we make some analogous defini-
tions.
∀x ∈ Bx0 (r′ ), f (x1 ) ≤ f (x), we say that f (x1 ) is a local maximum of f . A point
is called an extreme point if it is either a local minimum or local maximum. A
point t is called a critical point if f is differentiable at t and Dt (f ) = 0. A critical
point which is neither a maxima nor a minima is called a saddle point.
R2 → R
f :
(x, y) 7→ x2 + xy + y 2 + 2x + 3y
Solution: ◮ We have
" #
2x + y + 2
(∇f )(x, y) = ,
x + 2y + 3
2x + y + 2 = 0,
x + 2y + 3 = 0,
1 4
which yields x = − , y = − . Now,
3 3
" #
2 1
H(x,y) f = ,
1 2
" #
2 1
which is positive definite, since ∆1 = 2 > 0 and ∆2 = det =
1 2
1 4
3 > 0. Thus x0 = − , − is a relative minimum and we have
3 3
7 1 4
− = f − ,− ≤ f (x, y) = x2 + xy + y 2 + 2x + 3y.
3 3 3
◭
R3 → R
f : .
(x, y, z) 7→ x2 + y 2 + 3z 2 − xy + 2xz + yz
Solution: ◮ We have
2 3
2x − y + 2z
6 7
(∇f )(x, y, z) = 4 2y − x + z 5 ,
6z + 2x + y
Solution: ◮ We have
" # " #
3x2 + ay 0
(∇f )(x, y) = 2
= =⇒ 3x2 = −ay, 3y 2 = ax.
−3y + ax 0
a
If y = 0 then x = 0, so again (0, 0) is a critical point. If y = − then
3
3 a 2 a a a
x= · − = so ,− is a critical point.
a 3 3 3 3
Now,
" #
6x a
Hf (x,y) = =⇒ ∆1 = 6x, ∆2 = −36xy − a2 .
a −6y
2
a 0), a∆1 = 0, ∆2 = −a . If a 6= 0 thenthere
At (0,
a
is a saddle point.
a
2
At , − , ∆1 = 2a, ∆2 = 3a , whence ,− will be a local
3 3 3 3
minimum if a > 0 and a local maximum if a < 0. ◭
Homework
Problem 2.8.1 Determine the nature of the Problem 2.8.10 Find the extrema of
critical points of f (x, y) = y 2 −2x2 y+4x3 + f (x, y, z) = x2 y + y 2 z + 2x − z.
20x2 .
Problem 2.8.8 Determine the nature of the Problem 2.8.15 Find the minimum of
critical points of
√ 2
f (x, y, z) = 4x2 z − 2xy − 4x2 − z 2 + y. 144 − 16x2 p
F (x, y) = (x−y)2 + − 4 − y2 ,
3
Problem 2.8.9 Find the extrema of
for −3 ≤ x ≤ 3, −2 ≤ y ≤ 2.
f (x, y, z) = x2 + y 2 + z 2 + xyz.
Free to photocopy and distribute 128
Chapter 2
∇f (x0 ) = λ∇g(x0 ).
☞ The above theorem only locates extrema, it does not say anything concern-
ing the nature of the critical points found.
and
f (x, y) = 1 − y 2 − y 2 = 1 − 2y 2 ≤ 1.
Thus (±1, 0) are maximum points and (0, ±1) are minimum points. ◭
223 Example Find the maximum and the minimum points of f (x, y) = 4x + 3y,
subject to the constraint x2 + 4y 2 = 4, using Lagrange multipliers.
224 Example Let a > 0, b > 0, c > 0. Determine the maximum and minimum
x y z x2 y2 z2
values of f (x, y, z) = + + on the ellipsoid 2 + 2 + 2 = 1.
a b c a b c
x2 y2
Solution: ◮ We use Lagrange multipliers. Put g(x, y, z) = + +
a2 b2
2
z
− 1. Then
c2
2 3 2 3
1/a 2x/a2
6 7 6 7
∇f (x, y, z) = λ∇g(x, y, z) ⇐⇒ 4 1/b 5 = λ 4 2y/b2 5 .
1/c 2z/c2
a b c x2 y2
It follows that λ 6= 0. Hence x = ,y = ,z = . Since 2 + 2 +
2λ √ 2λ 2λ a b
z2 3 3
= 1, we deduce = 1 or λ = ± . Since a, b, c are positive, f
c2 4λ2 2
will have a maximum when all x, y, z are positive and a minimum when
all x, y, z are negative. Thus the maximum is when
a b c
x = √ ,y = √ ,z = √ ,
3 3 3
and
3 √
f (x, y, z) ≤ √ = 3
3
and the minimum is when
a b c
x = −√ , y = −√ , z = −√ ,
3 3 3
and
3 √
f (x, y, z) ≥ − √ = − 3.
3
225 Example Let a > 0, b > 0, c > 0. Determine the maximum volume of the
parallelepiped with sides parallel to the axes that can be enclosed inside the
x2 y2 z2
ellipsoid 2 + 2 + 2 = 1.
a b c
Solution: ◮ Let 2x, 2y, 2z, be the dimensions of the box. We must
x2
maximise f (x, y, z) = 8xyz subject to the constraint g(x, y, z) = 2 +
a
y2 z2
+ 2 − 1. Using Lagrange multipliers,
b2 c
2 3 2 3
8yz 2x/a2
6 7 6 7 x y z
∇f (x, y, z) = λ∇g(x, y, z) ⇐⇒ 48xz 5 = λ 4 2y/b2 5 =⇒ 4yz = λ , 4xz = λ , 4xy = λ .
2
a2 b2 c2
8xy 2z/c
Multiplying the first inequality by x, the second by y, the third by z, and
adding,
x2 y2 z2 x2 y2 z2
4xyz = λ , 4xyz = λ , 4xyz = λ , =⇒ 12xyz = λ + + = λ.
a2 b2 c2 a2 b2 c2
Hence
λ x2 y2 z2
=λ . =λ =λ
3 a2 b2 c2
6 0, then
If λ = 0, then 8xyz = 0, which minimises the volume. If λ =
a b c
x= √ , y= √ , z= √ ,
3 3 3
and the maximum value is
abc
8xyz ≤ 8 √ .
3 3
x2 y2 z2
2 2 2 1/3 + +
x y z 2 b2 c2 = 1 =⇒ 8xyz ≤ √
8
(x2 y 2 z 2 )1/3 = (abc)2/3 · · ≤ (abc)2/3 · a (abc).
a2 b2 a2 3 3 3 3
◭
Homework
Problem 2.9.1 A closed box (with six outer Problem 2.9.8 Find the axes of the ellipse
faces), has fixed surface area of S square
units. Find its maximum volume using La- 5x2 + 8xy + 5y 2 = 9.
grange multipliers. That is, subject to the
constraint 2ab + 2bc + 2ca = S, you must Problem 2.9.9 Optimise f (x, y, z) = x +
maximise abc. y + z subject to x2 + y 2 = 2, and x + z = 1.
Problem 2.9.2 Consider the problem of Problem 2.9.10 Let x, y be strictly positive
finding the closest point P ′ on the plane real numbers with x + y = 1. What is the
Π : ax + by + cz = d, a, b, c non-zero √
maximum value of x + xy?
constants with a + b + c 6= d to the point
P (1, 1, 1). In this problem, you will do this
in three essentially different ways. Problem 2.9.11 Let a, b be positive
real constants. Maximise f (x, y) =
1. Do this by a geometric argument, ar-
xa e−x y b e−y on the triangle in R2 bounded
guing the the point P ′ closest to P
by the lines x ≥ 0, y ≥ 0, x + y ≤ 1.
on Π is on the perpendicular passing
through P and P ′ .
2. Do this by means of Lagrange multi- Problem 2.9.12 Determine the extrema of
pliers, by minimising a suitable func- f (x, y) = cos2 x+cos2 y subject to the con-
π
tion f (x, y, z) subject to the con- straint x − y = .
4
straint g(x, y, z) = ax + by + cz − d.
3. Do this considering the uncon-
Problem 2.9.13 Determine the extrema of
strained
extrema of a suitable
func-
d − ax − by f (x, y, z) = x − 2y + 2z subject to the con-
tion h x, y, . straint x2 + y 2 + z 2 = 1.
c
Problem 2.9.3 Given that x, y are positive Problem 2.9.14 Find the points on the
real numbers such that x4 +81y 4 = 36 find curve determined by the equations
the maximum of x + 3y.
x2 + xy + y 2 − z 2 = 1, x2 + y 2 = 1
Problem 2.9.4 If x, y, z are positive real which are closest to the origin.
1
numbers such that x2 y 3 z = 2 , what is
6
the minimum value of f (x, y, z) = 2x + Problem 2.9.15 Does there exist a polyno-
3y + z? mial in two variables with real coefficients
p(x, y) such that p(x, y) > 0 for all x and
Problem 2.9.5 Find the maximum and the y and that for all real numbers c > 0 there
minimum values of f (x, y) = x2 + y 2 sub- exists (x0 , y0 ) ∈ R2 such that p(x0 , y0 ) =
ject to the constraint 5x2 + 6xy + 5y 2 = 8. c?
x1 , x2 , . . . , x n
In other words, dxj1 ∧ dxj2 ∧ · · · ∧ dxjk (a1 , a2 , . . . , ak ) is the xj1 xj2 . . . xjk compo-
nent of the signed k-volume of a k-parallelotope in Rn spanned by a1 , a2 , . . . , ak .
133
Differential Forms
☞ Anti-commutativity yields
da ∧ da = 0.
232 Example
g(x, y, z, w) = x + y 2 + z 3 + w4
is a 0-form in R4 .
X
n
∂f
df = dxi .
i=1
∂xi
Furthermore, if
ω = f (x1 , x2 , . . . , xn )dxj1 ∧ dxj2 ∧ · · · ∧ dxjk
is a k-form in Rn , the exterior derivative dω of ω is the (k + 1)-form
dω = df (x1 , x2 , . . . , xn ) ∧ dxj1 ∧ dxj2 ∧ · · · ∧ dxjk .
241 Example Consider the transformation of coordinates xyz into uvw coordi-
nates given by
z y+z
u = x + y + z, v = , w= .
y+z x+y+z
Then
du = dx + dy + dz,
z y
dv = − dy + dz,
(y + z)2 (y + z)2
y+z x x
dw = − dx + dy + dz.
(x + y + z)2 (x + y + z)2 (x + y + z)2
Multiplication gives
zx y(y + z)
du ∧ dv ∧ dw = − −
(y + z)2 (x + y + z)2 (y + z)2 (x + y + z)2
z(y + z) xy
+ − dx ∧ dy ∧ dz
(y + z)2 (x + y + z)2 (y + z)2 (x + y + z)2
2 2
z − y − zx − xy
= dx ∧ dy ∧ dz.
(y + z)2 (x + y + z)2
3.2 Zero-Manifolds
242 Definition A 0-dimensional oriented manifold of Rn is simply a point x ∈ Rn ,
with a choice of the + or − sign. A general oriented 0-manifold is a union of
oriented points.
244 Example Let M = −{(1, 0, 0)} ∪ +{(1, 2, 3)} ∪ −{(0, −2, 0)}2 be an oriented
0-manifold, and let ω = x + 2y + z 2 . Then
Z
ω = −ω((1, 0, 0)) + ω(1, 2, 3) − ω(0, 0, 3) = −(1) + (14) − (−4) = 17.
M
2
Do not confuse, say, −{(1, 0, 0)} with −(1, 0, 0) = (−1, 0, 0). The first one means that the
point (1, 0, 0) is given negative orientation, the second means that (−1, 0, 0) is the additive inverse
of (1, 0, 0).
3.3 One-Manifolds
245 Definition A 1-dimensional oriented manifold of Rn is simply an oriented smooth
curve Γ ∈ Rn , with a choice of a + orientation if the curve traverses in the di-
rection of increasing t, or with a choice of a − sign if the curve traverses in the
direction of decreasing t. A general oriented 1-manifold is a union of oriented
curves.
→
− dx
If f : R2 → R2 and if d→
−
r = , the classical way of writing this is
dy
Z
→
−
f • d→
−
r.
Γ
whence Z Z 2
ω = (3t3 + 2t2 )dt
Γ −1
2
2 3
= t3 + t4
3 4 −1
69
= .
4
What would happen if we had given the curve above a different parametri-
sation? First observe that the curve travels from (−1, 1) to (2, 4) on the
parabola
√ y = x2 .√These conditions are met with the parametrisation
x = t − 1, y = ( t − 1)2 , t ∈ [0; 9]. Then
√ √ √ √ √
xydx + (x + y)dy = ( t − 1)3 d( t − 1) + (( t − 1) + ( t − 1)2 )d( t − 1)2
√ √ √
= (3( t − 1)3 + 2( t − 1)2 )d( t − 1)
1 √ √
= √ (3( t − 1)3 + 2( t − 1)2 )dt,
2 t
whence
Z Z √
1 9 √
ω = √ (3( t − 1)3 + 2( t − 1)2 )dt
Γ
0 22 t 3/2 9
3t 7t 5t √
= − + − t
4 3 2 0
69
= ,
4
as before.
To solve this problem using Maple you may use the code below.
> with(Student[VectorCalculus]):
> LineInt( VectorField( <x*y,x+y> ), Path( <t,tˆ2>, t=-1..2));
◭
☞ It turns out that if two different parametrisations of the same curve have
the same orientation, then their integrals are equal. Hence, we only need to
worry about finding a suitable parametrisation.
where Γ is the line segment from (0, 0) to (1, 1) in the positive direction.
= 2π.
To solve this problem using Maple you may use the code below.
> with(Student[VectorCalculus]):
> LineInt( VectorField( <(x+y)/(xˆ2+yˆ2),(x-y)/(xˆ2+yˆ2)> ),
> LineSegments(<1,1>,<-1,1>,<-1,-1>,<1,-1>,<1,1>));
◭
To solve this problem using Maple you may use the code below.
> with(Student[VectorCalculus]):
> LineInt( VectorField( <yˆ2,xˆ2> ),Ellipse(9*xˆ2 + 4*yˆ2 -36));
◭
Solution: ◮ The lines passing through the given points have equations
−x − 4
LAB : y = , and LBC : y = −4x + 6. On LAB
3
Ê √
È
1 2 x 10dx
x||dx|| = x (dx)2 + (dy)2 = x 1 + − dx = ,
3 3
and on LBC
È È √
2
x||dx|| = x (dx)2 + (dy)2 = x( 1 + (−4) )dx = x 17dx.
Hence Z Z Z
x||dx|| = x||dx|| + x||dx||
Γ LAB √ LBC
Z 2 Z 1 √
x 10dx
= + x 17dx
√
−1 3√ 2
10 3 17
= − .
2 2
To solve this problem using Maple you may use the code below.
> with(Student[VectorCalculus]):
> PathInt( x, [x,y]=LineSegments( <-1,-1>, <2,-2>,<1,2> ) );
◭
Homework
2 b
Chapter 3
1
-3 -2 -1 1 2 3
b -1
-2 b
-3
Z
Problem 3.3.1 Consider xdx + ydy and ing at O(0, 0) going on a straight line to
Z C A 4 cos π6 , 4 sin π6 and continuing on an
xy||dx||. arc of a circle to B 4 cos π5 , 4 sin π5 .
B
C
Z b
and
X
n
∂f
d(dω) = d ∧ dxk
k=1
∂xk
Xn X
n
∂2f
= ∧ dxj ∧ dxk
k=1 j=1
∂xj ∂xk
X
n
∂2f ∂2f
= − dxj ∧ dxk
1≤j≤k≤n
∂xj ∂xk ∂xk ∂xj
= 0,
But
d da ∧ dxj1 ∧ dxj2 ∧ · · · dxjp = dda ∧ dxj1 ∧ dxj2 ∧ · · · dxjp
+da ∧ d dxj1 ∧ dxj2 ∧ · · · dxjp
= da ∧ d dxj1 ∧ dxj2 ∧ · · · dxjp ,
dF = ω.
xdx + ydy
is exact, since
1
xdx + ydy = d (x2 + y 2 ) .
2
ydx + xdy
is exact, since
ydx + xdy = d (xy) .
dF = ω.
We have a choice here of integrating either the first, or the second ex-
pression. Since integrating the second expression (with respect to y) is
easier, we find
ey
F (x, y) = + φ(x),
1 + x2
where φ(x) is a function depending only on x. To find it, we differentiate
the obtained expression for F with respect to x and find
∂F 2xey
=− + φ′ (x).
∂x (1 + x2 )2
∂F
Comparing this with our first expression for , we find
∂x
2x
φ′ (x) = ,
(1 + x2 )2
that is
1
φ(x) = − + c,
1 + x2
where c is a constant. We then take
ey − 1
F (x, y) = + c.
1 + x2
◭
dF = ω = y 2 z 3 dx + 2xyz 3 dy + 3xy 2 z 2 dz ?
Solution: ◮ We have
dω = (2yz 3 dy + 3y 2 z 2 dz) ∧ dx
+(2yz 3 dx + 2xz 3 dy + 6xyz 2 dz) ∧ dy
+(3y 2 z 2 dx + 6xyz 2 dy + 6xy 2 zdz) ∧ dz
= 0,
∂F ∂F ∂F
dF = dx + dy + dz = y 2 z 3 dx + 2xyz 3 dy + 3xy 2 z 2 dz.
∂x ∂y ∂z
Then
∂F
= y 2 z 3 =⇒ F = xy 2 z 3 + a(y, z),
∂x
∂F
= 2xyz 3 =⇒ F = xy 2 z 3 + b(x, z),
∂y
∂F
= 3xy 2 z 2 =⇒ F = xy 2 z 3 + c(x, y),
∂z
Comparing these three expressions for F , we obtain F (x, y, z) = xy 2 z 3 .
◭
where Γ is the closed polygon with vertices at A = (0, 0), B = (5, 0), C = (7, 2),
D = (3, 2), E = (1, 1), traversed in the order ABCDEA.
Solution: ◮ Since
∂ ∂
(x2 − y) = −1 = (y 2 − x),
∂y ∂x
the form is exact, and since this is a closed simple path, the integral is
0. ◭
3.5 Two-Manifolds
262 Definition A 2-dimensional oriented manifold of R2 is simply an open set (re-
gion) D ∈ R2 , where the + orientation is counter-clockwise and the − orientation
is clockwise. A general oriented 2-manifold is a union of open sets.
☞ In this section, unless otherwise noticed, we will choose the positive ori-
entation for the regions considered. This corresponds to using the area form
dxdy.
is the area of D.
263 Theorem (Fubini’s Theorem) Let D = [a; b] × [c; d], and let f : A → R be
continuous. Then
Z Z b Z d Z d Z b
f dA = f (x, y)dy dx = f (x, y)dx dy
a c c a
D
Fubini’s Theorem allows us to convert the double integral into iterated (single)
integrals.
264 Example Z
Z Z 1 3
xydA = xydy dx
3 !
[0;1]×[2;3] 0 2
Z
1 xy 2
= dx
0 2 2
Z 1
9x
= − 2x dx
2
1 0
5x2
=
4 0
5
= .
4
Notice that if we had integrated first with respect to x we would have obtained
the same result:
Z Z Z 1 !
3 1 3 x2 y
xydx dy = dy
2 0 2 2 0
Z 3 y
= dx
2 2 3
2
y
=
4 2
5
= .
4
Also, this integral is “factorable into x and y pieces” meaning that
Z Z 1 Z 3
xydA = xdx ydy
[0;1]×[2;3] 0 2
1 5
=
2 2
5
=
4
To solve this problem using Maple you may use the code below.
> with(Student[VectorCalculus]):
> int(x*y, [x,y]=Region(0..1,2..3));
Solution: ◮ The lines passing through the given points have equations
−x − 4 3x + 1
LAB : y = , LBC : y = −4x + 6, LCA : y = . Now, we
3 2
draw the region carefully. If we integrate first with respect to y, we
must divide the region as in figure 3.3, because there are two upper
lines which the upper value of y might be. The lower point of the dashed
line is (1, −5/3). The integral is thus
Z Z Z Z
1 (3x+1)/2 2 −4x+6 11
x y dy dx + x y dy dx = − .
−1 (−x−4)/3 1 (−x−4)/3 8
To solve this problem using Maple you may use the code below.
> with(Student[VectorCalculus]):
> int(x*y, [x,y]=Triangle(<-1,-1>,<2,-2>,<1,2>);
◭
267 Example Consider the region inside the parallelogram P with vertices at A :
(6, 3), B : (8, 4), C : (9, 6), D : (7, 5), as in figure 3.5. Find
Z
xy dxdy.
P
2 b 2 b 5 Chapter 3b b
b b
4
1 1
b
3
-3 -2 -1 1 2 3 -3 -2 -1 1 2 3 1
b -1 b -1 b
b
b b 1 2 3 4 5 6 7 8 9 10
-2 -2
-3 -3
x 3
LCD : y = + ,
2 2
LDA : y = 2x − 9.
The integral is thus
Z 4 Z 2y Z 5 Z (y+12)/2 Z 6 Z (y+12)/2 409
xy dxdy+ xy dxdy+ xy dxdy = .
3 (y+9)/2 4 (y+9)/2 5 2y−3 4
To solve this problem using Maple you may use the code below. Notice
that we have split the parallelogram into two triangles.
> with(Student[VectorCalculus]):
> int(x*y, [x,y]=Triangle(<6,3>,<8,4>,<7,5>))
> + int(x*y, [x,y]=Triangle(<8,4>,<9,6>,<7,5>));
◭
where
D = {(x, y) ∈ R2 |x ≥ 0, x2 + y 2 ≤ 1}.
Solution: ◮ We have
y √
0 ≤ y ≤ 4, ≤x≤ y =⇒ 0 ≤ x ≤ 2, x2 ≤ y ≤ 2x.
2
We then have
Z Z √ Z Z
4 y 2 2x
y/x
e dx dy = ey/x dy dx
0 y/2 x2
Z02
= xey/x 2x
x2 dx
Z02
= (xe2 − xex ) dx
0
= 2e2 − (2e2 − e2 + 1)
= e2 − 1
Z
√ √
271 Example Evaluate Tx2 + y 2 UdA, where R is the rectangle [0; 2] × [0; 2]
. R
Now the integrals can be computed by realising that they are areas of
quarter annuli, and so,
ZZ
1 πk
kdA = k · · π(k + 1 − k) = .
4 4
k≤x2 +y2 <k+1,x≥0,y≥0
Hence Z
π 3π
Tx2 + y 2 UdA = (1 + 2 + 3) = .
R 4 2
◭
Homework
Problem
Z 3 Z3.5.1 Evaluate the iterated inte- Problem 3.5.4 Find
x
1 Z
gral dydx.
1 0 x xydxdy
D
S = {(x, y) ∈ R2 : x ≥ 0, y ≥ 0, 1 ≤ x2 +y 2 ≤ 4}.
where D = [0; π]2 .
Z
Find x2 dA. Z 1 Z 1
Problem 3.5.6 Find min(x2 , y 2 )dxdy.
S
0 0
Z
Problem 3.5.7 Find xydxdy where Problem 3.5.12 Find
D Z
D = {(x, y) ∈ R2 : x > 0, y > 0, 9 < x2 +y 2 < 16, 1 < x2 −y 2
loge<
(1 16}.
+ x + y)dA
Z D
Z
b
Z 1 Z 1
2
Problem 3.5.9 Find 2ex dxdy
0 y
Z
Problem 3.5.10 Evaluate min(x, y 2 )dA.
[0;1]2
Z
Problem 3.5.11 Find xydA, where R
R
is the (unoriented) △OAB in figure 3.10 Figure 3.11: Problem 3.5.15.
with O(0, 0), A(3, 1), and B(4, 4).
Z
Problem 3.5.16 Evaluate xdA where R
B R
b
is the E-shaped figure in figure 3.12.
b
A
b
O b b
b b
b b
b b
where
D = {(x, y) ∈ R2 : x4 + y 4 + x2 − y 2 ≤ 1}.
Z 1 Z 1
ex/y
Problem 3.5.24 Evaluate dydx.
0
√
x y
where
Problem 3.5.21 Find
Z D = {(x, y) ∈ R2 : |x| ≤ 1, |y| ≤ 1}.
xydA
Z
D
Problem 3.5.26 Find (2x + 3y + 1) dA,
where D
where D is the triangle with vertices at
D = {(x, y) ∈ R2 : x ≥ 0, y ≥ 0, xy+y+x ≤ A(−1,
1}. −1), B(2, −4), and C(1, 3).
where Z 1 Z 1
2
D = {(x, y) ∈ R |x ≥ 0, y ≥ 0, x + y ≤ 1}. Problem 3.5.37 Find xy dxdy.
Z 0
1
0
x−1
Then demonstrate that dx =
Problem 3.5.29 Let f, g : [0; 1] → [0; 1] be 0 log x
continuous, with f increasing. Prove that log 2.
Z 1 Z 1 Z 1
(f ◦g)(x)dx ≤ f (x)dx+ g(x)dx. Problem 3.5.38 Evaluate
0 0 0
Z 4 Z √ 4−y p
Z
12x − x3 dxdy.
Problem 3.5.30 Compute (xy + y 2 )dA 0 0
S
where Z 2 Z 2 p
2 1/2 1/2 Problem 3.5.39 Evaluate y 1 + x3 dxdy.
S = {(x, y) ∈ R : |x| + |y| ≤ 1}. 0 y
D = {(x, y) ∈ R2 |x ≥ 1, y ≥ 1, x + y ≤ 4}.
Problem 3.5.33 A rectangle R on the plane
is the disjoint union R = ∪N
k=1 R k of rectan-
gles Rk . It is known that at least one side Problem 3.5.42 Prove that
of each of the rectangles Rk is an integer. Z +∞ Z +∞
xe−y sin y 1
Shew that at least one side of R is an inte- dydx = .
ger. 0 2x y2 16
D Z b
2 Z b
Z b
2 2
where (f (x)g(x))dx ≤ (f (x)) dx (g(x)) dx .
a a a
2
D = {(x, y) ∈ R |y ≥ 0, x−y+1 ≥ 0, x+2y−4 ≤ 0}.
272 Theorem Let (D, ∆) ∈ (Rn )2 be open, bounded sets in Rn with volume and
let g : ∆ → D be a continuously differentiable bijective mapping such that
1
det g ′ (u) 6= 0, and both | det g ′ (u)|, are bounded on ∆. For f : D → R
| det g ′ (u)|
bounded and integrable, f ◦ g| det g ′ (u)| is integrable on ∆ and
Z Z Z Z
··· f = ··· (f ◦ g)| det g ′ (u)|,
D ∆
that is
Z Z
··· f (x1 , x2 , . . . , xn )dx1 ∧ dx2 ∧ . . . ∧ dxn
D
Z Z
= ··· f (g(u1 , u2 , . . . , un ))| det g ′ (u)|du1 ∧ du2 ∧ . . . ∧ dun .
∆
One normally chooses changes of variables that map into rectangular regions, or
that simplify the integrand. Let us start with a rather trivial example.
7 7
b
6 6
b
5 5
b b b
4 4
b b b
3 3
2 2
1 1
1 2 3 4 5 6 7 8 9 10 1 2 3 4 5 6 7 8 9 10
Figure 3.15: Example 273. xy-plane. Figure 3.16: Example 273. uv-plane.
u = x + 2y =⇒ du = dx + 2dy,
v = 2x + y =⇒ dv = 2dx + dy,
giving
du ∧ dv = −3dx ∧ dy.
Now, " #" #
1 2 x
(u, v) =
2 1 y
is a linear transformation, and hence it maps quadrilaterals into quadri-
laterals. The corners of the rectangle in the area of integration in the xy-
plane are (0, 3), (1, 3), (1, 4), and (0, 4), (traversed counter-clockwise)
and they map into (6, 3), (7, 5), (9, 6), and (8, 4), respectively, in the
uv-plane (see figure 3.16). The form dx ∧ dy has opposite orientation to
du ∧ dv so we use
dv ∧ du = 3dx ∧ dy
du = 2xdx − 2ydy,
dv = 2ydx + 2xdy,
and so
du ∧ dv = (4x2 + 4y 2 )dx ∧ dy.
We now determine the region ∆ into which the square D = [0; 1]2 is
mapped. We use the fact that boundaries will be mapped into bound-
aries. Put
AB = {(x, 0) : 0 ≤ x ≤ 1},
BC = {(1, y) : 0 ≤ y ≤ 1},
CD = {(1 − x, 1) : 0 ≤ x ≤ 1},
DA = {(0, 1 − y) : 0 ≤ y ≤ 1}.
We deduce that
Z Z
4 4 1
(x − y )dA = (x4 − y 4 ) dudv
[0;1]2 ∆Z 4(x2 + y2 )
1
= (x2 − y 2 )dudv
4 Z∆
1
= ududv
4 Z ∆ Z
1 2 1−v 2 /4
= udu dv
4 0 v 2 /4−1
= 0,
as before. ◭
where
Solution: ◮ We have
dx = 2uvdu + u2 dv,
dy = v 2 du + 2uvdv,
dx ∧ dy = 3u2 v 2 du ∧ dv.
The region transforms into
π
b
2
1 b b
0 b b b b
π
0 1
2
Figure 3.17: Example 276. xy-plane. Figure 3.18: Example 276. uv-plane.
276 Example In this problem we will follow an argument of Calabi, Beukers, and
X 1
+∞
π2
Kock to prove that = .
n=1
n2 6
X
+∞
1 3 X
+∞
1
1. Prove that if S = 2
, then S = .
n=1
n 4 n=1
(2n − 1)2
X
+∞ Z 1 Z 1
1 dxdy
2. Prove that = .
n=1
(2n − 1)2 0 0 1 − x2 y 2
sin u sin v
3. Use the change of variables x = , y = in order to evaluate
Z Z cos v cos u
1 1 dxdy
.
0 0 1 − x2 y 2
Solution: ◮
X
+∞
1 1 X 1
+∞
1
= = S,
n=1
(2n)2 4 n=1
n2 4
a quarter of the sum, hence the sum of the odd terms must be three
3
quarters of the sum, S.
4
2. Observe that
Z 1 2 Z 1 Z 1 Z 1 Z 1
1 2n−2 1 2n−2 2n−2
= x dx =⇒ = x dx y dy = (xy)2n−2 dx
2n − 1 0 2n − 1 0 0 0 0
Thus
X
+∞
1 X Z 1Z 1
+∞ Z 1 Z 1 +∞
X Z 1Z 1
dxdy
2n−2 2n−2
= (xy) dxdy = (xy) dxdy = ,
n=1
(2n − 1)2 n=1 0 0 0 0 n=1 0 0 1 − x2 y 2
as claimed.3
sin u sin v
3. If x = ,y= , then
cos v cos u
dx = (cos u)(sec v)du+(sin u)(sec v)(tan v)dv, dy = (sec u)(tan u)(sin v)du+(sec u)(cos v)d
from where
dx∧dy = du∧dv−(tan2 u)(tan2 v)du∧dv = 1 − (tan2 u)(tan2 v) du∧dv.
Also,
sin2 v sin2 v
1 − x2 y 2 = 1 − · = 1 − (tan2 u)(tan2 v).
cos2 v cos2 u
This gives
dxdy
= dudv.
1 − x2 y 2
We now have to determine the region that the transformation x =
sin u sin v
,y= forms in the uv-plane. Observe that
cos v cos u
Ê s
1 − y2 1 − x2
u = arctan x 2
, v = arctan y .
1−x 1 − y2
This means that the square in the xy-plane in figure 3.17 is trans-
formed into the triangle in the uv-plane in figure 3.18.
We deduce,
Z Z Z Z Z
1 1 dxdy π/2 π/2−v π/2 π v2 π/2 π2 π2 π
= dudv = (π/2 − v) dv = v− = − =
0 0 1 − x2 y 2 0 0 0 2 2 0 4 8
Finally,
3 π2 π2
S= =⇒ S = .
4 8 6
◭
Homework
Problem 3.6.1 Let D ′ = {(u, v) ∈ R2 : u ≤ 1, −u ≤ v ≤ u}. Consider
R2 → R2
Φ:
u+v u−v .
(u, v) 7→ ,
2 2
3
This exchange of integral and sum needs justification. We will accept it for our purposes.
Z
Problem 3.6.4 Find f (x, y)dA where
D
Problem 3.6.5 Use the following steps (due to Tom Apostol) in order to prove that
X
∞
1 π2
2
= .
n 6
n=1
π2 u
➍ Finally, prove that the above integral is by using the substitution θ = arcsin .
6 2
whence
dx ∧ dy = (ρ cos2 θ + ρ sin2 θ)dρ ∧ dθ = ρdρ ∧ dθ.
where
D = {(x, y) ∈ R2 |x ≥ 0, y ≥ 0, y ≤ x, x2 + y 2 ≤ 1}.
Z
278 Example Evaluate xdA, where R is the region bounded by the circles
R
x2 + y 2 = 4 and x2 + y 2 = 2y.
b b
Figure 3.19: Example 277. Figure 3.20: Example 278. Figure 3.21: Example 279.
π
The angle clearly sweeps from 0 to . Thus the integral becomes
2
Z Z π/2 Z 2 sin θ
xdA = r 2 cos θdrdθ
R 0Z 2
1 π/2
= (8 cos θ − 8 cos θ sin3 θ)dθ
3 0
= 2.
◭
Z
2
−xy−y2
279 Example Find e−x dA, where
D
D = {(x, y) ∈ R2 : x2 + xy + y 2 ≤ 1}.
√ 2
2 2 y 2 3y
x + xy + y = x + + .
2 2
√
y 3y
Put U = x + ,V = . The integral becomes
2 2
Z Z
2
−xy−y2 2 2
+V 2 )
e−x dxdy = √ e−(U dU dV.
{x2 +xy+y2 ≤1} 3 {U 2 +V 2 ≤1}
Z ©
1
280 Example Evaluate dA over the region (x, y) ∈ R2 : x2 + y 2 ≤ 4, y ≥ 1
R (x2 + y 2 )3/2
(figure 3.22).
1
Solution: ◮ The radius sweeps from r = to r = 2. The desired
sin θ
integral is
Z Z 5π/6 Z 2
1 1
dA = drdθ
R (x2 + y 2 )3/2 r2
Zπ/6
5π/6
csc θ
1
= sin θ − dθ
π/6 2
√ π
= 3− .
3
◭
Z
281 Example Evaluate (x3 + y 3 )dA where R is the region bounded by the
R
x2 y2
ellipse + = 1 and the first quadrant, a > 0 and b > 0.
a2 b2
Solution: ◮ Put x = ar cos θ, y = br sin θ. Then
x2 y2 a2 r 2 cos2 θ b2 r 2 sin2 θ
+ = 1 =⇒ + = 1 =⇒ r = 1.
a2 b2 a2 b2
Thus the required integral is
Z Z π/2 Z 1
(x3 + y 3 )dA = abr 4 (cos3 θ + sin3 θ)drdθ
R 0 Z 0
1 Z π/2
= ab r 4 dr (a3 cos3 θ + b3 sin3 θ)dθ
0 0
1 2a3 + 2b3
= ab
5 3
2ab(a3 + b3 )
= .
15
Homework
Z
Problem 3.7.1 Evaluate xydA where R is the region
R
R = (x, y) ∈ R2 : x2 + y 2 ≤ 16, x ≥ 1, y ≥ 1 ,
as in the figure 3.23. Set up the integral in both Cartesian and polar coordinates.
Z
Problem 3.7.4 Find f dA where
D
where
R = {(x, y) ∈ R2 : 0 < x < 1, o < y < x2 }.
where
D = {(x, y) ∈ R2 |x ≥ 0, y ≥ 0, x2 + y 2 ≤ 1, x2 + y 2 − 2y ≥ 0}.
Z
Problem 3.7.9 Find f dA where
D
D = {(x, y) ∈ R2 |y ≥ 0, x2 + y 2 − 2x ≤ 0}
and f (x, y) = x2 y.
Z
Problem 3.7.10 Let D = {(x, y) ∈ R2 : x ≥ 1, x2 + y 2 ≤ 4}. Find xdA.
D
Z
Problem 3.7.11 Find f dA where
D
D = {(x, y) ∈ R2 |x ≥ 1, x2 + y 2 − 2x ≤ 0}
1
and f (x, y) = .
(x2 + y 2 )2
D = {(x, y) ∈ R2 : x2 + y 2 − y ≤ 0, x2 + y 2 − x ≤ 0}.
Problem 3.7.15 William Thompson (Lord Kelvin) is credited to have said: “A mathemati-
cian is someone to whom Z +∞ √
2 π
e−x dx =
0 2
is as obvious as twice two is four to you. Liouville was a mathematician.” Prove that
Z +∞ √
2 π
e−x dx =
0 2
➋ Let a > 0 be a real number and put ∆a = {(x, y) ∈ R2 ||x| ≤ a, |y| ≤ a}. Let
Z
2 2
Ja = e−(x +y ) dxdy.
∆a
Prove that
Ia ≤ Ja ≤ Ia √ 2 .
➌ Deduce that Z √
+∞
2 π
e−x dx = .
0 2
1
Problem 3.7.16 Let D = {(x, y) ∈ R2 : 4 ≤ x2 + y 2 ≤ 16} and f (x, y) = 2 .
Z x + xy + y 2
Find f (x, y)dA.
D
Problem 3.7.17 Prove that every closed convex region in the plane of area ≥ π has two
points which are two units apart.
Problem 3.7.18 In the xy-plane, if R is the set of points inside and on a convex polygon,
let D(x, y) be the distance from (x, y) to the nearest point R. Show that
Z +∞ Z +∞
e−D (x,y ) dxdy = 2π + L + A,
−∞ −∞
3.8 Three-Manifolds
282 Definition A 3-dimensional oriented manifold of R3 is simply an open set (body)
V ∈ R3 , where the + orientation is in the direction of the outward pointing nor-
mal to the body, and the − orientation is in the direction of the inward pointing
normal to the body. A general oriented 3-manifold is a union of open sets.
☞ In this section, unless otherwise noticed, we will choose the positive ori-
entation for the regions considered. This corresponds to using the volume form
dx ∧ dy ∧ dz.
= (ex − x − 1)dx
0
5
= e− .
2
◭
Z
284 Example Find z dV if
R
√ √ √
R = {(x, y, z) ∈ R3 |x ≥ 0, y ≥ 0, z ≥ 0, x + y + z ≤ 1}.
Solution: ◮ We have
Z Z c Z b−bz/c Z a−ay/b−az/c
xdxdydz = xdxdydz
0 0 0
V
Z Z
1 c b−bz/c ay az 2
= a− − dydz
2 Z0 0 b c
3
1 c a2 (−z + c) b
= dx
6 0 c3
a2 bc
=
24
◭
Z
286 Example Evaluate the integral xdV where S is the (unoriented) tetrahe-
S
dron with vertices (0, 0, 0), (3, 2, 0), (0, 3, 0), and (0, 0, 2). See figure 3.24.
b
C
b B Three-Manifolds
b A
b
O
b
b
C
x
B y
b
A
Figure 3.25: xy-projection.
Figure 3.24: Problem 286.
Solution: ◮ We must find the projections of the solid on the the coor-
dinate planes.
1. With the order dzdxdy, the limits of integration of z can only de-
pend, if at all, on x and y. Given an arbitrary point in the solid, its
2. With the order dxdydz, the limits of integration of x can only de-
pend, if at all, on y and z. Given an arbitrary point in the solid, x
sweeps √ from the plane to x = 2, so the limits for x are from x = y
to x = 4 − z. The projection of the solid on the yz-plane is the
area bounded by z = 4 − y 2 , and the z and y axes.
Z 4 Z √
4−z Z 2 Z 4 Z √
4−z
1
xyzdxdydz = (4y − y 3 )zdydz
0 0 y 2
Z 0 0
4 z3
= 2z − dz
0 8
= 8.
Homework
Z
Problem 3.8.1 Compute zdV where E is the region in the first octant bounded by the
E
planes y + z = 1 and x + z = 1.
Problem 3.8.2 Consider the solid S in the first octant, bounded by the
Z parabolic cylinder
x2 2
z = 2− and the planes z = 0, y = x, and y = 0. Prove that xyz = first by
2 S 3
integrating in the order dzdydx, and then by integrating in the order dydxdz.
Z Z
Problem 3.8.3 Evaluate the integrals 1dV and xdV , where R is the tetrahedron
R R
with vertices at (0, 0, 0), (1, 1, 1), (1, 0, 0), and (0, 0, 1).
Z
Problem 3.8.4 Compute xdV where E is the region in the first octant bounded by the
E
plane y = 3z and the cylinder x2 + y 2 = 9.
Z
dV
Problem 3.8.5 Find where
(1 + x2 z 2 )(1 + y 2 z 2 )
D
D = {(x, y, z) ∈ R3 : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, z ≥ 0}.
Z
289 Example Find z 2 dxdydz if
R
◭
Z
290 Example Evaluate (x2 + y 2 )dxdydz over the first octant region bounded
D
by the cylinders x2 + y 2
= 1 and x2 + y 2 = 4 and the planes z = 0, z = 1,
x = 0, x = y.
Z 1 Z π/2 Z 2 15π
ρ3 dρdθdz = .
0 π/4 1 16
291 Example Three long cylinders of radius R intersect at right angles. Find the
volume of their intersection.
D ′ = {(x, y, z) ∈ R3 : 0 ≤ y ≤ x, 0 ≤ z, x2 +y 2 ≤ R2 , y 2 +z 2 ≤ R2 , z 2 +x2 ≤ R2 }.
Now,
Z Z Z √R2 −ρ2 cos2 θ ! !
π/4 R
V′ = dz ρdρ dθ
0 0 0
Z π/4 Z R È
= ρ R2 − ρ2 cos2 θdρ dθ
Z0π/4 0
R
1
= − 2
(R2 − ρ2 cos2 θ)3/2 dθ
0 Z π/43 cos θ 0
R3 1 − sin3 θ
= dθ
3 0 cos2 θ √ !
Z 2
= R3 π/4 2 1 − u2
u = cos θ [tan θ]0 + du
3 1 u2
22
√
R3
= 1 − u−1 + u
3
√
1
2−1 3
= √ R .
2
Finally
√
V = 16V ′ = 8(2 − 2)R3 .
◭
We have
This gives
dx ∧ dy ∧ dz = −ρ2 sin φdρ ∧ dθ ∧ dφ.
From this derivation, the form dρ ∧ dθ ∧ dφ is negatively oriented, and so we
choose
dx ∧ dy ∧ dz = ρ2 sin φdρ ∧ dφ ∧ dθ
instead.
Z
292 Example Let (a, b, c) ∈]0; +∞[3 be fixed. Find xyz dV if
R
3 x2 y2 z2
R= (x, y, z) ∈ R : + + ≤ 1, x ≥ 0, y ≥ 0, z ≥ 0 .
a2 b2 c2
We have
dx ∧ dy ∧ dz = abcρ2 sin φdρ ∧ dφ ∧ dρ.
The integration region is mapped into
π π
∆ = [0; 1] × [0; ] × [0; ].
2 2
The integral becomes
Z Z Z
2
π/2 1
5
π/2
3 (abc)2
(abc) cos θ sin θ dθ ρ dρ cos φ sin φ dφ = .
0 0 0 48
Homework
p
Problem 3.9.1 Consider the region R below the cone z = x2 + y 2 and above the
2 2
paraboloid z = x + y for 0 ≤ z ≤ 1. Set up integrals for the volume of this region
in Cartesian, cylindrical and spherical coordinates. Also, find this volume.
Z
Problem 3.9.2 Consider the integral xdV , where R is the region above the paraboloid
R
z = x2 + y 2 and under the sphere x2 + y 2 + z 2 = 4. Set up integrals for the volume of
this region in Cartesian, cylindrical and spherical coordinates. Also, find this volume.
Problem 3.9.3 Consider the region R bounded by the sphere x2 + y 2 + z 2 = 4 and the
plane z = 1. Set up integrals for the volume of this region in Cartesian, cylindrical and
spherical coordinates. Also, find this volume.
x2 y2 z2
2
+ 2 + 2 =1
a b c
4πabc
is . Here a > 0, b > 0, c > 0.
3
Z
Problem 3.9.5 Compute ydV where E is the region between the cylinders x2 +y 2 = 1
E
and x2 + y 2 = 4, below the plane x − z = −2 and above the xy-plane.
Z
Problem 3.9.7 Compute y 2 z 2 dV where E is bounded by the paraboloid x = 1 − y 2 − z 2
E
and the plane x = 0.
Z
p
Problem 3.9.8 Compute z x2 + y 2 + z 2 dV where E is is the upper solid hemisphere
E
bounded by the xy-plane and the sphere of radius 1 about the origin.
where
R = {(x, y, z) ∈ R3 : x ≥ 0, y ≥ 0, z ≥ 0, x + y + z ≤ 1}.
☞ In this section, unless otherwise noticed, we will choose the positive ori-
entation for the regions considered. This corresponds to using the ordered basis
{dy ∧ dz, dz ∧ dx, dx ∧ dy}.
Here 2 È
d x = (dx ∧ dy)2 + (dz ∧ dx)2 + (dy ∧ dz)2
Now, Z Z
p
z d2 x = (u2 + v 2 ) 1 + 4u2 + 4v 2 dudv.
Σ D
297 Example Find the area of that part of the cylinder x2 + y 2 = 2y lying inside
the sphere x2 + y 2 + z 2 = 4.
Solution: ◮ We have
x2 + y 2 = 2y ⇐⇒ x2 + (y − 1)2 = 1.
whence
and so
2 È
d x = (dx ∧ dy)2 + (dz ∧ dx)2 + (dy ∧ dz)2
È
= cos2 u + sin2 u du ∧ dv
= du ∧ dv.
where Σ is the top side of the triangle with vertices at (2, 0, 0), (0, 2, 0), (0, 0, 4).
Solution: ◮ Observe that the plane passing through the three given
points has equation 2x + 2y + z = 4. We project this plane onto the
coordinate axes obtaining
Z Z 4 Z 2−z/2 8
xdydz = (2 − y − z/2)dydz = ,
0 0 3
Σ
Z Z 2 Z 4−2x
(z 2 − zx)dzdx = (z 2 − zx)dzdx = 8,
0 0
Σ
Z Z 2 Z 2−y 2
− xydxdy = − xydxdy = − ,
0 0 3
Σ
and hence
Z
xdydz + (z 2 − zx)dzdx − xydxdy = 10.
Σ
Homework
Z
Problem 3.10.1 Evaluate y d2 x where Σ is the surface z = x + y 2 , 0 ≤ x ≤ 1, 0 ≤
Σ
y ≤ 2.
p
Problem 3.10.2 Consider the cone z = x2 + y 2 . Find the surface area of the part of
the cone which lies between the planes z = 1 and z = 2.
Z
Problem 3.10.3 Evaluate x2 d2 x where Σ is the surface of the unit sphere x2 + y 2 +
Σ
z 2 = 1.
Z
p
Problem 3.10.4 Evaluate z d2 x over the conical surface z = x2 + y 2 between
S
z = 0 and z = 1.
Problem 3.10.5 You put a perfectly spherical egg through an egg slicer, resulting in n
slices of identical height, but you forgot to peel it first! Shew that the amount of egg shell
in any of the slices is the same. Your argument must use surface integrals.
where Σ is the top of the triangular region of the plane 2x + 2y + z = 6 bounded by the
first octant.
Solution: ◮ We will first use Green’s Theorem and then evaluate the
integral directly. We have
dω = d(x − y 3 ) ∧ dx + d(x3 ) ∧ dy
= (dx − 3y 2 dy) ∧ dx + (3x2 dx) ∧ dy
= (3y 2 + 3x2 )dx ∧ dy.
Z
301 Example Evaluate (x − y)dydz + zdzdx − ydxdy where S is the surface
S
of the sphere
x2 + y 2 + z 2 = 9
dω = (dx − dy) ∧ dy ∧ dz + dz ∧ dz ∧ dx − dy ∧ dx ∧ dy
= dx ∧ dy ∧ dz.
Also Z
zdzdx = 0,
Σ
In general, let
be a 2-form in R3 . Then
dω = df (x, y, z)dy ∧ dz + dg(x, y, z)dz ∧ dx + dh(x, y, z)dx ∧ dy
∂ ∂ ∂
= f (x, y, z)dx + f (x, y, z)dy + f (x, y, z)dz ∧ dy ∧ dz
∂x ∂y ∂z
∂ ∂ ∂
+ g(x, y, z)dx + g(x, y, z)dy + g(x, y, z)dz ∧ dz ∧ dx
∂x ∂y ∂z
∂ ∂ ∂
+ h(x, y, z)dx + h(x, y, z)dy + h(x, y, z)dz ∧ dx ∧ dy
∂x ∂y ∂z
∂ ∂ ∂
= f (x, y, z) + g(x, y, z) + h(x, y, z) dx ∧ dy ∧ dz,
∂x ∂y ∂z
which gives the classical Gauss’s Theorem
Z Z
∂ ∂ ∂
f (x, y, z)dydz+g(x, y, z)dzdx+h(x, y, z)dxdy = f (x, y, z) + g(x, y, z) + h(x, y, z) dx
∂x ∂y ∂z
∂M M
Solution: ◮ We have
dω = (dy) ∧ dx + (2dx − dz) ∧ dy + (dz − dx) ∧ dz
= −dx ∧ dy + 2dx ∧ dy + dy ∧ dz + dz ∧ dx
= dx ∧ dy + dy ∧ dz + dz ∧ dx.
Since on C, z = 1, the surface Σ on which we are integrating is the
inside of the circle x2 + y 2 + 1 = 4, i.e., x2 + y 2 = 3. Also, z = 1 implies
dz = 0 and so Z Z
dω = dxdy.
Σ Σ
Since this is just the area of the circular region x2 + y 2 ≤ 3, the integral
evaluates to Z
dxdy = 3π.
Σ
In general, let
be a 1-form in R3 . Then
then Z Z
→
−
(∇ × →
−
a)•dS = →
−
a • d→
−
r.
M ∂M
Homework
I
Problem 3.11.1 Evaluate x3 ydx + xydy where C is the square with vertices at (0, 0),
C
(2, 0), (2, 2) and (0, 2).
Problem 3.11.2 Consider the triangle △ with vertices A : (0, 0), B : (1, 1), C : (−2, 2).
➊ If LP Q denotes the equation of the line joining P and Q find LAB , LAC , and LB C .
➋ Evaluate I
y 2 dx + xdy.
△
➌ Find Z
(1 − 2y)dx ∧ dy
D
M = {(x, y, z) ∈ R3 : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, 0 ≤ z ≤ 2},
U = {(x, y, z) ∈ R3 : 0 ≤ x ≤ 1, 0 ≤ y ≤ 1, z = 2},
the boundary of the box without the upper top S = ∂M \ U , and the differential form
Z
4. Prove that the integral on the open box is (arctan y − x2 )dydz +(cos x sin z −
∂ M \U
y 3 )dzdx + (2zx + 6zy 2 )dxdy = −4.
1. Prove that the equation of the plane that contains the triangle T is 2x + y + 4z = 12.
2. Prove that dω = ydy ∧ dz + (2x − y) dz ∧ dx + zdx ∧ dy.
Z
y2
3. Prove that (2xz + arctan ex ) dx+(xz + (y + 1)y ) dy+ xy + + log(1 + z 2 ) dz =
2
Z 3 Z 12−4z ∂ T Z 6 Z 3−x/2
ydydz + 2xdzdx=108.
0 0 0 0
2 2
where Γ is the circle (x − 2) + y = 4. Also, prove this directly by using a path integral.
Problem 3.11.7 Let Γ denote the curve of intersection of the plane x + y = 2 and the
sphere x2 − 2x + y 2 − 2y + z 2 = 0, oriented clockwise when viewed from the origin. Use
Stoke’s Theorem to prove that
Z
√
ydx + zdy + xdz = −2π 2.
Γ
Prove this directly by parametrising the boundary of the surface and evaluating the path
integral.
where C is the positively oriented boundary of the region between the circles x2 + y 2 = 2
and x2 + y 2 = 4.
→
−
1.1.1 No. The zero vector 0 , has magnitude but no direction.
−→ −→ −→ −→ −→ −→ −→ −→ −→
1.1.2 We have 2BC = BE + EC. By Chasles’ Rule AC = AE + EC, and BD = BE + ED.
We deduce that
−→ −→ −→ −→ −→ −→ −→ −→
AC + BD = AE + EC + BE + ED = AD + BC.
−→ −→
But since ABCD is a parallelogram, AD = BC. Hence
−→ −→ −→ −→ −→
AC + BD = AD + BC = 2BC.
−
→ −
→ −
→ −
→ −→ −
→ −→
1.1.4 We have IA = −3IB ⇐⇒ IA = −3(IA + AB) = −3IA − 3AB. Thus we deduce
−
→ −
→ −→ −
→ −→
IA + 3IA = −3AB ⇐⇒ 4IA = −3AB
−
→ −→
⇐⇒ 4AI = 3AB
−
→ 3 −→
⇐⇒ AI = AB.
4
Similarly
−→ 1−→ −→ −
→
JA = − JB ⇐⇒ 3JA = −JB
3 −→ −→ −→
⇐⇒ 3JA = −JA − AB
−→ −→
⇐⇒ 4JA = −AB
−→ 1 −→
⇐⇒ AJ = AB
4
.
−
→ 3 −→ −
→ 1 −→
Thus we take I such that AI = AB and J such that AJ = AB.
4 4
Now
−−→ −−
→ −→ − → −
→
MA + 3MB = MI + IA + 3IB
−→ − → −
→
= 4MI + IA + 3IB
−→
= 4MI,
and
−−→ −−→ −→ −→ −→ − →
3MA + MB = 3MJ + 3JA + MJ + JB
−→ −→ − →
= 4MJ + 3JA + JB
−→
= 4MJ.
1.1.5
x + 1 = t = 2 − y =⇒ y = −x + 1.
4 3 4 2 1
1.1.6 α = ,β = , l = ,m = ,n = .
7 7 9 9 3
→
− →
− →
− →
− →
−
1.1.8 [A]. 0 , [B]. 0 , [C]. 0 , [D]. 0 , [E]. 2→
−
c (= 2 d )
186
Appendix A
1.3.2 Plainly,
a b − a −1 a+b 1
= + .
b 2 1 2 1
1
1.4.1 a = −3, b = − .
2
1.4.2 The desired transformations are in figures A.1 through A.3.
5 5
5
4 4
4
3 3
3
2 2
2
1 1
1
0 0
0
-1 -1
-1
-2 -2
-2
-3 -3
-3
-4 -4
-4
-4 -3 -2 -1 0 1 2 3 4 5 -4 -3 -2 -1 0 1 2 3 4 5
-4 -3 -2 -1 0 1 2 3 4 5
Figure A.1: Horizontal Figure A.2: Vertical Stretch. Figure A.3: Horizontal and
Stretch. Vertical Stretch.
1.4.3 The desired transformations are shewn in figures through A.4 A.7.
a b
1.4.9 , bc = −a2
c −a
3 3 3
2 2 2
Answers and Hints
1 1 1
0 0 0
-1 -1 -1
-2 -2 -2
-3 -3 -3
5
-4 -4 -4
4 -4 -3 -2 -1 0 1 2 3 4 5 -4 -3 -2 -1 0 1 2 3 4 5 -4 -3 -2 -1 0 1 2 3 4 5
3
2 Figure A.4: Levogyrate rota- Figure A.5: Levogyrate rota- Figure A.6: Dextrogyrate ro-
π π π
1 tion radians. tion radians. tation radians.
2 4 2
0
-1
-2
-3
-4
-4 5 -3 -2 -1 0 1 2 3 4 5 5 5
4 4 4
Figure
3 A.7: Dextrogyrate ro- 3 3
π
tation
2 radians. 2 2
4
1 1 1
0 0 0
-1 -1 -1
-2 -2 -2
-3 -3 -3
-4 -4 -4
-4 -3 -2 -1 0 1 2 3 4 5 -4 -3 -2 -1 0 1 2 3 4 5 -4 -3 -2 -1 0 1 2 3 4 5
Figure A.8: Reflexion about Figure A.9: Reflexion about Figure A.10: Reflexion about
the x-axis. the y-axis . the origin.
→
−
1.5.2 Since →
−
a • b = 0, we have
→
− →
− →
−
||→
−
a + b ||2 = (→
−a + b )•(→ −a + b)
→
− →
− →
− →
−
= a •→
−a + 2→ −a•b + b•b
→
− →
− →
−
= a •→
−a +0+ b•b
→
− →
−
= || a ||2 + || b ||2 ,
1.5.5 We have
||→
−
u +→
−v ||2 − ||→
−
u −→
−
v ||2 = (→
−u +→ −v )•(→
−
u +→−v ) − (→
−
u −→ −
v )•(→
−
u −→ −v)
= →
−u • u + 2 u • v + v • v − ( u • u − 2 u •→
→
− →
− →
− →
− →
− →
− →
− →
− −
v +→
−
v •→
−
v)
= →
−
4u•v, →
−
y−x
1.6.1 We have y − x = 4t =⇒ t = and so
4
y−x 3 y−x
x= −2
4 4
is the Cartesian equation sought.
1.6.3
1. ay − cx = ad − bc, this is a straight line with positive slope.
2. −1 ≤ x ≤ 1, y = 0, this is the line segment on the plane joining (−1, 0) to (1, 0).
x2 y2
3. 2
− 2 = 1, x > 0. This is one branch of a hyperbola.
a b
x2 y2
4. 2 − 2 = 1. This is a hyperbola.
a b
π
1.6.4 We may simply give the trivial parametrisation: x = t, y = log cos t, 0 ≤ t ≤ .
3
Then
(dx)2 + (dy)2 = (1 + tan2 t)(dt)2 = sec2 t(dt)2 .
Hence the arc length is
Z π /3 √
sec tdt = log(2 + 3).
0
1.6.5 Observe that y = 2x + 1, so the trace is part of this line. Since in the interval
[0; 4π], −1 ≤ sin t ≤ 1, we want the portion of the line y = 2x + 1 with −1 ≤ x ≤ 1
(and, thus −1 ≤ y ≤ 3). The curve starts at the middle point (0, 1) (at t = 0), reaches
π 3π
the high point (1, 3) at t = , reaches its low point (−1, 1) at t = , reaches its high
2 2
5π 7π
point (1, 3) again at t = , it goes to its low point (−1, 1) at t = , and finishes in the
2 2
middle point (0, 1) when t = 4π.
√ √
1.6.6 First observe that x + y = 1 demands x ∈ [0; 1] and [0; 1]. Again, one can give
many parametrisations. One is x = t2 , y = (1 − t)2 , t ∈ [0; 1]. This gives
È È É
p 2 1 2
2 2 2 2 2
(dx) + (dy) = 4t + (2 − 2t) dt = 2 2t − 2t + 1 dt = √ 4 t− + 1 dt.
2 2
To integrate
Z 1É
2 1 2
√ 4 t− + 1 dt,
2 0 2
we now use the trigonometric substitution
1 1
2 t− = tan θ =⇒ dt = sec2 θdθ.
2 2
1.6.10 Observe that the parametrisation traverses the curve once clockwise if t ∈ [0; 2π].
The area is given by
I I
1 x dx 1
det = xdy − ydx
2 Γ y dy 2
Z 0
4
= (sin3 t(− sin t(1 + sin2 t) + 2 sin t cos2 t)
2 π /2
− cos t(1 + sin2 t)(3 sin2 t cos t))dt
Z 0
= 2 (−3 sin2 t + sin4 t)dt
Z 0
π /2
9 1
= 2 − + cos 2t + cos 4t dt
π /2 8 8
9π
= .
8
☞ A shorter way of obtaining xdy − ydx would have been to argue that
y 9t2
xdy − ydx = x2 d = dt.
x (1 + t3 )2
1.6.12 See figure 1.58. Let θ be the angle (in radians) of rotation of the circle, and let C
be the centre of the circle. At θ = 0 the centre of the circle is at (0, ρ), and P = (0, ρ − d).
Suppose the circle is displaced towards the right, making the point P to rotate an angle
of θ radians. Then the centre of the circle has displaced rθ units horizontally, and so
is now located at (ρθ, ρ). The polar coordinates of the point P are (d sin θ; d cos θ), in
relation to the centre of the circle (notice that the circle moves clockwise). The point P
has moved x = ρθ − d sin θ horizontal units and y = ρ − d cos θ units. This is the desired
parametrisation.
1.6.14 We have
È È
dx = et (cos t−sin t)dt, dy = et (sin t+cos t)dt =⇒ (dx)2 + (dy)2 = et (cos t − sin t)2 + (sin t + cos t)2 dt =
1.6.15 Choose coordinates so that the origin is at the position of the gun, the y-axis is
vertical, and the airplane is on a point with coordinates (u, h) with u ≥ 0.
gt2
x = V t cos a; y = V t sin a − ,
2
where a is the angle of elevation, t is the time and g is the acceleration due to gravity.
Since we know that the shell strikes the plane, we must have
gt2
u = V t cos a; h = V t sin a − ,
2
whence 2
1
u2 + h + gt2 = V 2 t2 ,
2
and thus
g 2 t4
+ (gh − V 2 )2 t2 + h2 + u2 = 0.
4
The quadratic equation in t2 has a real root if
1.6.16 Suppose the parabolas have a point of contact P = (4px20 , 4px0 ). By symmetry,
the vertex V of the rolling parabola is the reflexion of the origin about the line tangent to
1
their point of contact. The slope of the tangent at P is , from where the equation of
2x0
the tangent is
x
y= + 2px0 .
2x0
The line normal to this line and passing through the origin is hence
y = −2xx0 ,
or
(x2 + y 2 )x + 2py 2 = 0,
giving the equation of the locus.
whence
É 2 È
→
− −
→ →
− 2 →
− →
− −
→
a − b = 2 a + 2 b − a + b = 2(13)2 + 2(19)2 − (24)2 = 22.
Ǒ
Ǒ 2 3
x 1 −2
6 7
1.7.2 y = 2 + t 4−15 .
z 3 0
(x − 1) − 2(y − 1) − (z − 1) = 0.
Ǒ
0
1.7.4 Observe that 0 (as 0 = 2(0) = 3(0)) is on the line, and hence on the plane.
0
Thus the vector 2 3 2 3
1−0 1
6 7 6 7
4−1 − 05 = 4−15
−1 − 0 −1
lies on the plane. Now, if x = 2y = 3z = t, then x = t, y = t/2, z = t/3. Hence, the
vectorial form of the equation of the line is
2 3 2 3 2 3
0 1 1
→
− 6 7 6 7 6 7
r = 405 + t 41/25 = t 41/25 .
0 1/3 1/3
2 3
1
6 7
This means that 41/25 also lies on the plane, and thus
1/3
2 3 2 3 2 3
1 1 1/6
6 7 6 7 6 7
4−15 × 41/25 = 4−4/35
−1 1/3 3/2
which is clearly impossible, and so the lines are skew. Let θ be the angle between them.
Then
2 · 2 + 1 · (−1) + 1 · 1 4 2
cos θ = p p = √ √ =⇒ θ = arccos .
(2)2 + (1)2 + (1)2 (2)2 + (−1)2 + (1)2 6 6 3
1
1.7.15 x + y + z = 1. .
6
→
− −
1.7.16 Assume contrariwise that →−
a, b, → c are three unit vectors in R3 such that the
2π →
− 1 →− − 1
angle between any two of them is > . Then →
−a • b < − , b •→ c < − , and →
−c •→
−
a <
3 2 2
1
− . Thus
2 2
→ →
− − − 2 →
− 2 − 2
− →
a + b + c = → a + b + →
c
→
− →
− −
+2→−
a • b + 2 b •→c + 2→ −
c •→
−
a
< 1+1+1−1−1−1
= 0,
which is impossible, since a norm of vectors is always ≥ 0.
→
− •→
t − s →− →
− →
−
1.7.17 Let projts = s be the projection of t over → −s 6= 0 . Let x0 be the point
(||s||)2
−→ →
−
on the plane that is nearest to b. Then bx0 = → −
x0 − b is orthogonal to the plane, and the
distance we seek is
→ →
− − →
− −
(−
r0 − b )•→ |(→
−
r0 − b )•→
→
− →
−
r0 − b n→−
n|
||projn || = → n = .
−
n
2
→−n
as we wanted to shew.
1.7.18 There are 7 vertices (V0 = (0, 0, 0), V1 = (11, 0, 0), V2 = (0, 9, 0), V3 = (0, 0, 8),
V4 = (0, 3, 8), V5 = (9, 0, 2), V6 = (4, 7, 0)) and 11 edges (V0 V1 , V0 V2 , V0 V3 , V1 V5 , V1 V6 ,
V2 V4 , V2 V6 , V3 V4 , V3 V5 , V4 V5 , and V4 V6 ).
x y
Figure A.11: Problem 1.7.18.
→
−
1.8.2 We have →
−
x ×→−
x = −→
−
x ×→
−
x by letting →
−
y = →
−
x in ??. Thus 2→
−
x ×→
−
x = 0 and
→
−
hence →
−
x ×→ −
x = 0.
as obtained before.
that is,
2ax + 3a2 y − az = a2 .
→
−v ×→−
w →
−v ×→−
w
1.8.6 Either of →
− →
− or − →− →
− will do. Now
|| v × w || || v × w ||
→
− →
− →
− →
− →
−
v ×→
−
w = (−a j + a k ) × ( i + a j )
→
− →
− →
− →
− →
− →
− →
− →
−
= −a( j × i ) − a2 ( j × j ) + a( k × i ) + a2 ( k × j )
→
− →
− →
−
= a k + a j − a2 i
Ǒ
−a2
= a ,
a
p p
and ||→
−
v ×→
−
w || = a4 + a2 + a2 = 2a2 + a4 . Hence we may take either
Ǒ
−a2
1
√ 2 a
2a + a4
a
or
Ǒ
−a2
1
−√ a .
2a2 + a4
a
→
− − 2 →
− 2 →
− →
−
||→
− x i − (→
x × i ||2 = → −
x • i )2 = 1 − (→
−
x • i )2 ,
→
− − 2 →
− 2 →
− →
−
x × k ||2 = →
||→
− x j − (→ −
x • j )2 = 1 − (→
−
x • j )2 ,
→
− − 2 →
− 2 →
− →
−
x × j ||2 = →
||→
− x k − (→ −
x • k )2 = 1 − (→
−
x • k )2 ,
→
− →
− →
− − 2
and since (→
−
x • i )2 +(→
−
x • j )2 +(→x • k )2 = →
− x = 1, the desired sum equals 3− 1 = 2.
1.8.9
→
− →
− →
− →
− − →
− →
− − → →
− − → →
− −
a ×(→
−
x × b ) = b ×(→
−
x ×→
−
a ) ⇐⇒ (→
−
a • b )→
x −(→
−
a •→
−
x ) b = ( b •→
a )−
x −( b •→
x )−
a ⇐⇒ →
−
a •→
−
x = b •→
x = 0.
→
−
The answer is thus {→
−
x :→
−
x ∈ R→
−
a × b }.
1.8.12
→
− − →
−
→
− (→
−
a • b )→
a + 6 b + 2→ −
a ×→−
c
x = →
12 + 2 a − 2
→
−
→
− (→
−
a •→
−c )→
−
a + 6→−c + 3→ −
a × b
y = →
18 + 3− a
2
This is so because both sides give the volume of the parallelogram spanned by →
−
x,→
−
y,→
−
z.
→
− →
− →
− →
− →
− →
− →
−
Now, putting x = a × b , y = c and z = d we gather that
→
− →
− →
− →
−
(→
−
a × b ) • (→
−
c × d) = ((→
−
a × b)×→
−
c)• d.
Now, again,
→
− →
− →
− − →
− →
− →
− −
(→
−
a × b )× →
−
c = −→
−
c × (→
−
a × b ) = −((→
−
c • b )→
a − (→
−
c •→
−
a ) b ) = (→
−
c •→
−
a ) b − (→
−
c • b )→
a.
This gives
→
− →
− →
− →
− − →
− →
− → − →
− − •→ −
((→
−
a × b )× →
−
c ) • d = ((→
−
c •→
−
a ) b − (→
−
c • b )→
a ) • d = (→
−
c •→
−
a )( b • d ) − (→
−
c • b )(→
a d ),
(→
−
x ×→
−
y ) • (→
−
u ×→
−
v ) = (→
−
x •→
−
u ) • (→
−
y •→
−
v ) − (→
−
x •→
−
v ) • (→
−
y •→
−
u ).
Adding these three equalities, and using the fact that the dot product is commutative, we
see that all the terms on the dextral side cancel out and we obtain 0, as required.
1.8.15
1. We have
2 3 2 3 2 3
6 0 12
−→ 6 7 −→ 6 7 −→ −→ →
− →
− →
− →
− →
− →
− →
− 6 7
CA = 4 0 5 , CB = 4 4 5 =⇒ CA×CB = (6 i −3 k )×(4 j −3 k ) = 24 k +18 j +12 i = 4185 .
−3 −3 24
2. We have p √ √
−→ −→
CA × CB = 122 + 182 + 242 = 1044 = 6 29.
5. From the preceding items, the line LC A is x = 6t, y = 0, z = 3 − 3t and the line
LD E is x = 3, y = 0, z = −3s. If the line intersect then 6t = 3,
0 = 0, 3 − 3t =
1 1 3
−3s gives t = and s = − . The point of intersection is thus 3, 0, . item The
2 2 2
area is √ √
1 −→ −→ 1
CA × CB = · 6 29 = 3 29.
2 2
6. Observe that
Ǒ
Ǒ
Ǒ
Ǒ
3 3 x 0
P = 0 , Q= y , R= 3 , S = 3 .
z 0 0 z
Now
An = AAn−1 = A(3n−2 A) = 3n−2 A2 = 3n−2 3A = 3n−1 A,
and so the assertion is proved by induction.
1.9.3 First, we will prove that
2 3
1 2 3 4 ... n−1 n
60 n−2 n − 17
6 1 2 3 ... 7
6 7
A2 = 6
6
0 0 1 2 ... n−3 n − 27 .
7
6 .. .. .. 7
4. . . ... . . . ... ... 5
0 0 0 0 ... 0 1
Observe that A = [aij ], where aij = 1 for i ≤ j and aij = 0 for i > j.
S1 = a, S2 = k.
k=i k=i
The first sum is trivial: there are j − i + 1 integers in the interval [i; j], and hence
X
j
X
j
S1 = a = S1 = a 1 = a(j − i + 1).
k=i k=i
For the second sum, we use Gauß trick: summing the sum forwards is the same as
summing the sum backwards, and so, adding the first two rows below,
(i + j)(j − i + 1)
which gives S2 = .
2
= (k − i + 1)
k=i
X
j
X
j
X
j
= k− i+ 1
k=i k=i k=i
(j + i)(j − i + 1)
= − i(j − i + 1) + (j − i + 1)
2
(j − i + 1)(j − i + 2)
= .
2
Assume now that i > j. Then
X
n
X
n
2 3
1 0 a+b
6 (a + b)2 7
= 6−(a + b) 1 − 7
4 2 5
0 0 1
= m(a + b)
For the second part, observe that using the preceding part of the problem,
2 3
1 0 0
6 02 7
m(a)m(−a) = m(a − a) = m(0) = 6 7
4−0 1 − 2 5 = I3 ,
0 0 1
a
1.11.3 √
3
1.12.1 Consider a right triangle △ABC rectangle at A with legs of length CA = h and
AB = r, as in figure A.12. The cone is generated when the triangle rotates about CA.
The gyrating curve is thep hypotenuse, whose centroid is its centre. The length of the
generating curve
r is thus r 2 + h2 and the length of curve described by the centre of
p
gravity is 2π = πr. The lateral area is thus πr r 2 + h2 .
2
b
C
b
r/2 b
′
b b
G
h/3 G B
b b
rh
To find the volume, we gyrate the whole right triangle, whose area is . We need
2
to find the centroid of the triangle. But from example 17, we know that the centroid G
of the triangle is is two thirds of the way from A to the midpoint of BC. If G′ is the
perpendicular projection of G onto [CA], then this means that G′ is at a vertical height
h 2 h GG′ h/3 r
of · = . By similar triangles = =⇒ GG′ = . Hence, the length of the
2 3 3 r h 3
2
curve described by the centre of gravity of the triangle is πr. The volume of the cone is
3
2 rh π
thus πr · = r 2 h.
3 2 3
−−−→ −−−→
1.14.1 Find a vector → −a mutually perpendicular to V1 V2 and V1 V3 and another vector
→
− −−−→ −−−→ 1
and a vector b mutually perpendicular to V1 V3 and V1 V4 . Then shew that cos θ = ,
3
→
−
where θ is the angle between → −a and b .
Ǒ
x
1.15.1 Let y be a point on S. If this point were on the xz plane, it would be on the
z
1p
ellipse, and its distance to the axis of rotation would be |x| = 1 − z 2 . Anywhere else,
2
Ǒ
Ǒ
x 0
the distance from y to the z-axis is the distance of this point to the point 0 :
z z
p
x2 + y 2 . This distance is the same as the length of the segment on the xz-plane going
from the z-axis. We thus have
p 1p
x2 + y 2 = 1 − z 2,
2
or
4x2 + 4y 2 + z 2 = 1.
Ǒ
x
1.15.2 Let y be a point on S. If this point were on the xy plane, it would be on the
z
1
line, and its distance to the axis of rotation would be |x| = |1 − 4y|. Anywhere else, the
Ǒ 3
Ǒ
Ǒ
x x 0
distance of y to the axis of rotation is the same as the distance of y to y ,
z z 0
p
that is x2 + z2. We must have
p 1
x2 + z 2 = |1 − 4y|,
3
which is to say
9x2 + 9z 2 − 16y 2 + 8y − 1 = 0.
1.15.6 The planes A : x + z = 0 and B : y = 0 are secant. The surface has equation of
2 2
the form f (A, B) = eA +B − A = 0, and it is thus a cylinder. The directrix has direction
→
− →
−
i − k.
1.15.7 Rearranging,
1
(x2 + y 2 + z 2 )2 − ((x + y + z)2 − (x2 + y 2 + z 2 )) − 1 = 0,
2
and so we may take A : x + y + z = 0, Σ : x2 + y 2 + z 2 = 0, shewing that the surface is
→
− →
− →
−
of revolution. Its axis is the line in the direction i + j + k .
1.15.9 Rearranging,
(x + y + z)2 − (x2 + y 2 + z 2 ) + 2(x + y + z) + 2 = 0,
so we may take A : x + y + z = 0, Σ : x2 + y 2 + z 2 = 0 as our plane and sphere. The
→
− →
− →
−
axis of revolution is then in the direction of i + j + k .
1.15.11 The largest circle has radius b. Parallel cross sections of the ellipsoid are sim-
ilar ellipses, hence we may increase the size of these by moving towards the centre of
the ellipse. Every plane through the origin which makes a circular cross section must
intersect the yz-plane, and the diameter of any such cross section must be a diameter of
y2 z2
the ellipse x = 0, 2 + 2 = 1. Therefore, the radius of the circle is at most b. Arguing
b c
similarly on the xy-plane shews that the radius of the circle is at least b. To shew that
circular cross section of radius b actually exist, one may verify that the two planes given
by a2 (b2 − c2 )z 2 = c2 (a2 − b2 )x2 give circular cross sections of radius b.
1.15.12 Any hyperboloid oriented like the one on the figure has an equation of the form
z2 x2 y2
= + − 1.
c2 a2 b2
When z = 0 we must have
1
4x2 + y 2 = 1 =⇒ a = , b = 1.
2
Thus
z2
= 4x2 + y 2 − 1.
c2
Hence, letting z = ±2,
4 1 y2 1 1 3
2
= 4x2 + y 2 − 1 =⇒ 2 = x2 + − = 1− = ,
c c 4 4 4 4
y2
since at z = ±2, x2 + = 1. The equation is thus
4
3z 2
= 4x2 + y 2 − 1.
4
1.16.1 The arc length element is
È p
(dx)2 + (dy)2 + (dz)2 = 4t2 + 36t + 81 dt = (2t + 9) dt.
We need t = 1 to t = 4. The desired length is
Z 4 4
(2t + 9) dt = (t2 + 9t) = (16 + 36) − (1 + 9) = 42.
1 1
1.16.2 Observe that z = 1 + x2 + y 2 is a paraboloid opening up, with vertex (lowest point)
at (0, 0, 1). On the other hand, z = 3 − x2 − y 2 is another paraboloid opening down, with
highest point at (0, 0, 3). Adding the equations we obtain
2z = z + z = (1 + x2 + y 2 ) + (3 − x2 − y 2 ) = 4 =⇒ 2z = 4 =⇒ z = 2,
(1 + x2 + y 2 ) − (3 − x2 − y 2 ) = z − z = 0 =⇒ 2x2 + 2y 2 − 2 = 0 =⇒ x2 + y 2 = 1,
1.16.3 Let →
−
r (t) lie on the plane ax + by + cz = d. Then we must have
t4 t3 t2
a +b +c = d =⇒ (at4 +bt3 +ct2 ) = d(1+t2 ) =⇒ at4 +bt3 +(c−d)t2 −d = 0,
1 + t2 1 + t2 1 + t2
which means that if → −r (t) is on the plane ax + by + cz = d, then t must satisfy the quartic
polynomial p(t) = at4 + bt3 + (c − d)t2 − d = 0. Hence, the tk are coplanar if and only
if they are roots of p(t). Since the coefficient of t in this polynomial is 0, then the sum of
the roots of p(t) taken three at a time is 0, that is,
t1 t2 t3 + t1 t2 t4 + t1 t3 t4 + t2 t3 t4 1 1 1 1
t1 t2 t3 +t1 t2 t4 +t1 t3 t4 +t2 t3 t4 = 0 =⇒ = 0 =⇒ + + + = 0,
t1 t2 t3 t4 t1 t2 t3 t4
as required.
1.16.4 Observe that in this problem you are only parametrising the ellipsoid! The tricky
part is to figure out the bounds in your parameters so that only the part above the plane
x + y + z = 0 is described. A common parametrisation found was:
The projection of the plane x + y + z = 0 onto the xy-plane is the line y = −x. To
be “above” this line, the angle2θ,3measured from the positive x-axis needs to be in the
1
π 3π 6 7
interval − ≤ θ ≤ . Since 415 is normal to the plane x + y + z = 0, the plane makes
4 4
1
π
an angle of with the z-axis. Recall that φ is measured from φ = 0 (positive z-axis) to
4
3π
φ = π (negative z-axis). Hence to be above the plane we need 0 ≤ φ ≤ .
4
1.16.6
2 3
−2a
6 7
1. 4a2 − 15
a2 + 1
2. x2 + y 2 = c2 + 1
Z 1
4π
3. π (c2 + 1)dc =
0 3
1.17.1
1.17.2 By CBS,
2
1 1 1 X
n
√ 1
(x1 + x2 + . . . + xn ) + +... + ≥ xi √
x1 x2 xn xi
i=1
= n2 .
1.17.3 By CBS,
(a + b + c + d)2 ≤ (1 + 1 + 1 + 1) a2 + b2 + c2 + d2 = 4 a 2 + b2 + c 2 + d 2 .
Hence,
16
(8 − e)2 ≤ 4 16 − e2 ⇐⇒ e (5e − 16) ≤ 0 ⇐⇒ 0 ≤ e ≤ .
5
16 6
The maximum value e = is reached when a = b = c = d = .
5 5
As there is equality,
(a1 , a2 , . . . , an ) = t(a31 , a32 , . . . , a3n )
for some real number t. Hence a1 = a2 = . . . = an = a, from where na = 96, na2 = 144
3
gives a = y n = 32.
2
1 +2 +··· + n n+1
n!1/n = (1 · 2 · · · n)1/n < = ,
n 2
with strict inequality for n > 1.
1.17.6 If x ∈ [−a; a], then a + x ≥ 0 and a − x ≥ 0, and thus we may use AM-GM with
a+x a−x
n = 8, a1 = a2 = · · · = a5 = and a6 = a7 = a8 = . We deduce that
5 3
a + x
a − x 8
5 3 5 +3 a 8
a+x a−x 5 3
≤ = ,
5 3 8 4
from where
55 33 a8
f (x) ≤ ,
48
a+x a−x
with equality if and only if = .
5 3
1 1 1
1, 1 + ,1 + ,...,1 + ,
n n n
has arithmetic mean
1
1+
n+1
and geometric mean
1 n/(n+1)
1+ .
n
Therefore,
1 1 n/(n+1)
1+ > 1+ ,
n+1 n
that is n+1
1 1 n
1+ > 1+ ,
n+1 n
which means
xn+1 > xn ,
giving the assertion.
2.1.2 Since polynomials are continuous functions and the image of a connected set is
connected for a continuous function, the image must be an interval of some sort. If the
image were a finite interval, then f (x, kx) would be bounded for every constant k, and so
the image would just be the point f (0, 0). The possibilities are thus
1. a single point (take for example, p(x, y) = 0),
2. a semi-infinite interval with an endpoint (take for example p(x, y) = x2 whose image
is [0; +∞[),
3. a semi-infinite interval with no endpoint (take for example p(x, y) = (xy − 1)2 + x2
whose image is ]0; +∞[),
4. all real numbers (take for example p(x, y) = x).
2.3.4 0
2.3.5 2
2.3.6 c = 0.
2.3.7 0
2.3.10 By AM-GM,
x2 y 2 z 2 (x2 + y 2 + z 2 )3 (x2 + y 2 + z 2 )2
≤ = →0
x2 +y +z2 2 2 2
27(x + y + z ) 2 27
as (x, y, z) → (0, 0, 0).
2.4.1 We have
→
− →
− →
−
F (→
−
x + h ) − F (→
−
x) = (→
−x + h ) × L(→
−x + h)−→ −x × L(→−
x)
→
− →
−
= ( x + h ) × (L( x ) + L( h )) − x × L(→
→
− →
− →
− −
x)
→
− →
− →
− →
− →
−
= x × L( h ) + h × L(→−x ) + h × L( h )
→
− →
− →
− →
− →
−
Now, we will prove that || h × L( h )|| = o h as h → 0 . For let
→
− X
n
h = hk →
−
e k,
k=1
where the →
−
e k are the standard basis for Rn . Then
→
− X
n
L( h ) = hk L(→
−
e k ),
k=1
≤ |hk | 2
||L(→
−
e k )||2
k=1 k=1
X n
→
−
= h ( ||L(→
−
e k )||2 )1/2 ,
k=1
as it was to be shewn.
→
− t
2.4.2 Assume that →
− 6
x = 0 . We use the fact that (1 + t)1/2 = 1 + + o (t) as t → 0. We
2
have
→
− →
− −
f (→
−
x + h ) − f (→
−
x ) = ||→ x + h || − →
− x
È −
→
− →
−
= (→
−x + h )•(→−x + h ) − → x
É
→ 2 →
− →
− 2 −
= −
x + 2→−x • h + h − → x
2
→
− →
−
2→
−
x • h + h
= É 2 .
→ 2 →
− →−
−x + 2→x • h + h + →
− x
−
→
− →
−
As h → 0 , É
→ 2 2 −
− →
− →
−
x + 2→
−
x • h + h + →
−
x → 2→
x .
2
→
− →
− →
− →
−
Since h = o h as h → 0 , we have
2
→
− →
−
2→
−
x • h + h →
− →
−
x•h
→
−
É 2 → + o h ,
→ →
− →
− →−
−
x + 2→
2 −x • h + h + →
x
− h
or
−
→
L( h )
1 − = o (1) .
→
−
h
→
− →
− →
− →
−
But the sinistral side → 1 as h → 0 , and the dextral side → 0 as h → 0 . This is a
→
−
contradiction, and so, such linear transformation L does not exist at the point 0 .
2.5.1 We have
2
∂f 2 r −r 2 /4t r 2 tn−2 2
= ntn−1 e−r /4t + tn e = nt n−1
+ e−r /4t ,
∂t 4t2 4
∂f 2r 2 1 2
r2 = r 2 − rtn e−r /4t = − r 3 tn−1 e−r /4t ,
∂r 4t 2
1 ∂ 2 ∂f 1 ∂ 1 3 n−1 −r 2 /4t
r = − r t e
r 2 ∂r ∂r
r 2 ∂r 2
1 3 2 n−1 −r 2 /4t 1 4 n−1 −r 2 /4t
= − r t e + r t e
r 2 2 4t
3 1 2
= − tn−1 + r 2 tn−2 e−r /4t .
2 4
We conclude that
r 2 tn−2 3 1 3
ntn−1 + = − tn−1 + r 2 tn−2 =⇒ n = − .
4 2 4 2
2.5.2 Observe that ¨
x if x ≤ y 2
f (x, y) = 2
y if x > y 2
Hence ¨
∂ 1 if x ≤ y 2
f (x, y) =
∂x 0 if x > y 2
and ¨
∂ 0 if x ≤ y 2
f (x, y) =
∂y 2y if x > y 2
and hence
′ 1 −1 2 ′ ′ 0 0
g (1, 0, 1) = , f (g(1, 0, 1)) = f (3, 0) = .
0 1 0 0 9
2.5.5 We have
∂f
(x, y) = 2xyg(x2 y),
∂x
and
∂f
(x, y) = x2 g(x2 y).
∂y
1 b
2.5.6 Differentiating both sides with respect to the parameter a, the integral is arctan +
2a3 a
b
2a2 (a2 + b2 )
2.5.10 We have
∂ ∂ ∂ ∂z ∂z
(x + z)2 + (y + z)2 = 8 =⇒ 2(1 + )(x + z) + 2 (y + z) = 0.
∂x ∂x ∂x ∂x ∂x
At (1, 1, 1) the last equation becomes
∂z ∂z ∂z 1
)+4
4(1 + = 0 =⇒ =− .
∂x ∂x ∂x 2
2 3 2 3
ey z e
6 7 6 7
2.6.1 ∇f (x, y, z) = 4xzey z 5 =⇒ (∇f )(2, 1, 1) = 42e5.
xyey z 2e
2 3 2 3
0 0
6 7 6 7
2.6.2 (∇ × f )(x, y, z) = 4 x 5 =⇒ (∇ × f )(2, 1, 1) = 4 2 5.
yexy e2
2 3
1
6 7
2.6.4 The vector 4−75 is perpendicular to the plane. Put f (x, y, z) = x2 + y 2 − 5xy +
0
2 3
2x − 5y + z
6 7
xz − yz + 3. Then (∇f )(x, y, z) = 42y − 5x − z 5. Observe that ∇f (x, y, z) is parallel
x−y
2 3
1
6 7
to the vector 4−75, and hence there exists a constant a such that
0
2 3 2 3
2x − 5y + z 1
6 7 6 7
42y − 5x − z 5 = a 4−75 =⇒ x = a, y = a, z = 4a.
x−y 0
x − 7y = −6 =⇒ a − 7a = −6 =⇒ a = 1.
Thus x = y = 1 and z = 4.
2.6.8 This is essentially the product rule: duv = udv +vdu, where ∇ acts the differential
operator and × is the product. Recall that when we defined the volume of a parallelepiped
→
− −
spanned by the vectors → −a, b,→ c , we saw that
→
− →
− →
−
a •(b ×→
−
c ) = (→
−
a × b)•→
−
c.
Treating ∇ = ∇→ →
− →
−
u + ∇→
− −v as a vector, first keeping v constant and then keeping u
constant we then see that
→
− →
− →
− →
− →
− →
− →
− →
− →
− →
− →
− →
−
∇→
u • ( u × v ) = (∇ × u ) • v ,
− ∇→
v • ( u × v ) = −∇ • ( v × u ) = −( ∇ × v ) • u .
−
Thus
→
− → − →
− → − →
− → − → − →
− → − → −
∇•(u×v) = (∇→ v )•(u×v) = ∇→
u +∇→
− − − v •( u × v ) = ( ∇ × u )• v −( ∇ × v )• u .
u •( u × v )+∇→
−
π π
2.6.11 An angle of with the x-axis and with the y-axis.
6 3
2.8.6 We have
2 3 2 3
2x − y + 1 0
6 7 6 7 2 1
∇f (x, y, z) = 4 2y − x 5 = 405 =⇒ x = − , y = − , z = 1,
3 3
2z − 2 0
2 1
whence − , − , 1 is the only critical point. The hessian is
3 3
2 3
2 −1 0
6 7
Hf (x, y, z) = 4−1 2 05 .
0 0 2
whence the hessian is positive definite at the critical point, from where we gather that the
critical point is a local minimum.
2.8.7 We have
4x3 − 4(x − y) 0
(∇f )(x, y) = = =⇒ 4x3 = 4(x − y) = −4y 3 =⇒ x = −y.
4y 3 + 4(x − y) 0
Hence
√ √
4x3 − 4(x − y) = 0 =⇒ 4x3 − 8x = 0 =⇒ 4x(x2 − 2) = 0 =⇒ x ∈ {− 2, 0, 2}.
√ √ √ √
Since x = −y, the critical points are thus (− 2, 2), (0, 0), ( 2, − 2). The Hessian is
now,
12x2 − 4 4
Hf (x, y) = ,
4 12y 2 − 4
and its principal minors are ∆1 = 12x2 − 4 and ∆2 = (12x2 − 4)(12y 2 − 4) − 16.
√ √ √ √
If (x, y) = (− 2, 2) or (x, y) = ( 2, − 2), then ∆1 = 20 > 0 and ∆2 = 384 > 0,
so the matrix is positive definite and we have a local minimum at each of these points.
f (x, x) = 2x4 > 0, f (x, −x) = 2x4 − 4x2 = 2x2 (x2 − 1).
If x is close enough to 0, = 2x2 (x2 −1) < 0, which means that the function both increases
and decreases to 0 in a neighbourhood of (0, 0), meaning that there is a saddle point there.
2.8.8 We have
2 3 2 3
8xz − 2y − 8x 0
6 7 6 7
∇f (x, y, z) = 4 −2x + 1 5 = 405 =⇒ x = 1/2; y = −1; z = 1/2.
2
4x − 2z 0
The hessian is 2 3
8z − 8 −2 8x
6 7
H = 4 −2 0 0 5.
8x 0 −2
The principal minors are 8z − 8; −4, and 8. At z = 1/2, the matrix is negative definite
and the critical point is thus a saddle point.
2.8.9 We have 2 3
2x + yz
6 7
(∇f )(x, y, z) = 42y + xz 5 ,
2z + xy
and 2 3
2 z y
6 7
Hr f = 4z 2 x5 .
y x 2
2 z
We see that ∆1 (x, y, z) = 2, ∆2 (x, y, z) = det = 4 − z 2 and ∆3 (x, y, z) =
z 2
det Hr f = 8 − 2x2 − 2 2y32
− 2z 2 + 2xyz.
0
6 7
If (∇f )(x, y, z) = 405 then we must have
0
2x = −yz,
2y = −xz,
2z = −xy,
and upon multiplication of the three equations,
8xyz = −x2 y 2 z 2 ,
that is,
xyz(xyz + 8) = 0.
Clearly, if xyz = 0, then we must have at least one of the variables equalling 0, in which
case, by virtue of the original three equations, all equal 0. Thus (0, 0, 0) is a critical
point. If xyz = −8, then none of the variables is 0, and solving for x, say, we must have
8
x = − , and substituting this into 2x + yz = 0 we gather (yz)2 = 16, meaning that
yz
either yz = 4, in which case x = −2, or zy = −4, in which case x = 2. It is easy to
see then that either exactly one of the variables is negative, or all three are negative. The
other critical points are therefore (−2, 2, 2), (2, −2, 2), (2, 2, −2), and (−2, −2, −2).
At (0, 0, 0), ∆1 (0, 0, 0) = 2 > 0, ∆2 (0, 0, 0) = 4 > 0, ∆1 (0, 0, 0) = 8 > 0, and thus
(0, 0, 0) is a minimum point. If x2 = y 2 = z 2 = 4, xyz = −8, then ∆2 (x, y, z) = 0, ∆3 =
−32, so these points are saddle points.
2.8.10 We have 2 3
2xy + 2
6 7
(∇f )(x, y, z) = 4x2 + 2yz 5 ,
y2 − 1
and 2 3
2y 2x 0
6 7
Hr f = 42x 2z 2y 5 .
0 2y 0
2y 2x
We see that ∆1 (x, y, z) = 2y, ∆2 (x, y, z) = det = 4yz − 4x2 and ∆3 (x, y, z) =
2x 2z
det Hr f = −8y 3 . 2 3
0
6 7
If (∇f )(x, y, z) = 405 then we must have
0
xy = −1,
x2 = −2yz,
y = ±1,
1 1 1
and hence (1, −1, ), and (−1, 1, − ) are the critical points. Now, ∆1 (1, −1, ) = −2,
2 2 2
1 1 1
∆2 (1, −1, ) = −6, and ∆3 (1, −1, ) = 8. Thus (1, −1, ) is a saddle point. Similarly,
2 2 2
1 1 1
∆1 (−1, 1, − ) = 2, ∆2 (−1, 1, − ) = −6, and ∆3 (−1, 1, − ) = −8, shewing that
2 2 2
1
(1, −1, ) is also a saddle point.
2
2.8.11 We find 2 3
4yz − 4x3
6 7
∇f (x, y, z) = 44xz − 4y 3 5 .
4xy − 4z 3
Assume ∇f (x, y, z) = 0. Then
Thus xyz ≥ 0, and if one of the variables is 0 so are the other two. Thus (0, 0, 0) is the
only critical point with at least one of the variables 0. Assume now that xyz 6= 0. Then
1
(xyz)3 = x4 y 4 z 4 = (xyz)4 =⇒ xyz = 1 =⇒ yz = =⇒ x4 = 1 =⇒ x = ±1.
x
Similarly, y = ±1, z = ±1. Since xyz = 1, exactly two of the variables can be negative.
Thus we find the following critical points:
(0, 0, 0), (1, 1, 1), (−1, −1, 1), (−1, 1, −1), (1, −1, −1).
The Hessian is 2 3
−12x2 4z 4y
6 7
Hx f = 4 4z −12y 2 4x 5.
4y 4x −12z 2
If 1 = xyz = x2 = y 2 = z 2 , we have ∆1 = −12x2 = −12 < 0, ∆2 = 144x2 y 2 − 16z 2 =
144 − 16 = 128 > 0, and
This means that for xyz 6= 0 the Hessian is negative definite and the function has a
local maximum at each of the four points (1, 1, 1), (−1, −1, 1), (−1, 1, −1), (1, −1, −1).
Observe that at these critical points f = 1. Now f (0, 0, 0) = 0 and f (−1, 1, 1) = −7.
4 − 2x − y − z = 0, 4 − x − 2y − z = 0, 4 − x − y − 2z = 0 =⇒ x = y = z = 1.
In this case 2 3
−2 −1 −1
6 7
Hf (1, 1, 1) = 4−1 −2 −15
−1 −1 −2
and the principal minors are ∆1 = −2 < 0, ∆2 = 3 > 0, and ∆3 = −4 < 0, so the
matrix is negative definite and we have a local maximum at (1, 1, 1).
Thus as x → 0+ then f (x, x, x) > 0 and f (x, −x, x) < 0, which means that in some
neighbourhood of (0, 0, 0) the function is both decreasing towards 0 and increasing to-
wards 0, which means that (0, 0, 0) is a saddle point.
2 2 2
2.8.13 To facilitate differentiation observe that g(x, y, z) = (xe−x )(ye−y )(ze−z ). Now
2 2 2 2
3
(1 − 2x2 )(yz)(e−x )(e−y )(e−z )
6 2 2 2 7
∇g(x, y, z) = 4(1 − 2y 2 )(xz)(e−x )(e−y )(e−z )5 .
2 2 2
(1 − 2z 2 )(xy)(e−x )(e−y )(e−z )
The function is 0 if any of the variables is 0. Since the function clearly assumes positive
and negative values, we can discard any point with a 0. If ∇(x, y, z) = 0, then x =
1 1 1
± √ ; y = ± √ z = ± √ . We find
2 2 2
2 3
(4x3 − 6x)(yz) (1 − 2x2 )(1 − 2y 2 )z (1 − 2x2 )(1 − 2z 2 )y
6 7
Hx g = t(x, y, z) 4(1 − 2y 2 )(1 − 2x2 )z (4y 3 − 6y)(xz) (1 − 2y 2 )(1 − 2z 2 )x5 ,
(1 − 2z 2 )(1 − 2x2 )y (1 − 2z 2 )(1 − 2y 2 )x (4z 3 − 6z)(xy)
2 2 2
with t(x, y, z) = (e−x )(e−y )(e−z ). Since at the critical points we have 1 − 2x2 =
1 − 2y 2 = 1 − 2z 2 = 0, the Hessian reduces to
2 3
(4x3 − 6x)(yz) 0 0
6 7
Hx g = (e−3/2 ) 4 0 (4y 3 − 6y)(xz) 0 5.
0 0 (4z 3 − 6z)(xy)
We have
∆1 = (4x3 − 6x)(yz)
∆2 = (4x3 − 6x)(4y 3 − 6y)(xyz 2 )
∆3 = (4x3 − 6x)(4y 3 − 6y)(4z 3 − 6z)(x2 y 2 z 2 ).
Also,
3 3
1 1 √ 1 1 √
4 √ −6 √ = −2 2 < 0, 4 −√ − 6 −√ = 2 2 > 0.
2 2 2 2
This means that if an even number of the variables is negative (0 or 2), then we the Hessian
is negative definite, and if an odd numbers of the variables is positive (1 or 3), the Hessian
is positive definite. We conclude that we have local maxima at
1 1 1 1 1 1 1 1 1 1 1 1
( √ , √ , √ ), (− √ , − √ , √ ), (− √ , √ , − √ ), ( √ , − √ , − √ )
2 2 2 2 2 2 2 2 2 2 2 2
Hence
∂f
(x, y) = 2xG′ (x2 + y) + G′ (y 2 − x) = 2xg(x2 + y) + g(y 2 − x);
∂x
∂f
(x, y) = G′ (x2 + y) − 2yG′ (y 2 − x) = g(x2 + y) − 2yg(y 2 − x).
∂y
This gives
∂f ∂f
(0, 0) = g(0) = (0, 0) = 0,
∂x ∂y
so (0, 0) is a critical point. Now, the Hessian of f is
2g(x2 + y) + 4x2 g ′ (x2 + y) − g ′ (y 2 − x) 2xg ′ (x2 + y) + 2yg ′ (y 2 − x)
Hf (x, y) = ,
2xg ′ (x2 − y) + 2yg ′ (y 2 − x) g ′ (x2 + y) − 2g(y 2 − x) − 4y 2 g ′ (y 2 − x)
and so ′
−g (0) 0
Hf (0, 0) = .
0 g ′ (0)
Regardless of the sign of g ′ (0), the determinant of of this last matrix is −(g ′ (0))2 < 0,
and so (0, 0) is a saddle point.
√
144 − 16x2
2.8.15 Since the coordinates (x, ), −3 ≤ x ≤ 3 describe an ellipse centred
3 p
at the origin and semi-axes 3 and 4, and the coordinates (y, 4 − y 2 ), −2 ≤ y ≤ 2
describe a circle centred at the origin with radius 2, the problem reduces to finding the
minimum between the boundaries of the circle and the ellipse. Geometrically this is easily
seen to be 1.
2.9.1 We have
2 3 2 3
bc 2b + 2c
6 7 6 7
∇(abc) = λ∇(2ab + 2bc + 2ca − S) =⇒ 4ca5 = λ 42a + 2c5
ab 2b + 2a
bc = 2λ(b + c)
=⇒ ca = 2λ(a + c)
ab = 2λ(b + a)
2.9.2
2 3
a
6 7
1. The vector 4 b 5 is perpendicular to Π. Hence, the equation of the perpendicular
c
passing through P is
Ǒ
Ǒ 2 3
x 1 a
6 7
y = 1 + t 4 b 5 =⇒ x = 1 + ta, y = 1 + tb, z = 1 + tc.
z 1 c
2. Let f (x, y, z) = (x − 1)2 + (y − 1)2 + (z − 1)2 be the square of the distance from
P to a point on the plane and let g(x, y, z) = ax + by + cz − d. Using Lagrange
multipliers,
2 3 2 3
2(x − 1) a
6 7 6 7
∇f (x, y, z) = λ∇g(x, y, z) =⇒ 42(y − 1)5 = λ 4 b 5 =⇒ 2(x−1) = λa, 2(y−1) = λb, 2(z−1) = λc.
2(z − 1) c
Since (1, 1, 1) is not on the plane and abc 6= 0, we gather that λ 6= 0. Now,
λa λb λc
x =1+ , y =1+ , z =1+ .
2 2 2
Putting these into the equation of the plane,
λa λb λc d−a−b−c
a 1+ +b 1+ +c 1+ = d =⇒ λ = 2 · 2 .
2 2 2 a + b2 + c 2
Then the coordinates of P ′ are
λa d−a−b−c λb d−a−b−c λc d−a−b−c
x = 1+ = 1+a· 2 , y = 1+ = 1+b· 2 , z = 1+ = 1+c· 2 ,
2 a + b2 + c 2 2 a + b2 + c 2 2 a + b2 + c 2
as before.
3. Consider the function
2
d − ax − by
t(x, y) = (x − 1)2 + (y − 1)2 + −1 ,
c
which is the square of the distance from a point (x, y, z) on the plane to the point
(1, 1, 1).
Now, 2 3
a d − ax − by
2(x − 1) − 2 −1 0
∇t(x, y) = 4 c c 5 =
b d − ax − by 0
2(y − 1) − 2 −1
c c
which implies
2.9.5 We put g(x, y) = 5x2 + 6xy + 5y 2 − 8 and argue using Lagrange multipliers. We
have
2x 10x + 6y
∇f (x, y) = λ∇g(x, y) =⇒ =λ .
2y 6x + 10y
This gives the three equations
The linear system (the first two equations) will have the unique solution (0, 0) as long as
25(λ−1)2 −9 6= 0, but this solution does not lie on the third equation. If 25(λ−1)2 −9 = 0,
then we deduce that x = ±y. Substituting this into the third equation we gather that
√ 1
10x2 ± 6x2 = 8, resulting in x = ± 2 or x = ± √ . Taking into account the third equa-
√ √ √ √2 √ √ √ √
tion, the feasible values are ( 2, − 2), (− 2, 2), (1/ 2, 1/ 2), (−1/ 2, −1/ 2)
The desired maximum is thus
√ √ √ √
f (− 2, 2) = f ( 2, − 2) = 4
8 6 8 6 x2 + y 2 5 8
5x2 +6xy +5y 2 = 8 =⇒ x2 +y 2 = − xy ≥ − · =⇒ x2 +y 2 ≥ · = 1.
5 5 5 5 2 8 5
2.9.6 Put g(x, y) = xp + y p − 1. We need a = pλxp−1 and b = pλy p−1 . Clearly then ,
λ 6= 0. We then have
1/(p−1) 1/(p−1)
a b
x= , y= .
λp λp
Thus p/(p−1) p/(p−1)
p p a b
1=x +y = + ,
λp λp
which gives
p/(p−1) p/(p−1) (p−1)/p
a b
λ= + .
p p
This gives
a1/(p−1) b1/(p−1)
x= , y= .
(a1/(p−1)
+b 1/(p−1) )1/p (a 1/(p−1) + b1/(p−1) )1/p
ap/(p−1) bp/(p−1)
ax + by ≤ + .
(a1/(p−1) + b1/(p−1) )1/p (a1/(p−1) + b1/(p−1) )1/p
from where √
14
λ =1± .
2
This gives the two points
2 4 6
(x, y, z) = 1 + √ ,2 + √ ,3 + √
14 14 14
and
2 4 6
(x, y, z) =1− √ ,2 − √ ,3 − √ .
14 14 14
√
12 14
The first point gives an absolute maximum of 18 + and the second an absolute
√ 7
12 14
minimum of 18 − .
7
2.9.8 Observe that the ellipse is symmetric about the origin. Now maximise and minimise
the distance between a point on the ellipse and the origin. If a and b are the semi-axes,
you will find that 2a = 2 and 2b = 6
2.9.10 One can use Lagrange multipliers here. But perhaps the easiest approach is to
put y = 1 − x and maximise
È
f (x) = x + x(1 − x).
Since
(1 − 2x)2 1
f ′′ (x) = − − p < 0,
4(x(1 − x))3/2 x(1 − x)
the value sought is a maximum. This maximum is thus
√ √
1 2 1 2
f + = + .
2 4 2 2
2.9.11 Claim: the function achieves its maximum on the boundary of the triangle. To
prove this claim we have to prove that there are no critical points strictly inside the trian-
gle. For this we compute the gradient and set it equal to the zero vector:
−axa−1 y b e−(x+y ) 0
(∇f )(x, y) = = =⇒ x = 0 or y = 0,
−bxa y b−1 e−(x+y ) 0
which means that the critical points occur on the boundary. Since the function is iden-
tically 0 for x = 0 or y = 0, we only need to look on the line x + y = 1 for the maxima.
Hence we maximise f subject to the constraint x + y = 1. Since x + y = 1, we can see
that f (x, y) = xa y b e−(x+y ) = xa y b e−1 on the line, so the problem reduces to maximising
h(x, y) = xa y b subject to the constraint x + y = 1. Using Lagrange multipliers,
a−1 b
ax y 1
(∇h)(x, y) = λ(∇g)(x, y) =⇒ a b−1
=λ ,
bx y 1
which in turn
b a
=⇒ axa−1 y b = λ = bxa y b−1 =⇒ ay = bx =⇒ ay = b(1−y) =⇒ y = , x= .
a+b a+b
Finally,
a b
a b
f (x, y) = xa y b e−(x+y ) ≤ xa y b e−1 ≤ e−1 .
a+b a+b
so f is doubly-periodic
of period π in
each argument. This means that we may restrict
ourselves to x ∈ 0; π and y ∈ 0; π . Put g(x, y) = x − y. Then
π π
If cos 2y = 0 then y = and so x = . In this case
4 2
π π 1
f , = .
2 4 2
If, however, − cos 2x 6= 0 then
π π
− sin 2x = − cos 2x =⇒ tan 2x = 1 =⇒ x = =⇒ y = − .
8 8
In this case we have π π π
f ,− = 2 cos2 ≈ 1.707106781.
8 8 8
If, however, cos 2y 6= 0 then
3π 5π
sin 2y = − cos 2y =⇒ tan 2y = −1 =⇒ y = =⇒ x = .
8 8
In this case we have
5π 3π 3π
f , = 2 cos2 ≈ .2928932190.
8 8 8
3.3.1
1. Let L1 : y = x + 1, L2 : −x + 1. Then
Z Z Z
xdx + ydy = xdx + ydy + xdx + ydy
C ZL11 L 2Z
1
= xdx(x + 1)dx + xdx − (−x + 1)dx
−1 0
= 0.
√
Also, both on L1 and on L2 we have ||dx|| = 2dx, thus
Z Z Z
xy||dx|| = xy||dx|| + xy||dx||
C L1 Z L2 Z 1
√ 1 √
= 2 x(x + 1)dx − 2 x(−x + 1)dx
−1 0
= 0.
π π
2. We put x = sin t, y = cos t, t ∈ − ; . Then
2 2
Z Z π /2
xdx + ydy = (sin t)(cos t)dt − (cos t)(sin t)dt
C −π /2
= 0.
È
Also, ||dx|| = (cos t)2 + (− sin t)2 dt = dt, and thus
Z Z π /2
xy||dx|| = (sin t)(cos t)dt
C −π /2
(sin t)2 π /2
=
2 −π /2
= 0.
√
3.3.2 Let Γ1 denote the straight line segment path from O to A = (2 3, 2) and Γ2 denote
π
the arc of the circle centred at (0, 0) and radius 4 going counterclockwise from θ = to
6
π
θ= .
5
←→ x
Observe that the Cartesian equation of the line OA is y = √ . Then on Γ1
3
x x 4
xdx + ydy = xdx + √ d √ = xdx.
3 3 3
Hence
Z Z 2
√
3
4
xdx + ydy = xdx = 8.
Γ1 0 3
π
On the arc of the circle we may put x = 4 cos θ, y = 4 sin θ and integrate from θ =
6
π
to θ = . Observe that there
5
xdx + ydy = (cos θ)d cos θ + (sin θ)d sin θ = − sin θ cos θdθ + sin θ cos θdθ = 0,
To solve this problem using Maple you may use the code below.
> with(Student[VectorCalculus]):
> LineInt( VectorField( <x,y> ), Line( <0,0>, <2*sqrt(3),2> ))
> +LineInt( VectorField( <x,y> ), Arc(Circle( <0,0>, 4), Pi/6, Pi/5) );
3.3.3 Using the parametrisations from the solution of problem 3.3.3, we find on Γ1 that
È É
1 2
x||dx|| = x (dx)2 + (dy)2 = x 1 + dx = √ xdx,
3 3
whence
Z Z 2
√
3
2 √
x||dx|| = √ xdx = 4 3.
Γ1 0 3
On Γ2 that
È È
x||dx|| = x (dx)2 + (dy)2 = 16 cos θ sin2 θ + cos2 θdθ = 16 cos θdθ,
whence
Z Z π /5
π π π
x||dx|| = 16 cos θdθ = 16 sin − 16 sin = 4 sin − 8.
Γ2 π /6 5 6 5
To solve this problem using Maple you may use the code below.
> with(Student[VectorCalculus]):
> PathInt(x, [x,y]= Line( <0,0>, <2*sqrt(3),2> ))
> +PathInt(x, [x,y]=Arc(Circle( <0,0>, 4), Pi/6, Pi/5) );
3π π
Maple gives 16 cos rather than our 16 sin . To check that these two are indeed
10 5
the same, use the code
> is(16*cos(3*Pi/10)=16*sin(Pi/5));
which returns true.
3.3.4 The curve lies on the sphere, and to parametrise this curve, we dispose of one of
the variables, y say, from where y = 1 − x and x2 + y 2 + z 2 = 1 give
x2 + (1 − x)2 + z 2 = 1 =⇒ 2x2 − 2x + z 2 = 0
1 2 1
=⇒ 2 x− + z2 =
2 2
1 2
=⇒ 4 x− + 2z 2 = 1.
2
So we now put
We must integrate on the 2 3side of the plane that can be viewed from the point (1, 1, 0)
1
6 7 1 2
(observe that the vector 415 is normal to the plane). On the zx-plane, 4 x − +2z 2 =
2
0
1 is an ellipse. To obtain a positive parametrisation we must integrate from t = 2π to t = 0
(this is because when you look at the ellipse from the point (1, 1, 0) the positive x-axis is
to your left, and not your right). Thus
I Z 0
sin t 1 cos t
zdx + xdy + ydz = √ d +
Z 0 2 2 2
Γ 2π
1 cos t 1 cos t
+ + d −
Z2π0 2 2
2 2
1 cos t sin t
+ − d √
2
Z 0 2π 2 2
sin t cos t cos t sin t 1
= + √ + − √ dt
2π 4 2 2 4 2 2
π
= √ .
2
3.5.1 2
1
3.5.2
3
15π
3.5.3
16
Z Z Z √
1 x
xydxdy = x y dy dx
0 x2
D Z 1
1
= x(x − x4 ) dx
0 2
1
= .
12
Z Z Z 1 Z y Z 1 Z 1
x2 dxdy + y 2 dxdy = x2 dxdy + y 2 dxdy
0 0 0 y
x≤y y ≤x
Z Z 1
y3
1
= dy + y 2 − y 3 dy
0 3 3 0 4 1
y 4 1 y y
= + −
12 0 3 4 0
1 1 1
= + −
12 3 4
1
= .
6
21
3.5.7
8
The integral is
Z 2
√
3 Z √16−x 2 Z 2 p√
3 √
2
3
√ xdydx = x 16 − x + x − 4 dx
0 − 33 x+4 0 3
√ 2√ 3
1 3 3
= − (16 − x2 )3/2 + x − 2x2
3 9 0
8
= .
3
3.5.9 e − 1
3.5.10 We have
Z Z Z
min(x, y 2 )dA = xdA + y 2 dA
[0;1]2 [0;1]2 [0;1]2
x≤y 2 y 2 <x
Z 1 Z y2 Z 1Z 1
= xdxdy + y 2 dxdy
0Z 0 Z 10 y 1 2
1 y 2
1
= x2 dy + y 2 x dy
2 Z0 0
Z 1 0 y2
1
1 4 2 4
= y dy + (y − y )dy
2 0 0
1 2
= +
10 15
7
= .
30
←→ x
3.5.11 Begin by finding the Cartesian equations of the various lines: for OA is y =
3
←→ ←→
(0 ≤ x ≤ 1), for AB is y = 3x − 8 (3 ≤ x ≤ 4), and for BO is y = x (0 ≤ x ≤ 4).
Z Z 3 Z x Z 4 Z x
xydA = xydydx + xydydx
R 0Z x/3 3Z 3x−8
3 x 4 x
1 1
= xy 2 dx + xy 2 dx
2 Z0 x/3
2 3 Z 3x−8
1
3
x 2
1
4
= x x2 − dx + x x2 − (3x − 8)2 dx
2 Z0 9Z 2 3
4
3
1
4
= x3 dx + −8x3 + 48x2 − 64x dx
9 0 2 3
= 9+9
= 18.
To solve this problem using Maple you may use the code below.
> with(Student[VectorCalculus]):
> int(x*y, [x,y]=Triangle(<0,0>,<3,1>,<4,4>));
Maple can also provide the limits of integration, but this command is limited, since
Maple is quite whimsical about which order of integration to choose. It also evaluates
expressions that it deems below its dignity to return unevaluated.
= dA + 2 dA + 3 dA + 4 dA + 5 dA
[0;2]2 [0;2]2 [0;2]2 [0;2]2 [0;2]2
1≤x+y 2 <2 2≤x+y 2 <3 3≤x+y 2 <4 4≤x+y 2 <5 5≤x+y 2 <6
By looking at the regions (as in figures A.13 through A.17 below) (I am omitting the details
of the integrations, relying on Maple for the evaluations), we obtain
Z Z 1 Z √ 2−x Z 2 Z √ 2−x
4 4√ 2 2 4√
dA = √ dydx + dydx = − + 2+ =− + 2.
0 1−x 1 0 3 3 3 3 3
[0;2]2
1≤x+y 2 <2
Z Z 2 Z √3−x
√ 8√ 4
2 dA = 2 √ dy dx = 4 3 − 2− .
0 2−x 3 3
[0;2]2
2≤x+y 2 <3
Z Z 2 Z √4−x
√ √
3 dA = 3 √ dy dx = 18 − 6 3 − 4 2.
0 3−x
[0;2]2
3≤x+y 2 <4
Z Z 1 Z 2 Z 2 Z √ 5−x
40 √ 64 √ 16 √
4 dA = 4 √ dy dx+4 √ dy dx = − +8 3+ −16 3+ 2.
0 4−x 1 4−x 3 3 3
[0;2]2
4≤x+y 2 <5
Z Z 2 Z 5
50 √
5 dA = 5 √ dy dx = − + 10 3.
1 5−x 3
[0;2]2
4≤x+y 2 <5
1 1 1 1Appendix A
0 0 0 0
2 0 1 2 0 1 2 0 1 2 0 1 2
1
Figure A.13: 1 ≤ x + Figure A.14: 2 ≤ x + Figure A.15: 3 ≤ x + Figure A.16: 4 ≤ x +
y2 < 2. y2 < 3. y2 < 4. y2 < 5.
0
0 1 2
Figure A.17: 5 ≤ x +
y2 < 6.
= 1dydx + 2dydx
1 1−x 1 2−x
= 4.
Z 1 Z √ 4−x 2 Z 2 Z √4−x 2
7
3.5.15 √ xdydx + xdydx = .
0 1−x 2 1 0 3
Z 1 Z 3 Z 2 Z −1 Z 2 Z 1 Z 2 Z 3
3.5.16 xdxdy + xdxdy + xdxdy + xdxdy = 7.
0 −2 1 −2 1 0 1 2
3.5.18 Upon splitting the domain of integration, we find that the integral equals
Z Z Z y 2
2 y2 2
πx 2y πx
sin dx dy = − cos dy
2y π 2y
1 y 1
Z 2
y
2y πy
= − cos −dy
π1 2
4(π + 2)
= ,
π3
upon integrating by parts.
3.5.19 The integral is 0. Observe that if (x, y) ∈ D then (−x, y) ∈ D. Also, f (−x, y) =
−f (x, y).
Z √ Z 1
√ √
2/ 5 4−y 2
2 8 5
3.5.20 √ dxdy = + 4 arcsin .
−2/
√
5 1
−2 4−y 2 5 5
3.5.25 Let
D1 = {(x, y) ∈ R2 | − 1 ≤ x ≤ 1, x ≤ y},
D2 = {(x, y) ∈ R2 | − 1 ≤ x ≤ 1, x > y}.
Then D = D1 ∪ D2 , D1 ∩ D2 = ∅ and so
Z Z Z
f (x, y)dxdy = f (x, y)dxdy + f (x, y)dxdy.
D1 D2
D
By symmetry, Z Z
f (x, y)dxdy = f (x, y)dxdy,
D1 D2
and so Z Z
f (x, y)dxdy = 2 f (x, y)dxdy
D ZD11 Z 1
= 2 (y − x) dy dx
Z 1−1 x
= (1 − 2x + x2 ) dx
−1
8
= .
3
3.5.26 The line joining A, and B has equation y = −x − 2, line joining B, and C has
equation y = −7x + 10, and line joining A, and C has equation y = 2x + 1. We split
the triangle along the vertical line x = 1, and integrate first with respect to y. The desired
integral is then
Z Z 1 Z 2x+1
(2x + 3y + 1) dxdy = (2x + 3y + 1)dy dx
D −1 Z 2
Z −7x+10
−x−2
+ (2x + 3y + 1)dy dx
Z 1 1 −x−2
21 2 3
= x + 9x − dx
−1 Z 2 2
2
+ 60x2 − 198x + 156 dx
1
= 4−1
= 3.
3.5.30 Put f (x, y) = xy + y 2 . If I, II, III, IV stand for the intersection of the region with
each quadrant, then
Z Z Z
f (x, y)dxdy = f (−x, y)d(−x)dy = − f (x, y)dxdy,
II II I
Z Z Z
f (x, y)dxdy = f (x, −y)dxd(−y) = − f (x, y)dxdy,
IV IV I
and Z Z Z
f (x, y)dxdy = f (−x, −y)d(−x)d(−y) = + f (x, y)dxdy.
III III I
Thus
Z Z Z Z Z
(xy + y 2 )dxdy = f (x, y)dxdy + f (x, y)dxdy + f (x, y)dxdy + f (x, y)dxdy
S Z
I IZI IZ
II ZV
I
3.5.31 We split the rectangle [0; a] × [0; b] into two triangles, depending on whether bx <
ay or bx ≥ ay. Hence
Z a Z b Z Z
2
x 2 ,a 2 y 2 ) 2
x 2 ,a 2 y 2 ) 2
x 2 ,a 2 y 2 )
emax(b dydx = emax(b dydx + emax(b dydx
0 0
Z
bx<ay Z bx≥ay
a2 y 2 b2 x2
= e dydx + e dydx
bx<ay bx≥ay
Z b Z ay /b Z a Z bx/a
2
y2 2 2
= ea
dxdy + ea y dydx
Z0 b 0 a 2 y 2 Z a 0
2 2
0
aye bxea y
= dy + dx
0 b 0 a
2 2 2 2
ea b − 1 ea b − 1
= +
2ab
2 2
2ab
ea b − 1
= .
ab
1
3.5.32 Observe that x ≥ (x + y)2 ≥ 0. Hence we may take the positive square root
√ 2 √
giving y ≤ 2x − x. Since y ≥ 0, we must have 2x − x ≥ 0 which means that x ≤ 2.
The integral equals
Z 2 Z √ 2x−x Z 2
√ 2 √ √
xydy dx = x( 2x − x)3/2 dx
0 0 3 0
Z √2
4 √
= u2 (u 2 − u2 )3/2 du
3 Z0
1
1
= (1 − v 2 )3/2 (1 + v)2 dv
6 −1
Z π /2
1
= cos4 θ(1 + sin2 θ)dθ
6 −π /2
7π
= .
96
3.5.33 Observe that
Z (
a 0 if a is an integer
sin 2πx dx = 1
0 (1 − cos 2πa) if a is not an integer
2π
Thus Z a
sin 2πx dx = 0 ⇐⇒ a is an integer.
0
Now
N Z
X
sin 2πx sin 2πy dxdy = 0
k=1 Rk
we deduce that at least one of the sides of R is an integer, finishing the proof.
3.5.34 We have
Z 1Z 1 Z 1 n Z
Y 1 Y
n
1 1
··· (x1 x2 · · · xn )dx1 dx2 . . . dxn = xk dxk = = n.
0 0 0 0 2 2
k=1 k=1
3.5.35 This is
Z 1Z 1 Z n Z Z Z
1 X
n
X 1 1 1
··· xk dx1 dx2 . . . dxn = ··· xk dx1 dx2 . . . dxn
0 0 0 k=1 k=1 0 0 0
Xn
1
=
2
k=1
n
= .
2
3.5.41 The integral equals
Z Z 3
Z 4−x
1 dy
dxdy = dy dx
(x + y)4 1 1 (x + y)4
D Z 3h i4−x
1
= − (x + y)−3 dx
1Z 3
3
1
1 1 1
= − dx
3 1 (1 + x)3 64
1
= .
48
3.5.45 The integral equals
Z Z 2/3
Z x+1
Z 4
Z 2− x
2
xdxdy = dy x dx + dy x dx
−1 0 2/3 0
D
Z Z x
2/3 4
= x(x + 1) dx + x 2− dx
−1 2/3 2
275
= .
54
3.5.46 Make the change of variables xk = 1 − yk . Then
Z 1Z 1 Z 1 π
I = ··· cos2 (x1 + x2 + · · · + xn ) dx1 dx2 . . . dxn
0 0 0 2n
equals
Z 1 Z 1 Z 1 π
··· sin2 (y1 + y2 + · · · + yn ) dy1 dy2 . . . dyn .
0 0 0 2n
2 2 1
Since sin t + cos t = 1, we have 2I = 1, and so I = .
2
u+v u−v
3.6.1 ➊ Put x = and y = . Then x + y = u and x − y. Observe that D ′
2 2
is the triangle in the uv plane bounded by the lines u = 0, u = 1, v = u, v = −u.
Its image under Φ is the triangle bounded by the equations x = 0, y = 0, x + y = 1.
Clearly also
1
dx ∧ dy = du ∧ dv.
2
➋ From the above
Z Z
2
−y 2 1
(x + y)2 ex dA = u2 euv dudv
D 2 ZD ′ Z
1 u
1
= u2 euv dudv
2 0 −u
Z 1
1 2 2
= u(eu − e−u )du
2 0
1
= (e + e−1 − 2).
4
dv = −2xdx + 2ydy.
Taking the wedge product of differential forms,
Hence
1
f (x, y)dx ∧ dy = (y 4 − x4 ) du ∧ dv
2(y 2 + x2 )
1 2
= (y − x2 )du ∧ dv
2
v
= du ∧ dv
2
The region transforms into
∆ = [a; b] × [0; 1].
The integral becomes
Z Z
f (x, y)dx ∧ dy = v du ∧ dv
D ∆ Z b Z 1
1
= du v dv
2 a 0
b−a
= .
4
3.6.5 ➊ Formally,
Z 1 Z 1 Z 1 Z 1
dxdy
= (1 + xy + x2 y 2 + x3 y 3 + · · · )dxdy
1 − xy
0 0 Z0 1 0 1
xy 2 x2 y 3 x3 y 4
= y+ + + +··· dx
Z0 1 2 3 4 0
2 3
x x x
= (1 + + + + · · · )dx
0 2 3 4
1 1 1
= 1 + 2 + 2 + 2 + ···
2 3 4
➋ This change of variables transforms the square [0; 1] × [0; 1] in the xy plane into
the square with vertices at (0, 0), (1, 1), (2, 0), and (1, −1) in the uv plane. We will
split this region of integration into two disjoint triangles: T1 with vertices at (0, 0),
(1, 1), (1, −1), and T2 with vertices at (1, −1), (1, 1), (2, 0). Observe that
1
dx ∧ dy = du ∧ dv,
2
Z 1 Z 1 Z
dxdy 1 du ∧ dv
= u 2 −v 2
0 0 1 − xy 2 1− 4
T 1 ∪T Z
Z 1
2
Z u Z 2 2−u
dv dv
= 2 du + 2 du,
0 −u 4 − u2 + v 2 1 u−2 4 − u2 + v 2
as desired.
Z t
dω
= arctan t.
0 1 + Ω2
Z √
15 Z √ 16−y 2 Z Z √ 15
1
xydxdy = 15y − y 3 dy
1 1 2 1
49
= .
2
3.7.12 Put
D ′ = {(x, y) ∈ R2 : y ≥ x, x2 + y 2 − y ≤ 0, x2 + y 2 − x ≤ 0}.
3.7.13 Observe that D = D2 \D1 where D2 is the disk limited by the equation x2 +y 2 = 1
and D1 is the disk limited by the equation x2 + y 2 = y. Hence
Z Z Z
dxdy dxdy dxdy
2 + y 2 )2
= 2 + y 2 )2
− .
D (1 + x D2 (1 + x D1 (1 + x2 + y 2 )2
and
Z Z π /2 Z sin θ Z π /2
dxdy ρ sin2 θdθ
= 2 dρdθ =
(1 + x2 + y 2 )2 2
(1 + ρ )2
1 + sin2 θ
D1 Z 0
+∞
0 0
√
dt dt π π 2
= − 2 = − .
0 t2 + 1 2t + 1 2 4
(We evaluated this last integral using t = tan θ) Finally, the integral equals
√ √
π π π 2 π 2
− − = .
2 2 4 4
3.7.14 We have
2xdx = cos θdρ − ρ sin θdθ, 2ydy = sin θdρ + ρ cos θdθ,
whence
4xydx ∧ dy = ρdρ ∧ dθ.
It follows that
p 1 2 2 p
x3 y 3 1 − x4 − y 4 dx ∧ dy = (x y )( 1 − x4 − y 4 )(4xy dx ∧ dy)
4 p
1 3
= (ρ cos θ sin θ 1 − ρ2 )dρ ∧ dθ
4
Observe that
x4 + y 4 ≤ 1 =⇒ ρ2 cos2 θ + ρ2 sin2 θ ≤ 1 =⇒ ρ ≤ 1.
Since the integration takes place on the first quadrant, we have 0 ≤ θ ≤ π/2. Hence the
integral becomes
Z π /2 Z 1 Z π /2 Z 1
1 3 p 1 p
(ρ cos θ sin θ 1 − ρ2 )dρdθ = cos θ sin θdθ ρ3 1 − ρ2 dρ
0 0 4 4 0 0
1 1 2
= · ·
4 2 15
1
= .
60
3.7.15 ➊ Using polar coordinates
Z 2π Z a
2 2
Ia = ρe−ρ dρ dθ = π(1 − e−a ).
0 0
3.7.16
Z Z 2π Z 4
1 r
dA = drdθ
x2 + xy + y 2 0 2 r 2 + r 2 sin θ cos θ
4≤x 2 +y 2 ≤16
Z 2π Z 4
1
= drdθ
Z
0 r(1 + sin θ cos
2 θ)
Z
2π 4
dθ dr
=
Z0 2π 1 + sin θ cos θ 2 r
dθ
= log 2
Z
0 1 + sin θ cosθ
2π
dθ
= 2 log 2
Z0 π 2 + sin 2θ
dθ
= 4 log 2
0 2 + sin 2θ
= 4I log 2,
Z π
dθ
so the problem reduces to evaluate I = . To find this integral, we now
0 2 + sin 2θ
use what has been dubbed as “the world’s sneakiest substitution”1 : we put tan θ = t.
1
by Michael Spivak, whose Calculus book I recommend greatly.
In so doing we have to pay attention to the fact that θ 7→ tan θ is not continuous on
π π
[0; π], so we split the interval of integration into two pieces, [0; π] = [0; ]∪] ; π]. Then
2 2
2t 1 − t2 dt
sin 2θ = , cos 2θ = , dθ = . Hence
1 + t2 1 + t2 1 + t2
Z π Z π /2 Z π
dθ dθ dθ
= +
0 2 + sin 2θ 0 2 + sin 2θ π /2 2 + sin 2θ
Z +∞ dt Z 0 dt
1+t 2 1+t 2
= 2t
+ 2t
2 + 1+t 2 −∞ 2 + 1+t 2
Z0 +∞ Z 0
dt dt
= +
2(t2 + t + 1) 2
−∞ 2(t + t + 1)
0
Z +∞ Z 0
2 dt 2 dt
= 2t 1
+ 2t 1
3 0 (√ + √ )2 + 1 3 −∞ ( √ + √ )2 + 1
√ +∞ 3 3 √ √ √ 30 3 √ √
3 2t 3 3 3 2t 3 3
= arctan + + arctan +
3
√ 0 √ 3 3 3 −∞ 3 3
3 π π 3 π π
= − + − −
3√ 2 6 3 6 2
π 3
= .
3
We conclude that Z √
1 4π 3 log 2
dA = .
x2 + xy + y 2 3
4≤x 2 +y 2 ≤16
3.7.17 Recall from formula 1.14 that the area enclosed by a simple closed curve Γ is given
by Z
1
xdy − ydx.
2 Γ
Using polar coordinates
xdy − ydx = (ρ cos θ)(sin θdρ + ρ cos θdθ) − (ρ sin θ)(cos θdρ − ρ sin θdθ)
= ρ2 dθ.
Parametrise the curve enclosing the region by polar coordinates so that the region is
tangent to the polar axis at the origin. Let the equation of the curve be ρ = f (θ). The
area of the region is then given by
Z π Z π Z π /2
1 1 1
ρ2 dθ = (f (θ))2 dθ = ((f (θ))2 + (f (θ + π/2))2 )dθ.
2 0 2 0 2 0
By the Pythagorean Theorem, the integral above is the integral of the square of the chord
in question. If no two points are farther than 2 units, their squares are no farther than 4
units, and so the area Z π /2
1
< 4dθ = π,
2 0
a contradiction.
3.7.18 Let I(S) denote the integral sought over a region S. Since D(x, y) = 0 inside R,
I(R) = A. Let L be a side of R with length l and let S(L ) be the half strip consisting of
the points of the plane having a point on L as nearest point of R. Set up coordinates uv
so that u is measured parallel to L and v is measured perpendicular to L. Then
Z l Z +∞
I(S(L )) = e−v dudv = l.
0 0
The sum of these integrals over all the vertices of R is 2π. Assembling all these integrals
we deduce the result.
3.8.1 We have
Z Z 1 Z 1−y Z 1−z
zdV = zdxdzdy
E Z0 1 Z0 1−y 0
= z − z 2 dzdy
Z0 1 0
(1 − y)2 (1 − y)3
= − dy
2 3
3 1
0
4
(1 − y) (1 − y)
= −
12 6 0
1
= = .
12
3.8.3 Let A = (1, 1, 1), B = (1, 0, 0), C = (0, 0, 1), and O = (0, 0, 0). We have four
planes passing through each triplet of points:
P1 : A, B, C, x−y+z =1
P2 : A, B, O z=y
P3 : A, C, O x=y
P4 : B, C, O y = 0.
Using the order of integration dzdxdy, z sweeps from P2 to P1 , so the limits are z = y to
z = 1 − x + y. The projection of the solid on the xy plane produces the region bounded
by the lines x = 0, x = 1 and x = y on the first quadrant of the xy-plane. Thus
Z 1 Z x Z 1−x+y Z 1 Z x
dzdydx = (1 − x)dydx
0 0 y Z0 1 0
= x − x2 dx
0 2
x x 3 1
= −
2 3 0
1
= .
6
3.9.1 Cartesian:
Z 1 Z √ 1−y 2 Z √ x 2 +y 2
√ dzdxdy.
−1 − 1−y 2 x 2 +y 2
Cylindrical: Z Z Z
1 2π r
rdzdθdr.
0 0 r2
Spherical:
Z π /2 Z 2π Z (cos φ)/(sin φ)2
r 2 sin φdrdθdφ.
π /4 0 0
π
The volume is .
3
3.9.2
p
1. Since x2 + y 2 ≤ z ≤ 4 − x2 − y 2 , we start our integration with the z-variable.
Observe that if (x, y, z) is on the intersection of the surfaces then
√
2 −1 ± 17
z + z = 4 =⇒ z = .
2
√
2 2 2 17 − 1
Since x + y + z = 4 =⇒ −2 ≤ z ≤ 2, we must have z =
only.
2
The projection of the circle of intersection of the paraboloid and the sphere onto the
xy-plane satisfies the equation
√
17 − 1
z 2 + z = 4 =⇒ x2 + y 2 + (x2 + y 2 )2 = 4 =⇒ x2 + y 2 = ,
2
r√
17 − 1
a circle of radius . The desired integral is thus
2
È√ È√ √
Z 17−1 Z 17−1
−x 2 Z 4−x 2 −y 2
2 2
È√ È√ xdzdydx.
17−1 17−1
− 2
− 2
−x 2 x 2 +y 2
3. Observe that
z = x2 +y 2 =⇒ r cos φ = r 2 (cos θ)2 (sin φ)2 +r 2 (sin θ)2 (sin φ)2 =⇒ r ∈ {0, (csc φ)(cot φ)}.
It is clear that the limits of the angle θ are from θ = 0 to θ = 2π. The angle φ starts
at φ = 0. Now,
√
17 − 1 √
2 17 − 1
z = r cos φ =⇒ cos φ = =⇒ φ = arccos
2 4
a conclusion that is easily reached, since the integrand is an odd function of x and the
domain of integration is symmetric about the origin in x.
3.9.3 Cartesian:
Z √
3 Z √3−y 2 Z √ 4−x 2 −y 2
√ √ dzdxdy.
− 3 − 3−y 2 1
Cylindrical:
Z √
3 Z 2π Z √ 4−r 2
rdzdθdr.
0 0 1
Spherical:
Z π /3 Z 2π Z 2
r 2 sin φdrdθdφ.
0 0 1/ cos φ
5π
The volume is .
3
3.9.5 We have Z Z Z Z
2π 2 2+r cos θ
ydV = r 2 sin θdzdrdθ = 0.
E 0 1 0
π
3.9.7
96
π
3.9.8
14
3.9.9 We put
Now,
dx = cos θ sin φ sin tdr − ρ sin θ sin φ sin tdθ + ρ cos θ cos φ sin tdφ + ρ cos θ sin φ cos tdt
dy = sin θ sin φ sin tdr + ρ cos θ sin φ sin tdθ + ρ sin θ cos φ sin tdφ + ρ sin θ sin φ cos tdt
.
du = cos φ sin tdr − ρ sin φ sin tdφ + ρ cos φ cos tdt
dv = cos tdr − ρ sin tdt
Therefore
Z ZZ Z Z π Z 2π Z π Z 1
2
+y 2 +u 2 +v 2 2
ex dxdydudv = r 3 er sin φ sin2 tdrdφdθdt
0 0 0 0
x 2 +y 2 +u 2 +v 2 ≤1
Z 1 Z 2π Z π Z π
2
= r 3 er dr dθ sin φdφ sin2 tdt
0 π 0 0 0
1
= (2π)(2)
2 2
= π 2.
u = x + y + z =⇒ du = dx + dy + dz,
Thus
0 ≤ uvw ≤ 1, 0 ≤ uv ≤ 1, 0 ≤ u ≤ 1,
which finally give
0 ≤ u ≤ 1, 0 ≤ v ≤ 1, 0 ≤ w ≤ 1.
The integral sought is then, using the fact that for positive integers m, n one has
Z 1
m!n!
xm (1 − x)n dx = ,
0 (m + n + 1)!
we deduce, Z Z Z
1 1 1
u20 v 18 w 8 (1 − u)4 (1 − v)(1 − w)9 dudvdw,
0 0 0
which in turn is
Z 1 Z 1 Z 1
1 1 1
u20 (1 − u)4 du v 18 (1 − v)dv w 8 (1 − w)9 dw = · ·
0 0 0 265650 380 437580
which is
1
= .
44172388260000
dx ∧ dy = du ∧ dv,
dy ∧ dz = −du ∧ dv,
dz ∧ dx = −2vdu ∧ dv,
and so 2 p
d x = 2 + 4v 2 du ∧ dv.
3.10.3 We use spherical coordinates, (x, y, z) = (cos θ sin φ, sin θ sin φ, cos φ). Here θ ∈
[0; 2π] is the latitude and φ ∈ [0; π] is the longitude. Observe that
whence
u v
dx ∧ dy = du ∧ dv, dy ∧ dz = − du ∧ dv, dz ∧ dx = − du ∧ dv,
z z
and so È
2
d x = (dx ∧ dy)2 + (dz ∧ dx)2 + (dy ∧ dz)2
É
u2 + v 2
= 1+ du ∧ dv
√ z2
= 2 du ∧ dv.
Hence
Z Z Z Z √
p √ √ 2π 1
2π 2
z d2 x = u2 + v 2 2 dudv = 2 ρ2 dρdθ = .
0 0 3
Σ u 2 +v 2 ≤1
3.10.5 If the egg has radius R, each slice will have height 2R/n. A slice can be parametrised
by 0 ≤ θ ≤ 2π, φ1 ≤ φ ≤ φ2 , with
R cos φ1 − R cos φ2 = 2R/n.
The area of the part of the surface of the sphere in slice is
Z 2π Z φ 2
R2 sin φdφdθ = 2πR2 (cos φ1 − cos φ2 ) = 4πR2 /n.
0 φ1
This means that each of the n slices has identical area 4πR2 /n.
3.10.6 We project this plane onto the coordinate axes obtaining
Z Z 6 Z 3−z/2
27
xydydz = (3 − y − z/2)ydydz = ,
0 0 4
Σ
Z Z 3 Z 6−2x
27
− x2 dzdx = − x2 dzdx = − ,
0 0 2
Σ
Z Z 3 Z 3−y
27
(x + z)dxdy = (6 − x − 2y)dxdy = ,
0 0 2
Σ
and hence Z
27
xydydz − x2 dzdx + (x + z)dxdy = .
4
Σ
3.11.1 Evaluating this directly would result in evaluating four path integrals, one for each
side of the square. We will use Green’s Theorem. We have
dω = d(x3 y) ∧ dx + d(xy) ∧ dy
= (3x2 ydx + x3 dy) ∧ dx + (ydx + xdy) ∧ dy
= (y − x3 )dx ∧ dy.
The region M is the area enclosed by the square. The integral equals
I Z 2Z 2
x3 ydx + xydy = (y − x3 )dxdy
C 0 0
= −4.
3.11.2 We have
1 4
➊ LAB is y = x; LAC is y = −x, and LB C is clearly y = − x + .
3 3
➋ We have
Z Z 1
5
y 2 dx + xdy = (x2 + x)dx =
6
ZAB Z0 −2
1 4 2 1 15
y 2 dx + xdy = − x+ − x dx = −
ZB C Z1 0 3 3 3 2
2 2 14
y dx + xdy = (x − x)dx =
CA −2 3
Adding these integrals we find
I
y 2 dx + xdy = −2.
△
➌ We have
Z Z 0
Z −x/3+4/3
(1 − 2y)dx ∧ dy = (1 − 2y)dy dx
−2 −x
D
Z 1
Z −x/3+4/3
+ (1 − 2y)dy dx
0 x
44 10
= − −
27 27
= −2.
(To evaluate the integrals you may resort to the fact that the area of the elliptical region
(x − x0 )2 (y − y0 )2
2
+ ≤ 1 is πab).
a b2
If we were to evaluate this integral directly, we would set
√
y = 1 + cos θ, z = 2 sin θ, x = 2 − y = 1 − cos θ.
The integral becomes
Z 2π
√ √
(1 + cos θ)d(1 − cos θ) + 2 sin θd(1 + cos θ) + (1 − cos θ)d( 2 sin θ)
0
which in turn
Z 2π √ √ √
= sin θ + sin θ cos θ − 2+ 2 cos θdθ = −2π 2.
0
Everyone is permitted to copy and distribute verbatim copies of this license document, but changing it is not allowed.
Preamble
The purpose of this License is to make a manual, textbook, or other functional and useful document “free” in the sense of freedom: to assure everyone the
effective freedom to copy and redistribute it, with or without modifying it, either commercially or noncommercially. Secondarily, this License preserves for the
author and publisher a way to get credit for their work, while not being considered responsible for modifications made by others.
This License is a kind of “copyleft”, which means that derivative works of the document must themselves be free in the same sense. It complements the
GNU General Public License, which is a copyleft license designed for free software.
We have designed this License in order to use it for manuals for free software, because free software needs free documentation: a free program should
come with manuals providing the same freedoms that the software does. But this License is not limited to software manuals; it can be used for any textual
work, regardless of subject matter or whether it is published as a printed book. We recommend this License principally for works whose purpose is instruction
or reference.
2. VERBATIM COPYING
You may copy and distribute the Document in any medium, either commercially or noncommercially, provided that this License, the copyright notices, and
the license notice saying this License applies to the Document are reproduced in all copies, and that you add no other conditions whatsoever to those of this
License. You may not use technical measures to obstruct or control the reading or further copying of the copies you make or distribute. However, you may accept
compensation in exchange for copies. If you distribute a large enough number of copies you must also follow the conditions in section 3.
You may also lend copies, under the same conditions stated above, and you may publicly display copies.
3. COPYING IN QUANTITY
250
Appendix A
If you publish printed copies (or copies in media that commonly have printed covers) of the Document, numbering more than 100, and the Document’s
license notice requires Cover Texts, you must enclose the copies in covers that carry, clearly and legibly, all these Cover Texts: Front-Cover Texts on the front
cover, and Back-Cover Texts on the back cover. Both covers must also clearly and legibly identify you as the publisher of these copies. The front cover must
present the full title with all words of the title equally prominent and visible. You may add other material on the covers in addition. Copying with changes limited
to the covers, as long as they preserve the title of the Document and satisfy these conditions, can be treated as verbatim copying in other respects.
If the required texts for either cover are too voluminous to fit legibly, you should put the first ones listed (as many as fit reasonably) on the actual cover,
and continue the rest onto adjacent pages.
If you publish or distribute Opaque copies of the Document numbering more than 100, you must either include a machine-readable Transparent copy
along with each Opaque copy, or state in or with each Opaque copy a computer-network location from which the general network-using public has access to
download using public-standard network protocols a complete Transparent copy of the Document, free of added material. If you use the latter option, you must
take reasonably prudent steps, when you begin distribution of Opaque copies in quantity, to ensure that this Transparent copy will remain thus accessible at the
stated location until at least one year after the last time you distribute an Opaque copy (directly or through your agents or retailers) of that edition to the public.
It is requested, but not required, that you contact the authors of the Document well before redistributing any large number of copies, to give them a chance
to provide you with an updated version of the Document.
4. MODIFICATIONS
You may copy and distribute a Modified Version of the Document under the conditions of sections 2 and 3 above, provided that you release the Modified
Version under precisely this License, with the Modified Version filling the role of the Document, thus licensing distribution and modification of the Modified
Version to whoever possesses a copy of it. In addition, you must do these things in the Modified Version:
A. Use in the Title Page (and on the covers, if any) a title distinct from that of the Document, and from those of previous versions (which should, if there
were any, be listed in the History section of the Document). You may use the same title as a previous version if the original publisher of that version
gives permission.
B. List on the Title Page, as authors, one or more persons or entities responsible for authorship of the modifications in the Modified Version, together
with at least five of the principal authors of the Document (all of its principal authors, if it has fewer than five), unless they release you from this
requirement.
C. State on the Title page the name of the publisher of the Modified Version, as the publisher.
E. Add an appropriate copyright notice for your modifications adjacent to the other copyright notices.
F. Include, immediately after the copyright notices, a license notice giving the public permission to use the Modified Version under the terms of this
License, in the form shown in the Addendum below.
G. Preserve in that license notice the full lists of Invariant Sections and required Cover Texts given in the Document’s license notice.
I. Preserve the section Entitled “History”, Preserve its Title, and add to it an item stating at least the title, year, new authors, and publisher of the Modified
Version as given on the Title Page. If there is no section Entitled “History” in the Document, create one stating the title, year, authors, and publisher of
the Document as given on its Title Page, then add an item describing the Modified Version as stated in the previous sentence.
J. Preserve the network location, if any, given in the Document for public access to a Transparent copy of the Document, and likewise the network
locations given in the Document for previous versions it was based on. These may be placed in the “History” section. You may omit a network location
for a work that was published at least four years before the Document itself, or if the original publisher of the version it refers to gives permission.
K. For any section Entitled “Acknowledgements” or “Dedications”, Preserve the Title of the section, and preserve in the section all the substance and tone
of each of the contributor acknowledgements and/or dedications given therein.
L. Preserve all the Invariant Sections of the Document, unaltered in their text and in their titles. Section numbers or the equivalent are not considered
part of the section titles.
M. Delete any section Entitled “Endorsements”. Such a section may not be included in the Modified Version.
N. Do not retitle any existing section to be Entitled “Endorsements” or to conflict in title with any Invariant Section.
If the Modified Version includes new front-matter sections or appendices that qualify as Secondary Sections and contain no material copied from the
Document, you may at your option designate some or all of these sections as invariant. To do this, add their titles to the list of Invariant Sections in the Modified
Version’s license notice. These titles must be distinct from any other section titles.
You may add a section Entitled “Endorsements”, provided it contains nothing but endorsements of your Modified Version by various parties–for example,
statements of peer review or that the text has been approved by an organization as the authoritative definition of a standard.
You may add a passage of up to five words as a Front-Cover Text, and a passage of up to 25 words as a Back-Cover Text, to the end of the list of Cover
Texts in the Modified Version. Only one passage of Front-Cover Text and one of Back-Cover Text may be added by (or through arrangements made by) any one
entity. If the Document already includes a cover text for the same cover, previously added by you or by arrangement made by the same entity you are acting on
behalf of, you may not add another; but you may replace the old one, on explicit permission from the previous publisher that added the old one.
The author(s) and publisher(s) of the Document do not by this License give permission to use their names for publicity for or to assert or imply endorsement
of any Modified Version.
5. COMBINING DOCUMENTS
You may combine the Document with other documents released under this License, under the terms defined in section 4 above for modified versions,
provided that you include in the combination all of the Invariant Sections of all of the original documents, unmodified, and list them all as Invariant Sections of
your combined work in its license notice, and that you preserve all their Warranty Disclaimers.
The combined work need only contain one copy of this License, and multiple identical Invariant Sections may be replaced with a single copy. If there are
multiple Invariant Sections with the same name but different contents, make the title of each such section unique by adding at the end of it, in parentheses, the
name of the original author or publisher of that section if known, or else a unique number. Make the same adjustment to the section titles in the list of Invariant
Sections in the license notice of the combined work.
In the combination, you must combine any sections Entitled “History” in the various original documents, forming one section Entitled “History”; likewise
combine any sections Entitled “Acknowledgements”, and any sections Entitled “Dedications”. You must delete all sections Entitled “Endorsements”.
6. COLLECTIONS OF DOCUMENTS
You may make a collection consisting of the Document and other documents released under this License, and replace the individual copies of this License
in the various documents with a single copy that is included in the collection, provided that you follow the rules of this License for verbatim copying of each of
the documents in all other respects.
You may extract a single document from such a collection, and distribute it individually under this License, provided you insert a copy of this License into
the extracted document, and follow this License in all other respects regarding verbatim copying of that document.
8. TRANSLATION
Translation is considered a kind of modification, so you may distribute translations of the Document under the terms of section 4. Replacing Invariant
Sections with translations requires special permission from their copyright holders, but you may include translations of some or all Invariant Sections in addition
to the original versions of these Invariant Sections. You may include a translation of this License, and all the license notices in the Document, and any Warranty
Disclaimers, provided that you also include the original English version of this License and the original versions of those notices and disclaimers. In case of a
disagreement between the translation and the original version of this License or a notice or disclaimer, the original version will prevail.
If a section in the Document is Entitled “Acknowledgements”, “Dedications”, or “History”, the requirement (section 4) to Preserve its Title (section 1) will
typically require changing the actual title.
9. TERMINATION
You may not copy, modify, sublicense, or distribute the Document except as expressly provided for under this License. Any other attempt to copy, modify,
sublicense or distribute the Document is void, and will automatically terminate your rights under this License. However, parties who have received copies, or
rights, from you under this License will not have their licenses terminated so long as such parties remain in full compliance.
Each version of the License is given a distinguishing version number. If the Document specifies that a particular numbered version of this License “or any
later version” applies to it, you have the option of following the terms and conditions either of that specified version or of any later version that has been published
(not as a draft) by the Free Software Foundation. If the Document does not specify a version number of this License, you may choose any version ever published
253