Vector-Calculus Class Note Vector-Calculus Class Note
Vector-Calculus Class Note Vector-Calculus Class Note
Vector-Calculus Class Note Vector-Calculus Class Note
Posenjit Das
Department of Mathematics,
Indian Institute of Space Science and Technology,
Valiamala P.O., Trivandrum 695 547, India
email: [email protected]
ii
Contents
4 Line Integration 27
4.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
4.2 Line Integration of Scalar Fields . . . . . . . . . . . . . . . . . . 27
4.3 Line Integration of Vector Fields . . . . . . . . . . . . . . . . . . 28
4.4 Conservative Vector Fields . . . . . . . . . . . . . . . . . . . . . . 30
4.5 1st Structure Theorem of CVFs . . . . . . . . . . . . . . . . . . . 32
4.6 Fundamental Theorem in R2 : Green’s Theorem . . . . . . . . . . 35
4.6.1 Application: 2nd Strure Theorem of CVFs in R2 . . . . . 36
4.6.2 Application: Extended Green’s Theorem . . . . . . . . . . 36
4.6.3 Application: Recognizing CVF in R2 . . . . . . . . . . . . 38
4.6.4 Application: Finding Area enclosed by a loop . . . . . . . 38
5 Surface Integration 41
iii
iv CONTENTS
Chapter 1
1.1 Introduction
Though the need of the conception of vectors was realized after the discovery of
physical quantities carrying multiple values together, e.g., velocity, force, e.t.c.,
significant development in the systematic study of vector is done in the recent
era, in the mid of 19th century. The primary development in the theory of
vectors is witnessed in the late 17th century, while studying the of structure
of complex numbers. By the passage of time the theory got enriched, as the
need grew, to understand the geometrical objects and physical phenomenons.
The modern theory vectors which is a quite sophisticated in nature is developed
between late 19th century and early 20th century.
In the early development, vectors were defined to be finite straight line seg-
ments with unique arrow heads, where the length of the straight line segments
represented the “magnitude” of the vectors and the arrow determined the “di-
rection” of the vector; but it failed to explain many properties of vectors which
should be satisfied naturally. Later, a corrected model of vectors was formalized
by defining vectors as pairs of points in the corresponding space; and subse-
quently more refined and formal definition represents vectors as points in the
ambient space.
(αx, αy), i.e., coordinate wise multiplication and by definition it is obvious that
“.” satisfies the properties (6)-(10) of Definition 1.2.1. This shows that V is a
space of vectors over R.
(II) In a similar fashion Rn , n ≥ 1 forms space of vectors over R with respect to
coordinate wise addition and multiplication.
(III) Suppose (x0 , y0 ) ∈ R2 . Let L(x0 ,y0 ) = {(λx0 , λy0 ) | λ ∈ R}. One should be
able to see that L(x0 ,y0 ) is a space of vectors over R. Note that the space of
vectors L(x0 ,y0 ) is a subset of the space of vectors R2 .
(IV) Let V = C[0, 1] be the set of all real valued continuous functions over [0, 1].
For f, g ∈ V, define f + g by (f + g)(x) = f(x) + g(x) for all x ∈ [0, 1], i.e.,
point wise addition; and for α ∈ R and f ∈ V define α.f by (α.f)(x) = αf(x) for
all x ∈ [0, 1]. It is easy to see that with respect these two operations V forms
a space of vectors over R. The additive identity of this space is the constant
function zero, i.e., the function 0(x) = 0 for all x ∈ [0, 1]; and for each element
f ∈ V, the additive inverse of f is the function −f given by (−f)(x) = −f(x) for
all x ∈ [0, 1].
(V) Let V = Mm×n (R) be the set of all m×n matrices with real entries. Let “+”
denote the usual matrix addition operation, i.e., for A = (aij ) and B = (bij ) in
V, A + B is defined by A + B = (aij + bij ). Define an operation . : R × V −→ V
by α.(aij ) = (αaij ) where α ∈ R and A = (aij ) ∈ V. It is easy to check that
with respect to these two operations V forms a space of vectors over R.
Definition 1.2.4. Let V be a space of vectors over R and S := {v̄1 , v̄2 , · · · , v̄n }
be a collection of vectors from V. By a linear combination of the vectors in S
we mean an element of the form α1 v̄1 + α2 v̄2 + · · · + αn v̄n where αi ∈ R are
constants.
Let < S >:= {α1 v̄1 + α2 v̄2 + · · · + αn v̄n | αi ∈ R} be the collection of all
possible linear combinations of vectors from S. It can be proved that < S > forms
a vector space “being in V”. < S > is called the linear span/hull of vectors in
S. Sometimes < S > is denoted by L(S).
Example 1.2.5. 1. Let v̄ ∈ R2 . Then < {v̄} > in R2 is a line passing though
origin (0, 0) and v̄.
2. Let v̄1 , v̄2 ∈ R3 , then < {v̄1 , v̄2 } > is either a straight line passing though
origin or a plane containing (0, 0), v̄1 and v̄2 in R3 .
3. Let v̄1 , v̄2 ∈ R2 , then < {v̄1 , v̄2 } > is either a straight line passing though
origin or the whole plane R2 .
4. Let e1 = (1, 0) and e2 = (0, 1). It can be checked that < {e1 , e2 } >= R2 .
5. Let e1 = (1, 0, 0), e2 = (0, 1, 0) and e3 = (0, 0, 1). It can be checked that
< {e1 , e2 , e3 } >= R3 .
Definition 1.2.6. A subset S of a vector space V over R is called a basis, if S
generates V, i.e., < S >= V; and any proper subset of S can not generate V
Example 1.2.7. 1. {e1 , e2 } is a basis of R2 .
2. {(1, 1), (2, 2)} is not a basis of R2 .
For any vector v̄ in V, kv̄k is called the norm or length of v̄. A space of
vectors endowed with a norm function is called a normed space. Geometrically
norm of a vector v̄ represents the length of line segment joining 0̄ and v̄.
Example 1.3.2.
(I) Let V = Rm where n is a natural number. We already have seen that V
is
q a vector space over R. For any v̄ = (v1 , v2 , · · · , vm ) in V define kv̄k =
v21 + v22 + · · · v2m . Then one can check that the function k.k : V −→ R is a
norm on V. This norm is called the Euclidean norm on Rm .
(II) Let V = C[0, 1]. Then V is a vector space over R. Define a function k.k :
V −→ R by kfk = sup {|f(x)| | x ∈ [0, 1]}. It can be seen that the function k.k is
a norm on V.
(III) For the same vector space as in the previous example, for any f ∈ C[0, 1]
Z 1 !1/2
define kfk = f2 . One can show that the defined function k.k : V −→ R
0
is a norm called the L2 -norm..
(IV) Let V = Mm×n (R). It has already been established that V is a vector space
over R. Define k.k : V −→ R by kAk = sup {kAx̄T k | x̄ ∈ Rn , kx̄k = 1} where A
is in V. The reader may take it as an exercise to show that the function k.k is
a norm on V.
d(f, g) = sup {|(f − g)|(x) | x ∈ [0, 1]} for all f, g ∈ C[0, 1].
This metric is called the supremum metric or the uniform metric on C[0, 1].
In high school Mathematics we have learned the measurement of angle be-
tween two vectors in R2 is done by dot product of vectors. For general space
of vectors we introduce a notion called inner product of vectors which exactly
catches the notion of angle between two vectors.
Definition 1.3.6. Let V be a space of vectors over R. A function h. , .i :
V × V −→ R is called an inner product if it satisfies the following properties:
IP-1 (Bilinearity): hα1 v̄1 + α2 v̄2 , v̄3 i = α1 hv̄1 , v̄3 i + α2 hv̄2 , v̄3 i for all α1 , α2 in
R and v̄1 , v̄2 , v̄3 in V.
IP-2 (Symmetricity): hv̄1 , v̄2 i = hv̄2 , v̄1 i for all v̄1 , v̄2 in V.
IP-3 (Positive-definiteness): hv̄, v̄i ≥ 0 for all v̄ in V; and hv̄, v̄i = 0 if and only
if v̄ = 0̄.
Example 1.3.7.
(I) Define h. , .i : R2 × R2 −→ R by h(x1 , y1 ) , (x2 , y2 )i = x1 x2 + y1 y2 for
all (x1 , y1 ), (x2 , y2 ) in R2 . This function is the well known dot product on R2 .
One can easily verify that the function h. , .i is an inner product on R2 .
(II) In more generality, the function h. , .i : Rm × Rm −→ R defined by hx̄ , ȳi =
x1 y1 + x2 y2 + · · · + xm ym for all x̄ = (x1 , x2 , · · · , xm ), ȳ = (y1 , y2 , · · · , ym ) in
Rm can be verified to be an inner product on Rm , which called the Euclidean
inner product.
(III) On the space of vectors V = C[0, 1], define h. , .i : V × V −→ R by hf , gi =
Z1
fg for all f, g in V. It is east to check that the function h. , .i is an inner
0
product on V, and is called the L2 inner product.
Note that thep property IP-1pof Definition 1.3.6 implies that hαf , αfi =
α2 hf , fi, i.e., hαf , αfi = |α| hf , fi for all f in V and for all α in R. Thus
looking at IP-1 and IP-3 of Definition 1.3.6 one may naturally ask whether inner
product leads to giving a norm in the corresponding space. We shall see that
the answer, in fact, is affirmative. To see this we shall observe some properties
of inner product.
Now we are ready to show that an inner product on a space of vectors defines
a norm on the corresponding space.
Proof. By the definition we see that kfk ≥ 0 for all f in V, and kαfk = |α|kfk
for all f in V. By corollary 1.3.9 we see that kf + gk ≤ kfk + kgk for all f, g in
V. This completes the proof.
Example 1.3.11.
(I) Let Rm be endowed with Euclidean inner product. Then it can be seen easily
that the corresponding norm is the Euclidean norm on Rm .
(II) The L2 inner product in C[0, 1] induces the L2 norm on the corresponding
space.
Example 1.3.12.
(I) Let x̄, ȳ ∈ R3 be such that the line segments 0̄, x̄ and 0̄, ȳ are perpendicular
to each other. Then one can check that the Euclidean inner product of x̄ and ȳ
is zero.
(II) Let x̄, ȳ ∈ R3 be such that the line segments 0̄, x̄ and 0̄, ȳ are parallel to
each other. Then one can check that the absolute value of the Euclidean inner
product of x̄ and ȳ is product of the norm of x̄ and ȳ.
In view of the above examples, for a general space of vectors we define the
following:
Definition 1.3.13. Let V be a space of vectors with an inner product h. , .i.
Two vectors v̄1 and v̄2 in V are called
(I) Orthogonal, if hv̄1 , v̄2 i = 0.
(II) Parallel, if |hv̄1 , v̄2 i| = kv̄1 kkv̄2 k.
In this chapter so far we have discussed on general space of vectors. But in the
rest of topics we shall mainly restrict our discussion in the space of vectors Rm ,
specially R2 and R3 , with Euclidean inner product. In the next section we shall
continue to explore the spatial structure of Rm .
Example 1.4.2.
(A) In R, the open interval (0, 1) is the open 1-ball with center 12 and radius 12 ;
i.e., (0, 1) = B( 12 , 21 ). In general, for a, b ∈ R where a 6= b we have (a, b) =
B( a+b a−b a+b a−b
2 , 2 ) and [a, b] = B( 2 , 2 ).
The unit open 1-ball, i.e, the unit open ball in R, is given by B(0) = (−1, 1);
and similarly B(0) = [−1, 1].
(B) The unit closed 2-ball is the circular disk { (x, y) ∈ R2 | x2 + y2 ≤ 1 } where
0̄ = (0, 0) ∈ R2 , the and unit open 2-ball B(0̄) is the “boundary-less” circular
disk { (x, y) ∈ R2 | x2 + y2 < 1 } where 0̄ = (0, 0) ∈ R2 .
(C) The closed 3-ball with center (x0 , y0 , z0 ) and radius r is the solid sphere
{ (x, y, z) ∈ R3 | (x − x0 )2 + (y − y0 )2 + (z − z0 )2 ≤ r2 }. The open 3-ball
B((x0 , y0 , z0 ), r) is the “boundary-less” solid sphere { (x, y, z) ∈ R3 | (x − x0 )2 +
(y − y0 )2 + (z − z0 )2 < r2 }.
Example 1.4.4.
(A) Any open ball B(x̄, r) in Rm is open, as for each ȳ ∈ B(x̄, r), we have
B(ȳ, r0 ) ⊆ B(x̄, r) where r0 = min {d(x̄, ȳ), r − d(x̄, ȳ)}.
(B) Any closed interval [a, b] is not open, for the points a and b in [a, b] there
exists no open 1-ball, i.e., open intervals, centered at a and centered at b, which
lies in [a, b]. Note that for a = b, we have [a, b] = {a}.
(C) The closed interval [a, b] is a closed set. To see this we need to observe
that R\[a, b] = (−∞, a) ∪ (b, ∞) is open. Suppose x ∈ R\[a, b]. Let r0 =
min {d(x, a), d(x, b)}. Then it is easy to show that B(x, r0 ) ⊆ R\[a, b]. Since
x ∈ R\[a, b] was arbitrary, we get that R\[a, b] is open.
Proposition 1.4.5.
1. Finite intersection of open sets is open.
2. Arbitrary union of open sets is open.
Proof. Exercise
Proposition 1.4.6.
1. Finite union of closed sets is closed.
2. Arbitrary intersection of closed sets is closed.
Example 1.4.8. Arbitrary intersection of open sets may not be open. Set
m
\1/n) ⊂ R . Then each Um is an open n-ball and hence an open
Un := B(0̄,
set. But Un = {0̄}, which is closed.
n≥1
Another approach would be to see that any kind of translation of the line
Lv̄ (0̄) look like P∗ + Lv̄ ((0, 0)) for some P∗ ∈ Rm , and which is a straight line.
Note that this straight line is not a vector space as it does not contain 0̄. It is
easy to see that this straight line is passing though P∗ and is parallel to v̄, i.e.,
to the straight line Lv̄ (0̄); and therefore the straight line passing though P0 and
parallel to v̄ will be P0 + Lv̄ ((0, 0)).
As Lv̄ (P0 ) = P0 + Lv̄ (0̄) = {P0 + αv̄ ∈ Rm | α ∈ R}, the parametric rep-
resentation of the line Lv̄ (P0 ), when m = 3, i.e., in R3 would be (x, y, z) =
(x0 + αv1 , y0 + αv2 , z0 + αv3 ) where P0 = (x0 , y0 , z0 ) and v̄ = (v1 , v2 , v3 ).
L(P1 , P2 ) = P1 + LP2 −P1 (0̄) = LP2 −P1 (P1 ) = {P1 + α(P2 − P1 ) | α ∈ R}.
Chapter 2
2.1 Introduction
Let G be an open set in R and f : G −→ R a real valued function of real variable.
f(x0 + h) − f(x0 )
We call f to be differentiable at x0 ∈ G if lim exists, and the
h→0 h
d ′ ′
limit is termed as dx (f(x))|x0 or f (x0 ). f (x0 ) represents the rate of change of
f along x-axis at the point x0 . We already have learned that if f is differentiable
at x0 , it also continuous at x0 .
11
One should be careful that the mere existence of the partial deriva-
tives at P0 does not ensure differentiability of the given function. Fur-
ther, if the partial derivatives of the function at the point P0 do not
exists, the function is not differentiable at P0 . One may check that
that if f is differentiable at P0 , the it is continuous at P0 (Prove it!).
x + y if x = 0 or y = 0
Exercise 2.1.1. 1. Let f(x, y) = . Show that f
1 otherwise
has partial derivatives at (0, 0), but f is not continuous at (0, 0).
2. Show that the function f(x, y) = |x| + |y| for all (x, y) ∈ R2 is continuous,
but partial derivatives pf f at (0, 0) does not exists, and therefore f is not
differentiable at (0, 0).
2. Under what condition the partial differentiability can replace the concept
of differentiability of a scalar field?
In this chapter we mainly try to get the answer to the above questions.
∂ f(x0 + h, y0 , z0 ) − f(x0 , y0 , z0 )
∂x (f)|P0 = lim
h→0 h
f((x0 , y0 , z0 ) + t(1, 0, 0)) − f((x0 , y0 , z0 ))
= lim
t→0 t
f(P0 + te1 ) − f(P0 )
= lim
t→0 t
Imitating this we define the rate of change of f along any line passing though
P0 and parallel to a unit vector v̄ (i.e., ||v̄|| = 1), i.e., along the line Lv̄ (P0 ) as
f(P0 + tv̄) − f(P0 )
lim , if the limit exists. We define this quantity the direc-
t→0 t
tional derivative of f at P0 along the unit vector v̄ and is denoted by
f(P0 + tv̄) − f(P0 )
Dv̄ (f)|P0 , i.e., Dv̄ (f)|P0 = lim , if the limit exists.
t→0 t
d
Exercise 2.2.1. 1. Show that Dv̄ (f)|P0 = dt (f(P0 + tv̄))|t=0 .
2. For any direction v̄, find the directional derivative Dv̄ (f)|(0,0) where f(x, y) =
x2 + y2 .
x + y if x = 0 or y = 0
3. Let f(x, y) = . Find directions v̄ such that
1 otherwise
Dv̄ (f)|(0,0) exists.
p
4. Let f(x, y) = x2 + y2 for all (x, y) ∈ R2 . Show that Dv̄ (f)|(x0 ,y0 ) exists
for all direction v̄ if (x0 , y0 ) 6= (0, 0).
It is obvious that if DP0 (f)|v̄ exists for all direction v̄ (unit vector), then all
∂
the partial derivatives of f exists at P0 , as De1 (f)|P0 = ∂x (f)|P0 , De2 (f)|P0 =
∂ ∂
∂y (f)| P0
, and D e3
(f)| P0
= ∂z (f)| P0
. So, we ask
Problem 2.2.2. If all the partial derivatives of f exists at P0 , does Dv̄ (f)|P0
exit for all direction v̄?
This question comes naturally as we know that any vector v̄ can be written
as linear combinations of e1 , e2 and e3 in R3 , i.e., in R3 we can reach any point
v̄ from 0̄ by moving along x-axis, y-axis and z-axis.
g(x) − g(a)
lim exists. If this limit exists we call this value to be g ′ (a),
x→a x−a
g(x) − g(a)
i.e., g ′ (a) = lim . Now define a function g1 : G −→ R by
x→a x−a
f(x)−f(a)
if x 6= a
g1 (x) = x−a
α if x = a
Note that if we set α = g ′ (a), then g1 is continuous at x = a, as g is
differentiable at x = a. Now suppose it is not known whether g is differentiable
at a. If we assume that g1 is continuous at a, then it is easy to see that g is
differentiable at a, and α = f ′ (a). Thus, we have
f(x) − f(a) d
Specifically, f1 (x) = for x ∈ G\{a}; and f1 (a) = dx (f)|a = f ′ (a).
x−a
It can be seen that the above result also holds true for functions of several
variables.
∂ ∂
Further, it can be seen that f1 (ā) = ∂x (f)|ā , f2 (ā) = ∂y (f)|ā , and, f3 (ā) =
∂
∂z (f)|ā .
Proof. Assignment!
Now, we assume that f : G −→ R be differentiable at P0 where G ⊂ R3
is open and P0 ∈ G. Then by Theorem 2.2.5 there exists functions f1 , f2 , f3 :
G −→ R continuous at P0 such that
f(x, y, z)−f(x0 , y0 , z0 ) = (x−x0 )f1 (x, y, z)+(y−y0 )f2 (x, y, z)+(z−z0 )f3 (x, y, z)
for all (x, y, z) ∈ G. Now, let v̄ be any unit vector in R3 and we choose the
points (x, y, z) ∈ G such that they lie on the straight line P0 + tv̄ where t ∈ R.
Then we get
f(P0 +tv̄)−f(P0 ) = (x0 +tv1 −x0 )f1 (P0 +tv̄)+(y0 +tv2 −y0 )f2 (P0 +tv̄)+(z0 +tv3 −z0 )f3 (P0 +tv̄)
f(P0 + tv̄) − f(P0 ) (x0 + tv1 − x0 )f1 (P0 + tv̄) (y0 + tv2 − y0 )f2 (P0 + tv̄) (z0 + tv3 − z0 )f3 (P0 + tv̄)
= + +
t t t t
for all t ∈ R\{0}.
∂ ∂ ∂
Dv̄ (f)|P0 = v1 (f)|P0 + v2 (f)|P0 + v3 (f)|P0 =< ∇f|P0 , v̄ > .
∂x ∂y ∂z
From this we conclude the following
Theorem 2.2.6. Let G ⊂ R3 be open, P0 ∈ G and f : G −→ R differentiable at
P0 , then
(I) Dv̄ (f)|P0 exists for all direction v̄, and
(II) Dv̄ (f)|P0 =< ∇f|P0 , v̄ > .
Now, if we follow the statements of Theorem 2.2.3 and Theorem 2.2.6, we
shall see that two different “hypotheses” gives the same result. So, naturally
one asks whether those two “hypotheses” in Theorem 2.2.3 and Theorem 2.2.6
have any relation. Clearly, mere differentiability of a function can not ensure
continuity of partial derivatives or even the existence of partial derivatives in
neighbourhood points. So, we should check whether the “conditions” of partial
differentiations in Theorem 2.2.3 ensure differentiability of the function in the
corresponding point. The answer is, in fact, “Yes”! From the technique applied
in Theorem 2.2.3, and the technique of forming the increment function f1 in
Theorem 2.2.4 we see a scope/ we sense our intuition of proving differentiability
of a scalar field under the condition the partial derivatives are continuous.
Let G ⊂ Rn be open and P0 = (x0 , y0 , z0 ) ∈ G. Suppose f : G −→ R a
scalar field. Then,
f(x, y, z) − f(x0 , y0 , z0 ) = f(x, y, z) − f(x, y, z0 )+
f(x, y, z0 ) − f(x, y0 , z0 )+
f(x, y0 , z0 ) − f(x, y, z)
Therefore, if we define
f(x,y,z)−f(x,y,z0 ) f(x,y,z0 )−f(x,y0 ,z0 )
x−x0 if z 6
= z0 y−y0 if y 6= y0
f1 (x, y, z) = ∂ f2 (x, y, z) = ∂
∂z (f)| (x,y,z0 ) if z = z0 ∂y (f)| (x,y ,z ) if y = y0
0 0
∂ ∂
If we further assume that ∂z (f)|(x,y,z) , and ∂y (f)|(x,y,z) are continuous at
P0 , then we see that the functions f1 , f2 and f3 are continuous at P0 . And, it
is easy to check that
f(x, y, z)−f(x0 , y0 , z0 ) = (z−z0 )f1 (x, y, z)+(y−y0 )f2 (x, y, z)+(x−x0 )f1 (x, y, z)
Thus, in view of Theorem 2.2.7 and Theorem 2.2.6, we see that Theorem
2.2.3 follows as a corollary; and combining those two results we can write the
following for G ⊂ R3 open, P0 ∈ G and f : G −→ R:
Remark 2.2.8.
A1 : all the partial derivatives of f exist in an open ball araound P0
A:
A2 : two of the three partail derivatibes of fare continuous at P0
⇓
B: f is differentiable at P0
⇓
C1 : Dv̄ (f)|P0 exists for all direction v̄
C:
C2 : Dv̄ (f)|P0 =< ∇f|P0 , v̄ >
4. Let f(x, y) = |x| + |y| for all (x, y) ∈ R2 . It can be seen that Dv̄ (f)|(0,0)
does not exist for any direction v̄, and therefore f is not differentiable at
(0, 0).
2 2
x y
if (x, y) 6= (0, 0)
Exercise 2.2.10. 1. f(x, y) = x4 +y2 . Show that f is
0 if (x, y) = (0, 0)
differentiable at (0, 0), though fy is not continuous at (0, 0).
2. Let f(x, y) = |xy| for all (x, y ∈ R2 ). Show that f is differentiable at (0, 0).
p
3. let f(x, y) = |xy| for all (x, y) ∈ R2 . Show that f is not differentiable at
(0, 0).
x3 y
if (x, y) 6= (0, 0)
4. f(x, y) = x4 +y2 . Show that f fails to be differen-
0 if (x, y) = (0, 0)
tiable at (0, 0), though Dv̄ (f)|(0,0) =< ∇(f)|(0,0) , v̄ > for any direction
v̄.
x2 y
if (x, y) 6= (0, 0)
5. f(x, y) = x2 +y2 . Show that f is continuous at (0, 0);
0 if (x, y) = (0, 0)
Dv̄ (f)|P0 exists for all direction v̄; but f is not differentiable at (0, 0).
Dv̄ (f)|P0 =< ∇(f)|P0 , v̄ >= ||∇(f)|P0 || ||v̄|| cos(θ) = ||∇(f)|P0 || cos(θ).
Where θ is the angle inscribed by ∇(f)|P0 and v̄ with 0̄. Then we have
∇(f)|P0
1. Dv̄ (f)|P0 is maximum if cos(θ) = 1, i.e., if v̄ = ||∇(f)|P0 || .
∇(f)|
2. Dv̄ (f)|P0 is mimimum if cos(θ) = −1, i.e., v̄ = − ||∇(f)|PP0 || .
0
∂ ∂
3. Dv̄ (f)|P0 = 0 if cos(θ) = π/2, i.e., v̄ = ±( ∂y (f)|P0 , − ∂x (f)|P0 ).
Example 2.3.1. 1. Let f(x, y) = x2 + y2 for all (x, y) ∈ R2 . Then for each
c > 0, the level curve Sc (f) is the circle x2 + y2 = c.
2. Let f(x, y, z) = x2 + y2 + z2 for all (x, y, z) ∈ R3 . Then for each c > 0,
the level surface Sc (f) is the sphere x2 + y2 + z2 = c.
Let n = 3, and suppose that the scalar field f is smooth. Consider the
level surface Sc (f). Since Sc (f) is a surface, it is union of curves. Let γ(t) =
(x(t), y(t), z(t)) for all t ∈ I ⊂ R be any smooth curve on Sc (f) passing though
the point P0 . Then f(γ(t)) = c. From which, on differentiating with respect to
t, we get
∂ d ∂ d ∂ d
(f)|P0 (x)|t0 + (f)|P0 (y)|t0 + (f)|P0 (z)|t0 = 0 , i.e.,
∂x dt ∂y dt ∂z dt
This shows that ∇(f)|P0 is normal to the curve γ at the point P0 . Since
P0 is arbitrary, we see that at any point, the gradient of a scalar field is
normal to its corresponding level set at that point, i.e., ∇(f)|P0 ⊥ Sc (f)
at P0 .
Exercise 2.3.2. 1. Draw the level surfaces of f(x, y, z) = x2 + y2 − z2 for
all (x, y, z) ∈ R3 .
2. Draw the level curves of f(x, y) = x2 − y2 for all (x, y) ∈ R2 .
2 2
3. Let h(x, y) = 2e−x + e−3y denote the height on a mountain at position
(x, y). In what direction from (1, 0) should one begin walking in order to
climb the fastest?
4. Find a unit normal vector to √ the√ following surfaces at the specified point
(i) x2 + y2 + z2 = 9 at (0, 3, 3), (ii) x3 y3 + y − z + 2 = 0 at (0, 0, 2),
(iii) z = 1/(x2 + y2 ) at (1, 1, 1/2).
5. Suppose that a √ particle is ejected from the surface x2 + y2 − z2 = −1 at
the point (1, 1, 3) along the normal directed toward the xy plane to the
surface at time t = 0 with a speed of 10 units per second. When and where
does it cross the xy plane?
6. Starting from (1, 1), in which direction should one travel in order to obtain
the most rapid rate of decrease of the function f : R2 −→ R defined by
f(x, y) := (x + y − 2)2 + (3x − y − 6)2 ?
p
7. About how much will the function f(x, y) := ln x2 + y2 change if the
point (x, y) is moved from (3, 4) a distance 0.1 unit straight toward (3, 6)?
Chapter 3
3.1 Introduction
If L is a straight line segment (finite) or a finite union of straight line segments in
R3 , it is well known to us how to measure the length of the segment or totality
of of the length of the segments. In this chapter we shall drive a tool to measure
the length of a segment of a parametric curve.
21
Now we shall show that a parametric curve has infinitely many parametrisa-
tion. Suppose [a, b] be an intervals in R. Define a maps Φ[a,b] : [a, b] −→ [0, 1],
and Ψ[a,b] : [0, 1] −→ [a, b] by
b−t
Φ[a,b] (t) = for all t ∈ [0, 1] (this function scales the interval [a, b] onto
b−a
[0, 1])
and
Ψ[a,b] (s) = bs + a(1 − s) for all s ∈ [0, 1] (this function scales the interval [0, 1]
onto [a, b]).
Note that Φ[a,b] (a) = 0, Φ[a,b] (b) = 1, Ψ[a,b] (0) = a, Ψ[a,b] (1) = b,
Φ[a,b] ◦ Ψ[a,b] (s) = s for all s ∈ [0, 1], Ψ[a,b] ◦ Φ[a,b] (t) = t for all t ∈ [a, b],
(i.e., Φ[a,b] = Ψ−1 −1
[a,b] , and Ψ[a,b] = Φ[a,b] ) and both the functions are infinitely
many differentiable . Now, for given any pair of interval [a, b] and [c, d], we
[a,b]
shall define a map T[c,d] : [c, d] −→ [a, b], using these properties of Φ and Ψ as
follows
Φ[c,d] Ψ[a,b]
[c, d] −→−→−→ [0, 1] −→−→−→ [a, b], i.e.,
[a,b]
T[c,d] (t) = Ψ[a,b] ◦ Φ[c,d] (t) for all t ∈ [c, d]. Now for any parametric
[a,b]
curve γ : [a, b] −→ R3 , consider the map γ ◦ T[c,d] : [c, d] −→ R3 . Check
[a,b] [a,b]
that {γ} = {γ ◦ T[c,d] }, and therefore both γ and γ ◦ T[c,d] represent the same
parametric curve. If we assume that {γ} = C, then the above discussion shows
that C has infinitely many parametrisation. The developed theory will show
that it is enough to consider one favorable parametrisation of a parametric
curve, and study the topics under discussion.
If lim ℓP (γ) exists we say that length of γ exists, and the value of the limit is
n→∞
defined as the length of γ, and is denoted by ℓ(γ); i.e., ℓ(γ) = lim ℓP (γ), if the
n→∞
limit exists.
If we let ||P|| := max |tn − tn−1 |, we see that lim ℓP (γ) = lim ℓP (γ).
n n→∞ ||P||→0
Let Peq ([a, b]) = {P ∈ P([a, b]) | points in P are equidistant}. As can be
seen that lim ℓP (γ) = lim {ℓP (γ) | P ∈ Peq ([a, b])}, for the context under
||P||→0 ||P||→0
discussion it is enough to take partitions of [a, b] from Peq ([a, b]). Note that
any partition in Peq ([a, b]) is uniquely determined by the number of points in
it; and therefore, instead of notation ℓP (γ) it is enough to use the notation
ℓn (γ) to indicate that we have taken a partition from Peq ([a, b]) having n + 1
equidistant points. Thus, we shall write ℓ(γ) = lim ℓn (γ), if the limit exists.
n→∞
Let P ∈ Peq ([a, b]) with n + 1 points. Assume that ||P|| = ∆t. Clearly, it
implies that for any two consecutive points ti and ti+1 from P we have ti+1 =
ti + ∆t. We see that lim ℓP (γ) = lim ℓn (γ) = lim ℓn (γ). If we write
n→∞ ||P||→0 ∆t→0
γ(t) = (x(t), y(t), z(t)), then we see that
n
X q
lim ℓn (γ) = lim (x(ti ) − x(ti−1 ))2 + (y(ti ) − y(ti−1 ))2 + (z(ti ) − z(ti−1 ))2
∆t→0 ∆t→0
i=1
Xn p
(x(ti ) − x(ti−1 ))2 + (y(ti ) − y(ti−1 ))2 + (z(ti ) − z(ti−1 ))2
= lim ∆t
∆t→0 ∆t
i=1
Xn
r
(x(ti ) − x(ti−1 ))2 + (y(ti ) − y(ti−1 ))2 + (z(ti ) − z(ti−1 ))2
= lim ∆t
∆t→0 ∆t2
i=1 s
Xn 2 2 2
x(ti ) − x(ti−1 ) y(ti ) − y(ti−1 ) z(ti ) − z(ti−1 )
= lim + + ∆t
∆t→0 ∆t ∆t ∆t
i=1 s
Xn 2 2 2
d d d
= lim (x(t))|ti−1 + (y(t))|ti−1 + (z(t))|ti−1 ∆t
∆t→0 dt dt dt
i=1
if γ is differentiable, i.e., if γ ′ (t) exists
Xn
= lim ||γ ′ (t)||ti−1 ∆t
∆t→0
Zb i=1
= ||γ ′ (t)||dt
a
if γ ′ is continuous, i.e., if γ is C1 -type
2.
defines a map s : [a, b] −→ R by s(x) = ||γ ′ ||. This function is called the arc
a
length function of the curve γ. We shall study the properties of this function
after we observe few examples.
Example 3.4.1. Let γ(t) = (cos(2πt), sin(2πt)) for all t ∈ [0, 1]. Clearly γ
is a C1 -type curve, and ||γ ′ (t)|| = ||(−2π sin(2πt), 2π cos(2πt))|| = 2π for all
t ∈ [0, 1]. Therefore, the arc length function of γ is given by
= 2πx
We know that the curve {γ} has another parametrisation δ : [0, 2π] −→ R3
defined by δ(t) = (cos(t), sin(t)) for all t ∈ [0, 2π]. For this parametrisation it
can be seen that s(x) = x for all x ∈ [0, 2π].
Remark 3.4.2. The above example clearly shows that the different parametri-
sation to the same curve gives different arc length function, i.e., So, we ask
for a given parametric curve C whether there is any parametrisation η defined
on an interval I for which s(x) = x for all x ∈ I. Note that in that case, for
that special paramerisation,
Z say η, we shall get ||η ′ (t)||
Z = 1 for all t ∈ I. Sup-
x
pose I = [c, d]. Since ||η ′ (t)||dt = ℓ(C), and s(x) = ||η ′ (t)||dt it is clear that
I c
c = 0 and d = ℓ(C). Since the arc length function s is a surjective function from
the domain of the parametrised curve C to the interval [0, ℓ(C)], we naturally
ask whether the function s−1 can be used to convert domain of the parametric
curve to [0, ℓ(C)]. The following properties will help us to get the answer to this
question.
Before going into the properties let us recall a fact: “Any real valued
bijective continuous function defined on a closed interval has contin-
uous inverse. Further, if the function is smooth, the inverse is also
smooth.” (Prove it!!!)
Theorem 3.4.3. Let s : [a, b] −→ R be the arc length function of the C1 -type
curve γ : [a, b] −→ R3 . Then the following holds
(I) s(x) ≥ 0 for all x ∈ [a, b], s(a) = 0, and s(b) = ℓ(γ).
(II) s is non-decreasing on [a, b].
(III) s 6≡ 0 if and only γ ′ (t) 6= 0̄ for some t ∈ [a, b], i.e., s ≡ 0 if and only if
γ ′ = 0̄ on [a, b].
(IV) If γ ′ (t) 6= 0̄ for all t ∈ [a, b] (i.e., γ is smooth), then s(x) = 0 if and only
if x = a.
(V) If γ ′ (t) 6= 0̄ for all t ∈ [a, b] (i.e., γ is smooth), then s is strictly increasing
on [a, b].
(VI) s is continuous on [a, b].
(VII) s is differentiable on [a, b].
d
(VIII) s is C1 -type on [a, b], and dx s(x) = ||γ ′ (x)|| for all x ∈ [a, b].
(IX) s is smooth if γ is smooth.
Proof. Assignment!
which does not look like 1 dx, rather it is weighted, i.e., instead of the fixed
a
unit weight 1 at each point t, a weight ||γ ′ (t)|| being carried by each point t. So,
one may ask what should replace “dx”, and which will look like “dx” so that
Zb
d
‘‘dx ′′ = length of the curve. Now, recall that if dx (f(x)) = g(x), we write it
a Z
in terms of language of integration/ symbol of integration as g(x)dx = f(x) or
Z Z Z
′
f (x)dx = f(x), or sometimes in short d(f(x)) = f(x), or simply df = f;
and therefore, in the light of (VIII) it is clear that “dx” should be replaced by
Z Zb
“ds” so that ds = s, and hence ds = s(b) − s(a) = length of the curve.
a
′
Further, if γ (t) 6= 0̄ for all t ∈ [a, b], clearly, by (V)th property we see that
s has its inverse s1 : [0, ℓ(γ)] −→ [a, b], and (VIII)th property ensures C1 -type-
ness of it. So, we consider the C1 -type curve γ ◦ s−1 : [0, ℓ(γ)] −→ R3 . It can be
d
checked that {γ} = {γ ◦ s−1 }, and || dt (γ ◦ s−1 (t))|| = 1, and therefore for γ ◦ s−1
Zx Zx
d −1
the arc length function is s(x) = || (γ ◦ s (t))||dt = 1 dt = x for all
0 dt 0
−1
x ∈ [0, ℓ(γ)]. The parametrisation γ◦s is called parametrisation of the
curve γ by its arc length function, or the arc-length parametrisation
of the curve γ.
Chapter 4
Line Integration
4.1 Introduction
In this chapter we shall study theory of integration of vector fields along a curve
aiming to classify a certain kind of vector fields called conservative vector fields.
2. If
Z φ : Z[c, d] −→ [a, Z b] is a smooth increasing function, then show that
f= f, i.e., f is independent of parametrisation of γ.
γ γ◦φ γ
27
i.e., F̄.ds̄ = < F̄(γ(t)), T̄ (t) > ds. Note that if F̄ is continuous, then the
γ γ
function < F̄(γ(t)), T̄ (t) > is also continuous, and therefore by Theorem 4.2.1,
we have
Z Z
F̄.ds̄ = < F̄(γ(t)), T̄ (t) > ds
γ γ
Zb
= < F̄(γ(t)), T̄ (t) > ||γ ′ (t)|| dt
a
Zb
γ ′ (t)
= < F̄(γ(t)), > ||γ ′ (t)|| dt
a ||γ ′ (t)||
Zb
< F̄(γ(t)), γ ′ (t) > ′
= ||γ (t)|| dt
a ||γ ′ (t)||
Zb
= < F̄(γ(t)), γ ′ (t) > dt, as γ ′ (t) 6= 0̄.
a
Thus we have the following theorem.
Theorem 4.3.1. Let G ⊂ R3 and F̄ : G −→ R3 a continuous vector field.
Z
Suppose γ : [a, b] −→ R3 a smooth curve such that Im(γ) ⊂ G. Then F̄
γ
Z Zb
exists, and F̄ = < F̄(γ(t)), γ ′ (t) > dt.
γ a
Example 4.3.4. 1. Let F̄(x, y) = (x2 , y2 ) for all (x, y) ∈ R2 , and γ(t) =
(r cos(t), r sin(t)) for all t ∈ [0, 2π]. Clearly F̄ is continuous and γ is
C1 -type, and ′
Z check that γ (t) 6= 0 for all t ∈ [0, 2π] (i.e., γ is smooth),
therefore, F̄ exists and,
γ
Z Z 2π
F̄ = < F̄(γ(t)), γ ′ (t) > dt
γ 0
Z 2π
= r3 < (cos(t)2 , sin(t)2 ), (− sin(t), cos(t)) > dt
0
Z 2π
= r3 < (cos(t)2 , sin(t)2 ), (− sin(t), cos(t)) > dt
0
=0
−y x
2. Let F̄(x, y) = , for all (x, y) ∈ R2 \(0, 0), and γ(t) =
x2 + y2 x2 + y2
(r cos(t), r sin(t)) for all
Z t ∈ [0, 2π]. Clearly F̄ is continuous and γ is
C1 -type, and therefore, F̄ exists and,
γ
Z Z 2π
F̄ = < F̄(γ(t)), γ ′ (t) > dt
γ 0
Z 2π
= < (− sin(t), cos(t)), (− sin(t), cos(t)) > dt
0
Z 2π
= 1 dt
0
= 2π
3. Let F̄(x, y) = (f1 (x), f2 (y)) for all (x, y) ∈ R2 where f1 , f2 are continuous
functions such that g1′ = f1 , and g2′ = f2 ; and γ(t) = (r cos(t), r sin(t))
1
Zfor all t ∈ [0, 2π]. Clearly F̄ is continuous and γ is C -type, and therefore,
F̄ exists and,
γ
Z Z 2π
F̄ = < F̄(γ(t)), γ ′ (t) > dt
γ 0
= 0(Check!!)
Since we know that for a C1 -type scalar field f, the vector field ∇(f) is
continuous, and directional derivatives Dv̄ (f)|P exists for all direction v̄, and
for all point P; and Dv̄ (f)|P =< ∇(f)|P , v̄ >. And since for a continuous
vector field F̄ in R3 and a smooth parametric curve γ : [a, b] −→ R3 , we
Z Zb
have F̄ = < F̄(γ(t)), γ ′ (t) > dt, one naturally ask what happens if
γ a
< F̄(γ(t)), γ ′ (t) > represents aZdirectional derivatives of a C1 scalar field, i.e.,
we are asking what happens to F̄, if F̄ = ∇(f) for some C1 -type scalar field f.
γ
Note that thought γ ′ (t) is not a unit vector, we can choose a parametrisation of
the given curve such that the derivative of it a unit vector (e.g., consider γ◦s−1 ).
Z Zb Zb
′
So, ∇(f) = < ∇(f)(γ(t)), γ (t) > dt = Dγ ′ (t) (f)|γ(t) dt
γ a a
Zb
= < ∇f(γ(t)), γ ′ (t) > dt
a
Zb
∂ d ∂ d ∂ d
= f|γ(t) x(t) + f|γ(t) y(t) + f|γ(t) z(t) dt
a ∂x dt ∂y dt ∂z dt
Zb
d
= f(γ(t))
a dt
= f(γ(b)) − f(γ(a)).
Zb Z
Thus, Dγ ′ (t) (f)|γ(t) dt = ∇(f) = f(γ(b)) − f(γ(a)). This leads to the
a γ
following theorem called the Fundamental Theorem of Line Integration.
Theorem 4.3.5. Let G ⊂ R3 and F̄ : G −→ R3 a continuous vector field. If
F̄ = ∇(f) for some scalar field f : G −→ R, then for Zany smooth parametric
curve γ with initial point A and final point B we have F̄ = f(B) − f(A).
γ
end points of the curve. If the end points of the curve are the same point, i.e.,
the curve is a loop, then clearly the value of the line integration is zero.
Definition 4.4.2. Let F̄ : G −→ R3 a continuousZ vector
Z field. We shall say that
line integration of F̄ is path independent, if F̄ = F̄ for any two smooth
γ1 γ2
curves γ1 and γ2 with same initial and final points.
Clearly, by Theorem 4.3.5, line integration of a continuous conser-
vative vector field is path independent. Now, suppose F̄ is a general
continuous vector field (non necessarily conservative). If line
Z integration
Z of F
is path independent, for any smooth loop γ, we shall have F̄ = F̄, as γ
γ γ−1
and γ−1 has same initial and final points,
Z and therefore,
Z sinceZ by the properties
Z
of line integration of vector fields F̄ = − F̄, we get F̄ = − F̄, i.e.,
Z γ−1 γ γ γ
F̄ = 0. So, naturally one asks whether the converse is true, i.e., suppose F̄ is a
γ Z
continuous vector field such that F̄ = 0 for all smooth loop γ, is then the line
γ
integration of F̄ path independent? The answer is, in fact, “yes”!. We describe
it in the theorem below.
Theorem 4.4.3. Let G ⊂ R3 and F̄ : G −→ R3 a continuous vector field. Then
the following are equivalent
(I) Line integration of F̄ is independent of path.
Z
(II) F̄ = 0 for all smooth loop γ in G.
γ
(II) =⇒ (I): Let γ1 and γ2 are two paths in G with same initial point A
Z Z
and final point B. We shall show that F̄ = F̄. Consider γ = γ1 ∗ γ−1
2 .
γ1 γ2
Clearly,
Z γ is a loop with initial and final point A, and therefore, by hypothesis,
F = 0. From which we see
γ Z
0 = F̄
γ
Z
= F̄
γ1 ∗γ−1
2
Z Z
= F̄ + F̄ .
γ1 γ−1
2
Z Z
= F̄ + − F̄
γ1 γ2
Z Z
= F̄ − F̄
γ1 γ2
Z Z
This shows that F̄ = F̄, which completes the proof.
γ1 γ2
From Theorem 4.3.5 and Theorem 4.4.3 we see the following implications
A: F̄ is a continuous CVF
⇓
B: Line integration of F is path independent
m
Z
C: F̄ = 0 for all smooth loop γ
γ
It can be seen that if G is open and path-connected, then every two points in
G can be connected by a piece-wise smooth path.
(III) F is a CVF.
Remark 4.5.3. Note that we can not prove a vector field to be conservative
using Theorem 4.5.2; rather it will be extremely helpful to disprove a VF is
conservative, if the VF satisfies all the hypotheses of the theorem.
−y x
Example 4.5.4. Let F̄(x, y) = , for all (x, y) ∈ R2 \(0, 0).
x2 + y2 x2 + y2
NoteZthat the domain of F̄ is path-connected, and F̄ is continuous. We have seen
that F̄ = 2π 6= 0 for any circle γ around (0, 0), and therefore by 1st Structure
γ
Theorem of CVF we get F̄ is not conservative!
One approach would be trying to find out a scalar field, on using some algo-
rithm, for which we can write the vector field as a gradient of the scalar field.
There are disadvantages of this approach, which will be discussing after we dis-
cuss the algorithm. Suppose F̄ : G −→ R3 a conservative vector field. We shall
try to find f : G −→ R such that F̄ = ∇(f) on G.
Step 3: Differentiating (2) with respect to y and comparing with (1b) we get
Z
∂ ∂
φ1 (y, z) = F2 − F1 ,
∂y ∂y x
and therefore,
Z Z
∂
φ1 (y, z) = F2 − F1 + φ2 (z) − −(3)
y ∂y x
Step 4: Again differentiating (4) with respect to z and comparing with (1c) we
get
Z Z Z Z
∂ ∂
φ2 (z) = F3 − F1 + F2 − F1 + c − −(5)
z ∂z x y ∂y x
Since we do not have, so far, a definitive tool which tells us a given vector
field is conservative, we need check whether ∇(f) = F̄ holds after finding f in
this process.
Z
So, f = sin(x)ex + φ1 (y, z) = sin(z)ex + φ( y, z).
x
Z
We know that ∂
∂y φ1 = F2 − ∂
∂y sin(z)ex = y2 , and therefore φ1 = y2 +
y
φ1 (y, z) = 1/3z3 + φ2 (z).
In this case we see that ∇(f) = F̄ on the domain of F̄, and therefore F̄ is a
conservative vector field.
2. Let F̄(x, y) = (ey , x). Following the above process we see that f = xey +
xy − ey + c, but it is easy to note that ∇(f) 6= F̄.
Let G ⊂ R2 be open and F̄ = (F1 , F2 ) : G −→ R2 a conservative vector
∂
field. Then there exists f : G −→ R such that ∇(f) = F̄ in G, i.e, ∂x f = F1
2 2
∂ ∂ ∂ ∂ ∂
and ∂y f = F2 ; and therefore ∂x∂y f = ∂x F2 , and ∂y∂x f = ∂y F1 . Thus, if
∂2 ∂2 ∂ ∂
∂x∂y f = ∂y∂x f, we shall have ∂x F2 = ∂y F1 .
This shows that line integration of F̄ along any loop vanishes and which,
from 1st Structure Theorem of CVF we see that F̄ is conservative.
Remark 4.6.4. 1. Note that one never achieves 1st structure theorem of
CVF from the 2nd structure theorem, and vice versa.
the holes. Further suppose C2 , C3 , · · · , Cn+1 are loops in G such that the region
of G entangled by Ci contains only one hole Hi for all i = 1, 2, · · · , n, and the
entangled region is inside the region entangled by C1 (i.e., Ci ’s are smaller loops
than C1 ). Orientation of all the loops are positive. Then
ZZ Z XZ
n+1
∂ ∂
F2 − F1 dx dy = F̄ − F̄.
G ∂x ∂y C1 Ci
i=2
∂ ∂
So, we orient C positively and choose F̄ for which ∂x F2 − ∂y F1 = 1. It
is to be noted that there are infinite number of such F̄, e.g., (0, x); (−y, 0);
(−y/2, x/2); (−y/n, (n − 1)x/n), etc. We choose F̄ = (0, x), which is C1 -type.
Since G is simply connected, C = ∂G is positively oriented, and F̄ is of C1 -type,
we get applying Green’s theorem
ZZ ZZ Z Z
∂ ∂
Area(G) = 1 dx dy = ( F2 − G2 ) dx dy = F̄ = (0, x).
G G ∂x ∂y ∂G C
Suppose C has a parametric form γ(t) = (x(t), y(t)) for all t ∈ [a, b]. Then,
Zb Zb
′ ′
Area(G) = < (0, x(t), (x (t), y (t))) > dt = x(t) y ′ (t) dt.
a a
x2 y2
Example 4.6.10. We shall find the area enclosed by C : 2 + 2 = 1, where
a b
a, b > 0 by Green’s theorem.
Chapter 5
Surface Integration
41