Vector Integral Calculus
Vector Integral Calculus
Vector Integral Calculus
Engineering Mathematics-I
Jamkhongam Touthang
Assistant Professor
Department of Applied Mathematics
Delhi Technological University
Delhi-110042
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Unit-I. Infinite series: Test for convergence of series (Comparison, Ratio, Root,
Integral Test, Raabe’s, Logarithmic), Alternating series, Absolute convergence,
Conditional convergence.
Unit-II. Differential and Integral Calculus of a single variable: Taylor’s &
Maclaurin’s expansion, Radius of curvature, Tracing of some standard curves,
Applications of definite integral to Area, Arc length, Surface Area and Volume
(in cartesian, parametric and polar coordinates).
Unit-III. Calculus of several variables: Partial differentiation, Euler’s theo-
rem, total differential, Taylor’s theorem, Maxima-Minima, Lagrange’s method
of multipliers, Application in estimation of error and approximation.
Unit-IV. Mutiple Integrals: Double integral (Cartesian and polar co-ordinates),
Change of order of integration, Triple integrals (Cartesian, Cylindrical and
Spherical co-ordinates), Beta and Gamma functions, Applications of multiple
integration in area and volume.
Unit-V. Vector Differential Calculus: Continuity and differentiability of vector
functions, Scalar and Vector point function, Gradient, Directional derivative,
Divergence, Curl and their applications.
Unit-VI. Vector Integral Calculus: Line integral, surface integral and volume
integral, applications to work done by the force, applications of Green’s, Stoke’s
and Gauss divergence theorems.
2
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Unit-VI
Mathematics-I
3
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Vector Integration
Let f~(t) and ~g (t) be two vector functions such that
d~g (t)
= f~(t),
dt
then ~g (t) is called an integral of f~(t) with respect to the scalar variable t and
we write
Z
f~(t) dt = ~g (t).
d~g (t) d
f~(t) = = [~g (t) + K],
dt dt
so that
Z
f~(t) dt = ~g (t) + K.
This is called the indefinite integral of f~(t) and its definite integral is
Z b
f~(t) dt = [~g (t) + K]ba = ~g (b) − ~g (a).
a
4
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Line Integral
Let f~(x, y, z) be a continuous vector function defined at all points of a curve C in
space. We partition the curve C into n parts by the points A = P0 , P1 , . . . , Pi−1 , Pi , . . . , Pn−1 , Pn =
B whose position vectors are respectively ~r0 , ~r1 , . . . , ~ri−1 , ~ri , . . . , ~rn−1 , ~rn . We
consider a point P on the arc Pi−1 Pi whose position vector is ~ri . Form the sum
n
X
S= f~(~ri ).δ~ri , where δ~ri = ~ri − ~ri−1 .
i=0
The limit of this sum as n → ∞ such that |δ~ri | → 0, provided it exists, is called
the line integral of f~ along the curve C and we write
Z Z
d~r
f~. d~r or f~. dt .
C C dt
I Z
If C is a closed curve, we write in place of .
Work Done
If f~ represents the force acting on a particle moving along an arc P Q, then the
work done during the small displacement δ~r = f~.δ~r.
The total work done by f~ in displacing a particle from P to Q is given by the
Z Q
line integral f. d~r .
P
Examples
Z
1. If f~ = (5xy − 6x2 )î + (2y − 4x)ĵ, evaluate f~. d~r along the curve C in the
C
xy-plane, y = x3 from the point (1, 1) to (2, 8).
In the xy-plane, ~r = x î + y ĵ so that d~r = dx î + dy ĵ.
5
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
= (64 + 32 − 48 − 16) − (1 + 1 − 3 − 2)
= 35 .
Z
2. If f = x2 y(1 + z), evaluate f d~r, where C is the curve x = t, y = t2 ,
C
z = (1 − t) from t = 0 to t = 1.
In 3D space, ~r = x î + y ĵ + z k̂ so that d~r = dx î + dy ĵ + dz k̂.
Therefore,
Z Z
f d~r = x2 y(1 + z)(dx î + dy ĵ + dz k̂)
C C
Z Z Z
2 2
= î x y(1 + z)dx + ĵ x y(1 + z)dy + k̂ x2 y(1 + z)dz
C C C
Z 1 Z 1 Z 1
= î t4 (2 − t)dt + ĵ t4 (2 − t)2tdt + k̂ t4 (2 − t)(−1)dt
0 0 0
7 8 7
= î + ĵ − k̂ .
30 21 30
Z
~ 2
3. If f = 3xy î − y ĵ, evaluate f~. d~r, where C is the curve in the xy-plane
C
y = 2x2 from the point (0, 0) to (1, 2).
7
Ans. −
6
Z
4. If f~ = (3x2 + 6y)î − 14yz ĵ + 20xz 2 k̂, evaluate f~. d~r, where C is the curve
C
x = t, y = t2 , z = t3 from the point (0, 0, 0) to (1, 1, 1).
Ans. 5
Z
~ 2
5. If f = xy î − z ĵ + x k̂, evaluate f~ × d~r, where C is the curve x = t2 ,
C
6
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
8. Find the work done when a force f~ = (x2 − y 2 + x)î − (2xy + y)ĵ moves a
particle in the xy-plane from (0, 0) to (1, 1) along the straight line y = x. Is the
work done different if the path is the parabola y 2 = x?
In the xy-plane, ~r = x î + y ĵ so that d~r = dx î + dy ĵ.
The work done by f~ is
Z Z
f~. d~r = [(x2 − y 2 + x)î − (2xy + y)ĵ].[dx î + dy ĵ]
C C
Z
= [(x2 − y 2 + x)dx − (2xy + y)dy] (1)
C
7
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
9. Find the work done by the force f~ = 3x2 î + (2xz − y)ĵ + z k̂ in displacing a
particle from x = 0 to x = 2 along the curve x2 = 4y, 3x3 = 8z.
Ans. 16
−y î + x ĵ
11. Find the work done by f~ = along the upper half of the circle
x2 + y 2
x2 + y 2 = 1.
Ans. π
8
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Surface
Equation of a surface, say S, is of the form f (x, y, z) = 0. Its parametric
representation is of the form ~r(x, y) = r1 (x, y)î + r2 (x, y)ĵ + r3 (x, y)k̂. The
continuous functions x = φ(t) and y = ψ(t) of a scalar variable t represent a
curve C on the surface S.
A surface S is said to be a smooth surface if it has a unique normal at each of
its points whose direction depends continuously on the points lying on it. S is
said to be piecewise smooth if it is not smooth but it can be partitioned into
finitely many smooth pieces. For example, a sphere is a smooth surface while a
cube is piecewise smooth.
If the positive normal direction at any point of a surface S can be continued
in a unique and continuous way to the whole of S, then we say that S is two
sided or orientable. Otherwise, S is said to be one sided or non-orientable. For
example, a cylinder is a two sided surface while the Möbius strip is one sided.
Surface Integral
Let S be a two sided surface and f~ be a vector point function defined at all
~ be
points of S. We partition S into a finite number of sub-surfaces and let δ S
an elementary surface area surrounding a point P (on S) whose position vector
~ can be thought of as a vector whose magnitude is the area and
is ~r. Here, δ S
whose direction that of the outward normal to the element. We form the sum
X
f~. δ S
~
where the summation extends over all the sub-surfaces. The limit of this sum
as the number of sub-surfaces approaches infinity in such a way that the area of
each sub-surface tends to zero, is called the surface integral of f~(~r) over S and
is denoted by
Z Z
f~. dS
~ or, f~. n̂ dS,
S S
Volume Integral
Let S be a surface enclosing a region V on which a vector function f~ is defined.
The region V is partitioned into a finite number of sub-regions V1 , V2 , . . . , Vn .
Let δvi be the volume of the sub-region Vi enclosing a point whose position
vector is ~ri .
n
X
Form the sum f~(~ri ) δvi .
i=1
The limit of this sum as n → ∞ in such a manner that δvi → 0 is called the
9
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Z
volume integral of f~ over V and is denoted by f~ dv.
V
Integral Theorems
In vector calculus, the three theorems-Green’s theorem, Stoke’s theorem and
Gauss divergence theorem, are important both from theoretical and practical
considerations.
Green’s theorem in a plane provides a relationship between a line integral
(around a closed curve) and a double integral. Stoke’s theorem, a general-
ization of Green’s theorem, provides a relationship between a line integral and
a surface integral. Gauss divergence theorem provides a relationship between a
surface integral and a volume integral.
Green’s Theorem
If φ(x, y) and ψ(x, y) are continuous functions of x and y having continuous
∂φ ∂ψ
partial derivatives and in a region R of the xy-plane, bounded by a
dy dx
simple closed curve C, then
I ZZ
∂ψ ∂φ
(φ dx + ψ dy) = − dx dy
C R dx dy
Stoke’s Theorem
Let S be an open surface bounded by a closed curve C and f~ be any vector
function having continuous first order partial derivatives on S, then
I ZZ
f~. d~r = curl f~. n̂ dS,
C S
10
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Then,
f~. d~r = f1 dx + f2 dy + f3 dz
~ ∂f3 ∂f2 ∂f1 ∂f3 ∂f2 ∂f1
and curl f = − î + − ĵ + − k̂ .
∂y ∂z ∂z ∂x ∂x ∂y
Therefore,
curl f~.dS
~ = curl f~.n̂ dS
∂f3 ∂f2 ∂f1 ∂f3
= − (n̂.î)dS + − (n̂.ĵ)dS
∂y ∂z ∂z ∂x
∂f2 ∂f1
+ − (n̂.k̂)dS
∂x ∂y
∂f3 ∂f2 ∂f1 ∂f3 ∂f2 ∂f1
= − dy dz + − dz dx + − dx dy
∂y ∂z ∂z ∂x ∂x ∂y
Substituting these values in Stoke’s theorem, we obtain
I Z Z
∂f3 ∂f2 ∂f1 ∂f3 ∂f2 ∂f1
(f1 dx + f2 dy + f3 dz) = − dydz + − dzdx + − dxdy
C S ∂y ∂z ∂z ∂x ∂x ∂y
which is Stoke’s theorem in the scalar form.
where n̂ is the outward drawn unit normal vector to the surface S. Consider
S to be a plane surface (region) R in the xy-plane bounded by a simple closed
curve C. Now,
and
î ĵ k̂
= ∂f2 − ∂f1 . k̂
∂ ∂ ∂
~
curl f =
∂x ∂y ∂z
∂x ∂y
f1 f2 0
11
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Normal at any point of the region R in the xy-plane is parallel to the z-axis.
Considering the positive direction of the z-axis as the positive (outward) direc-
tion of n̂, we get n̂ = k̂. Also, dS is a surface area element in the xy-plane and
thus dS = dx dy. Therefore, (3) becomes
I ZZ
∂f2 ∂f1
(f1 dx + f2 dy) = − dx dy
C R ∂x ∂y
which is Green’s theorem in a plane.
Stoke’s theorem problems will reduce to Green’s theorem problems if the surface
S is a plane surface (region).
Green’s Formulas
Set f~ = g∇h, where g and h are scalar functions, so that
div f~ = div (g∇h) = ∇.(g∇h) = g∇2 h + ∇g.∇h. Making these changes, Gauss
divergence theorem becomes
ZZ ZZZ
~
g∇h. dS = (g∇2 h + ∇g.∇h) dV (6)
S V
12
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Examples
I
1. Verify Green’s theorem in a plane for [(xy + y 2 ) dx + x2 dy], where C is
C
the boundary of the region bounded by y = x and y = x2 .
I I
Comparing the given line integral [(xy + y 2 ) dx + x2 dy] with (φ dx + ψ dy),
C C
we have
φ = xy + y 2 , ψ = x2
so that
∂φ ∂ψ
= x + 2y, = 2x .
dy dx
To verify Green’s theorem, we show that
I ZZ
∂ψ ∂φ
(φ dx + ψ dy) = − dx dy .
C R dx dy
Here,
I I
(φ dx + ψ dy) = [(xy + y 2 ) dx + x2 dy]
C C
I I
2 2
= [(xy + y ) dx + x dy] + [(xy + y 2 ) dx + x2 dy] (8)
C1 C2
Hence,
I
19 1
(φ dx + ψ dy) = −1=−
C 20 20
13
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Also
ZZ Z 1 Z x
∂ψ ∂φ
− dx dy = (x − 2y) dy dx
R dx dy 0 x2
Z 1
= (x4 − x3 ) dx
0
1
=− .
20
I ZZ
∂ψ ∂φ
This shows that (φ dx + ψ dy) = − dx dy .
C R dx dy
I
2. Use Green’s theorem to evaluate [(2x2 − y 2 ) dx + (x2 + y 2 ) dy], where C
C
is the boundary of the region enclosed by upper-half of x2 + y 2 = a2 and the
x-axis.
Let R be the region eclosed by C. By Green’s theorem
I ZZ
2 2 2 2 ∂ 2 2 ∂ 2 2
[(2x − y ) dx + (x + y ) dy] = (x + y ) − (2x − y ) dx dy
C dx dy
Z ZR
= (x + y) dx dy
R
Set x = r cos θ and y = r sin θ so that r varies from 0 to a and θ varies from 0
to π. Therefore,
I Z aZ π
[(2x2 − y 2 ) dx + (x2 + y 2 ) dy] = 2 r(cos θ + sin θ). r dθ dr
C 0 0
Z a Z π
2
=2 r dr (cos θ + sin θ) dθ
0 0
a3
= 2. .2
3
4a3
= .
3
3. Apply Green’s theorem to find the area between the parabolas y 2 = 4x and
x2 = 4y.
16
Ans.
3
I
1
4. Prove that the area enclosed by a plane curve is (x dy − y dx) by using
2 C
Green’s theorem.
I
5. Verify Green’s theorem in the xy-plane for [(xy 2 − 2xy) dx + (x2 y + 3) dy],
C
where C is the boundary of the region bounded by y 2 = 8x and x = 2.
14
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
6. Verify Stoke’s theorem for f~ = (x2 − y 2 ) î + 2xy ĵ, integrated round the
rectangle in the xy-plane (z = 0) bounded by the lines x = 0, y = 0, x = a,
y = b.
To verify Stoke’s theorem, we show that
I ZZ
f~. d~r = curl f~. n̂ dS
C S
a3 a3
I
Hence (9) yields f~. d~r = + ab2 + ab2 − + 0 = 2ab2 .
C 3 3
15
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Now,
curl f~ = ∇ × f~
∂ ∂ ∂
= î + ĵ + k̂ × [(x2 − y 2 )î + 2xy ĵ]
∂x ∂y ∂z
î ĵ k̂
∂ ∂ ∂
= = î(0) − ĵ(0) + k̂(2y + 2y) = 4y k̂
2 ∂x 2 ∂y ∂z
x − y 2xy 0
For the surface S, n̂ = k̂. Therefore, curl f~. n̂ = 4y k̂.k̂ = 4y and hence
ZZ ZZ
dx dy
curl f~. n̂ dS = 4y
S S |k̂.n̂|
Z aZ b
= 4y dy dx = 2ab2
0 0
I ZZ
This verifies Stoke’s theorem, that is, f~. d~r = curl f~. n̂ dS .
C S
I
7. Apply Stoke’s theorem to evaluate (yz dx + zx dy + xy dz), where C is the
C
curve x2 + y 2 = 1, z = y 2 .
We have
I I
f~. d~r = (yz î + zx ĵ + xy k̂).(dx î + dy ĵ + dz k̂)
C
IC
= (yz dx + zx dy + xy dz)
C
By Stoke’s theorem,
I ZZ
f~. d~r = curl f~. n̂ dS
C S
But
î ĵ k̂
∂ ∂ ∂
~
curl f = = î(0) − ĵ(0) + k̂(0) = 0̂
∂x ∂y ∂z
yz zx xy
ZZ
so that curl f~. n̂ dS = 0.
S
I
Hence, it follows from Stoke’s theorem, (yz dx + zx dy + xy dz) = 0 .
C
16
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
I
8. Apply Stoke’s theorem to evaluate [(x + y)dx + (2x − z)dy + (y + z)dz],
C
where C is the boundary of the triangle with vertices P (2, 0, 0), Q(0, 3, 0) and
R(0, 0, 4).
We have, f~ = (x + y) î + (2x − z) ĵ + (y + z) k̂ and thus
î ĵ k̂
∂ ∂ ∂
curl f~ =
= 2î + k̂
∂x ∂y ∂z
x+y 2x − z y+z
Equation of the plane through the points P (2, 0, 0), Q(0, 3, 0) and R(0, 0, 4) is
x y z
+ + =1 or 6x + 4y + 3z = 12 .
2 3 4
Vector normal to the plane is
∂ ∂ ∂
∇(6x + 4y + 3z − 12) = î + ĵ + k̂ (6x + 4y + 3z − 12)
∂x ∂y ∂z
= 6î + 2ĵ + 3k̂
and
ZZ ZZ
1
curl f~. n̂ dS = (2î + k̂). (6î + 2ĵ + 3k̂) dS
S S 7
ZZ
15
= dS
7 S
15
= × Area of the 4P QR
7
15
= × 244 (use Heron’s formula to find the area)
7
I
which gives the required value of [(x + y)dx + (2x − z)dy + (y + z)dz] .
C
17
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
I
9. Apply Stoke’s theorem to evaluate [(sin x − y)dx − cos xdy], where C is
C
the boundary of the triangle with vertices A(0, 0), B(π/2, 0) and C(π/2, 1).
2 π
Ans. +
π 4
10. Verify Stoke’s theorem for f~ = y î + z ĵ + x k̂, where C is the boundary of
the upper half surface S of the sphere x2 + y 2 + z 2 = 1 .
On S, z = 0 so that dz = 0 and x2 + y 2 = 1.
Therefore
I I I
f~. d~r = (y î + z ĵ + x k̂).(dx î + dy ĵ + dz k̂) = y dx
C C x2 +y 2 =1
Set x =Icos θ, y = sinZ θ2πso that dx = − sin θ dθ and θ varies from 0 to 2π.
Then f~. d~r = − sin2 θ dθ = −π .
C 0
Also
ZZ ZZ
curl f~. n̂ dS = (−î − ĵ − k̂). k̂ dS
S S
(since curl f~ = −î − ĵ − k̂ and n̂ = k̂ on S so that k̂.n̂ = 1)
ZZ
dx dy
=− = −(area of x2 + y 2 = 1) = −π(1)2 = −π
S |k̂.n̂|
I
= f~. d~r , which means Stoke’s theorem is verified
C
11. Verify Stoke’s theorem when f~ = (2xy − x2 )î − (x2 − y 2 )ĵ and C is the
boundary of the region enclosed by the parabolas x2 = y and y 2 = x.
I ZZ
~ 3
Ans. f . d~r = curl f~. n̂ dS = −
C S 5
12. Verify Stoke’s theorem when f~ = xy î−2yz ĵ−zx k̂ and S is the surface of the
rectangular parallelopiped bounded by the planes x = 0, x = 1, y = 0, y = 2,
and z = 3 above the xy-plane.
I ZZ
Ans. f~. d~r = curl f~. n̂ dS = −1
C S
18
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
S1 : OAC 0 B, S2 : CB 0 P A0 , S3 : OBA0 C
S4 : AC 0 P B 0 , S5 : OCB 0 A, S6 : BA0 P C 0 y
O B
On S1 , n̂ = −k̂ and z = 0 so that
A C0
f~ = x2 î + y 2 ĵ and hence
x
ZZ ZZ
f~. n̂ dS = (x2 î + y 2 ĵ + z 2 k̂).(−k̂) dS
S1 S1
=0
On S2 , n̂ = k̂ and z = 1 so that f~ = x2 î + y 2 ĵ + k̂ and thus
ZZ ZZ ZZ
dx dy
f~. n̂ dS = (x2 î + y 2 ĵ + k̂).(k̂) dS = (1)
S2 S2 S2 |k̂.n̂|
Z 1Z 1
= dx dy
0 0
=1
On S4 , n̂ = î and x = 1 so that f~ = î + y 2 ĵ + z 2 k̂ .
ZZ ZZ Z 1 Z 1
dy dz
f~. n̂ dS = (î + y 2 ĵ + z 2 k̂).(î) dS = (1)
S4 S4 0 0 |î.n̂|
Z 1Z 1
= dy dz = 1 .
0 0
19
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
On S6 , n̂ = ĵ and y = 1 so that f~ = x2 î + ĵ + z 2 k̂ .
ZZ ZZ Z 1 Z 1
dx dz
f~. n̂ dS = (x2 î + ĵ + z 2 k̂).(ĵ) dS = (1)
S6 S4 0 0 |ĵ.n̂|
Z 1Z 1
= dx dz = 1 .
0 0
ZZ
Hence, f~. n̂ dS = 0 + 1 + 0 + 1 + 0 + 1 = 3 .
S
Further,
ZZZ Z 1 Z 1 Z 1
( div f~ ) dV = ∇.(x2 î + y 2 ĵ + z 2 k̂) dx dy dz
V 0 0 0
Z 1Z 1Z 1
∂ ∂ ∂
= î + ĵ + k̂ .(x2 î + y 2 ĵ + z 2 k̂) dx dy dz
0 0 0 ∂x ∂y ∂z
Z 1Z 1Z 1
= (2x + 2y + 2z) dx dy dz
0 0 0
1 11 Z 1Z 1
x2
Z Z
1
=2 + xy + xz dy dz = 2 y+z+ dy dz
0 0 2 0 0 0 2
Z 1 2 1
y 1
=2 + yz + y dz
0 2 2 0
Z 1
=2 (1 + z) dz
0
= 3.
div f~ = ∇.f~
∂ ∂ ∂
= î + ĵ + k̂ .(4x î − 2y 2 ĵ + z 2 k̂)
∂x ∂y ∂z
= 4 − 4y + 2z.
20
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
S1 : OAC 0 B, S2 : CB 0 P A0 , S3 : OBA0 C
O B X
S4 : AC 0 P B 0 , S5 : OCB 0 A, S6 : BA0 P C 0
A C0
Z
21
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
22
MA101 Mathematics-I jamkhongamtouthang@dtu.ac.in
Hence, Gauss divergence theorem is verified for the given vector field.
Z
18. Apply Gauss divergence theorem to evaluate f~. n̂ dS , where f~ = y 2 z 2 î +
S
z 2 x2 ĵ + x2 y 2 k̂ and S is the upper part of the sphere x2 + y 2 + z 2 = a2 above
the xy-plane.
πa6
Ans.
12
19. Verify Gauss divergence theorem for f~ = y î + x ĵ + z 2 k̂ over the cylindrical
region bounded by x2 + y 2 = 9, z = 0 and z = 2.
Note: Unsolved problems would be either solved in the classes or left as exer-
cises. Use this material as supplementary to class lectures.
23