Engineering Mathematics I
Engineering Mathematics I
Engineering Mathematics I
Mathematics-I
Prof. C. Nahak,
Prof. J. Kumar
Prof. S. Kumar,
ENGINEERING MATHEMATICS – I (3+0)
Course Developer
Lesson 1
1.1 Introduction
In this lesson first we will state the Rolle’s theorems, mean value theorems and
study some of its applications.
5 www.AgriMoon.Com
Rolle’s Theorem, Lagrange’s Mean Value Theorem, Cauchy’s Mean Value Theorem
Hence,
when
inequalities we have, .
Note: Rolle’s theorem shows that b/w any two zero’s of a function there
6 www.AgriMoon.Com
Rolle’s Theorem, Lagrange’s Mean Value Theorem, Cauchy’s Mean Value Theorem
Solution:
Example 2: in .
Solution:
and , for .
Example 3: Show that the equation , has only one real root
Solution:
7 www.AgriMoon.Com
Rolle’s Theorem, Lagrange’s Mean Value Theorem, Cauchy’s Mean Value Theorem
But , a contradiction to
Rolle’s therorem. Hence the equation has only one real root.
, where .
then there is a point between and at which the tangent is parallel to the
Cauchy's mean value theorem, also known as the extended mean value theorem,
is the more general form of the mean value theorem.
and are both continuous on the closed interval , and differentiable on the
8 www.AgriMoon.Com
Rolle’s Theorem, Lagrange’s Mean Value Theorem, Cauchy’s Mean Value Theorem
that
Note 1: Cauchy's mean value theorem can be used to prove L'Hospital's rule.
The mean value theorem (Lagrange) is the special case of Cauchy's mean value
theorem when .
Note 2: The proof of Cauchy's mean value theorem is based on the same idea as
the proof of the mean value theorem
on .
9 www.AgriMoon.Com
Rolle’s Theorem, Lagrange’s Mean Value Theorem, Cauchy’s Mean Value Theorem
Solution:
Clearly iff
, . . Hence
i.e., implies , so .
such that .
10 www.AgriMoon.Com
Rolle’s Theorem, Lagrange’s Mean Value Theorem, Cauchy’s Mean Value Theorem
number where .
Solution:
But and .
11 www.AgriMoon.Com
Rolle’s Theorem, Lagrange’s Mean Value Theorem, Cauchy’s Mean Value Theorem
Solution:
Solution:
(a) f ( x) = x 2 − 3 x + 2 on [1,2]
2. The function f ( x) = 4 x 3 + x 2 − 4 x − 1 has roots 1 and -1. Find the root of the
derivative f ′( x) mentioned in Rolle’ s throrem.
12 www.AgriMoon.Com
Rolle’s Theorem, Lagrange’s Mean Value Theorem, Cauchy’s Mean Value Theorem
sin x
5. Using Cauchy’s mean value theorem show that lim =1
x →0 x
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
13 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 2
2.1 Introduction
14 www.AgriMoon.Com
Taylor's theorem / Taylor’s expansion, Maclaurin’s expansion
difference between the Taylor polynomial of degree n and the original function.
The remainder term depends on and is small if is close enough to .
Several expressions are available for it. The Lagrange form is given by
where
Notes
• In fact, the mean value theorem is used to prove Taylor's theorem with the
Lagrange remainder term.
15 www.AgriMoon.Com
Taylor's theorem / Taylor’s expansion, Maclaurin’s expansion
• However, for some functions , one can show that the remainder term
Solution:
Here , , So
16 www.AgriMoon.Com
Taylor's theorem / Taylor’s expansion, Maclaurin’s expansion
Solution:
for , we have
17 www.AgriMoon.Com
Taylor's theorem / Taylor’s expansion, Maclaurin’s expansion
Ans: .
5. Using the results of above problem, estimate the error of the approximate
1 1
equation 1 + x ≈ 1 + x − x 2 when x = 0.2.
2 8
6. Write down the Taylor’s expansion for the function f ( x) = sin x about the point
π
a= with n = 4.
4
x x2 x
7. Applying Taylor’s theorem with remainder prove that 1 + − < 1 + x < 1 + if
2 8 2
x > 0.
(i) ln (1 + x) (ii) (1 + x) .
m
18 www.AgriMoon.Com
Taylor's theorem / Taylor’s expansion, Maclaurin’s expansion
References
W. Thomas Finny (1998). Calculus and Analytic Geometry, 6th Edition, Publishers,
Narsa, India.
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
19 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 3
3.1 Introduction
In the first limit if we put we will get and in the second limit if we
will behave in the same fashion that it’s largest power behaves). Both of these
are called Indeterminate form.
First limit can be found by the factorizing the numerator cancelling the common
factor. That is
20 www.AgriMoon.Com
Indeterminate forms ; L’Hospital’s Rule
first is a indeterminate form, but we can’t factor this one. The second is an
indeterminate form, but we can’t just factor an out of the numerator. Does
there exists some method to evaluate the limits? The answer is yes. By
(L'Hospital's Rule).
or
where can be any real number, infinity or negative infinity. In these cases we
have,
21 www.AgriMoon.Com
Indeterminate forms ; L’Hospital’s Rule
and so
..................(1)
22 www.AgriMoon.Com
Indeterminate forms ; L’Hospital’s Rule
and, finally,
Note 3.1: The theorem also holds for the case where the functions and
23 www.AgriMoon.Com
Indeterminate forms ; L’Hospital’s Rule
infinity.
Example 3.1:
, we find that
24 www.AgriMoon.Com
Indeterminate forms ; L’Hospital’s Rule
Solution:
Taking derivative both numerator and denominator five times we obtain: Ans: 3
The other indeterminate forms reduce to the following cases. (a) (b) (c)
(d) (e) .
25 www.AgriMoon.Com
Indeterminate forms ; L’Hospital’s Rule
or If &
( )- form
Example 3.3
26 www.AgriMoon.Com
Indeterminate forms ; L’Hospital’s Rule
have
So .
Ans: 1
27 www.AgriMoon.Com
Indeterminate forms ; L’Hospital’s Rule
Ans: 1
1.
2.
3.
4.
5.
6.
7.
28 www.AgriMoon.Com
Indeterminate forms ; L’Hospital’s Rule
8.
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
29 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 4
4.1 Introduction
30 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
two variables and . For every and we put in, we get a number out. The
set of all we allowed to put into the function is called the domain of the
function. Usually the domain is unspecified, and then the domain is the set of all
we can put into the formula for and not get square roots of negatives, or
division by zero, or some such. i.e.,the domain is usually the set of all we
range is [0,7].
31 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
parallel to the -plane. That is, let for some number , and plot
in the -plane.
now onwards is the set of all points in the -plane for which is
32 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Solution:
curve.
As a result, we define an open region to be the set of all points inside of but not
including a closed curve, and we define a closed region to be the set of all
points inside of and including a closed curve.
33 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Solution:
square root. Thus, the domain of is the set of points that satisfy
or .
That is, the domain is the set of points inside and on the circle of radius 3
Moreover, the domain is a closed region of the -plane since it contains the
Functions of two variables are important for reasons other than that their graph
is a surface. In particular, a function of the form u(x,t) is often interpreted to be
a function of x at a given point in time. For example, let's place an xy-
coordinate system on a violin whose strings have a length of l, If u(x,t) is
considered the displacement of a string above or below a horizontal line at a
point and at a time , then y = u(x,t) is the shape of the string at a fixed time t.
34 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Likewise, u (x, t) might represent the temperature at a distance xfrom one end of
the rod at time .
Now we will extend the properties of limits and continuity from the familiar
function of one variable to the new territory of functions of two or more
variables.
In less formal language this means that, if the limit holds, then gets closer
Good job if you saw this as “limit does not exist” indicating a vertical
asymptote at .
35 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
This limit is indeterminate. With some algebraic manipulation, the zero factors
could cancel and reveal a real number as a limit. In this case, factoring leads
to……
The limit exists as approaches 2 even though the function does not exist. In
the first case, zero in the denominator led to a vertical asymptote; in the second
case the zeros cancelled out and the limit reveals a hole in the graph at .
36 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
whenever < .
Graphically for any point in the disc with radius , the value
For the limit of this function to exist at (-1,3), values of must get closer to 13
Solution:
. Let .
37 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Solution:
Now
Put , whenever .
Example 4.9.
Solution:
Now
. Set .
38 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Solution:
39 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
notation:
Give Definition
point in .
The following results are presented without proof. As was the case in functions
of one variable, continuity is “user friendly”. In other words, if is a real
number and and are continuous functions at then the functions below
40 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Finally, the following result asserts that the composition of continuous functions
are also continuous. If is continuous at and is continuous at ,
and
Example 4.12 Find the limit and discuss the continuity of the function
Solution:
continuous at origin.
41 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Solution:
. Take .
continuous?
Solution:
Note that
every where.
42 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Solution:
0.
bounded domain , then there will be at least one point in such that
43 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Corollary to property 2.
If a function is continuous in a closed and bounded domain and
assumes both positive and negative values, then there will be a point inside the
domain at which the vanishes.
1. Find , if it exists.
2. Show that
3. Prove that .
4. Find
44 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
5. Find , if it exists.
function at .
8. Find the
References
45 www.AgriMoon.Com
Limit, Continuity of Functions of Two Variables
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
46 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 5
5.1 Introduction
define as
and similarly
Solution:
, .
47 www.AgriMoon.Com
Partial and total derivatives
Solution:
means
48 www.AgriMoon.Com
Partial and total derivatives
or equivalently
49 www.AgriMoon.Com
Partial and total derivatives
approximated very well by its tangent line. This observation is the bases for
linear approximation.
50 www.AgriMoon.Com
Partial and total derivatives
The vector
or
51 www.AgriMoon.Com
Partial and total derivatives
false, as we know existence of partial derivative does not mean the function of
two variable is continuous. We might suspect that if is continuous at
and the first order partial derivatives exist there, then is differentiable at
Solution:
as and . That is
52 www.AgriMoon.Com
Partial and total derivatives
get .
The natural question to ask then is under what conditions can we conclude that
is differentiable at . The answer is contained in the following theorem.
all of the first order partial derivatives of exist in a open disk about
Solution:
The domain of is all of except for the origin. We shall show that has
continuous partial derivatives everywhere in its domain (that is, the function
53 www.AgriMoon.Com
Partial and total derivatives
continuous everywhere except the origin (where the denominators are zero).
Thus, is differentiable everywhere in its domain.
the origin.
Solution:
Since the limit does not exit so does not exit. Similarly we can show
54 www.AgriMoon.Com
Partial and total derivatives
In Example 5.3 the partial derivatives and did not exist at the origin and
Solution:
Now
i.e.,
That is
55 www.AgriMoon.Com
Partial and total derivatives
If this limit exists, we get the same value no matter how and approach 0.
But this limit does not exist, since small values for will make the fraction
arbitrarily large. Thus, this function is not differentiable at the origin, even
though the partial derivatives and exist.
Theorem 5.1 :
56 www.AgriMoon.Com
Partial and total derivatives
Here subscript 1 and 2 denote the partial derivative with respect to its first and
second argument, respectively. The proof is given in Lesson 7.
The differentials and are independent variables; that is, they can be given
any values. Then the differential , also called the total differential, is defined
by
Further, if changes from 2 to 2.05 and changes from 3 to 2.96, compare the
Solution:
By definition,
57 www.AgriMoon.Com
Partial and total derivatives
The increment of is
5.2 Total derivative: In the mathematical field of differential calculus, the term
total derivative has a number of closely related meanings.
respect to one of its input variables, e.g., , is different from the partial
assume that the other arguments are constant while varies; instead, it allows
the other arguments to depend on . The total derivative adds in these indirect
58 www.AgriMoon.Com
Partial and total derivatives
depends on , some of that change will be due to the partial derivative of with
respect to . However, some of that change will also be due to the partial
Solution:
, , .
59 www.AgriMoon.Com
Partial and total derivatives
3.
4. Find at where
References
60 www.AgriMoon.Com
Partial and total derivatives
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
61 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 6
6.1 Introduction
Definition 6.1
Example 6.2
62 www.AgriMoon.Com
Homogeneous Functions, Euler's Theorem
the axes.
Example 6.3 .
for all .
obtain
Chain rule :
Finally set .
63 www.AgriMoon.Com
Homogeneous Functions, Euler's Theorem
Proof. Let .
But ,
Hence
.............................
..................(6.2)
64 www.AgriMoon.Com
Homogeneous Functions, Euler's Theorem
..................(6.3)
65 www.AgriMoon.Com
Homogeneous Functions, Euler's Theorem
Ans.: ).
Solution:
fucntion of degree 2
i.e.,
or
Solution:
66 www.AgriMoon.Com
Homogeneous Functions, Euler's Theorem
i.e.,
1. If , show that
2. If , show that
3. If , prove that
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
67 www.AgriMoon.Com
Homogeneous Functions, Euler's Theorem
Suggested Readings
68 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 7
7.1 Introduction
The chain rule works for functions of more than one variable. Consider the
function where and , and and are
differentiable with respect to , then
69 www.AgriMoon.Com
Composite and Implicit Functions (Chain Rule) for Two Variables
It is noted that total increment is not equal to the sum of the partial increments,
. Let us assume that has continuous partial derivatives
at the point under consideration. Express in terms of partial
derivatives. To do this we have
70 www.AgriMoon.Com
Composite and Implicit Functions (Chain Rule) for Two Variables
Theorem 7.1:
Example 7.1:
71 www.AgriMoon.Com
Composite and Implicit Functions (Chain Rule) for Two Variables
Hence,
This formula is known as the formula for calculating the total derivative (in
Solution:
72 www.AgriMoon.Com
Composite and Implicit Functions (Chain Rule) for Two Variables
, , , , , .
So
Solution:
, , .
73 www.AgriMoon.Com
Composite and Implicit Functions (Chain Rule) for Two Variables
7.1.1 Let us find the the total differential of the composite function
and and , we know the total differential
Where and
Solution:
74 www.AgriMoon.Com
Composite and Implicit Functions (Chain Rule) for Two Variables
75 www.AgriMoon.Com
Composite and Implicit Functions (Chain Rule) for Two Variables
which implies
find .
Solution:
we obtain .
1. Find at where
2. If ,
, then show that
3. .
4. Find , when
76 www.AgriMoon.Com
Composite and Implicit Functions (Chain Rule) for Two Variables
constant, .
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
77 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 8
8.1 Introduction
functions along with partial derivatives are continous upto order 3. First let us
consider the higher order partial derivatives.
functions.
78 www.AgriMoon.Com
Derivative of Higher Order
Solution:
We have
79 www.AgriMoon.Com
Derivative of Higher Order
Hence
So
80 www.AgriMoon.Com
Derivative of Higher Order
i.e., .
then
Solution:
But is discontinuous at .
81 www.AgriMoon.Com
Derivative of Higher Order
Solution:
i.e., .
So
82 www.AgriMoon.Com
Derivative of Higher Order
Solution:
1. For , find ,
dy
4. Find the derivatives of the functions represented implicitly
dx
2 2 2
(i ) sin ( xy ) − e − x y = 0 (ii ) xe + ye − e = 0 (iii ) y = x (iv) x + y = a
xy 2 y x xy x y 3 3 3
5. If r = xφ ( x + y ) + yψ ( x + y ) , show that
83 www.AgriMoon.Com
Derivative of Higher Order
∂ 2r ∂ 2r ∂ 2r
− 2 + = 0.
∂x 2 ∂x∂y ∂y 2
1
6. If u = [φ (ax + y ) + φ (ax − y )] , show that
y
∂ 2u a 2 ∂ 2 ∂u
= ⋅ y .
∂x 2 y 2 ∂y ∂y
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
84 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 9
9.1 Introduction
Let us apply the Taylor formula for function of the variable assuming
to be constant.
85 www.AgriMoon.Com
Taylor's Expansion for Function of Two Variables
where ,
where ,
86 www.AgriMoon.Com
Taylor's Expansion for Function of Two Variables
87 www.AgriMoon.Com
Taylor's Expansion for Function of Two Variables
, becomes
Solution:
Where is given by
88 www.AgriMoon.Com
Taylor's Expansion for Function of Two Variables
π π
1. Expand z = sin x sin y in powers of ( x − ) and ( y − ) . Find the terms of
4 4
the first and second orders and R2 (the remainder of second order).
R2 .
89 www.AgriMoon.Com
Taylor's Expansion for Function of Two Variables
1
e ax sin by = by + abxy + [(a 3 x 3 − 3ab 2 xy 2 ) sin (bθy ) + (3a 2bx 2 y − b 3 y 3 ) cos(bθy )]e aθx .
6
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
90 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 10
10.1 Introduction
function.
point if
91 www.AgriMoon.Com
Maximum and Minimum of Function of Two Variables
The maximum and minimum of a function are called extrema of the function;
we say that a function has an extremum of a given point if it has a maximum or
minimum at the given points.
at .
Solution:
of .
92 www.AgriMoon.Com
Maximum and Minimum of Function of Two Variables
This result is not sufficient for investigating the extreme points, but permits
finding these values for cases in which we are sure of the existence of a
maximum or minimum. Otherwise more investigation is required.
Solution:
which (or does not exist) (or does not exist) are called critical
Let
93 www.AgriMoon.Com
Maximum and Minimum of Function of Two Variables
or
94 www.AgriMoon.Com
Maximum and Minimum of Function of Two Variables
Solution:
References
95 www.AgriMoon.Com
Maximum and Minimum of Function of Two Variables
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
96 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 11
11.1 Introduction
without a constraint is a well known problem in calculus. One would normally use
the gradient to find critical points (gradient ( ) vanishes). Then check all
Example 1.
97 www.AgriMoon.Com
Lagrange's Multiplier Rule / Constrained Optimization
Second Derivative Test: The Second derivative test determines the optimality of
stationary point according to the following rules:
∂2 f ∂2 f ∂2 f ∂f ∂f
Let= A,= B= , 2 C , and = = 0 at the point ( x, y ) , then
∂x 2
∂x∂y ∂y ∂x ∂y
98 www.AgriMoon.Com
Lagrange's Multiplier Rule / Constrained Optimization
3. If AC − B 2 < 0 at ( x, y ) , then ( x, y ) is a
saddle point of f .
4. If AC − B 2 =
0 , further investigation is
required.
is .
stationary points found above will not work. This new problem can be thought of
as finding extreme values of when the point is restricted to lie on the
99 www.AgriMoon.Com
Lagrange's Multiplier Rule / Constrained Optimization
surfaces touch each other,i.e , they have a common tangent for line.
This means that the surfaces, gradient vectors at that point are parallel, hence,
))
or . Then finding the gradient and Hessian as was done above will
100 www.AgriMoon.Com
Lagrange's Multiplier Rule / Constrained Optimization
subject to with
Example 11.3: Find the rectangle of parameter l which has maximum area i.e.,
Maximize subject to
101 www.AgriMoon.Com
Lagrange's Multiplier Rule / Constrained Optimization
Solution:
i.e., i.e., .
Example 11.4 Find the shortest distance from the point (1,0) to the parabola
, i.e., Minimize subject to .
102 www.AgriMoon.Com
Lagrange's Multiplier Rule / Constrained Optimization
Now
If then , from
Hence
Hence
i.e.,
103 www.AgriMoon.Com
Lagrange's Multiplier Rule / Constrained Optimization
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
104 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 12
12.1 Introduction
Definition 12.1: We say that the curve is convex downward bending up on the
interval if all points of the curve lie above the tangent at any point on the
interval. Or when the curve turns anti-clock wise we call it is convex downward
(concave upward) (see Fig. 1).
Definition: We say that a curve is convex upwards for bending down on the
interval if all points of the curve lie below the tangent at any point on the
105 www.AgriMoon.Com
Convexity, Concavity and Points of Inflexion
The curve has a point of inflexion at , at which the curve changes from convex
downward.
If is convex upward on .
If is convex doward on .
P •
106 www.AgriMoon.Com
Convexity, Concavity and Points of Inflexion
Example 12.1: Find the ranges of values of for which the curve
Solution:
107 www.AgriMoon.Com
Convexity, Concavity and Points of Inflexion
Example 12.2: Determine the intervals where the graph of the function is
convex downward and convex upward of
Solution:
Hence,
Example 12.3: Determine the intervals where the graph of the function is
convex downward and convex upward of ,
Solution:
108 www.AgriMoon.Com
Convexity, Concavity and Points of Inflexion
Solution:
, . if , or . i.e.,
Example 12.5: What conditions must the coefficients satisfy for the
Solution:
is positive. i.e. .
109 www.AgriMoon.Com
Convexity, Concavity and Points of Inflexion
graph of .
7.
8.
10.Pint of inflexion of .
References
110 www.AgriMoon.Com
Convexity, Concavity and Points of Inflexion
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
111 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 13
Curvature
13.1 Introduction
Curvature measures the extent to which a curve is not contained in a straight line.
It curvature measures how curved the curve is. We have heard the comparison of
bending or curvature of a road at two of its points. The curvature of a straight line
is zero. It also measures how fast the tangent vector turns as a point moves along
the curve.
Fig.1.
Let be a fixed point on the curve. Let arc , and arc , so that
arc . Let , be the angles which the tangents at and make
with some fixed line (say - axis). denotes the angle formed by these tangents.
The symbol also denotes the angle through which the tangent turns from and
through a distance . will be large or small, as compared with ,
depending the degree of the sharpness of the bend. This suggests the following
definitions:
112 www.AgriMoon.Com
Curvature
∆y
A
N
S
P ∆x
Fig. 2.
113 www.AgriMoon.Com
Curvature
As , we have
Hence , where ,
is +ve or -ve i.e., accordingly as the curve is convex downward or convex upward.
But we consider is +ve here. Curvature is zero at point of inflection. Since is
independent of the choice of -axis and -axis, interchanging and , we see that
, is given by
114 www.AgriMoon.Com
Curvature
Given , . .
and
We know
115 www.AgriMoon.Com
Curvature
numerator becomes
denominator becomes
Hence
116 www.AgriMoon.Com
Curvature
Solution:
Hence
* We know
numerator becomes
117 www.AgriMoon.Com
Curvature
denominator becomes
Hence
Ans.:
Ans.:
118 www.AgriMoon.Com
Curvature
and . So
When , .
Solution:
derivative at this point is positive, is the point which gives the smallest
radius of curvature.
Solution:
119 www.AgriMoon.Com
Curvature
implies
i.e. .
Solution:
which implies .
120 www.AgriMoon.Com
Curvature
Solution:
, , , . Using this
When ,
5. at the point
121 www.AgriMoon.Com
Curvature
7. at the point
8. Find the point of the curve at which the radius of curvature is minimum.
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
122 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 14
Asymptotes
14.1 Introduction
then we say that the function has an asymptote . A function can have more than
123 www.AgriMoon.Com
Asymptotes
M P ( x, y )
d δ
Fig. 1
Vertical Asymptotes
asymptotes one has to find values of such that when they are approached by
the function , the latter approaches infinity. Then the straight line is a
vertical asymptote.
124 www.AgriMoon.Com
Asymptotes
for as when .
, with .
equation is .
125 www.AgriMoon.Com
Asymptotes
, hence
So .
Solution:
126 www.AgriMoon.Com
Asymptotes
Next to find the asymptotes of the form , i.e., the inclined asymptote.
Solution:
127 www.AgriMoon.Com
Asymptotes
• When a function has an asymptote (and not all functions have them) the
function gets closer and closer to the asymptote as the input value to the
function approaches either a specific value a or positive or negative infinity.
Remember, the simplified rational function has cancelled any factors common to
both the numerator and denominator.
The first step is to cancel any factors common to both numerator and denominator.
In this case there are none.
The second step is to see where the denominator of the simplified function equals
0. implies .
The vertical line is the only vertical asymptote for the function. As the
input value to this function gets closer and closer to -1 the function itself looks
128 www.AgriMoon.Com
Asymptotes
Example 14.7
First simplify the function. Factor both numerator and denominator and cancel any
common factors.
The vertical line is the only vertical asymptote for this function. As the input
value to this function gets closer and closer to 3 the function itself looks more
Example 14.8 If
Factor both the numerator and denominator and cancel any common factors.
In this case there are no common factors to cancel.
129 www.AgriMoon.Com
Asymptotes
5. Horizontal Asymptotes
Horizontal asymptotes occur when either one of the following conditions is met
(you should notice that both conditions cannot be true for the same function).
• The degree of the numerator is less than the degree of the denominator. In this
case the asymptote is the horizontal line .
• The degree of the numerator is equal to the degree of the denominator. In this
case the asymptote is the horizontal line where is the leading
When the degree of the numerator is greater than the degree of the denominator
there is no horizontal asymptote.
Example 14.9
then there is a horizontal asymptote at the line because the degree of the
130 www.AgriMoon.Com
Asymptotes
This means that as gets larger and larger in both the positive and negative
directions ( ) and ( ) the function itself looks more and more like the
horizontal line
Find the vertical asymptotes, horizontal asymptotes and inclined asymptotes for
each of the following functions Problems:
Exercises:
1.
2.
3. Ans.
4. Ans.
131 www.AgriMoon.Com
Asymptotes
5. Ans. No asymptotes
6. Ans.
7. Ans.
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
132 www.AgriMoon.Com
Module 1: Differential Calculus
Lesson 15
Tracing of Curves
15.1 Introduction
2. Critical points
3. Regions of increase
4. Regions of decrease
Suppose we have a function f defined for all sufficiently larger numbers. Then
we get substantial information concerning our function by investigating how it
behaves as x becomes large.
For example, sin x oscillates between -1 and +1 no matter how large x is.
133 www.AgriMoon.Com
Tracing of Curves
y = ax2 + bx + c, with a ≠ 0.
There are two essential cases, when a > 0 or a < 0. We have the parabola which
looks like in the figure
y = ax2 + bx + c y = ax2 + bx + c
a>0 a<0
y = f(x) - 3x2 + 5x – 1
5 1
f ( x)= x 2 (−3 + − ),
x x2
when is large positive or negative, then x2 is large positive and the factor on
the right is close to -3. Hence f(x) is large negative. This means that the
parabola has the shape as shown in figure.
134 www.AgriMoon.Com
Tracing of Curves
5
We have f ′( x) =−6 x + 5 . Thus f ′( x) = 0 iff x = . There is exactly one
6
2
5 5 25
critical point. We have f =−3 + −1 > 0
6 6 6
The critical point is a maximum, because we have already seen that the parabola
bends down.
−3 x 2 + 5 x − 1 =0
−5 ± 25 − 12 5 ± 13
=x =
−6 6
5− 13 5 5+ 13
6 6 6
(i) Looking at what happens when x becomes large positive or negative tells us
whether the parabola bends up or down.
135 www.AgriMoon.Com
Tracing of Curves
f ′( x)= 2ax + b= 0
−b
So x =
2a
Knowing whether the parabola bends up or down tells us whether the critical
−b
point is maximum or minimum, and the value x = tells us exactly where this
2a
critical point lies.
(iii) The points where the parabola crosses the x-axis are determined by the
quadratic formula.
2 1
x3 1 + 2 − 3 and, when x → + ∞ means f ( x) → + ∞
x x
x3 + 2 x − 1
Q( x) = 3
2x − x + 1
1
Here if x → ± ∞ , then Q( x) → .
2
136 www.AgriMoon.Com
Tracing of Curves
x3 − 1
Example 15.4(b) Consider the quotient Q( x) = 2
x +5
The meaning of the above limit is that there is no number which is the limit of
Q(x) as x → + ∞ or x → − ∞ .
1. If x → + ∞ then f ( x) → + ∞
If x → − ∞ then f ( x) → − ∞
2. We have f ′(=
x) 3x 2 − 2
f ′( x) =0⇔ x=± 2
3
3. Let g (=
x) f ′(=
x) 3 x 2 − 2 . Then the graph of g is a parabola which is given
as
− 2 2
3 3
Graph of g ( x) = f ′( x)
137 www.AgriMoon.Com
Tracing of Curves
4. f ′′( x) = 6 x , and f ′′( x) > 0 iff x > 0 and f ′′( x) < 0 iff for x > 0, therefore f
is bending up ( convex downward ) for x > 0 and bending down ( convex
upward ) for x < 0. There is an inflection point at x = 0.
Putting all this together, we find that the graph of f looks like this
graph of f(x) = x3 – 2x + 1
2
− 3
y = - x3 + 3x - 5
1. When x = 0 , we have y = -5. With general polynomial for degree ≥ 3 there is
in general no simple formula for those x such that f(x) = 0, so we do not give
explicitly in the intersection of the graph with the x – axis.
2. The derivative is f ′( x) =
−3 x 2 + 3
138 www.AgriMoon.Com
Tracing of Curves
f ′( x) =0⇔ x=±1
The graph of f ′( x) is given by
Putting all this information together, we see that graph of f looks like this
-1 1
graph of f(x) = - x3 + 3x + 5
139 www.AgriMoon.Com
Tracing of Curves
Solution:
Here we have f ′( x=
) 12 x 2 > 0 ∀ x ≠ 0 . There is only one critical point, when
x = 0. Hence the function is strictly increasing for all x, and its graph looks like
f ′′( x) =
24 x > 0 for all x > 0
Convex downward
convexup
Solution:
f ′( x)= 3 x 2 + 4 > 0 ∀ x
f ′′( x) =
6 x > 0 for x > 0
140 www.AgriMoon.Com
Tracing of Curves
Convex downward
Convex upward
x −1
=y f=
( x)
x +1
1. When x = 0, we have f(x) = 1. When x = 1, f(x) = 0.
2
2. The derivative is f ′( x) =
( x + 1)
2
141 www.AgriMoon.Com
Tracing of Curves
f(x)
x=0
f(x)=1
x=1
f(x)=0
−4
6. The second derivative is f ′′( x) = .
( x + 1)
3
There is no inflection point since f ′′( x) ≠ 0 for all x where the function is
defined. If x < -1, (x + 1)3 < 0, and f ′′( x) > 0 , f(x) is bending up or convex
downward. If x > -1, then x+1 > 0 ⇒ (x+1)3 > 0. So f ′′( x) < 0 i.e., f(x) is
bending down (convex upward).
1
x 1 −
x −1
( x) = lim =
x
7. As x → ∞ , f(x) → 1 f= 1
x + 1 x→∞ 1
x 1 +
x
when x → −∞ , f(x) → 1
142 www.AgriMoon.Com
Tracing of Curves
EXERCISES
Sketch the following curves, indicating all the information stated in the
examples etc.
x2 + 2
1. y =
x −3
x −3
2. y =
x2 + 1
3. =
y x4 + 4 x
4. =
y x8 + x
5. f ( x) = x 4 + 3 x 3 − x 2 + 5
x2 − 1
6. y =
x
References
143 www.AgriMoon.Com
Tracing of Curves
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
144 www.AgriMoon.Com
Module 2: Integral calculus
Lesson 16
Improper Integral
16.1 Introduction
Integral with infinite limits. Let a function f(x) be defined, positive and continuous
for all values of x such that a ≤ x < ∞ . Consider the integral
b
I (b) = ∫ f ( x) dx
a
a b x
O
Fig. 1
This integral is meaningful for b > a. This integral varies with b and is continuous
function of b. Let us consider the behavior of this integral when b → +∞ (Fig. 1).
Definition 16.1 if there exists a finite limit
b
lim ∫ f ( x) dx
b→∞ a
145 www.AgriMoon.Com
Improper Integral
Then this limit is called the improper integral of the function f(x) on the interval
[ a, +∞ ] and is denoted by the symbol
+∞
∫a
f ( x) dx
+∞ b
∫ a
f ( x) dx = lim ∫ f ( x) dx
b→∞ a
b
In this case it is said that the improper integral exists or converges. If ∫a
f ( x) dx
+∞
as b → +∞ does not have a finite limit, one say that ∫a
f ( x) dx does not exist or
diverges.
a a
∫−∞
f ( x) dx = lim
α →−∞ α∫ f ( x) dx
146 www.AgriMoon.Com
Improper Integral
+∞ c c
∫=
f ( x) dx
−∞ ∫−∞
f ( x) dx + ∫
−∞
f ( x) dx
The latter equation should be understood as if each of the improper integrals on the
right exists, then, by definition, the integral on the left also exists (converges).
+∞ dx
Example 16.1: Evaluate the integral ∫ 0 1 + x2
Solution:
+∞
dx b dx b
π
∫0 1 + x 2 b→∞ ∫0 1 + x 2 b→∞
−1
= =
lim =
lim tan x
0 2
Note that this integral expresses the area of an infinite curvilinear trapezoid crosses
x –axis as x → ∞ .
+∞ dx
Example 16.2: Evaluate ∫−∞ 1 + x 2
Solution:
+∞ dx dx
0 +∞ dx
∫−∞ 1 + x 2
= ∫−∞ 1 + x 2 ∫0 1 + x 2
+
π
The 2nd integral is equal to (see example 1)
2
147 www.AgriMoon.Com
Improper Integral
π
lim tan x= lim ( tan −1 0 − tan −1 =
b)
0
0dx 0 dx
∫−∞ 1 + = lim ∫ −1
=
x 2 b→−∞ b 1 + x 2 b→−∞ b b→−∞ 2
+∞dx π π
Hence, ∫−∞ 1 + x 2
= + =π
2 2
Theorem 16.1: Let f and g be continuous function on the interval [a, ∞) with
o ≤ f ( x) ≤ g ( x) ∀ a ≤ x < ∞ .
+∞ +∞
If ∫a
g ( x) dx converges then ∫
a
f ( x) dx also converges, and
+∞ +∞
∫a
f ( x) dx ≤ ∫
a
g ( x) dx
follows:
148 www.AgriMoon.Com
Improper Integral
c c −ε
∫
a
f ( x) dx = lim ∫
ε →0 + a
f ( x) dx
If the limit on the right exists, the integral is called an improper convergent
integral, otherwise it is divergent. If the function f(x) is discontinuous at x = a of
the interval [a,c] then by definition ,
c c
∫ a
f ( x) dx = lim ∫
ε →0 a +ε
f ( x) dx
If the function f(x) is discontinuous at some point x = x0 inside the interval [a,c] ,
we put
c x0 c
∫=
a
f ( x) dx ∫
a
f ( x) dx + ∫ f ( x) dx
x0
If both the improper integrals on the right hand side of the equation exist.
1 dx
Example 16.3 Evaluate ∫0
1− x
1 dx 1−ε dx
Solution: ∫0
= lim ∫
1 − x ε →0 0 1 − x
dx
1−ε
(1 − x ) 2 dx
−1
= lim ∫
ε →0 0
1−ε
(1 − x )
1
2
= −
− 12 + 1
0
149 www.AgriMoon.Com
Improper Integral
1−ε
=
−2 1 − x
0
=
lim − 2
ε →0
( ε −1 = 2 )
1dx
Example 16.4: Evaluate the integral ∫−1 x 2 .
Solution:
Since inside the interval of integration there exist a point x = 0, at which the
integrand is not continuous, we express the integration as:
1 dx 0−ε dx 1 dx
∫−1 x 2 ε →0 ∫−1 x 2 ε →0 ∫0+ε x 2
= lim + lim
−ε dx 1 dx
= lim ∫
x 2 ε →0 ∫ε x 2
+ lim
ε →0 −1
−ε 1
1 1
=lim − − lim
ε →0 x ε →0 x
−1 ε
1 1
=− lim − + 1 − lim 1 −
ε →0
ε ε →0 ε
1 1
But − lim − + 1 =∞ and − lim 1 − =
∞ i.e., the integral diverges on [-1,0]
ε →0
ε ε →0
ε
as well as on [0,1] .
Hence the given integral diverges on the entire interval [-1, 1].
150 www.AgriMoon.Com
Improper Integral
It should be noted that if we had evaluated the given integral without paying
attention to the discontinuity of the integrand at point x = 0, the result would have
dx −1
1
1 1
1
been wrong as ∫ 2 = =
− − = −2
−1 x x −1 1 −1
c dx c dx
∫a (c − x ) p
, also ∫
a ( x − a) p
c dx
It is easy to verify that ∫a (c − x ) p
converges for p < 1 and diverges for p ≥ 1.
EXERCISES
1 dx
1. ∫
0
1 − x2
∞
2. ∫
0
e − x dx
∞ dx
3. ∫
0 a + x2
2
151 www.AgriMoon.Com
Improper Integral
1 dx
4. ∫
0
1 − x2
1
5. ∫ ln x dx
0
π π
Ans.: 1. 1, 2. 1, 3. (a > 0) , 4. & 5. 1
2a 2
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
152 www.AgriMoon.Com
Module 2: Integral Calculus
Lesson 17
17.1 Introduction
b
lim ∫ f ( x) dx
b→∞ a
Then this limit is called the value of the improper integral of the function f(x)
on the interval [ a, +∞ ] and is denoted by the symbol
+∞
∫a
f ( x) dx
+∞ b
∫a
f ( x) dx = lim ∫ f ( x) dx
b→∞ a
153 www.AgriMoon.Com
Tests for Convergence
f(x), the x – axis and the ordinates x = a, x = b, it is natural to consider that the
+∞
improper integral ∫a
f ( x) dx expresses the area of an unbounded ( infinite )
a a
∫−∞
f ( x) dx = lim ∫
α →−∞ α
f ( x) dx
+∞ c c
∫=
f ( x) dx
−∞ ∫−∞
f ( x) dx + ∫
−∞
f ( x) dx
+∞ dx
Example 16.1 Find out at which p the integral ∫1 xp
converges and at which
it diverges.
a b
Solution:
Since (when p ≠ 1)
154 www.AgriMoon.Com
Tests for Convergence
b
dx 1 1− p 1
( b1− p − 1)
b
∫1 x p 1 − p
= = x
1− p
1
We have
∫1 x p b→+∞ 1 − p ( b − 1)
+∞ dx 1 1− p
= lim
Consequently, with respect to like this integral we conclude that if p > 1, then
+∞ dx 1
∫1
=
x p p −1
, and the integral converges.
+∞ dx
If p < 1, then ∫ 1 xp
= ∞ and integral diverges.
+∞ dx +∞
When p = 1, ∫1 x p
= ln x 1 = ∞ , and the integral diverges.
+∞ dx
Note: We call the p-integral ∫1 xp
converges for p > 1, and diverges for p ≤ 1
Theorem 17.1. Let f and g be continuous function on the interval [a, ∞) with
o ≤ f ( x) ≤ g ( x) ∀ a ≤ x < ∞ .
155 www.AgriMoon.Com
Tests for Convergence
+∞ +∞
If ∫ a
g ( x) dx converges then ∫ a
f ( x) dx also converges, and
+∞ +∞
∫a
f ( x) dx ≤ ∫
a
g ( x) dx
+∞ dx
Example 17.2 Investigate the integral ∫
1 x 2 (1 + e x )
for convergence.
Solution:
1 1
< 2
x (1 + e ) x
2 x
+∞
+∞ dx 1
And ∫1 x2
=
−
x1
=
1
+∞ dx
Consequently, ∫ 1 x 2 (1 + e x )
converges, and its value is less than 1. Hence
+∞ dx
∫1 x 2 (1 + e x )
converges.
+∞
Theorem 17.2. If for all x( x ≥ a ),0 ≤ g ( x) ≤ f ( x) holds true and ∫a
g ( x) dx
+∞
diverges, then the integral ∫
a
f ( x) dx also diverges.
Example 17.3 Find out whether the following integral converges or diverges.
+∞ x +1
∫1
x 3
dx
Solution:
156 www.AgriMoon.Com
Tests for Convergence
x +1 x 1
We note that > =
x3 x3 x
+∞ dx 1
But ∫
1
x
1
2
= ∞ as p=
2
< 2 . Hence the given integral is divergent.
+∞
Theorem17.3. If the integral ∫a
f ( x) dx converges, then the integral
+∞
∫a
f ( x)dx also converges.
+∞
Definition 17.1: An integral ∫
a
f ( x)dx converges conditionally if and only if
+∞ +∞
∫a
f ( x)dx converges but ∫ a
f ( x) dx is not convergent.
+∞ sin x
Example 17.3 Investigate the convergence of the integral ∫
1 x3
dx .
Solution:
sin x 1 +∞ 1
Here,
x3
≤
x3
. But ∫
1 x3
dx convergent as p = 3.
∞ sin x
Therefore, the integral ∫
1 x3
dx also converges.
157 www.AgriMoon.Com
Tests for Convergence
A function f(x) is defined and continuous when a ≤ x < c , and either not defined
or discontinuous when x = c. In this case, one cannot speak of the integral
c
∫a
f ( x) dx as limit of integral sums, because f(x) is not continuous on [a, c] and
c
The integral ∫a
f ( x) dx of the function f(x) discontinuous at a point c is defined
as follows:
c c −ε
∫
a
f ( x) dx = lim ∫
ε →0 a
f ( x) dx
If the limit on the right exists, the integral is called an improper convergent
integral, otherwise it is divergent. If the function f(x) is discontinuous at x = a of
the interval [a, c] then by definition,
c c
∫ a
f ( x) dx = lim ∫
ε →0 a +ε
f ( x) dx
If the function f(x) is discontinuous at some point x = x0 inside the integral [a, c]
, we put
c x0 c
∫=
f ( x) dx
a ∫a
f ( x) dx + ∫ f ( x) dx
x0
If both the improper integrals on the right hand side of the equation exist.
1 dx
Example 17.4 Test the convergence of the integral ∫−1 x 2 .
158 www.AgriMoon.Com
Tests for Convergence
Solution:
Since inside the interval of integration there exist a point x = 0, at which the
integrand is not continuous, we express the integration as:
1 dx 0−ε dx 1 dx
∫−1 x 2 ε →0 ∫−1 x 2 ε →0 ∫0+ε x 2
= lim + lim
−ε dx 1 dx
= lim ∫
x 2 ε →0 ∫ε x 2
+ lim
ε →0 −1
−ε 1
1 1
=lim − − lim
ε →0 x ε →0 x
−1 ε
1 1
=− lim − + 1 − lim 1 −
ε →0
ε ε →0 ε
1 1
But − lim − + 1 =∞ and − lim 1 − =
∞ i.e., the integral diverges on [-
ε →0
ε ε →0
ε
1, 0] as well as on [0,1] .
Hence the given integral diverges on the entire interval [-1, 1].
It should be noted that if we had evaluated the given integral without paying
attention to the discontinuity of the integrand at point x = 0, the result would
dx −1
1
1 1 1
have been wrong as ∫ 2 = =
− − = −2
−1 x x −1 1 −1
x=0
1
y= 2
x
159 www.AgriMoon.Com
Tests for Convergence
(-1, 1) (1, 1)
y=0
-1 -1
Fig. 3
Note: If the function f(x), defined on the interval [a, b], and has finite number
of discontinuity points a1, a2,…., an within the interval,
b a1 a2 b
Then ∫a
f (=
x) dx ∫
a
f ( x) dx + ∫ f ( x) dx + ..... + ∫ f ( x) dx
a1 an
If each of the improper integrals on the right side of the equation converges then
b
∫a
f ( x) dx is called convergent but if even one of these integrals diverges, then
b
∫a
f ( x) dx too is called divergent.
Theorem 17.3. Let f(x) and g(x) be continuous functions in [a,c] except at x = c
and at all points of this interval the inequalities g(x) ≥ f(x) hold and
c c
∫a
g ( x) dx converges, then ∫a
f ( x) dx also converges.
160 www.AgriMoon.Com
Tests for Convergence
Theorem 17.4. Let f(x) and g(x) be continuous functions on [a,c] except at x =
c and at all points of this interval the inequalities f(x) ≥ g(x) ≥ 0 hold and
c c
∫a
g ( x) dx diverges, then ∫ a
f ( x) dx also diverges.
c dx c dx
∫a (c − x ) p
, also ∫
a ( x − a) p
c dx
It is easy to verify that ∫ a (c − x ) p
converges for p < 1 and diverges for p ≥ 1.
1 dx
Example17.5 Does the integral ∫0
x + 4x 3
converge?
Solution:
1 1
Now ≤
x + 4 x3 x
dx 1 1 dx
∫
1
The improper integral ∫ 0
x
1
2
as
2
< 1 exists and hence
0
x + 4 x3
also
exists.
161 www.AgriMoon.Com
Tests for Convergence
EXERCISES
∞ dx
2. ∫
1
x
∞ dx
3. ∫−∞ x 2 + 2 x + 2
1 dx
4. ∫ 0
1
x3
2 dx
5. ∫ 0 x3
6. Let b > 2. Find the area under the curve y = e −2 x between 2 and b. Does this
area approach a limit when b → ∞ . If so what limit?
∞
7. Can an improper integral ∫
a
f ( x) dx ever be transformed onto a proper
3
Ans.: 1. The integral diverges, 2. The integral diverges, 3. π , 4. , 5. The
2
1 1 −4 1 1 1
integral diverges, 6. − e −2b + e , yes e −4 & 7. Yes, =
f ( x) = , x .
2 2 2 x2 t
References
162 www.AgriMoon.Com
Tests for Convergence
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
163 www.AgriMoon.Com
Module 2: Integral Calculus
Lesson 18
Rectification
18.1 Introduction
The method of finding the length of the arc of the curve of is called the
rectification. Let y = f ( x ) be a differentiable function defined on [a, b] with
a < b and assume that its derivative is continuous. Our aim is to determine the
length of the curve described by the graph. The main idea behind this is to
approximate the curve by small line segments and add these up.
(x1,f(x1)
(x2,f(x2) (x4,f(x4)
(x3,f(x3)
a = x0 x1 x2 x3 x4 = b
Fig .1
For each xi we have on the curve ( xi , f ( xi ) ) . We draw the line segments between
two successive points. The length of such a segments the length of the line
between
164 www.AgriMoon.Com
Rectification
(1)
( f ( xi+1)− f ( xi ) ) =( xi+1− xi ) f ′(ci )
n −1
165 www.AgriMoon.Com
Rectification
i.e., H (ci ) lies between the minimum and the maximum of on the
interval [ xi , xi +1 ] . Thus the sum we have written down lies between a lower sum
and an upper sum for the function H . We call such sums as Riemann sums.
This is true for every partition of the interval.
We know from basic integration theory that there is exactly one number lying
between every upper sum and every lower sum, and that number is the definite
interval. Therefore it is reasonable to define:
2
b dy b
∫a dx ∫a 1 + f ′( x) dx ---------(18.4)
= + = 2
1 dx
Similarly for x = φ ( y) and φ ′( y) are continuous on [a, b] , then the length of our
curve between a and b is
2
b dx b
∫a dy ∫a 1 + φ ′( y) dy
= + = 2
1 dy
Solution:
166 www.AgriMoon.Com
Rectification
3
As f , f ′( x) = x are both continuous on [0, 4], the length of the arc or length
1
2
2
2
3 1
4 4
9
of curve L = ∫
0 2
∫
1 + x 2 dx =
0
1 + x dx ,
4
9
Let 1 + x =
t , when x=0, t=1,
4
x=4, t=10
4 10 10
9 4 1 4 2 3 8 32
∫
0
1 + x dx = ∫ t 2 dt = × × t 2 =
4 91 9 3 1 27
10 − 1
Example 18.2 Find the length of the curve y=x2 between x =0 and x =1.
Solution:
1 1
0 0
1
12
Hence ∫ 1 + 4 x 2 dx = ∫ 1 + u 2 du ---------------------------- (18.5),
0 2 0
b
We can find the integral ∫ 1 + x 2 dx for b > 0, as
0
1 1
( ) ( ) ( )
2 −2
1
b + b + 1 + 2ln b + b 2 + 1 − b + b 2 + 1
2
4 2 2
167 www.AgriMoon.Com
Rectification
1 1
( ) ( ) 12 ( 2 + 5 )
2 2 −2
So ∫ 1 + u 2 du= 2+ 5 + 2ln 2 + 5 −
0 4 2
1 1
( ) ( ) ( )
2 1 −2
Hence (18.5) becomes: 2+ 5 + 2ln 2 + 5 − 2+ 5
8 2 2
There is one other way in which we can describe a curve. Suppose that we look
at a point which moves in the plane. Its coordinates can be given as a function
of time t. Thus, we get two functions of t, say
x = f(t), y = g(t),
We may view these as describing a point moving along a curve. The functions f
and g give the coordinates of the point as function of t.
=
Example 18.3 cosθ , y r sin θ . Then
Let, x r=
( cos θ ,sin θ )
θ
r =1
168 www.AgriMoon.Com
Rectification
This describes the motion of a bug around the circle with angular speed ω . Note
that the parametric representation of a curve is not unique. For
example x = r sin θ , y = r cosθ also represents a point on the circle.
and assume that both f, g have continuous derivatives. With eqn (18.4) it is very
reasonable to define the length of our curve (in parametric form) to be
b
= ∫ f ′(t ) 2 + g ′(t ) 2 dt .
b
l a
a
Observe that when a curve is given in usual form y = f(x) we can let
This shows how to view the usual form as a special case of the parametric form.
In that case g ′(t ) = 1 and the formula for the length in parametric form is seen to
be the same as the formula we obtained before for a curve y = f(x) .It is also
convenient to put the formula in the other standard notation for the derivative.
We have
169 www.AgriMoon.Com
Rectification
dx dy
= f ′(t ) and = y′(t )
dt dt
2 2
dx dy
b
= ∫ + dt
b
l a
a dt dt
2 2
dx dy
t
=s (t ) ∫
a
+ dt
dt dt
This gives
2 2
ds dx dy
= + = f ′(t ) 2 + g ′(t ) 2
dt dt dt
170 www.AgriMoon.Com
Rectification
2
ds dy
= 1+
dx dx
Example 18.4 Find the length of the curve x = cos t , y = sin t between
=t 0,=t π
Solution:
∫ ( sin t ) + ( cos t ) dt
2 2
=π
Example 18.5 Find the length of the curve x = et cos t , y = et sin t between t
=1 and t = 2.
Solution:
2 2
t ′ t ′
2
= ∫ ( e cos t ) + ( e sin t ) dt
2
l 1
1
2
=∫ (−e sin t + e cos t ) + (e cos t + e sin t ) dt
t t 2 t t 2
2
= ∫
1
(e 2t sin 2 t + e 2t cos 2 t + e 2t cos 2 t + e 2t sin 2 t ) dt
171 www.AgriMoon.Com
Rectification
2
= ∫ e dt
2= 2 e 2 − e
t
Solution:
dx
= 3cos 2 θ (− sin θ )
We have
dθ
dy
= 3sin 2 θ cosθ
dθ
Hence,
π
2
π
=l0 2
∫
0
9cos 4 θ + 9sin 4 θ cos 2 θ dθ
π
2
= 3∫ cos 2 θ sin 2 θ dθ
0
π
2
Hence
π π π
2
π 32 3 2
3
3∫ sin θ cosθ dθ =
2 ∫0
l0 =
2
sin 2θ dθ =
− cos 2θ =
0
4 0 2
EXERCISES
172 www.AgriMoon.Com
Rectification
1
=
1. y ln x , ≤x≤2,
2
2. y = 4 − x 2 , − 2 ≤ x ≤ 2 ,
1 x
=
3. y (e + e − x ) between x = 1 and x = -1
2
=
4. y ln cos x , 0 ≤ x ≤ π3 ,
9. Using exercise (9), find the length of the curve r = sin 2 θ2 from 0 to π .
4+2 5 17 + 4 4
Ans.: 1.
5
+ ln , 2. 2 17 + ln
1
(
, 3. e − , 4. ln 2 + 3 , )
2 1+ 5 17 − 4 e
5. 2π r , 6. 3, 7. 5, 8. 8 & 9. 2
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
173 www.AgriMoon.Com
Rectification
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
174 www.AgriMoon.Com
Module 2: Integral Calculus
Lesson 19
19.1 Introduction
xi xi +1
175 www.AgriMoon.Com
Volume and Surface of Revolution
n −1 n −1
∑ π f (ci )2 ( xi +1 − xi ) ≤ V ≤ ∑ π f (di )2 ( xi +1 − xi )
=i 0 =i 0
b
It is therefore reasonable to define this volume to be V = ∫ π f ( x)2 dx
a
a t1
x = a and x = b respectively.
Solution:
1
1
x3 2
∫ π (1 − x2
) dx =π ( x − ) =π
0 3 0 3
2 4
So the volume of full sphere is 2 × π = π
3 3
Example 19.2: Find the volume obtained by rotating the region between y = x3
and y = x in the first quadrant around the x − axis .
176 www.AgriMoon.Com
Volume and Surface of Revolution
y = x2
y = x2
0 0
1 1
= π ∫ x 2 dx −π ∫ x6 dx
0 0
π π
= −
3 7
1
Example 19.3: (Volume of Chimneys) . Consider the function f ( x) = .
x
a 1 b
a
177 www.AgriMoon.Com
Volume and Surface of Revolution
1
Let 0 < a < 1 . The volume of revolution of the curve y = between x = a and
x
1
dx
x = 1 is given by ∫ π = π ln x a = −π ln a ,
1
a
x
But we are interested to find finite volume for the infinite chimney.
1
Example 19.4: Compute the volume of revolution of the curve y = between
x4
a and 1. Find the limit as a → 0
Solution:
1
The volume of revolution of the curve y = between x = a and x = 1
x4
1
1
1 −1
1 1
is given by the integral ∫ π 1
π∫
dx =x 2 dx =×2 x2π 2π 1 − a
=
a a
x2 a
178 www.AgriMoon.Com
Volume and Surface of Revolution
Example 19.5 Find the volume of a cone whose base has a radius r , and a
height h , by rotating a straight line passing through the origin around the
x − axis
Solution:
y r
y= x
h r
x
h
r 1
The equation of the straight line is y = x . Slant height is y = . Hence the
h x2
2
h
r r2 h 2 π r 2 h3 1 2
volume of the cone is ∫ π =
x dx π= 2 ∫
x dx × =π r h
2
0 h h 0 h 3 3
a b
179 www.AgriMoon.Com
Volume and Surface of Revolution
S = ∫ 2π y 1 + dx
b
dy
a dx
The idea again is to approximate the curve by line segments. We use a partition
a = x0 ≤ x1 ≤ x2 ≤ x3....... ≤ xn = b
f ( xi +1 )
f ( xi )
fx( xi ) xfi+1
( xi +1 )
i
Li
Then Li = ( xi +1 − xi ) + ( f ( xi +1 )2 − f ( xi ) )
2 2
The length of a circle of radius y is 2π y . If we rotate the line segment about the
then the x− axis area of the surface of rotation will be between 2π f (ti ) Li and
2π f (si ) Li where f (ti ) and f (si ) are the minimum and maximum of f ,
respectively on the interval [ xi , xi +1 ] . This is illustrated on Fig 1.
180 www.AgriMoon.Com
Volume and Surface of Revolution
f ′(ci )( xi +1 − xi ) , ci ∈ ( xi , xi +1 )
f ( xi +1 ) − f ( xi ) =
Hence Li = ( xi +1 − xi ) + f (ci )2 ( xi +1 − xi )
2 2
1 + f ′(ci )2 ( xi +1 − xi )
=
n −1
Now take the sum ∑ 2π f (c ) i 1 + f ′(ci )2 ( xi +1 − xi )
i =0
This is a Riemann sum, between the upper and lower sums for the integral
b
=S ∫ 2π f ( x)
a
1 + f ′( x) 2 dx
Thus it is reasonable that the surface area should be defined by this integral, as
was to be shown.
Suppose that
181 www.AgriMoon.Com
Volume and Surface of Revolution
=x f (t=
), y g (t ), a ≤ t ≤ b
Then the length of Li between ( f (ti ), g (ti ) ) and ( f (ti +1 ), g (ti +1 ) ) is given by
= f ′(ci )2 + g ′(di )2 ( ti +1 − ti )
2 2
b dx dy
=S ∫ 2π y + dt
a dt dt
when t = x , this coincides with the formula found previously. It is also useful to
write this formula symbolically S = ∫ 2π yds
182 www.AgriMoon.Com
Volume and Surface of Revolution
= ( dx ) + ( dy )
2 2
ds
When using this symbolic notation, we don not put limits of integration. Only
when we use explicit parameter over an interval a ≤ t ≤ b we explicitly write the
surface area as
b
ds
S = ∫ 2π y dt
a dt
Example 19.6 We wish to find the area of a sphere for radius r > 0 .
Solution: we can view the sphere as the area of revolution of a circle for radius
r , and to express the circle in parametric form,
= θ , y r sin θ ,0 ≤ θ ≤ π
x r cos=
π
S = ∫ 2π r sin θ r 2 sin θ + r 2 cos θ dθ
0
π 2
= ∫ 2π r sin θ dθ
0
π
= 2π r 2 (− cos θ ) 0
= 4π r 2
183 www.AgriMoon.Com
Volume and Surface of Revolution
Exercises
limit ?
4 π2 π 2.54 π π π π
Ans.: 1. π r 3 , 2. π , 3. − , 4. , 5. π (e − 2) & 6. − 3 , yes:
3 8 4 3 24 3b 24
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
184 www.AgriMoon.Com
Volume and Surface of Revolution
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
185 www.AgriMoon.Com
Module 2: Integral Calculus
Lesson 20
Double Integration
20.1 Introduction
many other important interpretations (distance, volume, arc length, surface area,
moment of inertia, mass, hydrostatic pressure, work) depending on the nature and
interpretation of f .
186 www.AgriMoon.Com
Double Integration
In this Lesson we shall see that integrals of functions of two or more variables
which are called multiple integrals and defined I much the same way as integrals of
functions of single variable.
Double Integrals: Here we define the integral of a function f(x, y) of two variables
over a rectangular region in xy-plane. We then show how such an integral is
evaluated and generalize the definition to include bounded regions of a more
general nature.
c
a x
b
187 www.AgriMoon.Com
Double Integration
∆A =∆x∆y
Choose a point (xk, yk) in each piece of ∆Ak and from the sum
n
=Sn ∑ f ( x , y )∆A
k =1
k k k ------------------- (20.2)
n
=
Thus ∫∫
f ( x, y )dA lim
∆A→0
∑ f ( x , y )∆A
1
k k k ---------- (20.3)
R
As with functions of a single variable, the sums approach this limit no matter how
the interval [a, b] and [c, d] that determine R are subdivide, along as the lengths of
the subdivisions both go to zero. The limit (20.3) is independent of the order in
which the area ∆Ak are numbered, and independent of the choice of ( xk , yk )
within each ∆Ak . The continuity of f sufficient condition or the existence of the
double integral, but not a necessary one, and limit question exists for many
discontinuous functions also.
188 www.AgriMoon.Com
Double Integration
Like “single” integrals, we have the following properties for double integrals of
continuous functions which are useful in computations and applications.
(iv) ∫∫ f ( x, y)dA ≥ 0 if
R
f ( x, y ) ≥ 0 on R
∫∫=
f ( x, y )dA ∫∫ f ( x, y )dA + ∫∫ f ( x, y )dA
R1 R2 R1 R2
Volume: When f(x, y) > 0, we may interpret ∫∫ f ( x, y)dA as the volume of the
R
189 www.AgriMoon.Com
Double Integration
n
=Sn ∑ f ( x , y )∆A
k =1
k k k is the volume of a vertical rectangular prism y that
z
f ( xk , yk )
c d
• •
a• ∆Ak
b
•
approximate the volume of the portion of the solid that stands above the box -
∆Ak . The sum Sn thus approximates what we call the total volume of the solid,
and we define this volume to be
d b b d
∫∫ f ( x, y)dA
=
R
∫=
∫ f ( x, y)dxdy ∫ ∫ f ( x, y)dydx
c a a c
190 www.AgriMoon.Com
Double Integration
Fubbin’s theorem shows that double integrals over rectangles can be calculated as
iterated integrals. This means that we can evaluate a double integral by integrating
one variable at a time, using the integration techniques we already know for
function of a single variable.
Fubin’s theorem also says that we may calculate the double integral by integrating
in either order (a genuine convenience). In particular, when we calculate a volume
by slicing, we may use either planes perpendicular to the x-axis or planes
perpendicular to y-axis. We get same answer either way.
Even more important is the fact that Fubin’s theorem holds for any continuous
function f(x, y). In particular it may have negative values as well as positive values
on R, and the integrals we calculate with Fubin’s theorem may represent other
things besides volumes.
Example 20.1: Suppose we wish to calculate the volume under the plane z = 4-x-y
over the region R : 0 ≤ x ≤ 2, 0 ≤ y ≤ 1 in the xy – plane.
1 2
∫∫ f ( x, y)dA
R
= ∫ ∫ (4 − x − y )dx dy
0 0
191 www.AgriMoon.Com
Double Integration
2
1
x2
= ∫ 4 x − − xy dy
0 0
2
1
= ∫ (8 − 2 − 2 y ) dy
0
1
= ∫ ( 6 − 2 y ) dy
0
1
=6 y − y 2 =5
0
f ( x, y ) = 1 − 6 x 2 y and R : 0 ≤ x ≤ 2, − 1 ≤ y ≤ 1
1 2
∫∫ f ( x, y=
)dA ∫ ∫ (1 − 6 x 2
y )dx dy
R −1 0
1 2
= ∫ ( x − 2x y )
3
dy
−1 −1
1
= ∫ ( 2 − 16 y ) dy
−1
1
= 2 y − 8 y2
−1
= (2 − 8) − (−2 − 8) = 4
192 www.AgriMoon.Com
Double Integration
1 2 2 1
∫ ∫ (1 − 6 x ∫ y − 3x y
y )dy dx =
2 2 2
dx
−1 0 0 −1
2
= ∫ (1 − 3x ) − (−1 − 3 x 2 ) dx
2
2
= ∫ 1 − 3 y + 1 + 3 x 2 dx
2
= 2=
2
x 0 4.
The difficult part of evaluating a double integral can be finding the limits of
integration. But there is a procedure to follow:
If we want to evaluate over a region R, integrating first with respect to y and then
with respect to x, we take the following steps:
2. We integrate from the y-value where L enters R to the y-value where L leaves R
3. We choose x-limits that include all the vertical lines that pass through R
x =1 y= 1− x 2
193 www.AgriMoon.Com
Double Integration
y
Leaves where
y = 1- x2
Leaves where
1 y = 1- x
x
0 x 1
Smallest Biggest
x=0 x=1
1 =x 1− y 2
∫ ∫
0 x = 1− y
f ( x, y )dx dy
sin x
Example 20.4 Calculate ∫∫A x dA where A is the triangle in the xy-plane
bounded by the x-axis, the line y = x and the line y = 1.
x sin x 1
Solution: ∫ ∫ dy dx
00
x
sin x y = x
1
= ∫ y dx
x
0 y =0
∫ sin x dx =
= − cos x 0 =
− cos1+ .46
1
194 www.AgriMoon.Com
Double Integration
1 1
sin x sin x
∫∫
0 y
x
dx dy , we can’t evaluate it because we can’t express ∫ x
in terms of
elementary functions.
PROBLEM
Evaluate the following integrals and sketch the region over which each integration
takes place.
3 2
1. ∫∫
0 0
(4 − y 2 )dy dx
3 0
2. ∫∫
0 −2
( x 2 y − 2 xy )dy dx
π x
3. ∫ ∫
0 0
x sin y dy dx
π sin x
4. ∫ ∫
0 0
y dy dx
1 1
∫ ∫ ye xy dx dy
y
5. Find the value of the integral
10 0
2 2x
6. Sketch the region of integration of ∫ ∫ f ( x, y)dy dx and express the integral as
0 x2
Ans.: 1. 16, 2. 0, 3.
( 4 + π ) , 4. π , 5. 9 − 9e & 6.
2
2 4
195 www.AgriMoon.Com
Double Integration
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
196 www.AgriMoon.Com
Module 2: Integral Calculus
Lesson-21
Triple Integration
21.1 Introduction
The cells have dimensions ∆x by ∆y by ∆z . We number the cells that lie inside
n
=Sn ∑ F ( x , y , z )∆V
k =1
k k k k ------------- (21.1)
lim Sn = ∫∫∫ F ( x, y, z ) dV
D
We call this limit the triple integral of F over D. The limit also exists for some
discontinuous functions.
197 www.AgriMoon.Com
Triple Integration
Triple integrals share many algebraic properties with double and single
integrals. Writing F by F ( x, y, z ) and G for G ( x, y, z ) , we have the following
3. ∫∫∫ F dV ≥ 0 if F ≥ 0 in D
D
4. ∫∫∫ F dV ≥ ∫∫∫ G dV if F ≥ G on D
D D
5. ∫∫∫ F=
dV ∫∫∫ F dV + ∫∫∫ F dV + ..... + ∫∫∫ F dV
D D1 D2 Dn
The triple integral Evaluation is hardly evaluated directly from its definition as a
limit. Instead, one applies a three-dimensional version of Fubin’s theorem to
evaluate the integral by repeated single integrations.
198 www.AgriMoon.Com
Triple Integration
z = f2(x, y)
z = f1(x,y y)
C
x
R
Fig.12
f2 ( x, y )
∫∫∫ F ( x , y , z ) dV = ∫∫ ∫
F ( x , y , z ) dz dy dx
D R f1 ( x , y )
f2 ( x, y )
or ∫∫∫ F ( x, y, z )dV = ∫∫ ∫ F ( x, y, z )dz dy dx -------------- (21.1)
D R f1 ( x , y )
If we omit the parenthesis .The z-limits of integration indicate that for every
(x, y) in the region R, z may extend from the lower surface z = f1(x, y) to the
upper surface z = f2(x, y). The y – and x - limits of integration have not given
explicitly in Eq (21.1) but are to be determined in the usual way from the
boundaries of R.
199 www.AgriMoon.Com
Triple Integration
Leaves D at z = f2(x, y)
Enters D at z = f1(x, y)
R
(x, y)
Fig. 12
Example 21.1 Find the volume enclosed between the two surfaces z = x2+3y2
and
z = 8-x2-y2.
x2+3y2 = 8-x2-y2
or x2+2y2 = 4
which is elliptic .
200 www.AgriMoon.Com
Triple Integration
2 (4− x 2 )/2 8− x 2 − y 2
V= ∫ ∫ ∫ dz dy dx
−2 − (4− x 2 )/2 x +3 y 2 2
2 (4− x 2 )/2
= ∫ ∫ (8 − x − y − x − 3 y ) dy dx
2 2 2 2
−2 − (4− x 2 )/2
2 (4− x 2 )/2
= ∫ ∫ (8 − 2 x − 4 y ) dy dx
2 2
−2 − (4− x 2 )/2
2 (4− x 2 )/2
= ∫ ∫ (
2 8 − 2 x 2 − 4 y 2 dy dx )
−2 0
2 (4− x 2 )/2
4
= ∫ 2 (8 − 2 x 2 ) y − y 3 dx
−2
3 0
( ) dx
2 3
8
∫ 2(8 − 2 x )
(4− x 2 ) (4− x 2 )
= −
2 2
2 2
−2
3
3
2
( )
2
4 2
= ∫
3 −2
4 − x 2 dx
3
2
( )
2
8 2
=
3 0∫ 4 − x 2 dx
= 8π 2
As we know, there are sometimes two different orders in which the single
integrations that evaluate a double integral may be worked (but not always). For
triple integral there are sometimes (but not always) as many as six workable
orders of integration. The next example shows an extreme case in which all six
are possible.
Example 21.2 Each of the following integrals gives the volume of the solid
shown
201 www.AgriMoon.Com
Triple Integration
in Fig 3.
z
y + z =1
y
Fig 3.
1 1 2 1 1− y 2
(a ) ∫∫∫ dx dy dz (b) ∫ ∫ ∫ dx dz dy
x
0 0 0 0 0 0
1 1 1− z 2 1 1− z
(c) ∫∫ ∫ dy dx dz (d ) ∫∫ ∫ dy dz dx
0 0 0 0 0 0
1 2 1− y 2 1 1− y
(e) ∫∫ ∫ dz dx dy ( f ) ∫∫ ∫ dz dy dx
0 0 0 0 0 0
EXERCISES
1. Write six different iterated triple integrals for the volume of the rectangular
solid in the first octant bounded by the co-ordinate planes and the planes x = 1,
y = 2,
z = 3. Evaluate one of the integrals.
2. Write six different intersected triple integrals of the volume in the first octant
enclosed by the cylinder x2 + z2 = 4 and the plane y = 3. Evaluate one of the
integrals.
202 www.AgriMoon.Com
Triple Integration
3. Write an iterated triple integrals in the order dz dy dx for the volume of the
region bounded below by the xy-plane and above by the paraboloid z = x2+ y2
and lying inside the cylinder x2 + y2 = 4.
1 1 1− y
4. Rewrite the integral ∫ ∫ ∫ dz dy dx as an equivalent integrated integral in the
−1 x 2 0
order.
a) dy dz dx b) dy dx dz c) dx dy dz d) dx dz dy e) dz dx dy
1 2 3 21 3 3 21 31 2 2 31
Ans.: 1. ∫∫∫ dz dy dx , ∫∫∫ dz dx dy , ∫∫∫ dx dy dz , ∫∫∫ dy dx dz , ∫∫∫ dx dz dy ,
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
∫∫ ∫ dz dx dy , ∫∫ ∫ dz dy dx , ∫ ∫ ∫ dy dz dx ,
0 0 0 0 0 0 0 0 0
2 1− z 2 3 3 2 4− z 2 2 3 1− x 2
12 π .
2 1− x 2 x2 + y 2
3. 4 ∫ ∫ ∫ dz dy dx & 4.
0 0 0
References
203 www.AgriMoon.Com
Triple Integration
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
204 www.AgriMoon.Com
Module 2: Integral Calculus
Lesson 22
Area & Volume using Double and Triple Integration
22.1 Introduction
We have seen if we take f (x, y) = 1 in the definition of the double integral over
a region in Eqn (20.2), is the partial sum reduce to
n n
Sn = ∑ f ( xk , yk )∆Ak = ∑ ∆Ak ,
=k 1=k 1
and give area of the region as n → ∞ . In that case ∆x, ∆y approach zero. In
this case we define the area on a rectangular region R to be the limit
= lim ∑ ∆A=
Area k ∫∫ dA (22.1)
R
Example 22.1 Find the area of the region R bounded by y = x and y = x2 in the
first quadrant.
1 x 1
∫ ∫ dy =
dx ∫ ( x − x 2
)dx
2
0 x 0
6
x 2 x3 1
= − =
2 3 0 6
Example 22.2 Find the area of the region R enclosed by the parabola y = x2 and
the line y = x + 2.
205 www.AgriMoon.Com
Area & Volume using Double and Triple Integration
Solution: x 2 = x + 2 ⇒ x 2 − x − 2 = 0
x 2 − 2 x=
+ x − 2 0 i.e., x( x − 2) +=
1( x − 2) 0
( x + 1)( x − 2) =
0
x = −1, 2
2 x+2
Hence the area A = ∫ ∫ dy dx
−1 x 2
2 x+2
=∫ y dx
−1 x2
2
= ∫ −1
( x + 2 − x 2 )dx
2
x2 x3
= + 2x −
2 3 −1
8 1 1
= 2+ 4− − −2+
3 2 3
8 1 1
=2+4− − − +2
3 2 3
16 + 3 + 2
= 8−
6
7 9
=8 − =
2 2
Solution:
y
For order of integration reversed, draw a horizontal lin L2. It enters at x = ,
2
leaves at x = y . To include all such lines we let y to n from y = 0 to y=
4. The integral is
206 www.AgriMoon.Com
Area & Volume using Double and Triple Integration
4 y
∫∫
0
y
2
f ( x, y )dx dy
If F (r ,θ ) ≡ 1 the constant function whose value is one, then the value over R is
the areas of R (which agrees our earlier definition). Thus
Area of R = ∫∫ r dr dθ
R
π 4
-π 4
207 www.AgriMoon.Com
Area & Volume using Double and Triple Integration
π
2 a cos 2 2θ
π r = 2 a cos 2 2θ
4 4
r2
∫ ∫ r dr dθ = ∫ 2
dθ
− π4 0 − π4 r =0
π
4
= ∫a cos 2θ dθ
2
π
−4
π
2 4
a
= sin 2θ
2 − π4
a2
= [1 − (−1)]
2
= a2
As ∆ x , ∆ y , ∆ z all approaches zero, the cells ∆Vk become smaller and we need
n
Volume of D = lim ∑ ∆Vk =
∫∫∫ dV . k =1 D
208 www.AgriMoon.Com
Area & Volume using Double and Triple Integration
The triple integral Evaluation is hardly evaluated directly from its definition as a
limit. Instead, one applies a three-dimensional version of Fubin’s theorem to
evaluate the integral by repeated single integrations.
Z
. (0,0,c)
. Y
(a,0,0) . (0,b,0)
X
Fig. 4
Fig. 4 shows a system of mutually orthogonal coordinates axes OX, OY, OZ.
The Cartesian coordinates of a point P(x, y, z) in the space may be read from the
coordinates axes by passing planes through P perpendicular to each axis. The
points on the x-axis have their y- and z- ordinates both zero. Points in a plane
perpendicular to the z-axis, say, all have the same z - coordinate. Thus of the
points in the plane perpendicular to the z- axis and 5 units above the xy-plane
all have coordinates of the form (x,y,5) . We can write z = 5 as an evaluation for
this plane. The three planes x = 2, y = 3, z = 5 intersect in the point P (2, 3, 5).
The points of the yz- plane are obtained setting x = 0. The three coordinates
planes x = 0, y = 0, z = 0 divide the space into eight cells, called octants. The
octant in which all three coordinates are positive is the first octant, but there is
no conventional numbering of the remaining seven octants.
209 www.AgriMoon.Com
Area & Volume using Double and Triple Integration
Solution: The points all horizontal plane z = 3, and in this plane they lie in this
cirle x2 + y2 = 4.Thus we may describe the set of the circle in the plane x2 + y2 =
4 in the plane z = 3.
x = r cosθ =
r x2 + y 2
y
y = r sin θ tan θ =
x
z=z
Spherical coordinates are useful when there is a center of symmetry that we can
take as the origin. The spherical coordinates ( ρ ,ϕ ,θ ) are shown the first
φ ρ P (ρ ,θ ,φ )
y
θ
x
210 www.AgriMoon.Com
Area & Volume using Double and Triple Integration
The second spherical coordinate φ , is the angle measured down from the z-axis
to the line OP. The equation ρ = constant describes cone with vertex at O, axis
OZ and generating angle φ , provide we broaden our interpretation of the word
π
“cone” to include the xy- plane for which φ = and cones the generation
2
π
angles greater than .
2
z=z θ =θ z = ρ cosθ
∫∫∫ dx dy dz
Volume : = ∫∫∫
= dz rdr dθ ∫∫∫ ρ sin θ d ρ dφ dθ
2
211 www.AgriMoon.Com
Area & Volume using Double and Triple Integration
Exercises
1. Find the area of the region R enclosed by the parabola y = x2 and the line y =
x +1
3. Find the volume of the solid in the first octant bounded by the paraboloid.
z = 36 – 4x2 – 9y2
4. Find the volume of the solid enclosed between the surfaces x2 + y2 = 92 and
x2 + z2 = 92 .
x y z
5. The volumes of the tetrahedron bounded by the plane + + =1 and the
a b c
coordinate planes.
7. The volume of the wedge cut from the cylinder x2 + y2 = 1 and the plane z =
y above and plane below.
8. The volume of the region in the first octant bounded by the coordinate
planes, above by the cylinder x2 + z = 1 and on the right by the paraboloid y =
x2 + z2
16a 3 1 2 2 2
Ans.: 1. ,2. , 3. 27π , 4. , 5. abc , 6. , 7. & 8.
3 6 3 3 7
212 www.AgriMoon.Com
Area & Volume using Double and Triple Integration
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
213 www.AgriMoon.Com
Module 2: Integral Calculus
Lesson 23
Gamma Function
Definition 23.1 The Gamma Function: The gamma function is defined by the
improper integral
+∞
Γ ( λ + 1) =
∫ e t dt ------------- (23.1)
−t λ
0
To deduce some of the properties of the gamma function, let us integrate Eq.
(23.1) by parts:
+∞ R
∫0
e − t t λ dt = lim
R →+∞ 0 ∫ e − t t λ dt
= lim −e − t t λ + λ ∫ e − t t λ dt
R R
R →+∞
0 0
− Rλ +∞
= lim R + 0 + λ ∫ e − t t λ −1 dt
R →+∞
e 0
+∞
= λ ∫ e − t t λ −1 dt
0
214 www.AgriMoon.Com
Gamma Function
∞ ∞
λ ∫ e − t dt =
Γ(1) = −e − t
0
=
1
0
Γ(2) =1.Γ(1) =1
Γ(3) =2.Γ(2) =2!
Γ(4) =3.Γ(3) =3! ------------ (23.3)
Γ(λ + 1) =
λ ! ------------ (23.4)
215 www.AgriMoon.Com
Gamma Function
x
-4 -3 -2 -1 1 2 3 4
+∞
∫ e − x dx = π
2
1
Theorem 23.1. 2
0
∫∫ <∫∫ < ∫∫
D1 S D2
(24.1.5)
y (R,R)
S
D1
D2
O (R,0) (R 2 ,0) x
Fig 2.
216 www.AgriMoon.Com
Gamma Function
Now evaluate these integrals by iterate integrals, the centre one in rectangular
coordinates, and other two in polar coordinates:
R π R R R 2 π
∫ e − r r dr ∫ dθ < ∫ e − x dx ∫ e − y dy < ∫ e − r r dr ∫ dθ
2 2 2 2 2 2
0 0 0 0 0 0
π
( ) (∫ ) π
( )
R 2
− R2 − x2
1− e < < 1 − e −2 R
2
e dx
4 0 4
(∫ ) π
+∞ 2
e − x dx =
2
0 4
+∞ π
∫ e − x dx =
2
i.e.,
0 2
1
Theorem 23.2. Γ =π
2
1 +∞ +∞
∫ ∫
1 −1 −1
Now, =
Γ e − t=
t 2 dt e − t t 2 dt
2 0 0
+∞
∫ e dy π set t = y2
−y
= 2=
2
+∞
∫
5
Example 23.1 Evaluate the integral x 4 e − x dx
0
+∞ +∞
∫ x 4 e − x dx = 2 ∫ t 2 e − t dt = 2Γ( 92 )
5 7
0 0
217 www.AgriMoon.Com
Gamma Function
7 5 3 3 105 π
Γ(=
9
2) . . Γ(=
2)
2 2 2 8 2
+∞ 105 π 105 π
∫
5
x 4 e − x dx =×
2 × =
0 8 2 8
Solution: We have
r r r r
f=
(r ) + 1 + 2 ..... + (n − 1)h h n
h h h h
Γ ( hr + 1) Γ ( hr + 2 ) Γ ( hr + n )
=h n
....
Γ ( hr ) Γ ( hr + 1) Γ ( hr + n − 1)
Γ ( hr + n ) n
= .h
Γ ( hr )
r
obtained by the recursion Eq. 2 with λ =
h
Some special cases of the result of Example 2 are interesting. For particular
case, set r = 1 and h = 2. Then
218 www.AgriMoon.Com
Gamma Function
2n Γ(n + 12 )
1.3.5....(2n − 1) = 1
Γ( 2 )
1 1 π
But Γ( 2 ) =
Γ( 32 ) = .
2 2
Hence
2n Γ(n + 12 )
1.3.5....(2n − 1) =
π
However,
2.4.6.....2n
=
1.3.5....(2 n − 1) 1.3.5....(2n − 1)
2.4.6.....2n
(2n)!
=
2n n !
(2n)! π (2n)!
Γ ( n + 12=
) × = π
n
2 n! 2 n 2 n
2 n!
for n = 1, 2, ……
+∞
23.3. Γ ( x ) r
Theorem = x
∫e
− rt x −1
t dt , r > 0, x > 0
0
219 www.AgriMoon.Com
Gamma Function
+∞
Theorem 23.4. Γ ( x ) = ∫
2 e − t t 2 x−1 dt
2
Proof: Set t2 = y
Extension of definition
Γ ( x + n)
Γ ( x) = , −n < x < −n + 1
x( x + 1)( x + 2).....( x + n − 1)
Thus we have defined Γ ( x ) for all x except x = 0,-1,-2,…. Observe that when
n =1 the right hand side of (6) depends on the values of Γ ( x ) in the interval
0<x<1. It is clear that Γ ( x ) has been defined for negative x in such a way that
equation
Example 4. Compute Γ ( 12 )
Γ ( − 12 + 1) =− 12 Γ ( − 12 )
220 www.AgriMoon.Com
Gamma Function
i.e., Γ ( 12 ) =− 12 Γ ( − 12 )
i.e., Γ ( − 12 ) =−2 π
Exercise
+∞
∫ x 2 e − x dx
2
2.
0
+∞
∫x
−4
3. e − x dx
0
+∞
∫ (1 − x) e − x dx
3
4.
0
+∞
∫x e − x dx
3
5.
0
+∞
∫e
−t x
6. Show that the improper integral t dt converges for x > −1 and diverges
0
for x ≤ −1 .
1
dx
7. Compute ∫
0 x ln ( 1x )
∫
8. Evaluate 2−9 x dx using gamma function ( Hint : 2−9 x = e −9 x ln 2 )
2 2 2
221 www.AgriMoon.Com
Gamma Function
π
Ans.: 1. Γ ( 32 ) or , 2. 6, 3. ∞ , 4. -9394, 5. 2 × 7!, 6. ,7. 2π & 8.
2
1 π
6 ln 2
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
222 www.AgriMoon.Com
Module 2: Integral Calculus
Lesson 24
24.1 Introduction
Definition 2.2
For x, y positive we define the Beta function by
Theorem 24.1.
To evaluate the Beta function we usually use the Gamma function. To find their
relationship, one has to do a rather complicated calculation involving change of
variables (from rectangular into tricky polar) in a double integral.
When x and y are positive integers, it follows from the definition of the gamma
function that:
223 www.AgriMoon.Com
The Beta Function
∞ t x −1
β ( x, y ) = ∫0+ x + y dt
( )
1 + t
224 www.AgriMoon.Com
The Beta Function
Proof : When x and y are arbitrary positive numbers, the proof proceeds as
follows. From the double integral of the nonnegative function t 2 x −1u 2 y −1e − t −u 2
2
over the three regions D1 , D2 and S of figure 1 of Lesson 23. Now, however, t
and u are the variables, however, t and u are the variables x and y positive
constants. We have relation (23.5) of Lesson 23 as before. Again we evaluate
the central double integral by iteration in rectangular coordinates: the other two,
in polar coordinates:
π 2 x + 2 y −1 2 y −1
π 2 x + 2 y −1
< ∫0 cos θ sin θ dθ ∫0 e r
2 2 x −1 2 y −1 R 2 −r2
dr
Now, if we let R become infinite and use Theorems 23.4 and 24.3, we obtain
1 1 Γ( x) Γ( y )
B( y, x) =Γ( x + y ) ,0 < x , 0 < y
2 2 2 2
∫
1
Example 24.1 Evaluate 0
x 4 (1 − x)3dx
Γ(5)Γ(4) 1
∫0 x (1 − x) dx = ∫0 x (1 − x=
) 4−1dx B= =
1 4 3 5 −1 1
Solution: (5,4)
Γ(9) 280
1
∫
1
Example 24.2 Evaluate dx
0 3
x 2 (1 − x)
225 www.AgriMoon.Com
The Beta Function
1 −1 2
Γ( 13 )Γ( 23 )
∫0
−1
β
1
∫ − = = =Γ( 13 )Γ( 23 )
1
dx =
1 3
Solution: x (1 3
x ) dx ( ,
1 2
)
Γ(1)
3 3
0 3
x 2 (1 − x)
∫
1
Example 24.3 Evaluate 0
x .(1 − x)dx
3
Γ ( 32 ) Γ ( 2 )
∫0 x (1 − x) dx = β ( 2 ,2 ) =
−1
∫
1 1
x (1 − x)dx= 2 −1
2
3
Solution: 0
Γ ( 72 )
Γ ( 32 ) =
1
2 π
Γ ( 52 ) =
3
4 π
Γ ( 72 ) =
15
8 π
4
∫
1
Thus 0
x .(1 − x)dx =
15
Proof: We know,
for
= Γ(q )Γ(1 − q )
Solution: Let
226 www.AgriMoon.Com
The Beta Function
Exercises
1.
2.
3.
4.
5.
6.
7.
8.
9.
References
Widder, D.V. (2002). Advance Calculus 2nd Edition, Publishers, PHI, India.
227 www.AgriMoon.Com
The Beta Function
Piskunov, N. (1996). Differential and Integral Calculus Vol I, & II, Publishers,
CBS, India.
Suggested Readings
228 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 25
Introduction
where x(n) stands for the nth derivative of unknown function x(t) with respect to the
independent variable t. For example
5
d4 x d2 x dx
4
+ 2 + = et (25.2)
dt dt dt
dx
= x + sin x. (25.3)
dt
The order of a differential equation is referred to the highest order derivative involved in
the differential equation. For example, the order of the differential Equation (25.2) is four.
The degree of a differential equation is the degree of the highest order derivative which
occurs in it; after the differential equation has been made free from radicals and fractions
as far as derivatives are concerned, e.g. in differential Equation (25.2), the degree is one.
229 www.AgriMoon.Com
Introduction
A differential equation is called linear if (a) every dependent variable and every derivative
involved occurs in first degree only, and (b) no product of dependent variables and/or
derivatives occur. A differential is not linear is called nonlinear. For examples, Equation
(25.2) is linear and (25.3) is nonlinear.
Any relation between the dependent and independent variables, when substituted in the
differential equation, reduces it to an identity is called a solution of differential equation.
For example, y = e2x is a solution of y ′ = 2y .
25.2.1 Example
y = A/x + B. (25.5)
y ′ = −A/x2 . (25.6)
y ′′ = 2A/x3 . (25.7)
230 2 www.AgriMoon.Com
Introduction
Let
An n-parameter family of curves is a set of relations of the form {(x, y) : f (x, y, c1, c2 , ..., cn ) =
0}, where f is real valued function of x, y, c1, c2 , ..., cn and each ci (i = 1, 2, ...n) ranges over
an interval of real values.
Suppose we are given a family of curves containing n arbitrary constants. Then by differ-
entiating it successively n times and eliminating all arbitrary constants from the (n + 1)
equations we obtain an nth order differential equation whose solution is the given family
of curves. We now illustrate the procedure of forming differential equations with the help
of some examples.
25.4.1 Problem 1
Find the differential equation of the family of curves y = emx , where m is an arbitrary
constant.
Solution: We have the family of curves
y = emx . (25.9)
231 3 www.AgriMoon.Com
Introduction
y ′ = memx . (25.10)
Now, we eliminate m from (25.9) and (25.10) and using m = loge y , we obtain the required
differential equation as
y ′ = y loge y.
25.4.2 Problem 2
x2 + y 2 = a2 . (25.11)
x + yy ′ = 0,
25.4.3 Problem 3
Obtain the differential equation satisfied by xy = aex + be−x + x2 , where a and b are an
arbitrary constant.
Solution: Given family of curves
xy ′′ + 2y ′ = (xy − x2 ) + 2, (25.14)
232 4 www.AgriMoon.Com
Introduction
Remark: From the above examples we observed that the number of arbitrary con-
stants in a solution of a differential equation depends upon the order of the differential
equation and is the same as its order. Hence a general solution of an nth order differential
equation will contain n arbitrary constant.
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
Piskunov, N. (1996). Differential and Integral Calculus (Volume - 2). First Edition. CBS
Publisher, Moscow.
233 5 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 26
In this lesson we present solution techniques of differential equations of first order and
first degree. We shall mainly discuss differential equation of variable separable form,
homogeneous equations and equations reducible to homogeneous form.
There are two standard forms of differential equations of first order and first degree,
namely,
dy
= f (x, y) or Mdx + Ndy = 0
dx
Here M and N are functions of x and y , or constants. We discuss here some special forms
of these equations where exact solution can easily be obtained.
If in a differential equation, it is possible to get all the functions x and dx to one side and
all the functions of y and dy to the other, the variables are said to be separable. In other
words if a differential equation can be written in the form F (x)dx + G(y)dy = 0, we say
variables are separable and its solution is obtained by integrating the equation as
Z Z
F (x)dx + G(y)dy = c,
26.2.1 Problem 1
dy
Solve = ex+y + x2 ey .
dx
Solution: For separating variables, we rewrite the given equation as
234 www.AgriMoon.Com
Differential Equation of First Order
−e−y = ex + x3 /3 + c,
26.2.2 Problem 2
26.3.1 Example
dy
Solve = sec(x + y).
dx
Solution: Let, x + y = v so that
dy dv
= − 1. (26.1)
dx dx
Using (26.1), the given differential equation becomes
dv
= sec v + 1. (26.2)
dx
235 2 www.AgriMoon.Com
Differential Equation of First Order
1 2 cos2 (v/2) − 1
dx = dv ⇒ dx = dv
sec v + 1 1 + 2 cos2 (v/2) − 1
A differential equation of first order and first degree is said to be homogeneous if it can
be put in the form
dy
= f (y/x).
dx
These equations can be solved by letting y/x = v and differentiating with respect to x as
dv dv
v+x = f (v) ⇒ x = f (v) − v.
dx dx
Then, separating variables, we have
dx dv
=
x f (v) − v
26.4.1 Example
236 3 www.AgriMoon.Com
Differential Equation of First Order
Solution: Since the right hand side of the given equation is function of y/x alone, the
given problem is homogeneous equation. Substituting y/x = v so that
dy dv
=v+x (26.3)
dx dx
the given equation becomes
dv dx cos v
v+x = v + tan v → = dv
dx x sin v
Integrating and substituting the value of v , we get the solution as
y
cx = sin ,
x
where c is an arbitrary constant.
where X, Y are new variables and h, k are constants to be chosen so that the resulting
equation in X, Y becomes homogeneous. From above we have dx = dX , and dy = dY , so
that dy/dx = dY /dX . Now the given differential equation in new variables becomes
dX aX + bY + (ah + bk + c)
= ′ (26.5)
dY a X + b′ Y + (a′ h + b′ k + c′ )
In order to make (26.5) homogeneous, the constant h and k must satisfy the following
algebraic equations
ah + bk + c = 0 , a′ h + b′ k + c′ = 0 (26.6)
237 4 www.AgriMoon.Com
Differential Equation of First Order
X = x − h, Y = y − k. (26.8)
26.5.1 Example
238 5 www.AgriMoon.Com
Differential Equation of First Order
Thus, we get
Suggested Readings
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
Weir, M.D., Hass, J. and Giordano, F.R. (2005). Thomas Calculus. Eleventh Edition.
Pearson Education. New Delhi.
McQuarrie, D.A. (2009). Mathematical Methods for Scientist and Engineers. First Indian
Edition. Viva Books Pvt. Ltd. New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
Piskunov, N. (1996). Differential and Integral Calculus (Volume - 2). First Edition. CBS
Publisher, Moscow.
239 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 27
In this lesson we shall learn linear differential equations of first order. Such equations
are very often used in applications. Solution strategies of solving such equations will be
discussed. Further a another special form of differential equation which can be reduced
to linear differential equation of first order will be studied.
A first order differential equation is called linear if it can be written in the form
dy
+ P (x)y = Q(x) (27.1)
dx
where P and Q are constants or function of x only.
A method of solving (27.1) relies on multiplying the equation by a function called inte-
grating function so that the left hand side of the differential equation can be brought under
a common derivative. Suppose R(x) is an integrating factor of the (27.1). Multiplying the
(27.1) by R(x), we obtain
dy
R(x) + P (x)R(x)y = Q(x)R(x) (27.2)
dx
Suppose, we wish that the L.H.S of (27.2) is the differential coefficient of some product.
dy
Clearly, the term R(x) dx can only be obtained by differentiating the product R(x)y(x). In
other words, we wish to have
dy d
R(x) + P (x)R(x)y(x) = (R(x)y(x)). (27.3)
dx dx
This implies
dy dy dR
R(x) + P (x)R(x)y(x) = R(x) + y(x) .
dx dx dx
On cancelling the first term on both the sides we obtain
dR dR
P (x)R(x)y(x) = y(x) ⇒ = Rdx.
dx R
240 www.AgriMoon.Com
Linear Differential Equation of First Order
R
Integrating the above equation, we get log R = P dx. Note that the constant of integration
is not important here because the integrating factor will be used to multiplying both the
sides of the differential equation and therefore it will be cancelled. Thus, an integrating
factor (I.F.) of the differential Equation (27.1) is
R
P dx
R=e (27.4)
or
R Z R
P dx P dx
ye = Qe dx + c,
which is required solution of given differential equation. Here C is the constant of inte-
gration.
27.2.1 Problem 1
dy
Solve x cos x + y(x sin x + cos x) = 1, 0 < x < π/2.
dx
Solution: We rewrite the given equation as
dy 1 sec x
+ (tan x + ) = .
dx x x
An I.F. of the given differential equation is
R
(tan x+ x1 )dx
e = elog x sec x = x sec x.
or
yx sec x = tan x + c,
where, c is an arbitrary constant.
241 2 www.AgriMoon.Com
Linear Differential Equation of First Order
27.2.2 Problem 2
dy
Solve (1 + x2 ) = x(1 − y).
dx
Solution: Rewriting the given differential equation in standard form
dy x x
+ 2
y= .
dx 1 + x 1 + x2
The I.F. is R x
dx 1 2 p
I.F. = e 1+x2 = e 2 ln(1+x ) = 1 + x2 .
The solution is
x
p Z
2
y 1+x = √ +c ⇒ y = 1 + c (1 + x2 )−1/2
1 + x2
Here c is an arbitrary constant.
dv/dx + P v = Q, (27.6)
which is linear in v and x and its solution can be obtained with the help of I.F. as before.
R
Thus, we have an I.F.= e pdx and the solution is
R Z R
pdx P dx
ve = Qe dx + c.
27.4.1 Problem 1
dy
Solve cos y + 2x sin y = x.
dx
242 3 www.AgriMoon.Com
Linear Differential Equation of First Order
dy dv
Solution: Substitution sin y = v which implies cos y dx = dx reduces the given differential
equation to
dv
+ 2xv = x.
dx
R 2
This is a linear differential equation of first order and its I.F. is e 2xdx = ex . The solution
of the equation in v is given by
1
Z
x2 2 2
ve = xex dx + c ⇒ v = + ce−x .
2
Replacing v by sin y we get the required solution as
−1 1 −x2
y = sin + ce .
2
27.4.2 Problem 2
dy
Solve + x sin 2y = x3 cos2 y
dx
Solution: Dividing the given differential equation by cos2 y , we obtain
dy
sec2 y + 2x tan y = x3 .
dx
dy dv
Putting tan y = v so that sec2 y dx = dx . Hence the above equation becomes
dv
+ 2xv = x3 ,
dx
2
which is linear. Its I.F. is ex and its solution is given as follows
Z
x2 2
ve = ex x3 dx + c,
2 1 2
vex = (x2 − 1)ex + c.
2
Replacing v by tan y we obtain the required solution.
dy/dx + P y = Qy n (27.7)
243 4 www.AgriMoon.Com
Linear Differential Equation of First Order
where P and Q are constants or function of x only and n is constant except 0 and 1 is
called Bernoulli differential equation. This equation can easily be solved by multiplying
both sides by y −n as
dy 1 dv
Setting y 1−n = v , so that y −n dx = (1−n) dx , the Equation (27.8) becomes
27.5.1 Example
dy
Solve x + y = y 2 ln x .
dx
dv 1
− v = x−1 ln x. (27.10)
dx x
R 1
x dx 1
The I.F. of the differential Equation (27.10) is e− = x, and hence the solution be-
comes
1
Z
v =− x−2 log xdx + c
x
or by replacing v by y −1 we get
y −1 = 1 + ln x + cx,
244 5 www.AgriMoon.Com
Linear Differential Equation of First Order
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
Weir, M.D., Hass, J. and Giordano, F.R. (2005). Thomas Calculus. Eleventh Edition.
Pearson Education. New Delhi.
McQuarrie, D.A. (2009). Mathematical Methods for Scientist and Engineers. First Indian
Edition. Viva Books Pvt. Ltd. New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
Piskunov, N. (1996). Differential and Integral Calculus (Volume - 2). First Edition. CBS
Publisher, Moscow.
245 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 28
If M and N are functions of x and y , the equation Mdx + Ndy = 0 is called exact when
there exists a function f (x, y) such that
or equivalently
∂f ∂f
dx + dy = Mdx + Ndy.
∂y ∂x
28.1.1 Theorem
to be exact is
∂M ∂N
= . (28.2)
∂y ∂x
Proof: First we proof that the condition (28.2) is necessary. To prove we let the Equation
(28.1) to be exact. Then, by definition, there exists f (x, y) such that
∂f ∂f
dx + dy = Mdx + Ndy. (28.3)
∂y ∂x
246 www.AgriMoon.Com
Exact Differential Equation of First Order
247 2 www.AgriMoon.Com
Exact Differential Equation of First Order
28.2.1 Problem 1
Therefore
∂M ∂N
= −4x − 4y =
∂y ∂x
Hence, the given equation is exact and hence there exists a function f (x, y) such that
∂f ∂f
d(f (x, y)) = dx + dy = Mdx + Ndy
∂x ∂y
which implies
∂f ∂f
= M(x, y) and = N(x, y)
∂x ∂y
Integration of the first of above equations with respect to x gives
1
f = x3 − 2x2 y − 2y 2 x + c1 (y)
3
where c1 (y) is an arbitrary function of y only. Differentiating the above f with respect to
∂f
y and using = N(x, y) we get
∂y
∂f
= −2x2 − 4xy + c′1 (y) = +y 2 − 4xy − 2x2
∂y
This implies
y3
c′1 (y) = y 2 ⇒ c1 (y) = + c2
3
Hence the solution is given by
f (x, y) = c3 ⇒ x3 − 6xy(x + y) + y 3 = c
248 3 www.AgriMoon.Com
Exact Differential Equation of First Order
28.2.2 Problem 2
Determine whether the differential equation (x + sin y)dx + (x cos y − 2y)dy = 0 is exact
and solve it.
Solution: For given equation we have
Now we check
∂M ∂N
= cos y =
∂y ∂x
Hence the given differential equation is exact. For the solution we seek a function f (x, y)
so that
∂f ∂f
= (x + sin y) and = (x cos y − 2y)
∂x ∂y
From the first relation we get
x2
f (x, y) = + x sin y + c1 (y)
2
Differentiating w.r.t. y and using the second relation of (28.6) we get
Therefore, we have
x2
f (x, y) = + x sin y − y 2 + c2
2
Then the solution of the given differential equation
x2
f (x, y) = c3 ⇒ + x sin y − y 2 = c.
2
28.2.3 Problem 3
249 4 www.AgriMoon.Com
Exact Differential Equation of First Order
So the equation is exact. Then, there exists a function f (x, y) such that
∂f ∂f
= (2y 2x − 2y 3 ) and = (4y 3 − 6y 2x + 2yx2 )
∂x ∂y
This gives
∂f
f (x, y) = (y 2 x2 − 2xy 3 ) + c1 (y) ⇒ = (2yx2 − 6xy 2) + c′1 (y)
∂y
This implies
c′1 (y) = 4y 3 ⇒ c1 (y) = y 4 + c2
28.2.4 Problem 4
Solve that the differential equation (3xy + y 2)dx + (x2 + xy)dy = 0. is not exact and hence
it cannot be solve by the method discussed above.
Solution: For the given differential equation we have
∂M ∂N
= 3x + 2y, and = 2x + y;
∂y ∂x
∂M ∂N
Since = , the given equation is not exact.
∂y ∂x
Now we see that it cannot be solved by the procedure described previously where we seek
a function f such that
∂f ∂f
= 3xy + y 2 and = x2 + xy (28.7)
∂x ∂y
250 5 www.AgriMoon.Com
Exact Differential Equation of First Order
This provides
1
c′1 (y) = − x2 − xy
2
Since the right side of the above depends on x as well as on y , it is impossible to solve this
equation for c1 (y). Thus there is no f (x, y) exists and hence the given differential equation
cannot be solved in this way.
Suggested Readings
McQuarrie, D.A. (2009). Mathematical Methods for Scientist and Engineers. First Indian
Edition. Viva Books Pvt. Ltd. New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
Piskunov, N. (1996). Differential and Integral Calculus (Volume - 2). First Edition. CBS
Publisher, Moscow.
251 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 29
In general, equations of the type M(x, y)dx + N(x, y)dy = 0 are not exact. However, it
is sometimes possible to transform the equation into an exact differential equation multi-
plying it by a suitable function I(x, y). That is, if I(x, y) is an integrating factor then the
differential equation
becomes exact. A solution to the above equation is obtained by solving the exact differ-
ential equation as in the previous lesson. Note that the given equation may have several
integrating factors. This is exactly the procedure we have used for solving linear differ-
ential equations in earlier lesson. Here we deal with more general differential equation.
There is not much theory behind finding integrating factor by inspection. This method
works based on recognition of some standard exact differentials that occur frequently in
practice. The following list of exact differentials would be quite useful in solving exact
differential equations:
29.1.1 Example
252 www.AgriMoon.Com
Exact Differential Equations: Integrating Factors
xy 2 + x = cy
1
If the equation Mdx+Ndy = 0 is homogeneous and Mx+Ny 6= 0, then I(x, y) =
(Mx + Ny)
is an integrating factor. In order to prove the result, we need to show that
Mdx + Ndy
= d (some function x and y)
Mx + Ny
Given that M(x, y) and N(x, y) are homogeneous functions of some degree n, i.e., M(tx, ty) =
tn M(x, y) and N(x, y) = tn N(x, y). Then
x 1 1 1 x
M ,1 =M x, y = n M(x, y) ⇒ M(x, y) = y n M ,1
y y y y y
253 2 www.AgriMoon.Com
Exact Differential Equations: Integrating Factors
Similarly, we get
n x
N(x, y) = y N ,1
y
Now consider
x x x x x
(Mx − Ny) y nxM − y,1 y nyN y,1 yM y,1
x −N y,1
= = =f
(Mx + Ny) n x n x
y xM y , 1 + y yN y , 1 x x x y
yM y,1 + N y,1
Rewriting f (x/y) = f (exp(ln(x/y))) and defining g(x) := f (exp(x)), the above equation
becomes
Mdx + Ndy 1 x
= d (ln(xy)) + g (ln(x/y)) d ln
Mx + Ny 2 y
1
Thus is an integrating factor of the homogenous differential equation Mdx +
Mx + Ny
Ndy = 0.
29.2.1 Example
Solve the differential equation (x2 y − 2xy 2 )dx − (x3 − 3x2 y)dy = 0
Solution: The given equation is a homogeneous differential equation. Comparing it with
Mdx + Ndy = 0, we have M = x2 y − 2xy 2 and N = −(x3 − 3x2 y). Since
254 3 www.AgriMoon.Com
Exact Differential Equations: Integrating Factors
x − 2y ln x + 3y ln y = cy
1
If the equation Mdx + Ndy = 0 is of the form f1 (xy)ydx + f2(xy)xdy = 0, then
(Mx − Ny)
is an integrating factor provided Mx − Ny 6= 0. Similar to rule II we now show that
Mdx + Ndy
= d (some function x and y)
Mx − Ny
f1 (xy) + f2 (xy)
Let f (xy) := and g(x) := f (exp(x)), the above equation reduces to
f1 (xy) − f2 (xy)
Mdx + Ndy 1 x 1 x
= f (xy)d (ln xy) + d ln = g(ln xy)d (ln xy) + d ln
Mx − Ny 2 y 2 y
255 4 www.AgriMoon.Com
Exact Differential Equations: Integrating Factors
29.3.1 Example
f1 (xy)ydx + f2 (xy)xdy = 0
and we have
Mx − Ny = xy(x2 y 2 + 2) − xy(2 − 2x2 y 2 ) = 3x3 y 3 6= 0
This is an exact differential equation which can be solved with the technique discussed in
previous lesson.
Now we discuss the most general approach of finding integrating function. The idea is to
multiply the given differential equation
by a function I(x, y) and then try to choose I(x, y) so that the resulting equation
If a function I(x, y) satisfying the partial differential Equation (29.4) can be found, then
(29.3) will be exact. Unfortunately, solving Equation (29.4), is as difficult to solve as the
original Equation (29.2) by some other methods. Therefore, while in principle integrating
factors are powerful tools for solving differential equations, in practice they can be found
256 5 www.AgriMoon.Com
Exact Differential Equations: Integrating Factors
only in special cases. The cases we will consider are: (i) an integrating factor I that is
either as function of x only, or (ii) a function of y only.
Let us determine necessary conditions on M and N so that (29.2) has an integrating factor
I that depends on x only. Assuming that I is a function of x only, then Equation (29.4)
reduces to
dI dI IMy − INx
IMy = INx + N ⇒ = (29.5)
dx dx N
If (My − Nx )/N is a function of x only, say f(x), then there is an integrating factor I
R
that also depends only on x which can be found by solving (29.5) as I(x) = e f (x)dx . A
similar procedure can be used to determine a condition under which Equation (29.2) has
an integrating factor depending only on y . To conclude, we have:
R
1 ∂M ∂N
If − is function of x alone say f (x), then I(x) = e f (x)dx is an I.F.
N ∂y ∂x
1 ∂N ∂M R
If − is function of y alone say f (y), then I(y) = e f (y)dy is an I.F.
M ∂x ∂y
29.5.1 Problem 1
M = (x2 + y 2 + x) and N = xy
29.5.2 Problem 2
257 6 www.AgriMoon.Com
Exact Differential Equations: Integrating Factors
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
McQuarrie, D.A. (2009). Mathematical Methods for Scientist and Engineers. First Indian
Edition. Viva Books Pvt. Ltd. New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Piskunov, N. (1996). Differential and Integral Calculus (Volume - 2). First Edition. CBS
Publisher, Moscow.
258 7 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 30
In this lesson we discuss linear differential equation of higher order with constant coeffi-
cients. In particular, we shall learn about the techniques of finding solutions of homoge-
nous equations. Different cases will be considered with the help of several examples.
In a linear differential equation, the dependent variable and its differential coefficients
occur only in the first degree and are not multiplied together. The general form of the
equation is
dn y dn−1 y dn−2 y
+ a1 (x) + a2 (x) + . . . + an (x)y = F (x), (30.1)
dxn dxn−1 dxn−2
where a1 , a2 , . . . , an and F are either constants or functions of x only. If the right hand
side, i.e. F (x), is identically zero, the equation is said to be homogeneous; otherwise it is
called nonhomogeneous. Before we discuss some particular cases of the above equation
we state two facts about the solution of a linear homogeneous differential equation. The
first says that if we know n solutions y1 , y2 , . . . , yn of the linear homogeneous equation,
then any linear combination y = c1 y1 + c2 y2 + . . . cn yn is also a solution for any constants
c1 , c2 , . . . , cn . This can easily be proved by substituting y = c1 y1 + c2 y2 + . . . cn yn into
the equation and using linearity of the equation. The second important result concerns
about the general solutions (solution containing all solutions) to the linear homogeneous
equation. This result says that any solution is some linear combination of y1 , y2 , . . . , yn for
some suitable values of constants c1 , c2 , . . . , cn . However, this is not true for any combina-
tion of solutions but is true if the solutions y1 , y2 , . . . , yn are linearly independent.
259 www.AgriMoon.Com
Linear Differential Equation of Higher Order
where a1 , a2 , . . . , an are constants, is called linear differential equation with constant coef-
dn
ficients. Using the symbols Dn := , the Equation (30.2) becomes
dxn
(D n + a1 D n−1 + a2 D n−2 + ... + an )y = F (x), (30.3)
Further defining f (D) := Dn + a1 Dn−1 + a2 Dn−2 + ... + an , we can rewrite the given
differential equation in a more compact form as f (D)y = F (x). Here f (D) acts as oper-
ator on y to yield F (x). The general solution of (30.2) can be written as the sum of the
general solution of the corresponding homogeneous equation, refereed as complimentary
function (C.F.), and a particular solution or sometimes called particular integral (P.I) of
nonhomogeneous equation. Thus
y = C.F. + P.I. (30.4)
Note that the C.F. involves n arbitrary constants and P.I. does not involve any arbitrary
constant. It is readily evident that y in (30.4) is the general solution of the given non-
homogeneous differential equation because it satisfies the given differential equation as
f (D) (C.F. + P.I.) = f (D) (C.F.) + f (D) (P.I.) = 0 + F (x) and it has n arbitrary constants.
To solve Equation (30.5), we seek a function which satisfies the above equation. One
intelligent guess of such a function is the exponential function emx , where m is a constant.
Differentiations of this exponential function are just constant multiples of the original
exponential. If we substitute this function into the Equation (30.5), we obtain
(mn + a1 mn−1 + a2 mn−2 + ... + an )emx = 0 (30.6)
Since the exponential function is never zero, we can divide this last equation by emx .
Thus, y = emx is a solution to Equation (30.5) if and only if m is a solution to the algebraic
equation
mn + a1 mn−1 + a2 mn−2 + ... + an = 0 (30.7)
Equation (30.7) is called the auxiliary equation (A.E.) or characteristic equation (C.E.) of
the differential Equation (30.5).
260 2 www.AgriMoon.Com
Linear Differential Equation of Higher Order
If m1 , m2 , m3 , ..., mn be real and distinct then the solutions em1 x , emm x , . . . , emn x are lin-
early independent and the general solution of the given homogeneous differential equation
becomes
y = c1 em1 x + c2 em2 x + c3 em3 x + ... + cn emn x ,
30.4.1 Example
Find the general solution of the differential equation (D3 + 6D2 + 11D + 6)y = 0.
Solution: The A.E. is (m3 + 6m2 + 11m + 6) = 0. The roots are m = −1, −2, −3. Hence
the required solution is y = c1 e−x + c2 e−2x + c3 e−3x .
Let m1 = m2 are repeated roots of the A.E. Then, we have n − 1 linearly independent
solutions. It can be shown that a simple choice y = xem1 x is also a solution which is
independent to the rest n − 1 solutions. Thus, the general solution of the given differential
equation is given by
y = (c1 + c2 x)em1 x + c3 em3 x + ... + cn emn x
The above idea can be further extended by taking solutions xem1 x , x2 em1 x , . . . , xl−1 em1 x . . .
if the root m1 is repeating l−times.
30.5.1 Example
The roots of the A.E. are m = 1, 1, −2, −2. The required solution is y = (c1 + c2 x)ex + (c3 +
c4 x)e−2x .
261 3 www.AgriMoon.Com
Linear Differential Equation of Higher Order
y = c′1 eαx (cos βx + i sin βx) + c′2 eαx (cos βx − i sin βx) + c3 em3 x + ... + cn emn x .
Defining new constants c1 = c′1 + c′2 and c2 = i(c1 − c2 ), the general solution becomes
Similar to the case II, the solution for repeated complex roots can be found, see example
below.
30.6.1 Example
This is the case of repeated complex root, so case II and case III can be combined to give
the desired solution as y = (c1 + c2 x) cos x + (c3 + c4 x) sin x.
30.7.1 Problem 1
Find the general solution of the differential equation (D3 + 3D2 + 3D + 1)y = 0.
Solution: The A.E. and its root are given by (m + 1)3 = 0 and m = −1, −1, −1. Therefore,
the required solution is y = (c1 + c2 x + c3 x2 )e−x .
262 4 www.AgriMoon.Com
Linear Differential Equation of Higher Order
30.7.2 Problem 2
30.7.3 Problem 3
30.7.4 Problem 4
30.7.5 Problem 5
Find the general solution of the differential equation (D2 + 1)3 (D2 + D + 1)2 y = 0.
Solution: The A.E. of given equation is
(D 2 + 1)3 (D 2 + D + 1)2 = 0,
√ √
1 3 1 3
The roots are m = ±i, ±i, ±i, − ± i , − ± i . Therefore, the desired solution is
2 2 2 2
263 5 www.AgriMoon.Com
Linear Differential Equation of Higher Order
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
McQuarrie, D.A. (2009). Mathematical Methods for Scientist and Engineers. First Indian
Edition. Viva Books Pvt. Ltd. New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
Piskunov, N. (1996). Differential and Integral Calculus (Volume - 2). First Edition. CBS
Publisher, Moscow.
264 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 31
As we have seen in the earlier lesson that a general nonhomogeneous linear differential
equations with constant coefficients can be written in operator form as f (D)y = F (x).
The operator, 1/f (D) is called inverse operator which gives a particular integral when
operated on both the sides of the given differential equation. Hence, a particular integral
1
of the given differential equation is given as f (D) F (x). First we give a rather general idea
of getting a particular integral with this method and then state some other useful direct
results. Note that the operator f (D) can be expressed as (D − α1 )(D − α2 ) . . . (D − αn ) and
thus a particular integral is given as
1 1 1 1
F (x) = ... F (x) (31.1)
f (D) D − α1 D − α2 D − αn
1
We give a general idea of evaluating an expression of the type F (x). This procedure
D−α
can be repeatedly applied to find a particular integral (31.1). However, applicability of this
method depends upon the form of F (x).
We give a general theorem that can be applied to any problem for finding particular inte-
gral of a differential equation.
31.1.1 Theorem 1
265 www.AgriMoon.Com
Linear Differential Equation of Higher Order
Since our interest is finding a particular integrals, the constant of integration is dropped.
Thus,
1
Z
F (x) = eαx F (x)e−αx dx.
D−α
Now we state some useful result those will be used to find P.I. of certain special forms of
F (x).
31.1.2 Theorem 2
31.1.3 Theorem 3
If α is a constant and g(x) is any function, then f (D) (eαx g(x)) = eαx f (D + α)g(x)
Proof: We know that D (eαx g(x)) = αeαx g(x) + eαx Dg(x) = eαx (α + D)g(x). Similar to the
proof of previous theorem we can prove with induction that Dn eαx g(x) = eαx (α + D)n g(x)
for any natural number n. This proves the result f (D) (eαx g(x)) = eαx f (D + α)g(x). This
result is known as shifting property of operator f (D).
266 2 www.AgriMoon.Com
Linear Differential Equation of Higher Order
31.1.4 Theorem 4
Proof: It can easily be verified that D2 sin(αx + β) = −α2 sin(αx + β) and D2 cos(αx + β) =
−α2 cos(αx + β). In other words, we can replace D 2 by −α2 and this proves the desired
result.
Now we describe the method for some special form of F (x).
We know from Theorem 31.1.2 that f (D)eαx = f (α)eαx . Operating on both sides by
1/f (D) we get
1 1 αx
eαx = f (α)eαx ⇒ eαx = f (α) e
f (D) f (D)
This implies that
1 αx 1 αx
e = e , provided f (α) 6= 0
f (D) f (α)
If f (α) = 0, then (D − α) is a factor of f (D), say f (D) = (D − α)g(D). Then
1 αx 1 1 αx 1 1 αx
e = e = e provided g(α) 6= 0
f (D) (D − α) g(D) (D − α) g(α)
In case g(α) = 0 then , say f (D) = (D − α)2 h(D). In this case we get
1 αx 1 1 αx 1 x2 αx
e = e = e provided h(α) 6= 0
f (D) h(α) (D − α)2 g(α) 2!
Again, if h(α) = 0, the same procedure can be repeated. To conclude, we have the follow-
ing results:
1 αx 1 αx
(i) e = e , where f (α) 6= 0
f (D) f (α)
267 3 www.AgriMoon.Com
Linear Differential Equation of Higher Order
(ii) If f (α) = 0, then f (D) must posses a factor of the type (D − α)r , say f (D) =
(D − α)r g(D) where g(α) 6= 0. Then the following formula is applicable
1 αx xr αx
e = e
(D − α)r r!
.
31.3.1 Problem 1
Find the general solution of the differential equation (D2 − 3D + 2)y = e3x .
Solution: The auxiliary equation is
(m2 − 3m + 2) = 0 ⇒ (m − 1)(m − 2) = 0 ⇒ m = 1, 2.
31.3.2 Problem 2
268 4 www.AgriMoon.Com
Linear Differential Equation of Higher Order
We express f (D) as a function of D2 , say f (D) = φ(D2 ). From Theorem 31.1.4 we know
that φ(D2 ) sin(αx + β) = φ(−α2 ) sin(αx + β). Applying [φ(D2 )]−1 both sides we obtain
1
sin(αx + β) = 2
φ(−α2 ) sin(αx + β)
φ(D )
31.5.1 Problem 1
(m2 + 1) = 0 ⇒ m = ±i
Hence, C.F. = (c1 cos x + c2 sin x). The particular integral is given by
1 1 1
P.I. = cos 2x = cos 2x = cos 2x.
D2 +1 2
(−2 + 1) −3
1
The required solution is: y = (c1 cos x + c2 sin x) − cos 2x.
3
269 5 www.AgriMoon.Com
Linear Differential Equation of Higher Order
31.5.2 Problem 2
Suggested Readings
McQuarrie, D.A. (2009). Mathematical Methods for Scientist and Engineers. First Indian
Edition. Viva Books Pvt. Ltd. New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
270 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 32
Here we continue discussion for solving linear equation of the form f (D)y = F (x). In the
last lesson, we have found particular integral for two different types of functions F (x).
In this lesson we shall continue discussing various other situations for finding particular
integral.
Take out the lowest degree term from f (D), so as to reduce it in the form [1 ± f (D)]n . Take
it to numerator, i.e., [1 ± f (D)]−n and expand it in ascending powers of D with the help of
Binomial series:
32.1.1 Example
(m2 + m) = 0 ⇒ m = 0, −1.
Taking 1 + D into numerator and expending this into an infinite series we get
1 1
P.I. = (1 − D + D 2 − D 3 + ...)(x2 + 2x + 4) = (x2 + 2x + 4 − 2x − 2 + 2)
D D
271 www.AgriMoon.Com
Linear Differential Equation of Higher Order (Cont.)
1
Operating 1/D on each term, we obtain P.I. = (x2 + 4) = (x3 /3 + 4x). The desired gen-
D
eral solution is 3
x
y = c1 + c2 e−x + + 4x .
3
32.2 Rule IV: F (x) is of the form eαx V , where V is any function of x
Using shift property of the operator discussed in the last lesson we can easily prove that
1 αx 1
e V = eαx V.
f (D) f (D + α)
32.2.1 Example
m2 − 2m + 1 = 0, and m = 1, 1.
x4
The required solution is y = (c1 + c2 x)ex + ex .
12
272 2 www.AgriMoon.Com
Linear Differential Equation of Higher Order (Cont.)
where V is a function of x. We start with the fact that for a given function g(x) we have
In general, by the method of induction for any natural number n we can show that
n n d n
D (xg(x)) = xD (g(x)) + D g(x)
dD
273 3 www.AgriMoon.Com
Linear Differential Equation of Higher Order (Cont.)
This implies
1 d 1
xV (x) = f (D) x V (x) − f (D) V (x)
f (D) dD f (D)
32.3.1 Example
274 4 www.AgriMoon.Com
Linear Differential Equation of Higher Order (Cont.)
In this case Rule IV or Rule V can be applied. For the application of rule IV we should
note that
1 m 1 m iαx
1. x sin αx = Im x e
f (D) f (D)
1 m 1 m iαx
2. x cos αx = Re x e .
f (D) f (D)
32.4.1 Example
275 5 www.AgriMoon.Com
Linear Differential Equation of Higher Order (Cont.)
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
McQuarrie, D.A. (2009). Mathematical Methods for Scientist and Engineers. First Indian
Edition. Viva Books Pvt. Ltd. New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
Piskunov, N. (1996). Differential and Integral Calculus (Volume - 2). First Edition. CBS
Publisher, Moscow.
276 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 33
In the last lesson we have discussed operator method of finding particular integral. In
this lesson we lean method of undetermined coefficients for finding particular integral of
non-homogeneous differential equations. This method is relatively easier to apply once a
possible form of a particular integral is known. This method is mainly applicable to linear
differential equations with constant coefficients.
277 www.AgriMoon.Com
Method of Undetermined Coefficients
33.2.1 Problem 1
y ′′ + 5y ′ + 6y = 2x + 1 (33.1)
To find particular integral, the trick is to somehow to guess one particular solution to
Equation (33.1). Note that 2x + 1 is a polynomial, and the left hand side of the equation
will be a polynomial if we let y be a polynomial of the same degree. Let us try
yp = Ax + B
278 2 www.AgriMoon.Com
Method of Undetermined Coefficients
The general solution must satisfy the given initial conditions. First find
1
y ′ = −2C1 e−2x − 3C2 e−3x +
3
Then
1 1 1
0 = y(0) = C1 + C2 − , = y ′ (0) = −2C1 − 3C2 + .
9 3 3
We solve to get C1 = 1/3 and C2 = −2/9. The particular solution we want is
1 2 3x − 1 3e−2x − 2e−3x + 3x − 1
y(x) = e−2x − e−3x + = .
3 9 9 9
33.2.2 Problem 2
y ′′ + 2y ′ + 2y = cos(2x).
Solution: We start by guessing the solution that includes some multiple of cos(2x). We
may have to also add a multiple of sin(2x) to our guess since derivatives of cosine are
sines. We try
yp = A cos(2x) + B sin(2x).
The left hand side must equal to right hand side. We group terms and get −4A+4B +2A =
1 and −4B − 4A + 2B = 0. So −2A + 4B = 1 and 2A + B = 0 and hence A = −1 10 and
1
B = 5 . Hence a particular solution is
− cos(2x) + 2 sin(2x)
yp = A cos(2x) + B sin(2x) = .
10
Remark 1: If the right hand side contains exponentials we try exponentials. For
example, for
Ly = e3x ,
279 3 www.AgriMoon.Com
Method of Undetermined Coefficients
Remark 2: If the right hand side is a multiple of sines, cosines, exponentials, and
polynomials, we can use the product rule for differentiation to come up with a guess. We
need to guess a form for yp such that Lyp is of the same form, and has all the terms needed
to for the right hand side. For example,
We will plug in and then hopefully get equations that we can solve for A, B, C, D, E, and
F.
In this case we find u that solves Lu = e2x and v that solves Lv = cos x (that is, do each
term separately). Then note that if y = u + v , then Ly = e2x + cos x. This is because L is
linear; we have Ly = L(u + v) = Lu + Lv = e2x + cos x.
33.2.3 Problem 3
Solution: By splitting up the right side of the given differential equation, we obtain the
three differential equations
Solutions of these three equations can be found with appropriate guess of the particular
integral discussed above. Finally, a particular solution is their sum, namely,
1 3 5 10 2
Y (t) = e2t + cos t sin t + et cos2t + et sin 2t.
2 17 17 13 13
280 4 www.AgriMoon.Com
Method of Undetermined Coefficients
The procedure illustrated in these examples enables us to solve a large class of problems
in a reasonably efficient manner. However, there is one difficulty that sometimes occurs.
It could be that our guess actually solves the associated homogeneous equation. The next
example illustrates how it arises.
33.2.4 Problem 4
y ′′ − 9y = e3x
Solution: In order to find a particular integral an intelligent guess would be y = Ae3x , but
if we plug this into the left hand side of the equation we get
There is no way we can choose A to make the left hand side be e3x because our guess sat-
isfies homogeneous equation. Note that the general solution of the homogeneous equation
is
C.F. = C1 e−3x + C2 e3x
Thus our assumed particular solution is actually a solution of the corresponding homo-
geneous equation; consequently, it cannot possibly be a solution of the nonhomogeneous
equation. To find a particular solution we must therefore consider functions of a some-
what different form. We modify our guess to y = Axe3x and notice there is no difficulty
anymore. Note that y ′ = Ae3x + 3Axe3x and y ′′ = 6Ae3x + 9Axe3x . So
Thus 6Ae3x is supposed to equal e3x . Hence, 6A = 1 and so A = 61 . We can now write the
general solution as
1 3x
y = yc + yp = C1 e−3x + C2 e3x + xe .
6
33.2.5 Problem 5
281 5 www.AgriMoon.Com
Method of Undetermined Coefficients
Then, upon calculating yp′ and Yp′′ , substituting them into the given differential equation,
we find that
4A sin 2t + 4B cos 2t = 3cos2t
Therefore A = 0 and B = 3/4, so a particular solution of the given differential equation is
3
yp (t) = t sin 2t
4
Remark 4: It is also possible that multiplying by x does not get rid of the problem we
had faced in last two examples. For example,
y ′′ − 6y ′ + 9y = e3x .
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
282 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 34
In the last lesson we have discussed method of undetermined coefficients for finding par-
ticular integral. In this lesson we lean another rather general method, called method of
variation of parameters, of finding particular integral of non-homogeneous differential
equation. In contrast to the method of undetermined coefficients, this method is also ap-
plicable for solving linear equations with variable coefficients. For the sake of simplicity
we restrict ourselves for second order linear differential equations. However the method
is also applicable for higher order linear differential equations.
y ′′ + P y ′ + Qy = R (34.1)
y ′′ + P y ′ + Qy = 0 (34.2)
y = au + bv (34.3)
where a, b are two arbitrary constants and u, v are functions of x. Since u and v are solu-
tions of (34.2), we have
u′′ + P u′ + Qu = 0, v ′′ + P v ′ + Qv = 0. (34.4)
yp = Au + Bv (34.5)
283 www.AgriMoon.Com
Method of Variation of Parameters
is a particular integral of (34.1). To determine A(x) and B(x) we need to have two equa-
tions. These are obtained as follows. First we compute
In order to avoid second order derivatives of A and B and to simplify the above expression
we take
A′ u + B ′ v = 0, (34.7)
Using the values of y, y ′ and y ′′ given by (34.5), (34.8) and (34.9) into the Equation (34.1),
we get
A′ u′ + B ′ v ′ = R (34.10)
284 2 www.AgriMoon.Com
Method of Variation of Parameters
where Z Z
vR uR
f (x) = − dx, g(x) = dx
W W
Using (34.12) into (34.5), we have
yp = uf (x) + vg(x).
yp = au + bv + uf (x) + vg(x).
34.2.1 Problem 1
In order to find a particular integral, we have u = cos nx and v = sin nx and R = sec nx.
The Wronskian is given as
cos nx sin nx
W = = n 6= 0.
−n sin nx n cos nx
where
sin nx sec nx 1
Z Z
vR
f (x) = − dx = − dx = 2 ln(cos nx)
W n n
285 3 www.AgriMoon.Com
Method of Variation of Parameters
and
cos nx sec nx
Z Z
uR x
g(x) = dx = − dx =
W n n
Hence, the required solution is
1 x
y = (c1 cos nx + c2 sin nx) + cos nx 2
ln(cos nx) + sin nx
n n
34.2.2 Problem 2
To find particular integral we have u = cos nx, v = sin nx, R = sec nx. The Wronskian is
given by
cos nx sin nx
W = = n 6= 0.
−n sin nx n cos nx
where
sin nx tan nx 1
Z Z
vR
f (x) = − dx = − dx = 2 [sin nx − ln(sec nx + tan nx)]
W n n
and
cos nx tan nx 1
Z Z
uR
g(x) = dx = − dx = − 2 cos nx
W n n
The desired general solution is
cos nx 1
y = (c1 cos nx + c2 sin nx) + 2
[sin nx − ln(sec nx + tan nx)] − 2 sin nx cos nx
n n
34.2.3 Problem 3
d2 y
Solve the differential equation 2
+ n2 y = cot nx.
dx
286 4 www.AgriMoon.Com
Method of Variation of Parameters
In this case, we have u = cos nx, v = sin nx and R = cot nx. The Wronskian is given by
cos nx sin nx
W = = n 6= 0
−n sin nx n cos nx
where
sin nx cot nx 1
Z Z
vR
f (x) = − dx = − dx = − 2 sin nx
W n n
and
cos nx cot nx 1 h
Z Z
uR nx i
g(x) = dx = − dx = 2 cos nx + ln tan
W n n 2
The required solution is
1 1 h nx i
y = (c1 cos nx + c2 sin nx) − cos nx sin nx + cos nx + ln tan sin nx
n2 n2 2
34.2.4 Problem 4
Using the method of variation of parameters, find the general solution of the differential
equation
d2 y
+ y = sec2 x
dx2
287 5 www.AgriMoon.Com
Method of Variation of Parameters
where
Z Z Z
vR 2
f (x) = − dx = − sin x sec xdx = − sec x tan xdx = − sec x
W
and Z Z Z
uR 2
g(x) = dx = cos x sec xdx = sec xdx = ln[sec x + tan x]
W
The required solution is
Suggested Readings
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
Weir, M.D., Hass, J. and Giordano, F.R. (2005). Thomas Calculus. Eleventh Edition.
Pearson Education. New Delhi.
McQuarrie, D.A. (2009). Mathematical Methods for Scientist and Engineers. First Indian
Edition. Viva Books Pvt. Ltd. New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
288 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 35
In this lesson we shall study two special forms of linear equations with variable coeffi-
cients which can be reduced to linear differential equations with constant coefficients by
a suitable substitution. Those special forms which we study here are called Cauchy-Euler
homogeneous linear differential equations and Legendre’s homogeneous linear differen-
tial equations.
The above equation can be reduced to linear differential equation with constant coeffi-
cients by substituting
dz 1
x = ez , or ln x = z, so that = (35.3)
dx x
Using chain rule for differentiation we obtain
dy dy dz 1 dy
= =
dx dz dx x dz
289 www.AgriMoon.Com
Equations Reducible to Linear Differential Equations with Constant Coefficients
d
Defining =: D1 , we have
dz
dy dy
x = ⇔ xDy = D1 y
dx dz
Similarly, for the second order derivative
d2 y
ddy d 1 dy 1 dy 1 d dy
2
= = =− 2 +
dx dxdx dx x dz x dz x dx dz
1 d2 y
1 dy 1 d dy dz 1 dy
=− 2 + =− 2 + 2 2
x dz x dz dz dx x dz x dz
Thus, we have
d2 y d2 y dy
x2 = − ⇒ x2 D 2 y = D1 (D1 − 1)y.
dx2 dz 2 dz
Similarly, x3 D3 y = D1 (D1 − 1)(D1 − 2)y and so on. In general, we have the relationship
The Equation (35.4) is a linear differential equation with constant coefficients which can
solved with the methods discussed in previous lessons. Finally, by replacing z by ln x we
obtain the desired solution of the given differential equation.
35.2.1 Problem 1
290 2 www.AgriMoon.Com
Equations Reducible to Linear Differential Equations with Constant Coefficients
The roots of the corresponding characteristic equation are m = 2, −2. The required solu-
tion of the transformed equation is
y = c1 e2z + c2 e−2z
y = c1 x2 + c2 x−2 .
35.2.2 Problem 2
Find the general solution of the differential equation (x2 D2 + y)y = 3x2 .
Solution: Substituting x = ez , the given equation reduces to
Substituting z = ln x, we get
√ h √ √ i
C.F. = x c1 cos ln x 3/2 + c1 sin ln x 3/2
291 3 www.AgriMoon.Com
Equations Reducible to Linear Differential Equations with Constant Coefficients
This implies
dy dy
(a + bx) =b ⇔ (a + bx)Dy = bD1 y
dx dz
Similarly for the second order derivative we get
d2 y
d dy d b dy
2
= =
dx dx dx dx (a + bx) dz
d2 y b2
dy b d dy
=− +
dx2 (a + bx)2 dz (a + bx) dx dz
b2
dy b d dy dz
=− +
(a + bx)2 dz (a + bx) dz dz dx
d2 y b2 dy b2 d2 y
= − +
dx2 (a + bx)2 dz (a + bx)2 dz 2
This gives us
2y
d2 y dy
2d 2
(a + bx) =b − ⇔ (a + bx)2 D 2 y = b2 D1 (D1 − 1)y
dx2 dz 2 dz
292 4 www.AgriMoon.Com
Equations Reducible to Linear Differential Equations with Constant Coefficients
In general, we have
The methods of solving this transformed equation are same as discussed in previous sec-
tion.
35.3.1 Example
d
Solution: Using D = and dividing both sides by (x+1), the given differential equation
dx
can be rewritten as
This is the Legendre’s homogeneous linear equation which can be solved by substituting
(1 + x) = ez ⇔ ln(1 + x) = z
or
D13 − D12 − D1 + 1 y = e−2z
293 5 www.AgriMoon.Com
Equations Reducible to Linear Differential Equations with Constant Coefficients
m3 − m2 − m + 1 y = 0
The roots of the characteristics equations are m = 1, 1, −1. Hence the complimentary
function of the transformed differential equation is
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
McQuarrie, D.A. (2009). Mathematical Methods for Scientist and Engineers. First Indian
Edition. Viva Books Pvt. Ltd. New Delhi.
Raisinghania, M.D. (2005). Ordinary & Partial Differential Equation. Eighth Edition. S.
Chand & Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
294 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 36
Let x and y be the dependent and t be the independent variable. Thus, in such equations
there occur differential coefficients of x, y with respect to t. Let D = d/dt, then such
equations can be put into the form
where T1 and T2 are functions of the independent variable t and f1 (D), f2 (D), g1 (D), and
g2 (D) are all rational integral functions of D with constant coefficients. In general, the
number of equations will be equal to the number of dependent variables, i.e., if there are
n dependent variables there will be n equations.
In order to eliminate y between equations (36.1) and (36.2), operating on both sides of
(36.1) by g2 (D) and on both sides of (36.2) by f2 (D) and subtracting, we get
This is a linear differential equation with constant coefficients in x and t and can be solved
to give the value of x in terms of t. Substituting this value of x in either (36.1) or (36.2),
we get the value of y in terms of t.
295 www.AgriMoon.Com
Methods for Solving Simultaneous Ordinary Differential Equations
Remark 1: The above Equations (36.1) and (36.2) can be also solved by first elimi-
nating x between them and solving the resulting equation to get y in terms of t. Substitut-
ing this value of y in either (36.1) or (36.2), we get the value of x in terms of t.
Remark 2: In the general solutions of (36.1) and (36.2) the number of arbitrary
constants will be equal to the sum of the orders of the equations (36.1) and (36.2).
36.3.1 Problem 1
Solution: Writing D for d/dt, the given equations can be rewritten in the following
symbolic form as
(D − 7)x + y = 0 (36.6)
−2x + (D − 5)y = 0 (36.7)
or
(D 2 − 12D + 37)y = 0
296 2 www.AgriMoon.Com
Methods for Solving Simultaneous Ordinary Differential Equations
The roots of the auxiliary equation are m = 6 ± i. Therefore, we get the general solution
for the variable y as
where c1 and c2 being arbitrary constants. We now find x by using Equation (36.7). Now
from (36.10), differentiating w.r.t. t, we get
or on simplifications we obtain
Thus, equations (36.10) and (36.12) give the desired general solution.
36.3.2 Problem 2
D 2 y − y + 5Dv = x (36.13)
2Dy − D 2 v + 4v = 2 (36.14)
Solution: Multiplying (36.13) by 2D and (36.14) by (D2 − 1) and then subtracting (36.14)
from the Equation (36.13) we obtain
or
(D 4 + 5D 2 + 4)v = 4 (36.15)
297 3 www.AgriMoon.Com
Methods for Solving Simultaneous Ordinary Differential Equations
This is a linear differential equations with constant coefficients whose solution can easily
be found. The characteristic equation of the corresponding homogeneous equation is
Now we find an equation giving y in terms of v . This can be done by eliminating from the
equations (36.13) and (36.14) those terms which involve derivatives of y. So multiplying
Equation (36.13) by 2 and Equation (36.14) by D we get
or
1
y = −x + D 3 v + 3Dv (36.20)
2
Substitute v from (36.16) into the Equation (36.21) to obtain the expression for y as
5 5
y = −x − c1 cos x + c2 cos x + 2c4 cos 4x − 2c3 sin 2x, (36.21)
2 2
298 4 www.AgriMoon.Com
Methods for Solving Simultaneous Ordinary Differential Equations
36.4.1 Example
Dx − y = t (36.22)
x + Dy = 1 (36.23)
D 2 x − Dy = 1 (36.24)
(D 2 + 1)x = 2 (36.25)
The auxiliary equation of the above differential equation is m2 + 1 = 0 and therefore the
general solution of the homogeneous equation is
299 5 www.AgriMoon.Com
Methods for Solving Simultaneous Ordinary Differential Equations
Suggested Readings
300 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 37
37.1 Introduction
x2 y ′′ + xy ′ + (x2 a2 )y = 0,
where c is a constant.
Here we provide some definitions which will be very useful for finding series solution of
the differential equations.
301 www.AgriMoon.Com
Series Solutions about an Ordinary Point
exists and converges to f (x) for all x in the interval of convergence of (37.2).
where p(x) = Q(x)/P (x) and q(x) = R(x)/P (x). The Equation (37.3) is called equivalent
normalized form of the Equation (37.1).
If the point x = x0 is not an ordinary point of the differential Equation (37.1) or (37.3),
then it is called a singular point of the differential equation of (37.3). There are two types
of singular points: (i) regular singular points, and (ii) irregular singular points. A singular
point x = x0 of the differential Equation (37.3) is called a regular singular point of the
differential Equation (37.3) if both
are analytic at x = x0 . A singular point, which is not regular is called an irregular singular
point.
302 2 www.AgriMoon.Com
Series Solutions about an Ordinary Point
37.3.1 Problem 1
d2 y x dy 1
2
+ − y = 0. (37.4)
dx (x − 1)(x + 1) dx (x − 1)(x + 1)
Since both p(x) and q(x) are analytic at x = 0, the point x = 0 is an ordinary point of the
given Equation (37.4). Further note that both p(x) and q(x) are not analytic at x = 1, thus
x = 1 is not an ordinary point and so x = 1 is a singular point. Also
show that both (x − 1)P (x) and (x − 1)2 Q(x) are analytic at x = 1. Therefore x = 1 is a
regular singular point.
37.3.2 Problem 2
Determine whether the point x = 0 is an ordinary point or regular point of the differential
equation
xy ′′ + sin(x)y + x2 y = 0
303 3 www.AgriMoon.Com
Series Solutions about an Ordinary Point
37.3.3 Problem 3
is not analytic at x = 2. Hence both x = 0 and x = 2 are singular point of the differential
equations. At x = 0 we have
1 3x2
xp(x) = and x2 q(x) =
(x(x − 2)) (x − 2)2
Note that x2 q(x) is non-singular at x = 0 but xp(x) is not analytic at this point. Hence
x = 0 is an irregular singular point. At x = 2 we have
1
(x − 2)p(x) = and (x − 2)2 q(x) = 2
x2
Both functions are analytic at x = 2 and hence x = 2 is a regular singular point.
where x is a variable and cn are constants, called coefficients of the series. There are
three possibilities about the convergence of a power series. The series may converge only
at x = 0 or it may converge for all values of x. If this is not the case then a definite
positive number R exists such that the given series converges for every |x − x0 | < R and
304 4 www.AgriMoon.Com
Series Solutions about an Ordinary Point
diverges for every |x − x0 | > R. Such a number is known as the radius of convergence and
]x0 − R, x0 + R[, the interval of convergence, of the given series.
Among several formulas for determining convergence of the power series, ratio test is
most common and simple to use. Given a power series ∞
P
n=0 cn (x − x0 ) we compute
n
1 cn+1
= lim ,
R n→∞ cn
then the series is convergence for |x − x0 | < R and divergent |x − x0 | > R.
37.4.1 Example
If x = x0 is an ordinary point of (37.5), then (37.5) has two non-trivial linearly indepen-
dent power series solutions of the form
∞
X
Cn (x − x0 )n (37.6)
n=0
and these power series converge in some interval of convergence |x − x0 | < R, (where R
is the radius of convergence of (37.6)) about x0 .
305 5 www.AgriMoon.Com
Series Solutions about an Ordinary Point
and then to find out coefficients Cn we need to differentiate (37.7) and plug in the deriva-
tives into the Equation (37.6). Once we have the appropriate coefficients, we call (37.7)
the series solution to (37.5) near x = x0 . More precisely, differentiating twice, the Equa-
tion (37.7) yields
∞
X ∞
X
′ ′′
y = nCn (x − x0 ) n−1
and y = n(n − 1)Cn (x − x0 )n−2 (37.8)
n=0 n=0
A0 + A1 (x − x0 ) + A2 (x − x0 )2 + . . . + An (x − x0 )n + . . . = 0, (37.9)
where the coefficients A0 , A1 , A2 . . . etc. are now some functions of the coefficients C0 , C1 , C2 , . . .
etc. Since the Equation (37.9) is an identity, all the coefficients A0 , A1 , A2 . . . of (37.9)
must be zero, i.e.,
A0 = 0, A1 = 0, A2 = 0, . . . , An = 0 (37.10)
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Raisinghania, M.D. (2009). Advanced Differential Equations. Twelfth Edition. S. Chand
& Company Ltd., New Delhi.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
306 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 38
In the last lesson we have discussed series solution of the homogeneous differential equa-
tions. In this lesson we demonstrate the method by using a couple of basic examples. For
demonstration we take first example of a differential equation with constant coefficients
and then some more involved examples will be discussed.
38.1.1 Problem 1
Differentiating y , we have
∞
X ∞
X
′ ′′
y (x) = ncn xn−1 and y (x) = n(n − 1)cn xn−2
n=1 n=2
This implies
∞
X
[(n + 2)(n + 1)cn+2 − cn xn ] xn = 0
n=0
307 www.AgriMoon.Com
Series Solutions about an Ordinary Point (Cont.)
Since the series is always equal to 0 then each coefficient must be zero. Thus we have
Putting n = 0, 1, 2 . . ., we get
c0 c1 c0 c1
c2 = , c3 = , c4 = , c5 = ,...
2! 3! 4! 5!
In general, we have
c0 c1
c2k = , c2k+1 = . . . for k = 1, 2, . . . .
(2k)! (2k + 1)!
Putting these values into the series and collecting the c0 and c1 terms we get
x2 x2k x3 x2k+1
y(x) = c0 1+ + ...+ + . . . + c1 x + + ...+ + ...
2! (2k)! 3! (2k + 1)!
This is the desired series solution. It should be noted that this series solution can be
rewritten into the form of well known solution y(x) = c1 ex + c2 e−x of the given differential
equation as
x2 x2
x −x
c1 e + c2 e = c1 1+x+ + . . . + c2 1 − x + + ...
2! 2!
308 2 www.AgriMoon.Com
Series Solutions about an Ordinary Point (Cont.)
38.1.2 Problem 2
which leads to
∞
X ∞
X ∞
X ∞
X
n−2 n−1 n
n(n − 1)cn x −2 n(n − 1)cn x + n(n − 1)cn x − 2 cn xn = 0
n=2 n=2 n=2 n=0
In order to write the series in terms the coefficients of xn we shift the summation index as
∞
X ∞
X ∞
X ∞
X
n n n
(n + 2)(n + 1)cn+2 x − 2 n(n + 1)cn+1 x + n(n − 1)cn x − 2 cn xn = 0
n=0 n=1 n=2 n=0
The sum in second and third series can also start from 0 without changing the series. This
leads to
∞
X
[(n + 2)(n + 1)cn+2 − 2n(n + 1)cn+1 + n(n − 1)cn − 2cn ] xn = 0
n=0
Putting n = 0, 1, 2, . . . we get
1
c2 = c0 , c3 = (2c0 + c1 ) =: c, c4 = c, c5 = c . . .
3
Hence the series solution becomes
∞
X
2
y = c0 + c1 x + c0 x + c xn .
n=3
309 3 www.AgriMoon.Com
Series Solutions about an Ordinary Point (Cont.)
38.1.3 Problem 3
Find the power series solution of the equation (x2 + 1)y ′′ + xy ′ − xy = 0 in powers of x
(i.e. about x = 0).
Solution: Clearly x = 0 is an ordinary point of the given differential equation. Therefore,
to find the series solution, we take power series
∞
X
2 3
y = c0 + c1 x + c2 x + c3 x + . . . = cn xn . (38.3)
n=0
Putting the above value of y, y ′ and y ′′ in the given differential equation, we obtain
∞
X ∞
X ∞
X
(x2 + 1) n(n − 1)cn xn−2 + x ncn xn−1 − x ncn xn = 0
n=2 n=1 n=0
∞
X ∞
X ∞
X ∞
X
n n−2 n
⇒ n(n − 1)cn x + n(n − 1)cn x − ncn x − cn xn+1 = 0
n=2 n=2 n=1 n=0
This leads to
∞
X ∞
X ∞
X ∞
X
n n n
n(n − 1)cn x + (n + 2)(n + 1)cn+2 x + ncn x − cn−1 xn = 0
n=2 n=0 n=1 n=1
Equating the constant term and the coefficients of various powers of x, we get
310 4 www.AgriMoon.Com
Series Solutions about an Ordinary Point (Cont.)
y = c0 + c1 x + c2 x2 + c3 x3 + c4 x4 + c5 x5 + . . . ∞
38.1.4 Problem 4
Find the power series solution of the initial value problem xy ′′ + y ′ + 2y = 0, y(1) = 1,
y ′ (1) = 2 in powers of (x − 1).
This leads to
∞
X ∞
X ∞
X ∞
X
n−1 n−2 n−1
n(n − 1)cn (x − 1) + n(n − 1)cn (x − 1) + ncn (x − 1) +2 cn (x − 1)n = 0
n=2 n=2 n=1 n=0
311 5 www.AgriMoon.Com
Series Solutions about an Ordinary Point (Cont.)
Suggested Readings
312 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 39
39.1 Introduction
In this lesson we discuss series solution about a singular point. In particular, the power
series method discussed in last lessons will be generalized. The generalized power series
method is also known as Frobenius method.
Let us consider a simple first order differential equation 2xy ′ − y = 0 and try to apply the
power series method discussed in the last lessons. Note that x = 0 is a singular point. If
we plug in
∞
X
y= ak xk ,
k=0
Consider the differential equation of the form y ′′ + p(x)y ′ + q(x)y = 0. Note that xp(x) and
x2 q(x) are analytic at x = 0. We try a series solution of the from
∞
X
r
y=x cn xn = xr (c0 + c1 x + c2 x2 + . . .), where c0 6= 0
n=0
313 www.AgriMoon.Com
Series Solutions about a Regular Singular Point
∞
X
′′
y = (n + r)(n + r − 1)cn xn+r−2
n=0
xp(x) ′ x2 q(x)
y ′′ + y + y=0
x x2
Substituting all values of y , y ′, y ′′, xp(x) and x2 q(x) series into the above differential
equation we get
∞
X ∞
X ∞
X ∞
X ∞
X
(n+r)(n+r−1)cn xn+r−2 + an xn−1 × (n+r)cn xn+r−1 + bn xn−2 × cn xn+r = 0
n=0 n=0 n=0 n=0 n=0
Multiplying by x2 we get
∞
X ∞
X ∞
X ∞
X ∞
X
n+r n n+r n
(n + r)(n + r − 1)cn x + an x × (n + r)cn x + bn x × cn xn+r = 0
n=0 n=0 n=0 n=0 n=0
We can now equate coefficients of various powers of x to zero to form a system of equa-
tions involving unknown coefficients cn . Equating the coefficient of xr we obtain
[r(r − 1) + a0 r + b0 ]c0 = 0
Since c0 6= 0, we obtain
The above quadratic equation is known as the indicial equation of the given differential
equation. The general solution of the given differential equation depends on the roots of
the indicial equation. There are three possible general cases:
314 2 www.AgriMoon.Com
Series Solutions about a Regular Singular Point
39.2.1 Case I: The indicial equation has two real roots which do not differ by an
integer
Let r1 and r2 are the roots of the indicial equation. Then the two linearly independent
solution will follow from
∞
X ∞
X
r1 n r2
y1 (x) = x cn z y2 (x) = x cn z n
n=0 n=0
If the indicial equation has a doubled root r, then we find one solution
∞
X
y1 = xr ak xk ,
k=0
39.2.3 Case III: The indicial equation has two real roots which differ by an integer
If the indicial equation has two real roots such that r1 − r2 is an integer, then one solution
is ∞ X
y1 = xr1 ak xk ,
k=0
where we plug y2 into the given equation and solve for the constants bk and C .
315 3 www.AgriMoon.Com
Series Solutions about a Regular Singular Point
Remark 1: Note that the case-I also includes complex numbers because in that case
r1 − r2 will be a complex number which cannot be equal to a real integer.
Remark 2: Note that the mai idea is to find at least one Frobenius-type solution. If
we are lucky and find two, we are done. If we only get one, we either use the ideas above
or the method of variation of parameters to obtain a second solution.
39.3.1 Example
4xy ′′ + 2y ′ + y = 0
316 4 www.AgriMoon.Com
Series Solutions about a Regular Singular Point
Multiplying by x2 we obtain
∞ ∞ ∞
X
n+r 1X 1X
(n + r)(n + r − 1)cn x + (n + r)cn xn+r + cn xn+r+1 = 0 (39.3)
2 4
n=0 n=0 n=0
317 5 www.AgriMoon.Com
Series Solutions about a Regular Singular Point
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Raisinghania, M.D. (2009). Advanced Differential Equations. Twelfth Edition. S. Chand
& Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
318 6 www.AgriMoon.Com
Module 3: Ordinary Differential Equations
Lesson 40
In this lesson we continue series solution about a singular point. We shall demonstrate the
method with some useful differential equations.
40.1.1 Problem 1
where r is a real number, not necessarily an integer. Again if such a solution exists, it may
only exist for positive x. First let us find the derivatives
∞
X
′
y = (k + r) ak xk+r−1 ,
k=0
X∞
y ′′ = (k + r) (k + r − 1) ak xk+r−2 .
k=0
319 www.AgriMoon.Com
Series Solutions about a Regular Singular Point (Cont...)
Re-indexing leads to
∞
X ∞
X ∞
X ∞
X
k+r k+r k+r
4(k + r) (k + r − 1) ak x − 4(k + r − 1) ak−1 x + ak x − 2ak−1 xk+r = 0
k=0 k=1 k=0 k=1
4r(r − 1) + 1 = 0
It has a double root at r = 21 . All other coefficients of xk+r also have to be zero so
4(k + r) (k + r − 1) + 1 ak − 4(k + r − 1) + 2 ak−1 = 0.
1
If we plug in r = 2 and solve for ak , we get
4(k + 12 − 1) + 2 1
ak = 1 1 ak−1 = ak−1 .
4(k + 2 ) (k + 2 − 1) + 1 k
In other words,
∞ ∞ ∞
X
k+r
X 1 k+1/2 √ X 1 k √ x
y= ak x = x = x x = xe .
k! k!
k=0 k=0 k=0
So we have one solution of the given differential equation. Here we have written the series
in terms of elementary functions. However this is not always possible.
320 2 www.AgriMoon.Com
Series Solutions about a Regular Singular Point (Cont...)
40.1.2 Problem 2
x2 y ′′ + xy ′ + x2 − p2 y = 0.
(40.1)
which implies
∞
X ∞
X
y′ = cm (k + m)xk+m−1 , y ′′ = cm (k + m)(k + m − 1)xk+m−2
m=0 m=0
So the roots of indicial equation are k = p, −p. Next equating to zero the coefficient of
xk+1 in (40.3) gives
321 3 www.AgriMoon.Com
Series Solutions about a Regular Singular Point (Cont...)
cm (k + m + p)(k + m − p) + cm−2 = 0
1
⇒ cm = cm−2 .
(k + m + p)(p − k − m)
1
⇒ cm = cm−2 . (40.4)
(k + m + p)(p − k − m)
c1 = c3 = c5 = c7 = . . . = 0.
1 1
c4 = c2 = c0
(k + 4 + p)(p − k − 4) (k + 4 + p)(p − k − 4)(k + 2 + p)(p − k − 2)
and so on. Putting these values in (40.2) and also replacing c0 by 1, we get
x2 x4
y = 1+ + + ...
(k + 2 + p)(p − k − 2) (k + 4 + p)(p − k − 4)(k + 2 + p)(p − k − 2)
(−1)k x2k
∞
x2
X
p
y1 = x 1− + . . . = xp
4(1 + p) 22k k!(k + p)(k − 1 + p) · · · (2 + p)(1 + p)
k=0
(−1)k x2k
∞
x2
X
−p
y2 = x 1− + . . . = x−p
4(1 − p) 22k k!(k − p)(k − 1 − p) · · · (2 − p)(1 − p)
k=0
Therefore when 2p is not an integer, we have the general solution to Bessel’s equation of
order p
y = c1 y1 (x) + c2 y2 (x),
322 4 www.AgriMoon.Com
Series Solutions about a Regular Singular Point (Cont...)
Remark: We define the Bessel functions of the first kind Bessel function of the first
kind of order p and −p as
(−1)k
∞
1 X x 2k+p
Jp (x) = p y1 = ,
2 Γ(1 + p) k!Γ(k + p + 1) 2
k=0
(−1)k
∞
1 X x 2k−p
J−p (x) = −p y2 = .
2 Γ(1 − p) k!Γ(k − p + 1) 2
k=0
As these are constant multiples of the solutions we found above, these are both solutions to
Bessel’s equation of order p. When p is not an integer, Jp and J−p are linearly independent.
When 2p is an integer we obtain
(−1)k x 2k+p
∞
X
Jp (x) = .
k!(k + p)! 2
k=0
40.1.3 Problem 3
r(r − 1) + r = r 2 = 0.
(n + r)2 an = −an−1 , n ≥ 1.
323 5 www.AgriMoon.Com
Series Solutions about a Regular Singular Point (Cont...)
Suggested Readings
Boyce, W.E. and DiPrima, R.C. (2001). Elementary Differential Equations and Boundary
Value Problems. Seventh Edition, John Willey & Sons, Inc., New York.
Raisinghania, M.D. (2009). Advanced Differential Equations. Twelfth Edition. S. Chand
& Company Ltd., New Delhi.
Kreyszig, E. (1993). Advanced Engineering Mathematics. Seventh Edition, John Willey
& Sons, Inc., New York.
Arfken, G.B. (2001). Mathematical Methods for Physicists. Fifth Edition, Harcourt Aca-
demic Press, San Diego.
Grewal, B.S. (2007). Higher Engineering Mathematics. Fourteenth Edition. Khanna
Publishilers, New Delhi.
Dubey, R. (2010). Mathematics for Engineers (Volume II). Narosa Publishing House.
New Delhi.
Edwards, C.H., Penney, D.E. (2007). Elementary Differential Equations with Boundary
Value Problems. Sixth Edition. Pearson Higher Ed, USA.
324 6 www.AgriMoon.Com
Module-IV: Vector Calculus
Lesson 41
Introduction
A scalar function 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) is a function defined at each point in a certain domain 𝐷𝐷 in space. It
takes real values. It depends on the specific point 𝑃𝑃(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) in space, but not on any particular
coordinate system which may be used. For every point (𝑥𝑥, 𝑦𝑦, 𝑧𝑧) ∈ 𝐷𝐷, 𝑓𝑓 takes a real value. We
say the a scalar field 𝑓𝑓 is defined in 𝐷𝐷. For example, The distance function in the three
dimensional space taken as the Euclidean distance between the points 𝑃𝑃(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) and
𝑃𝑃0 (𝑥𝑥0 , 𝑦𝑦0, 𝑧𝑧0 )
𝑓𝑓(𝑃𝑃) = 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = �(𝑥𝑥 − 𝑥𝑥0 )2 + (𝑦𝑦 − 𝑦𝑦0 )2 + (𝑧𝑧 − 𝑧𝑧0 )2
V =V ( P ) =v1i + v2 j + v3 k
and we say that a vector field is defined in 𝐷𝐷. In Cartesian system of coordinates, it can be
written as
𝑉𝑉 = 𝑣𝑣1 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑖𝑖 + 𝑣𝑣2 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑗𝑗 + 𝑣𝑣3 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑘𝑘.
An example of a vector field is the velocity field 𝑉𝑉(𝑃𝑃) defined at any point 𝑃𝑃 on a rotating body.
Let 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) be a single valued continuous scalar function defined at every point 𝑃𝑃 ∈ 𝐷𝐷. Then
an equation of a surface is defined by 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 𝑐𝑐, a constant. It is called a level surface of the
function.
41.1.4 Example : We determine the level surface of the scalar field in space, defined by the
following function f ( x, y, z ) = x + y + z.
325 www.AgriMoon.Com
Introduction
We find that 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 𝑐𝑐 gives x + y + z =c which is equation of a plane. For different c they
define parallel planes. Therefore, the level surfaces are parallel planes.
41.1.5 Example: Determine the level surface of the scalar field in space, defined by the
function f ( x, y, z ) =x 2 + 9 y 2 + 16 z 2 .
The parametric representation of a curve 𝐶𝐶 in the two dimensional Cartesian plane is given by
=x x ( t=
) , y y (t ) , a ≤ t ≤ b . Using this the position vector of a point 𝑃𝑃 on the curve 𝐶𝐶 can be
r ( t ) x ( t ) i + y (t ) j .
written as =
Therefore, the position vector of a point on a curve defines a vector function. Similarly a three
dimensional curve or a space curve or a space curve 𝐶𝐶 can be parameterized as
The parametric form of a line passing through a point with position vector a and with the
direction of vector b is given by
𝑟𝑟(𝑡𝑡) = 𝑎𝑎 + 𝑡𝑡𝑡𝑡 = (𝑎𝑎1 + 𝑡𝑡𝑏𝑏1 )𝑖𝑖 + (𝑎𝑎2 + 𝑡𝑡𝑏𝑏2 )𝑗𝑗 + (𝑎𝑎3 + 𝑡𝑡𝑡𝑡3 )𝑘𝑘
2
326 www.AgriMoon.Com
Introduction
x2 y 2
+ =
1
a 2 b2
is given by
Let us consider the parabola 𝑦𝑦 2 = 4𝑎𝑎𝑥𝑥. Now take 𝑦𝑦 = 𝑡𝑡 as one parameter and then we can write
the parametric form of the parabola as
𝑡𝑡 2
𝑟𝑟(𝑡𝑡) = �4𝑎𝑎 � 𝑖𝑖 + 𝑡𝑡𝑡𝑡
We can give parametric representation of surfaces can be done using two parameters. Let
𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 𝑐𝑐 or 𝑔𝑔(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 0 be the equation of a surface. Let an explicit representation of the
surface be written as 𝑧𝑧 = ℎ(𝑥𝑥, 𝑦𝑦). Then, if we substitute 𝑢𝑢 = 𝑥𝑥, 𝑦𝑦 = 𝑣𝑣, the parametric form of
the surface can be reduced to
41.2.7 Example
41.2.8 Example
𝑟𝑟(𝑢𝑢, 𝑣𝑣) = 𝑎𝑎 cos 𝑢𝑢 cos 𝑣𝑣 𝑖𝑖 + 𝑎𝑎 sin 𝑢𝑢 cos 𝑣𝑣 𝑗𝑗 + 𝑎𝑎 sin 𝑣𝑣 𝑘𝑘, 0 ≤ 𝑢𝑢 ≤ 2𝜋𝜋, −𝜋𝜋/2 ≤ 𝑣𝑣 ≤ 𝜋𝜋/2
41.2.9 Example
𝑥𝑥 2 𝑦𝑦 2 𝑧𝑧 2
+ 𝑏𝑏 2 + 𝑐𝑐 2 = 1
𝑎𝑎 2
is given by
3
327 www.AgriMoon.Com
Introduction
The vector function 𝑣𝑣(𝑡𝑡) is said to have the limit p as 𝑡𝑡 → 𝑙𝑙 if 𝑣𝑣(𝑡𝑡) is defined in some
neighbourhood of l , except possibly at 𝑡𝑡 = 𝑙𝑙, and
We write 𝑝𝑝. In the Cartesian system, this implies that limits of the component functions
𝑣𝑣1 (𝑡𝑡), 𝑣𝑣2 (𝑡𝑡) and 𝑣𝑣3 (𝑡𝑡) exist as 𝑡𝑡 → 𝑙𝑙 and
lim𝑡𝑡→𝑙𝑙 𝑣𝑣1 (𝑡𝑡) = 𝑝𝑝1 , lim𝑡𝑡→𝑙𝑙 𝑣𝑣2 (𝑡𝑡) = 𝑝𝑝2 , lim𝑡𝑡→𝑙𝑙 𝑣𝑣3 (𝑡𝑡) = 𝑝𝑝3
41.3.2 Continuity
(𝑖𝑖) 𝑣𝑣(𝑡𝑡) is defined in some neighbourhood of l , (ii)lim𝑡𝑡→𝑙𝑙 𝑣𝑣(𝑡𝑡) exists, and (iii)lim𝑡𝑡→𝑙𝑙 𝑣𝑣(𝑡𝑡) =
𝑣𝑣(𝑙𝑙).
In Cartesian system, this implies that 𝑣𝑣(𝑡𝑡)is continuous at 𝑡𝑡 = 𝑙𝑙, if and only if the component
functions 𝑣𝑣1 (𝑡𝑡), 𝑣𝑣2 (𝑡𝑡) and 𝑣𝑣3 (𝑡𝑡) are continuous at 𝑡𝑡 = 𝑙𝑙.
41.3.3 Differentiability
𝑑𝑑𝑑𝑑
exists. If the limit exists, then we write it as 𝑣𝑣 ′ (𝑡𝑡) or as 𝑑𝑑𝑑𝑑
.
In Cartesian system, this implies that the component functions 𝑣𝑣1 (𝑡𝑡), 𝑣𝑣2 (𝑡𝑡) and 𝑣𝑣3 (𝑡𝑡) are
differentiable at a point 𝑡𝑡, and the limits
𝑣𝑣𝑖𝑖 (𝑡𝑡+∆𝑡𝑡)−𝑣𝑣𝑖𝑖 (𝑡𝑡)
lim ∆𝑡𝑡
, 𝑖𝑖 = 1,2,3 exist.
∆𝑡𝑡→0
Let 𝑣𝑣(𝑡𝑡) = 𝑟𝑟(𝑡𝑡) = 𝑥𝑥(𝑡𝑡)𝑖𝑖 + 𝑦𝑦(𝑡𝑡)𝑗𝑗 + 𝑧𝑧(𝑡𝑡)𝑘𝑘 be the parametric representation of a curve 𝐶𝐶.
4
328 www.AgriMoon.Com
Introduction
41.3.4 Example
𝜋𝜋 𝜋𝜋
Assume x = sin t . Then 𝑦𝑦 = 1 − 2𝑠𝑠𝑠𝑠𝑠𝑠2 𝑡𝑡 = 𝑐𝑐𝑐𝑐𝑐𝑐2𝑡𝑡 , The range of t is − 2 ≤ 𝑡𝑡 ≤ 2 . So
𝜋𝜋 𝜋𝜋
𝑟𝑟(𝑡𝑡) = 𝑠𝑠𝑠𝑠𝑠𝑠𝑠𝑠 𝑖𝑖 + 𝑐𝑐𝑐𝑐𝑐𝑐2𝑡𝑡 𝑗𝑗, − 2 ≤ 𝑡𝑡 ≤ 2
41.3.5 Example
We find the tangent vector to the curve with parametric representation given by
𝑡𝑡+1
𝑥𝑥 = 𝑡𝑡 3 , 𝑦𝑦 = 𝑡𝑡
, 𝑧𝑧 = 𝑡𝑡 2 + 1, at the point 𝑡𝑡 = 2.
First note that the position vector of a point on the given curve is
1
𝑟𝑟(𝑡𝑡) = 𝑡𝑡 3 𝑖𝑖 + �1 + 𝑡𝑡 � 𝑗𝑗 + (𝑡𝑡 2 + 1)𝑘𝑘, 𝑡𝑡 ≠ 0.
1
and 𝑟𝑟 ′ (2) = 12𝑖𝑖 − 4 𝑗𝑗 + 4𝑘𝑘 .
3
The position vector of the point at which 𝑟𝑟 ′ (2) is the tangent is 𝑟𝑟(2) = 8𝑖𝑖 + 2 𝑗𝑗 + 5𝑘𝑘.
Therefore we require the position vector of a point on the line passing through the point whose
position vector is 𝑟𝑟(2) and has the direction of 𝑟𝑟 ′ (2). Hence, parametric form of the line is given
by
3 𝑡𝑡 3 1
𝑥𝑥 = 8 + 12𝑡𝑡, 𝑦𝑦 = 2 − 4 , 𝑧𝑧 = 5 + 4𝑡𝑡 or 𝑟𝑟 ′ (𝑡𝑡) = �8, 2 , 5� + 𝑡𝑡(12, − 4 , 4).
5
329 www.AgriMoon.Com
Introduction
where . and × represent the dot and cross products, respectively. It must be mentioned that the
cross product of two vectors is not commutative.
41.3.7 Example
(𝑖𝑖) 𝑣𝑣(𝑡𝑡) = (cos 𝑡𝑡 + 𝑡𝑡 2 )(𝑡𝑡𝑡𝑡 + 𝑗𝑗 + 2𝑘𝑘) (𝑖𝑖𝑖𝑖) 𝑣𝑣(𝑡𝑡) = (3𝑡𝑡𝑡𝑡 + 5𝑡𝑡 2 𝑗𝑗 + 6𝑘𝑘). (𝑡𝑡 2 𝑖𝑖 − 2𝑡𝑡𝑡𝑡 + 𝑡𝑡𝑡𝑡)
Solution
(𝒊𝒊𝒊𝒊) 𝑣𝑣 ′ (𝑡𝑡) = (3𝑡𝑡𝑡𝑡 + 5𝑡𝑡 2 𝑗𝑗 + 6𝑘𝑘)′ . (𝑡𝑡 2 𝑖𝑖 − 2𝑡𝑡𝑡𝑡 + 𝑡𝑡𝑡𝑡) + (3𝑡𝑡𝑡𝑡 + 5𝑡𝑡 2 𝑗𝑗 + 6𝑘𝑘). (𝑡𝑡 2 𝑖𝑖 −
2𝑡𝑡𝑡𝑡 + 𝑡𝑡𝑡𝑡)′
= (3𝑖𝑖 + 10𝑡𝑡 𝑗𝑗)(𝑡𝑡 2 𝑖𝑖 + 2𝑡𝑡𝑡𝑡 + 𝑡𝑡𝑡𝑡) + (3𝑡𝑡𝑡𝑡 + 5𝑡𝑡 2 𝑗𝑗 + 6𝑘𝑘). (2𝑡𝑡𝑡𝑡 − 2𝑗𝑗 + 𝑘𝑘)
= 6 − 21𝑡𝑡 2
Let the curve 𝐶𝐶 represented in parametric form as 𝑟𝑟 = 𝑟𝑟(𝑡𝑡), 𝑎𝑎 ≤ 𝑡𝑡 ≤ 𝑏𝑏. In Cartesian system, we
have 𝑟𝑟(𝑡𝑡) = 𝑥𝑥(𝑡𝑡)𝑖𝑖 + 𝑦𝑦(𝑡𝑡)𝑗𝑗 + 𝑧𝑧(𝑡𝑡)𝑘𝑘. Then, the length of the curve is given by
6
330 www.AgriMoon.Com
Introduction
𝑏𝑏 2 2 2 𝑏𝑏
′ (𝑡𝑡)� ′ (𝑡𝑡)� ′ (𝑡𝑡)�
𝑙𝑙 = � [�𝑥𝑥 + �𝑦𝑦 + �𝑧𝑧 1/2
] 𝑑𝑑𝑑𝑑 = � [𝑟𝑟 ′ (𝑡𝑡). 𝑟𝑟 ′ (𝑡𝑡)]1/2
𝑎𝑎 𝑎𝑎
Sometimes the notation |𝑟𝑟 ′ (𝑡𝑡)| is also used instead of ||𝑟𝑟 ′ (𝑡𝑡)||.
𝑡𝑡 2 2 2 𝑡𝑡
𝑠𝑠(𝑡𝑡) = ∫𝑎𝑎 [�𝑥𝑥 ′ (𝜉𝜉)� + �𝑦𝑦 ′ (𝜉𝜉)� + �𝑧𝑧 ′ (𝜉𝜉)� ]1/2 𝑑𝑑𝑑𝑑 = ∫𝑎𝑎 ||𝑟𝑟 ′ (𝜉𝜉)|| 𝑑𝑑𝑑𝑑 (41.3.1)
Then, 𝑠𝑠(𝑡𝑡) is the arc length of the curve from its initial point (𝑥𝑥(𝑎𝑎), 𝑦𝑦(𝑎𝑎), 𝑧𝑧(𝑎𝑎)) to an arbitrary
point (𝑥𝑥(𝑡𝑡), 𝑦𝑦(𝑡𝑡), 𝑧𝑧(𝑡𝑡)) on the curve 𝐶𝐶. Therefore, 𝑠𝑠(𝑡𝑡) is the length function. Using relation
(41.3.1), it is possible to solve for 𝑡𝑡 as a function of 𝑠𝑠 , that is 𝑡𝑡 = 𝑠𝑠(𝑡𝑡). Then the curve 𝐶𝐶 can be
parameterised in terms of the arc length 𝑠𝑠 as
41.3.9 Example
=x (t ) a=
cos t , y ( t ) a=
sin t , z ( t ) ct .
Therefore , we have
2𝜋𝜋
𝑠𝑠 = arc length = ∫0 [𝑎𝑎2 𝑠𝑠𝑠𝑠𝑠𝑠2 𝑡𝑡 + 𝑎𝑎2 𝑐𝑐𝑐𝑐𝑐𝑐 2 𝑡𝑡 + 𝑐𝑐 2 ]1/2 𝑑𝑑𝑑𝑑 = (2𝜋𝜋)(𝑎𝑎2 + 𝑐𝑐 2 )1/2
7
331 www.AgriMoon.Com
Introduction
Suggested Readings
Courant, R. and John, F. (1989), Introduction to Calculus and Analysis, Vol. II, Springer-Verlag,
New York.
Jain, R.K. and Iyengar, S.R.K. (2002) Advanced Engineering Mathematics, Narosa Publishing
House, New Delhi.
Jordan, D.W. and Smith, P. (2002) Mathematical Techniques, Oxford University Press, Oxford.
Piskunov, N. (1974) Differentail and Integral Calculus, Vol. II, MIR Publishers, Moscow.
Wylie, C. R. and Barrett, L.C. (2003) Advanced Engineering Mathematics, Tata McGraw-Hill,
New Delhi.
8
332 www.AgriMoon.Com
Module-IV: Vector Calculus
Lesson 42
The gradient of a scalar field 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧), denoted by ∇𝑓𝑓 or grad (𝑓𝑓) is defined as
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
∇f = 𝑖𝑖 𝜕𝜕𝜕𝜕 + 𝑗𝑗 𝜕𝜕𝜕𝜕 + 𝑘𝑘 𝜕𝜕𝜕𝜕
Note that the del operator ∇ operates on a scalar field and produces a vector field.
42.1. 1 Example
Solution
𝜕𝜕 𝜕𝜕
∇f(x, y) = �𝑖𝑖 𝜕𝜕𝜕𝜕 + 𝑗𝑗 𝜕𝜕𝜕𝜕 � (𝑦𝑦 2 − 4𝑥𝑥𝑥𝑥) = −4𝑦𝑦 𝑖𝑖 + (2𝑦𝑦 − 4𝑥𝑥)𝑗𝑗
42.1. 2 Example
1
𝒓𝒓 = 𝑥𝑥𝑥𝑥 + 𝑦𝑦𝑦𝑦 + 𝑧𝑧𝑧𝑧, |𝒓𝒓| = 𝑟𝑟 and 𝑟𝑟̂ = 𝒓𝒓/𝑟𝑟, then show that grad�𝑟𝑟 � = −𝑟𝑟̂ /𝑟𝑟 2 .
Solution
1 𝜕𝜕 𝜕𝜕 𝜕𝜕 1 1 𝜕𝜕𝜕𝜕 1 𝜕𝜕𝜕𝜕 1 𝜕𝜕𝜕𝜕 1 𝑥𝑥 𝑦𝑦
Grad�𝑟𝑟 � = �𝑖𝑖 𝜕𝜕𝜕𝜕 + 𝑗𝑗 𝜕𝜕𝜕𝜕 + 𝑘𝑘 𝜕𝜕𝜕𝜕 � �𝑟𝑟 � = 𝑖𝑖 �− 𝑟𝑟 2 𝜕𝜕𝜕𝜕 � + 𝑗𝑗 �− 𝑟𝑟 2 𝜕𝜕𝜕𝜕 � + 𝑘𝑘 �− 𝑟𝑟 2 𝜕𝜕𝜕𝜕 � = − 𝑟𝑟 2 �𝑟𝑟 𝑖𝑖 + 𝑟𝑟 𝑗𝑗 +
𝑧𝑧
𝑟𝑟
𝑘𝑘�
1 𝒓𝒓 𝑟𝑟̂
= − 𝑟𝑟 2 �𝑟𝑟 � = − 𝑟𝑟 2
333 www.AgriMoon.Com
Gradient and Directional Derivative
Let 𝑓𝑓(𝑃𝑃) = 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) be a differentiable scalar field. Let 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 𝑘𝑘 be a level surface and
𝑃𝑃0 (𝑥𝑥0 , 𝑦𝑦0 , 𝑧𝑧0 ) be a point on it. There are infinite number of smooth curves on the surface passing
through the point 𝑃𝑃0 . Each of these curves has a tangent at 𝑃𝑃0 . The totality of these tangent lines
form a tangent plane to the surface at a point 𝑃𝑃0 . A vector normal to this plane at 𝑃𝑃0 is called the
normal vector to the surface at this point.
Consider now a smooth curve 𝐶𝐶 on the surface passing through a point 𝑃𝑃 on the surface. Let
𝑥𝑥 = 𝑥𝑥(𝑡𝑡), 𝑦𝑦 = 𝑦𝑦(𝑡𝑡), 𝑧𝑧 = 𝑧𝑧(𝑡𝑡) be the parametric representation of the curve 𝐶𝐶. Any point 𝑃𝑃 on 𝐶𝐶
has the position vector 𝑟𝑟(𝑡𝑡) = 𝑥𝑥(𝑡𝑡)𝑖𝑖 + 𝑦𝑦(𝑡𝑡)𝑗𝑗 + 𝑧𝑧(𝑡𝑡)𝑘𝑘. Since the curve lies on the surface, we
have
or ∇𝑓𝑓. 𝑟𝑟 ′ (𝑡𝑡) = 0
Let ∇𝑓𝑓(𝑃𝑃) ≠ 0 and 𝑟𝑟 ′ (𝑡𝑡) ≠ 0. Now 𝑟𝑟 ′ (𝑡𝑡) is a tangent to 𝐶𝐶 at the point 𝑃𝑃 and lies in the tangent
plane to the surface at . Hence ∇𝑓𝑓(𝑃𝑃) is orthogonal to every tangent vector at 𝑃𝑃. Therefore,
∇𝑓𝑓(𝑃𝑃) is the vector normal to the surface 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 𝑘𝑘 at the point 𝑃𝑃.
42.1. 4 Example
We will find a unit normal vector to the surface 𝑥𝑥𝑦𝑦 2 + 2𝑦𝑦𝑦𝑦 = 8 at the point (3, −2,1).
Therefore
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
∇f = 𝑖𝑖 𝜕𝜕𝜕𝜕 + 𝑗𝑗 𝜕𝜕𝜕𝜕 + 𝑘𝑘 𝜕𝜕𝜕𝜕 = 𝑦𝑦 2 𝑖𝑖 + (2𝑥𝑥𝑥𝑥 + 2𝑧𝑧)𝑗𝑗 + 2𝑦𝑦𝑦𝑦
At (3, −2,1), we obtain the normal vector as ∇f(3, −2,1) = 4i − 10j − 4k. The unit normal
vector at (3, −2,1) is given by
4𝑖𝑖−10𝑗𝑗 −4𝑘𝑘 2𝑖𝑖−5𝑗𝑗 −2𝑘𝑘
= .
√16+100+16 √33
2
334 www.AgriMoon.Com
Gradient and Directional Derivative
42.1. 5 Example
Here we will find the angle between the two surfaces 𝑥𝑥𝑥𝑥𝑥𝑥𝑥𝑥 𝑧𝑧 = 𝑦𝑦 2 − 1 and 𝑥𝑥 2 𝑦𝑦 = 2 − 𝑧𝑧 at the
given point (1,1,1).
First note that the angle between two surfaces at a common point is the angle between their
normals at that point. Now we have
𝑓𝑓1 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 𝑥𝑥𝑥𝑥𝑥𝑥𝑥𝑥 𝑧𝑧 − 𝑦𝑦 2 + 1 = 0, ∆𝑓𝑓1 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = (log 𝑧𝑧)𝑖𝑖 − 2𝑦𝑦𝑦𝑦 + (𝑥𝑥/𝑧𝑧)𝑘𝑘
Let 𝑓𝑓 and 𝑔𝑔 be any two differentiable scalar fields. The gradient satidfies the following algebraic
properties,
∆(𝑐𝑐1 𝑓𝑓 + 𝑐𝑐2 𝑔𝑔) = 𝑐𝑐1 ∆𝑓𝑓 + 𝑐𝑐2 ∆𝑔𝑔, where 𝑐𝑐1 , 𝑐𝑐2 are arbitrary constants
Let 𝑏𝑏� = 𝑏𝑏1 𝑖𝑖 + 𝑏𝑏2 𝑗𝑗 + 𝑏𝑏3 𝑘𝑘 be any unit vector. Let 𝑃𝑃0 be any point 𝑃𝑃0 : 𝑎𝑎 = 𝑎𝑎1 𝑖𝑖 + 𝑎𝑎2 𝑗𝑗 + 𝑎𝑎3 𝑘𝑘.
Then, the position vector of any point 𝑄𝑄 on the line passing through 𝑃𝑃0 and in the direction of 𝑏𝑏�
is given by
3
335 www.AgriMoon.Com
Gradient and Directional Derivative
𝑟𝑟 = 𝑎𝑎 + 𝑡𝑡𝑏𝑏� = (𝑎𝑎1 + 𝑡𝑡𝑏𝑏1 )𝑖𝑖 + (𝑎𝑎2 + 𝑡𝑡𝑏𝑏2 )𝑗𝑗 + (𝑎𝑎3 + 𝑡𝑡𝑡𝑡3 )𝑘𝑘 = 𝑥𝑥(𝑡𝑡)𝑖𝑖 + 𝑦𝑦(𝑡𝑡)𝑗𝑗 + 𝑧𝑧(𝑡𝑡)𝑘𝑘
This is, the point 𝑄𝑄(𝑎𝑎1 + 𝑡𝑡𝑏𝑏1 , 𝑎𝑎2 + 𝑡𝑡𝑏𝑏2 , 𝑎𝑎3 + 𝑡𝑡𝑡𝑡3 ) is on this line. Now, the vector formthe point
𝑃𝑃0 to 𝑄𝑄 is given by 𝑡𝑡𝑏𝑏�. Since | 𝑏𝑏�|=1, the distance from 𝑃𝑃0 to 𝑄𝑄 is 𝑡𝑡. Then
𝜕𝜕𝜕𝜕 𝑓𝑓(𝑄𝑄)−𝑓𝑓(𝑃𝑃)
𝜕𝜕𝜕𝜕
= lim𝑡𝑡→0 𝑡𝑡
if it exists, is called the directional derivative of 𝑓𝑓 at the point 𝑃𝑃0 in the direction to 𝑏𝑏� .
𝜕𝜕
Therefore 𝜕𝜕𝜕𝜕 𝑓𝑓(𝑥𝑥(𝑡𝑡), 𝑦𝑦(𝑡𝑡), 𝑧𝑧(𝑡𝑡)) is rate of change of 𝑓𝑓 with respect to the distance 𝑡𝑡.
We have
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝑑𝑑𝑑𝑑 𝜕𝜕𝜕𝜕 𝑑𝑑𝑑𝑑 𝜕𝜕𝜕𝜕 𝑑𝑑𝑑𝑑
𝜕𝜕𝜕𝜕
= 𝜕𝜕𝜕𝜕 𝑑𝑑𝑑𝑑
+ 𝜕𝜕𝜕𝜕 𝑑𝑑𝑑𝑑
+ 𝜕𝜕𝜕𝜕 𝑑𝑑𝑑𝑑
where
𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
,
𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
, are evaluated at 𝑡𝑡 = 0 .
We write
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑
𝜕𝜕𝜕𝜕
= �𝑖𝑖 𝜕𝜕𝜕𝜕 + 𝑗𝑗 𝜕𝜕𝜕𝜕 + 𝑘𝑘 𝜕𝜕𝜕𝜕 � . �𝑖𝑖 𝑑𝑑𝑑𝑑 + 𝑗𝑗 𝑑𝑑𝑑𝑑 + 𝑘𝑘 𝑑𝑑𝑑𝑑 � = ∇𝑓𝑓. 𝑑𝑑𝑑𝑑
𝑑𝑑𝑑𝑑
But 𝑑𝑑𝑑𝑑
= 𝑏𝑏�(a unit vector). Therefore, the directional derivative of 𝑓𝑓 in the direction of 𝑏𝑏� in given
by
which is denoted by 𝐷𝐷𝑏𝑏 (𝑓𝑓). Note that 𝑏𝑏� is a unit vector. If the direction is specified by a vector
𝑢𝑢, then 𝑏𝑏� = 𝑢𝑢/|𝑢𝑢|.
42.2.1 Example
We will determine the directional derivative of 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 𝑥𝑥𝑦𝑦 2 + 4𝑥𝑥𝑥𝑥𝑥𝑥 + 𝑧𝑧 2 at the point
(1,2,3) in the direction of 3𝑖𝑖 + 4𝑗𝑗 − 5𝑘𝑘.
Consider
At the point (1,2,3), we have ∇𝑓𝑓 = 28𝑖𝑖 + 16𝑗𝑗 + 14𝑘𝑘. The unit vector in the given direction is
𝑏𝑏� = (3𝑖𝑖 + 4𝑗𝑗 − 5𝑘𝑘)/5√2.
Therefore
4
336 www.AgriMoon.Com
Gradient and Directional Derivative
1 78
𝐷𝐷𝑏𝑏 (1,2,3) = 5√2 (28𝑖𝑖 + 16𝑗𝑗 + 14𝑘𝑘). (3𝑖𝑖 + 4𝑗𝑗 − 5𝑘𝑘) = 5√2
Suggested Readings
Courant, R. and John, F. (1989) Introduction to Calculus and Analysis, Vol. II, Springer-Verlag,
New York.
Jain, R.K. and Iyengar, S.R.K. (2002) Advanced Engineering Mathematics, Narosa Publishing
House, New Delhi.
Jordan, D.W. and Smith, P. (2002) Mathematical Techniques, Oxford University Press, Oxford.
Piskunov, N. (1974) Differentail and Integral Calculus, Vol. II, MIR Publishers, Moscow.
Wylie, C. R. and Barrett, L.C. (2003) Advanced Engineering Mathematics, Tata McGraw-Hill,
New Delhi.
5
337 www.AgriMoon.Com
Module-IV: Vector Calculus
Lesson 43
43.1.1 Example
Here we will find the divergence of the vector field 𝑣𝑣 = (𝑥𝑥 2 𝑦𝑦 2 − 𝑧𝑧 3 )𝑖𝑖 + 2𝑥𝑥𝑥𝑥𝑥𝑥𝑥𝑥 + 𝑒𝑒 𝑥𝑥𝑥𝑥𝑥𝑥 𝑘𝑘.
𝑖𝑖 𝑗𝑗 𝑘𝑘
𝜕𝜕 𝜕𝜕 𝜕𝜕
Curl 𝑣𝑣 = ∇ × 𝑣𝑣 = � �.
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝑣𝑣1 𝑣𝑣2 𝑣𝑣3
43.2.1 Example
Find the curl of the vector field 𝑣𝑣 = (𝑥𝑥 2 𝑦𝑦 2 − 𝑧𝑧 3 )𝑖𝑖 + 2𝑥𝑥𝑥𝑥𝑥𝑥𝑥𝑥 + 𝑒𝑒 𝑥𝑥𝑥𝑥𝑥𝑥 𝑘𝑘
Solution
338 www.AgriMoon.Com
Divergence and Curl
𝑖𝑖 𝑗𝑗 𝑘𝑘
𝜕𝜕 𝜕𝜕 𝜕𝜕
Curl 𝑣𝑣 = � 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
�
𝑥𝑥 𝑦𝑦 − 𝑧𝑧 3
2 2
2𝑥𝑥𝑥𝑥𝑥𝑥 𝑒𝑒 𝑥𝑥𝑥𝑥𝑥𝑥
Curl(grad𝑓𝑓)=0 or ∇ × ∇𝑓𝑓 = 0
𝑖𝑖 𝑗𝑗 𝑘𝑘
𝜕𝜕 𝜕𝜕 𝜕𝜕
∇ × ∇𝑓𝑓 = ��𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 ��
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝜕𝜕 2 𝑓𝑓 𝜕𝜕 2 𝑓𝑓 𝜕𝜕 2 𝑓𝑓 𝜕𝜕 2 𝑓𝑓 𝜕𝜕 2 𝑓𝑓 𝜕𝜕 2 𝑓𝑓
= 𝑖𝑖 � − � + 𝑗𝑗 � − � + 𝑘𝑘 � − �=0
𝜕𝜕𝜕𝜕𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕𝜕𝜕𝜕𝜕
43.2.4 Example
Solution
𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜕𝜕 𝜕𝜕
∇. (𝑓𝑓𝑓𝑓) = �𝑖𝑖 𝜕𝜕𝜕𝜕 + 𝑗𝑗 𝜕𝜕𝜕𝜕 + 𝑘𝑘 𝜕𝜕𝜕𝜕 � . (𝑓𝑓𝑓𝑓1 𝑖𝑖 + 𝑓𝑓𝑣𝑣2 𝑗𝑗 + 𝑓𝑓𝑣𝑣3 𝑘𝑘) = 𝜕𝜕𝜕𝜕 (𝑓𝑓𝑓𝑓1 ) + 𝜕𝜕𝜕𝜕 (𝑓𝑓𝑣𝑣2 ) + 𝜕𝜕𝜕𝜕 (𝑓𝑓𝑣𝑣3 )
2
339 www.AgriMoon.Com
Divergence and Curl
43.2.5 Example
If 𝒓𝒓 = 𝑥𝑥𝑥𝑥 + 𝑦𝑦𝑦𝑦 + 𝑧𝑧𝑧𝑧, |𝒓𝒓| = 𝑟𝑟, show that div (𝒓𝒓/𝑟𝑟 3 )=0
Solution
𝜕𝜕 𝑥𝑥
𝒓𝒓 𝜕𝜕 𝜕𝜕 𝜕𝜕 𝑥𝑥 𝑦𝑦 𝑧𝑧 � �
∆. �𝑟𝑟 3 � = �𝑖𝑖 𝜕𝜕𝜕𝜕 + 𝑗𝑗 𝜕𝜕𝜕𝜕 + 𝑘𝑘 𝜕𝜕𝜕𝜕 � . �𝑖𝑖 𝑟𝑟 3 + 𝑗𝑗 𝑟𝑟 3 + 𝑘𝑘 𝑟𝑟 3 � = ∑ 𝜕𝜕𝜕𝜕 𝑟𝑟 3
Since 𝑟𝑟 2 = 𝑥𝑥 2 +𝑦𝑦 2 + 𝑧𝑧 2
𝒓𝒓 3 3
Therefore, ∆. �𝑟𝑟 3 � = 𝑟𝑟 3 − 𝑟𝑟 3 = 𝟎𝟎
43.2.6 Example
Solution
𝜕𝜕 𝜕𝜕
(i) curl (𝑓𝑓𝑓𝑓) = ∇ × (𝑓𝑓𝑓𝑓) = ∇ × (𝑓𝑓𝑓𝑓1 𝑖𝑖 + 𝑓𝑓𝑣𝑣2 𝑗𝑗 + 𝑓𝑓𝑣𝑣3 𝑘𝑘) = ∑[𝜕𝜕𝜕𝜕 (𝑓𝑓𝑣𝑣3 ) − 𝜕𝜕𝜕𝜕 (𝑓𝑓𝑣𝑣2 )]
𝜕𝜕 𝜕𝜕 𝜕𝜕
= 𝑓𝑓(curl 𝑣𝑣) + �𝑖𝑖 𝜕𝜕𝜕𝜕 + 𝑗𝑗 𝜕𝜕𝜕𝜕 + 𝑘𝑘 𝜕𝜕𝜕𝜕 � × (𝑣𝑣1 𝑖𝑖 + 𝑣𝑣2 𝑗𝑗 + 𝑣𝑣3 𝑘𝑘)
3
340 www.AgriMoon.Com
Divergence and Curl
𝜕𝜕 𝜕𝜕𝑣𝑣 𝜕𝜕𝑣𝑣2
(iii) grad(div 𝑣𝑣)= ∇ × (∇ × 𝑣𝑣) = �∑ 𝑖𝑖 𝜕𝜕𝜕𝜕 � × �∑ 𝑖𝑖 � 𝜕𝜕𝜕𝜕3 − ��
𝜕𝜕𝜕𝜕
𝜕𝜕 𝜕𝜕 2 𝜕𝜕 2 𝜕𝜕 2
= �∑ 𝑖𝑖 𝜕𝜕𝜕𝜕 � (∇. 𝑣𝑣) − �𝜕𝜕𝑥𝑥 2 + 𝜕𝜕𝑦𝑦 2 + 𝜕𝜕𝑧𝑧 2 � (∑ 𝑖𝑖𝑣𝑣1 )
= ∇(∇. 𝑣𝑣) − ∇2 𝑣𝑣
Suggested Readings
Courant, R. and John, F. (1989), Introduction to Calculus and Analysis, Vol. II, Springer-Verlag,
New York.
Jain, R.K. and Iyengar, S.R.K. (2002) Advanced Engineering Mathematics, Narosa Publishing
House, New Delhi.
Jordan, D.W. and Smith, P. (2002) Mathematical Techniques, Oxford University Press, Oxford.
Piskunov, N. (1974) Differentail and Integral Calculus, Vol. II, MIR Publishers, Moscow.
Wylie, C. R. and Barrett, L.C. (2003) Advanced Engineering Mathematics, Tata McGraw-Hill,
New Delhi.
4
341 www.AgriMoon.Com
Module-IV: Vector Calculus
Lesson 44
Line Integral
44.1 Introduction
Let 𝐶𝐶 be a simple curve. Let the parametric representation of 𝐶𝐶 be written as
since
44.2.1 Example
1
Evaluate ∫𝐶𝐶 (𝑥𝑥 2 + 𝑦𝑦𝑦𝑦)𝑑𝑑𝑑𝑑, where 𝐶𝐶 is the curve defined by 𝑥𝑥 = 4𝑦𝑦, 𝑧𝑧 = 3 form (2, 2 , 3) to
(4,1,3).
Solution
√17 4 3 139√17
Hence ∫𝐶𝐶 (𝑥𝑥 2 + 𝑦𝑦𝑦𝑦)𝑑𝑑𝑑𝑑 = ∫ �𝑡𝑡 2 + 4 𝑡𝑡� 𝑑𝑑𝑑𝑑 = .
4 2 24
Let 𝐶𝐶 be a smooth curve whose parametric representation is given in Eqs. (44.1.1) and
(44.1.2). Let
342 www.AgriMoon.Com
Line Integral
𝑣𝑣(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 𝑣𝑣1 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑖𝑖 + 𝑣𝑣2 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑗𝑗 + 𝑣𝑣3 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑘𝑘
be a vector field that is continuous on 𝐶𝐶. Then, the line integral of 𝑣𝑣 over 𝐶𝐶 is defined by
∫𝐶𝐶 𝑣𝑣. 𝑑𝑑𝑑𝑑 = ∫𝐶𝐶 𝑣𝑣1 𝑑𝑑𝑑𝑑 + 𝑣𝑣2 𝑑𝑑𝑑𝑑 + 𝑣𝑣3 𝑑𝑑𝑑𝑑
𝑑𝑑𝑑𝑑
= ∫𝐶𝐶 𝑣𝑣�𝑥𝑥(𝑡𝑡), 𝑦𝑦(𝑡𝑡), 𝑧𝑧(𝑡𝑡)� . 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 (44.2.1)
Similarly, if 𝑣𝑣 = 𝑣𝑣2 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑗𝑗 or 𝑣𝑣 = 𝑣𝑣3 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑘𝑘, we respectively obtained
𝑑𝑑𝑑𝑑
∫𝐶𝐶 𝑣𝑣. 𝑑𝑑𝑑𝑑 = ∫𝐶𝐶 𝑣𝑣2 𝑑𝑑𝑑𝑑 = ∫𝐶𝐶 𝑣𝑣2 (𝑥𝑥(𝑡𝑡), 𝑦𝑦(𝑡𝑡), 𝑧𝑧(𝑡𝑡)) 𝑑𝑑𝑑𝑑
𝑑𝑑𝑑𝑑
𝑑𝑑𝑑𝑑
and ∫𝐶𝐶 𝑣𝑣. 𝑑𝑑𝑑𝑑 = ∫𝐶𝐶 𝑣𝑣3 𝑑𝑑𝑑𝑑 = ∫𝐶𝐶 𝑣𝑣3 (𝑥𝑥(𝑡𝑡), 𝑦𝑦(𝑡𝑡), 𝑧𝑧(𝑡𝑡)) 𝑑𝑑𝑑𝑑
𝑑𝑑𝑑𝑑.
44.2.2 Example
Evaluate the line integral of 𝑣𝑣 = 𝑥𝑥𝑥𝑥𝑥𝑥 + 𝑦𝑦 2 𝑗𝑗 + 𝑒𝑒 𝑧𝑧 𝑘𝑘 over the curve 𝐶𝐶 whose parametric
representation is given by 𝑥𝑥 = 𝑡𝑡 2 , 𝑦𝑦 = 2𝑡𝑡, 0 ≤ 𝑡𝑡 ≤ 1.
Solution:
1 37
= ∫0 (4𝑡𝑡 4 + 8𝑡𝑡 2 + 𝑒𝑒 𝑡𝑡 )𝑑𝑑𝑑𝑑 = 15 + 𝑒𝑒
44.2.3 Example
Evaluate the integral ∫𝑐𝑐 (𝑥𝑥 2 + 𝑦𝑦𝑦𝑦)𝑑𝑑𝑑𝑑, where 𝐶𝐶 is given by 𝑥𝑥 = 𝑡𝑡, 𝑦𝑦 = 𝑡𝑡 2 , 𝑧𝑧 = 3𝑡𝑡, 1 ≤ 𝑡𝑡 ≤ 2.
Solution:
2 163
We have ∫𝑐𝑐 (𝑥𝑥 2 + 𝑦𝑦𝑦𝑦)𝑑𝑑𝑑𝑑 = 2 ∫1 (𝑡𝑡 2 + 3𝑡𝑡 3 ) 𝑑𝑑𝑑𝑑 = 4
∫𝐶𝐶 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑑𝑑𝑑𝑑 + 𝑔𝑔(𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑑𝑑𝑑𝑑 + ℎ(𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑑𝑑𝑑𝑑
2
343 www.AgriMoon.Com
Line Integral
44.3.1 Example
Evaluate ∫𝐶𝐶 (𝑥𝑥 + 𝑦𝑦)𝑑𝑑𝑑𝑑 − 𝑥𝑥 2 𝑑𝑑𝑑𝑑 + (𝑦𝑦 + 𝑧𝑧)𝑑𝑑𝑑𝑑 , where 𝐶𝐶 is 𝑥𝑥 2 = 4𝑦𝑦, 𝑧𝑧 = 𝑥𝑥, 0 ≤ 𝑡𝑡 ≤ 2.
Solution
𝑡𝑡 2
First we consider parametric form of 𝐶𝐶 as 𝑥𝑥 = 𝑡𝑡, 𝑦𝑦 = , 𝑧𝑧 = 2, 0 ≤ 𝑡𝑡 ≤ 2.
4
Therefore,
2 𝑡𝑡 2 𝑡𝑡 𝑡𝑡 2 10
∫𝐶𝐶 (𝑥𝑥 + 𝑦𝑦)𝑑𝑑𝑑𝑑 − 𝑥𝑥 2 𝑑𝑑𝑑𝑑 + (𝑦𝑦 + 𝑧𝑧)𝑑𝑑𝑑𝑑 = ∫0 ��𝑡𝑡 + 4
� − 𝑡𝑡 2 �2� + � 4 + 𝑡𝑡�� 𝑑𝑑𝑑𝑑 = 3
Let 𝑣𝑣(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) = 𝑣𝑣1 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑖𝑖 + 𝑣𝑣2 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑗𝑗 + 𝑣𝑣3 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑘𝑘 be a vector function defined and
continuous at every point on 𝐶𝐶. Then the line integral of tangential component of 𝑣𝑣 along the
curve 𝐶𝐶 from a point 𝑃𝑃 to the point 𝑄𝑄 is given by
𝑄𝑄
∫𝑃𝑃 𝑣𝑣. 𝑑𝑑𝑑𝑑 = ∫𝐶𝐶 𝑣𝑣. 𝑑𝑑𝑑𝑑 = ∫𝑐𝑐 𝑣𝑣1 𝑑𝑑𝑑𝑑 + 𝑣𝑣2 𝑑𝑑𝑑𝑑 + 𝑣𝑣3 𝑑𝑑𝑑𝑑
Let now 𝑣𝑣 = 𝐹𝐹, a variable force acting on a particle which moves along a curve 𝐶𝐶. Then, the
work 𝑊𝑊 done by the force 𝐹𝐹 in displacing the particle from the point 𝑃𝑃 to the point 𝑃𝑃 along
the curve 𝐶𝐶 is given by
𝑄𝑄
𝑊𝑊 = ∫𝑃𝑃 𝐹𝐹. 𝑑𝑑𝑑𝑑 = ∫𝐶𝐶 ∗ 𝐹𝐹. 𝑑𝑑𝑑𝑑
where 𝐶𝐶 ∗ is the part of𝐶𝐶 , whose initial and terminal point are 𝑃𝑃 and 𝑄𝑄.
Suppose that 𝐹𝐹is a conservative vector field . Then 𝐹𝐹 can be written as 𝐹𝐹 = grad(𝑓𝑓), where 𝑓𝑓
is a scalar potential(field). Then, the work done
3
344 www.AgriMoon.Com
Line Integral
44.4.1 Example
Find the work done by the force 𝐹𝐹 = −𝑥𝑥𝑥𝑥𝑥𝑥 + 𝑦𝑦 2 𝑗𝑗 + 𝑧𝑧𝑧𝑧 in moving a particle over the circular
path𝑥𝑥 2 + 𝑦𝑦 2 = 4, 𝑧𝑧 = 0 form (2,0,0) to (0,2,0).
Solution
The parametric representation of the given curve is 𝑥𝑥 = 2 cot 𝑡𝑡, 𝑦𝑦 = 2 sin 𝑡𝑡, 𝑧𝑧 = 0, 0 ≤ 𝑡𝑡 ≤
𝜋𝜋2 . Therefore, work done 𝑊𝑊 is given by
𝜋𝜋/2
16
� [−4 sin 𝑡𝑡 cos 𝑡𝑡 (−2 sin 𝑡𝑡 ) + 4 𝑠𝑠𝑠𝑠𝑠𝑠2 𝑡𝑡(2𝑐𝑐𝑐𝑐𝑐𝑐) ] 𝑑𝑑𝑑𝑑 =
0 13
44.4.2 Circulation
A line integral of a vector field 𝑣𝑣 around a simple closed curve 𝐶𝐶 is defined as the
circulation of 𝑣𝑣 around 𝐶𝐶.
𝑑𝑑𝑑𝑑
Circulation = ∮𝐶𝐶 𝑣𝑣. 𝑑𝑑𝑑𝑑 = ∮𝐶𝐶 𝑣𝑣. 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 = ∮𝑐𝑐 𝑣𝑣. 𝑇𝑇𝑇𝑇𝑇𝑇,
where 𝑇𝑇 is the tangent vector to 𝐶𝐶. For example, in fluid mechanics, let 𝑣𝑣 represents the
velocity field of a fluid and 𝐶𝐶 be a closed curve in its domain. Then, circulation gives the
amount by which the fluid tends to turn the curve rotating or circulating around 𝐶𝐶. If
∮𝑐𝑐 𝑣𝑣. 𝑇𝑇𝑇𝑇𝑇𝑇 > 0 then the fluid tends to rotate 𝐶𝐶in the anti-clockwise direction, while if
∮𝑐𝑐 𝑣𝑣. 𝑇𝑇𝑇𝑇𝑇𝑇 < 0 , then the fluid tends to rotate 𝐶𝐶 in the clockwise direction perpendicular to𝑇𝑇
at every point on 𝐶𝐶, then ∮𝑐𝑐 𝑣𝑣. 𝑇𝑇𝑇𝑇𝑇𝑇 = 0, that is the curve does not move at all.
Therefore, a differential expression expre 𝑑𝑑𝑑𝑑 = 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑑𝑑𝑥𝑥 + 𝑔𝑔(𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑑𝑑𝑑𝑑 + ℎ(𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑑𝑑𝑑𝑑 is
an exact differential, if there exists a scalar function 𝜙𝜙(𝑥𝑥, 𝑦𝑦, 𝑧𝑧) such that
𝑑𝑑𝑑𝑑 = 𝑓𝑓(𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑑𝑑𝑑𝑑 + 𝑔𝑔(𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑑𝑑𝑑𝑑 + ℎ(𝑥𝑥, 𝑦𝑦, 𝑧𝑧)𝑑𝑑𝑑𝑑.
We now present the result on the independence of the path of a line integral
4
345 www.AgriMoon.Com
Line Integral
44.5.1 Theorem
Let 𝐶𝐶 be a curve in simply connected domain 𝐷𝐷 in space. Let 𝑓𝑓, 𝑔𝑔 and ℎ be continuous
function having continuous first partial derivatives in 𝐷𝐷. Then ∫𝐶𝐶 𝑓𝑓𝑓𝑓𝑓𝑓 + 𝑔𝑔𝑔𝑔𝑔𝑔 + ℎ𝑑𝑑𝑑𝑑 is
independent of path 𝐶𝐶 if and only if the integrand is exact differential in 𝐷𝐷.
44.5.2 Example
𝑥𝑥𝑥𝑥𝑥𝑥 +𝑦𝑦𝑦𝑦𝑦𝑦
Show that ∫𝐶𝐶 is independent of path of integration which does not pass through the
�𝑥𝑥 2 +𝑦𝑦 2
origin. Find the value of the integral from the point 𝑃𝑃(−1,2) to the point 𝑄𝑄(2,3).
Solution
𝑥𝑥 𝑦𝑦
We have 𝑓𝑓(𝑥𝑥, 𝑦𝑦) = and 𝑔𝑔(𝑥𝑥, 𝑦𝑦) =
�𝑥𝑥 2 +𝑦𝑦 2 �𝑥𝑥 2 +𝑦𝑦 2
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
Now 𝜕𝜕𝜕𝜕 = −𝑥𝑥𝑥𝑥/(𝑥𝑥 2 + 𝑦𝑦 2 )3/2 and 𝜕𝜕𝜕𝜕
= −𝑥𝑥𝑥𝑥/(𝑥𝑥 2 + 𝑦𝑦 2 )3/2
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
Since 𝜕𝜕𝜕𝜕
= 𝜕𝜕𝜕𝜕 , the integral is independent of any path of integration which does not pass
through the origin. Also, the integrand is an exact differential. Therefore, there exists a
function 𝜙𝜙(𝑥𝑥, 𝑦𝑦) such that
𝜕𝜕𝜕𝜕 𝑥𝑥 𝜕𝜕𝜕𝜕 𝑦𝑦
= 𝑓𝑓(𝑥𝑥, 𝑦𝑦) = and = 𝑔𝑔(𝑥𝑥, 𝑦𝑦) =
𝜕𝜕𝜕𝜕 �𝑥𝑥 2 +𝑦𝑦 2 𝜕𝜕𝜕𝜕 �𝑥𝑥 2 +𝑦𝑦 2
Integrating the first equation with respect to 𝑥𝑥, we get 𝜙𝜙(𝑥𝑥, 𝑦𝑦) = �𝑥𝑥 2 + 𝑦𝑦 2 + ℎ(𝑦𝑦).
𝜕𝜕𝜕𝜕 𝑦𝑦 𝑦𝑦 𝑑𝑑ℎ 𝑑𝑑ℎ
Substituting in 𝜕𝜕𝜕𝜕 = = + 𝑑𝑑𝑑𝑑 or 𝑑𝑑𝑑𝑑 = 0 or ℎ(𝑦𝑦) = 𝑘𝑘, constant.
�𝑥𝑥 2 +𝑦𝑦 2 �𝑥𝑥 2 +𝑦𝑦 2
Suggested Readings
Courant, R. and John, F. (1989), Introduction to Calculus and Analysis, Vol. II, Springer-
Verlag, New York.
Jain, R.K. and Iyengar, S.R.K. (2002) Advanced Engineering Mathematics, Narosa
Publishing House, New Delhi.
Jordan, D.W. and Smith, P. (2002) Mathematical Techniques, Oxford University Press,
Oxford.
5
346 www.AgriMoon.Com
Line Integral
Piskunov, N. (1974) Differentail and Integral Calculus, Vol. II, MIR Publishers, Moscow.
Wylie, C. R. and Barrett, L.C. (2003) Advanced Engineering Mathematics, Tata McGraw-
Hill, New Delhi.
6
347 www.AgriMoon.Com
Module-IV: Vector Calculus
Lesson 45
45.1 Introduction
The theorem provides a relationship between a double integral over a region and the line
integral over the closed curve C bounding R. Green’s theorem is also called the first
fundamental theorem of integral vector calculus.
∂g ∂f
∫ f ( x, y)dx + g ( x, y)dy =
C
∫∫ ∂x − ∂y dxdy
R
The integration being carried in the positive direction (counter clockwise direction) of C.
Proof: We shall prove Green’s theorem for a particular case of the region R.
R : u1 ( x ) ≤ y ≤ u2 ( x ) , a ≤ x ≤ b
R : v1 ( x ) ≤ x ≤ v2 ( x ) , c ≤ y ≤ d
We obtain
∂g
d v2 ( y ) ∂g d
∫∫R ∂x dxdy
= ∫c v=
∫( x ) ∂x dx dy ∫ [ g (v ( y), y) − g (v ( y), y)]dy
2 1
1 c
d c
= ∫ g (v2 ( y), y)dy + ∫ g (v1 ( y), y)dy =
c d
∫ g ( x, y)dy
C
We obtain
∂f
b u2 ( y ) ∂f b
∫∫R ∂x dxdy
= ∫a u=
∫( x ) ∂y dy dy ∫ [ f ( x, u ( x)) − f ( x, u ( x))]dx
2 1
1 a
b a
∫ f ( x, u2 ( x))dx + ∫ g ( x, u1 ( x))dx =
=
a b
− ∫ f ( x, y )dx
C
348 www.AgriMoon.Com
Green’s Theorem in the Plane
∂g ∂f
∫ f ( x, y)dx + g ( x, y)dy =
C
∫∫ ∂x − ∂y dxdy .
R
∫ (x + y 2 )dx + ( y + 2 x)dy, where C is the boundary of the region in the first quadrant that is
2
Solution: The curves intersect at (0,0) and (1,1). The bounding curve is C. We have
f ( x, y=
) x 2 + y 2 and g ( x, y )= y + 2 x .
C R
1 x 1
x
∫ ∫ (2 − 2 y)dydx =
= ∫ (2 y − y ) |
2
2
dx
0 x2 0 x
1
= ∫ (2
0
x − x − 2 x 2 + x 4=
)dx 11/ 30
C C
The closed path C bounds the region R. Using the Green’s theorem, we obtain
∫ (x − y 3 )dx + ( x + y )dy= ∫∫ (1 + 3 y
2 2
)dxdy.
C R
It is convenient to use polar coordinates to evaluate the integral. The region R is given by
= θ , y r sin θ , 0 ≤ r ≤ 4, −π / 4 ≤ θ ≤ π / 4.
R : x r cos=
Therefore,
2
349 www.AgriMoon.Com
Green’s Theorem in the Plane
π /4 4 π /4
r2 3 4 2 4
∫∫ (1 + 3 y )dxdy =∫ ∫ + θ θ =∫ 2 + 4 r sin θ 0 dθ
2 2 2
(1 3r sin ) rdrd
R −π /4 0 −π /4
π /4 π /4
= ∫ (8 + 192sin θ )dθ = ∫ [8 + 96(1 − cos 2θ )]dθ
2
−π /4 −π /4
π /4
2[104θ − 48sin 2θ ]
= 52π − 96.
=
0
−x
=
45.2.4 Example: Verify the Green’s theorem for f ( x, y ) e= sin y, g ( x, y ) e − x cos y and C
is the square with vertices at (0,0), (π/2,0), (π/2,π/2), (0,π/2).
∫
C
fdx + gdy = ∫ + ∫ + ∫ + ∫
C1 C2 C3 C4
( fdx + gdy )
∫e
−x
(sin ydx + cos ydy ) =
0,
C1
: x π / 2, 0 ≤ y ≤ π / 2 and
along C2=
π /2
∫e ∫e
−π /2
−x
(sin ydx + cos =
ydy ) ydy e −π /2 ,
cos=
C2 0
C3 : y π / 2, π / 2 ≤ x ≤ 0 and
along=
∫e ∫ e dx =
e −π /2 − 1,
−x −x
(sin ydx + cos ydy ) =
C3 π /2
: x 0, π / 2 ≤ y ≤ 0 and
along C4=
∫e ∫ cos ydy =
−x
(sin ydx + cos ydy ) = −1.
C4 π /2
Therefore,
π /2 π /2
∫ fdx + gdy =∫∫ (−2e cos y )dxdy =∫ ∫ (−2e cos ydxdy ) =2(e −π /2 − 1).
−x −x
C R 0 0
3
350 www.AgriMoon.Com
Green’s Theorem in the Plane
∂u
∫ ∂n ds= ∫∫ ∇ udxdy,
2
C R
Solution:
dr dx dy
=
T = i+ j
ds ds ds
dy dx
=
n i− j.
ds ds
∂u
∫ ∂n ds= ∫ ∇u.nds
C C
∂u ∂u ∂y ∂u ∂x ∂u ∂u
∫ ∂n ds =
C
∫ ∂x ∂s − ∂y ∂s ds =
C
∫ − ∂y dx + ∂x dy
C
∂ u ∂ u
2 2
∫∫R ∂x 2 + ∂y 2 dxdy =
= ∫∫R ∇ udxdy.
2
Suggested Readings
Courant, R. and John, F. (1989), Introduction to Calculus and Analysis, Vol. II, Springer-
Verlag, New York.
Jain, R.K. and Iyengar, S.R.K. (2002) Advanced Engineering Mathematics, Narosa
Publishing House, New Delhi.
Jordan, D.W. and Smith, P. (2002) Mathematical Techniques, Oxford University Press,
Oxford.
4
351 www.AgriMoon.Com
Green’s Theorem in the Plane
Piskunov, N. (1974) Differentail and Integral Calculus, Vol. II, MIR Publishers, Moscow.
Wylie, C. R. and Barrett, L.C. (2003) Advanced Engineering Mathematics, Tata McGraw-
Hill, New Delhi.
5
352 www.AgriMoon.Com
Module-IV: Vector Calculus
Lesson 46
Surface Integral
46.1 Introduction
b
The double and triple integrals are the generalizations of the definite integral ∫ f ( x)dx to two
a
and three dimensions respectively. The surface area integral is a generalization of the arc
length integral
b
∫
a
1 + ( y ') 2 dx.
We shall now present a generalization of the line integral ∫ f ( x, y)ds to three dimensions.
C
Let g ( x, y, z ) be a given function defined in the three dimensional space and let S be surface
which is the graph of a function z = f ( x, y ), or y = h1 ( x, z ), or x = h( y, z ). We assume that
(i) g ( x, y, z ) is continuous at all points on S, (ii) S is smooth and bounded and (iii) the
projection R of the surface S on x-y plane, x-z plane, or y-z plane respectively expressed in
the forms as assumed in the proof of the Green’s theorem. For example, the projection R on
the x-y plane can be expressed in the forms
R : u1 ( x ) ≤ y ≤ u2 ( x ) , a ≤ x ≤ b
or R : v1 ( x ) ≤ x ≤ v2 ( x ) , c ≤ y ≤ d .
The surface integral can be defined in a similar way as the double integral is defined.
Subdivide S into n parts S1 , S 2 ,..., S n of areas ∆A1 , ∆A2 ,..., ∆An . The projection R of S is
therefore partitioned into n rectangles R1 , R2 ,..., Rn . We choose an arbitrary point
Pk ( xk , yk , zk ) on each element of the surface area S k and form the sum
n
=In ∑ g ( x , y , z )∆A .
k =1
k k k k
Let n → ∞ , such that the largest element of the surface area shrinks to a point. This implies
that as n → ∞ ,the length of the longest diagonal of the projected rectangles tends to zero. In
the limit as n → ∞ , the sequence {I n } has a limiting value which is independent of the way S
is subdivided and the choice of Pk on S k . This limiting value is called the surface integral of
g ( x, y, z ) over S.
353 www.AgriMoon.Com
Surface Integral
∫∫ g ( x, y, z )dA
= ∫∫ g[ x(u, v), y(u, v), z(u, v)] | r × r | dudv u v
S R*
∫∫ g[ x(u, v), y(u, v), z (u, v)][r rv 2 − (ru .rv ) 2 ]1/2 dudv
2
u
R*
∫∫=
g ( x, y, z )dA ∫∫ g[ x, h ( x, z ), z ][1 + f + f y 2 ]1/2 dxdy
2
1 x
S R
∫∫ g=
( x, y, z )dA ∫∫ g[h( y, z ), y, z ][1 + (h ) + (h1 ) 2 y ]1/2 dxdz
2
1 x
S R
∫∫=
g ( x, y, z )dA ∫∫ g[h( y, z ), y, z ][1 + h + hz 2 ]1/2 dydz
2
y
S R
∫∫ g=
S
( x, y, z )dA ∫∫ g ( x, y, z )dA + ∫∫ g ( x, y, z ) dA +... + ∫∫ g ( x, y, z ) dA .
S1 S2 Sk
2
354 www.AgriMoon.Com
Surface Integral
m = ∫∫ ρ ( x, y, z )dA.
S
I = ∫∫ ρ ( x, y, z )d 2 dA.
S
where d is the distance of the point (x,y,z) from the reference axis l. If the surface is
homogeneous, then ρ ( x, y, z ) = constant and ρ ( x, y, z ) =m/A, where A is the surface area of
S. Then,
m
I= ∫∫
A S
d 2 dA
46.3.1 Example: Find the mass of the surface of the cone z = 2 + x 2 + y 2 , 2 ≤ z ≤ 7, in the
first octant, if the density ρ ( x, y, z ) at any point of the surface is proportional to its distance
from the x-y plane.
x y
z=f ( x, y ) =+
2 x2 + y 2 , f x = , fy =
x2 + y 2 x2 + y 2
x2 y2
dA = 1 + f x + f y dxdy = 1 + 2 + dxdy = 2dxdy.
2
2
x + y 2 x2 + y 2
m= ∫∫ czdA = ∫∫ c[2 +
S R
x 2 + y 2 ] 2dxdy
= c 2 ∫∫ [2 + x 2 + y 2 ]dxdy
R
=
Substituting θ , y r sin θ , 0 ≤ θ ≤ π / 2, we obtain
x r cos=
5 π /2 π /2 5
2 r3
m= c 2 ∫ ∫ (2 + r )rdrdθ= c 2 ∫ r + dθ
0 0 0 3 0
125 π 100 2
= c 2 25 + = π c.
3 2 3
3
355 www.AgriMoon.Com
Surface Integral
hy =
−2 y, hz =
0, (1 + hy 2 + hz 2 )1/2 =
(1 + 4 y 2 )1/2 .
The projection of S on the y-z plane is the rectangle OABC with sides=
y 0,=
y =
6, z 0
and z = 8. Therefore,
6 8
∫∫ ydA =
∫∫ y(1 + 4 y ) dydz =
∫ ∫ y(1 + 4 y ) dydz
2 1/2 2 1/2
S R 0 0
6
(1 + 4 y 2 )3/2 2 248
= 8 = =
[(25)3/2 − 1] .
8(3 / 2) 0 3 3
gradf 1
grad f = 2i + 3 j + 4k , n = = (2i + 3 j + 4k ).
| gradf | 29
Consider the projection of S on the x-y plane. The projection of the portion of the plane ABC
in the first octant is the rectangle bounded by=
x 0,=y 0 and 2 x + 3 y = 12. We have
dxdy dxdy
=
dA = .
n.k 4 / 29
1
Therefore, ∫∫=
S
F .ndA ∫∫
S 29
(12 z + 18 + 12 y )dA.
4
356 www.AgriMoon.Com
Surface Integral
1 1
∫∫ F .ndA
S
= ∫∫
S 29
(54 − 6 x + 3 y )=
dA
4 ∫∫
R
(54 − 6 x + 3 y )dxdy
1
6 (12 − 2 x )/3 6
1
= ∫ ∫ (54 − 6 x + 3 y )dy = dx ∫ (360 − 102 x + 7 x 2 )dx
y 0
4 x 0=
= 60
1 7 36
= 360 x − 51x 2
+ x = 138.
6 3 0
Suggested Readings
Courant, R. and John, F. (1989), Introduction to Calculus and Analysis, Vol. II, Springer-
Verlag, New York.
Jain, R.K. and Iyengar, S.R.K. (2002) Advanced Engineering Mathematics, Narosa
Publishing House, New Delhi.
Jordan, D.W. and Smith, P. (2002) Mathematical Techniques, Oxford University Press,
Oxford.
Piskunov, N. (1974) Differentail and Integral Calculus, Vol. II, MIR Publishers, Moscow.
Wylie, C. R. and Barrett, L.C. (2003) Advanced Engineering Mathematics, Tata McGraw-
Hill, New Delhi.
5
357 www.AgriMoon.Com
Module-IV: Vector Calculus
Lesson 47
Stokes’s Theorem
47.1 Introduction
Let C be a curve in two dimensions which is written in the parametric form r = r ( s ) . Then,
the unit tangent vector to C is given by
dx dy
=
T i+ j
ds ds
dx dy dx dy
Then v.T =
( gi − fj ). i + j =g −f .
ds ds ds ds
This result can be considered as a particular case of the Stokes’s theorem. Extension of the
Green’s theorem to three dimensions can be done under the following generalizations.
(i) The closed curve C enclosing R in the plane → the closed curve C bounding an
open smooth orientable surface S (open two sided surface).
(ii) The unit normal n to C → the unit outward or inward normal n to S.
(iii) Counter clockwise direction of C → the direction of C is governed by the direction
of the normal n to S. If n is taken as outward normal, then C is oriented as right
handed screw and if n is taken as inward normal, then C is oriented as left handed
screw.
47.2.1 Theorem (Stokes’s Theorem): Let S be a piecewise smooth orient able surface
bounded by a piecewise smooth simple closed curve C. Let
v( x, y, z ) =v1 ( x, y, z )i + v2 ( x, y, z ) j + v3 ( x, y, z ) k be a vector function which is continuous and
has continuous first order partial derivatives in a domain which contains S. If C is traversed
in the positive direction, then
358 www.AgriMoon.Com
Stokes’s Theorem
47.2.2 Remark: As in divergence theorem, the theorem holds if the given surface S can be
subdivided into finitely many special surfaces such that each of these surfaces can be
described in the required manner.
47.2.3 Remark: To prove the Stokes’s theorem, it is not necessary that the equation of the
surface should be simultaneously written in the =
forms z f=
( x, y ), y g ( x, z ) and x = h( y, z )
. For example, if we take the question of the surface as z = f ( x, y ) and assume that f ( x, y )
has continuous second order partial derivatives then the theorem can be easily proved.
We know that in rigid body rotation, if v denotes the tangential (linear) velocity of a point on
it, then curl v represents the angular velocity of the uniformly rotating body. We also know
that a line integral of a vector field v around a simple closed curve C defines the circulation
of v around C. For example, if v denotes the velocity of a fluid, then circulation gives the
amount by which the fluid tends to turn the curve by rotating or circulating around C.
Therefore, circulation (line integral) is closely related to curl of the vector field. To see this,
let Cr be a small circle with centre at P* ( x* , y* , z * ) . Then, by Stokes’s theorem, we have
∫ v.dr = ∫∫ curlv.ndA
Cr Sr
where S r is a small surface whose bounding curve is Cr . Let P ( x, y, z ) be any arbitrary point
on Cr . We approximate curlv( P ) ≈ curlv( P* ). Then, we have
= [curlv( P* ).n( P )] Ar
*
where Ar is the surface area of S r . Let the radius r of Cr tend to zero. Then, the
approximation curlv( P ) ≈ curlv( P* ) becomes more accurate and in the limit as r → 0, we get
1
curlv( P* ).n( P ) = lim ∫ v.dr.
*
r →0 Ar Cr
The left hand side of the above equation is the normal component of curl v. The right hand
side of equation is circulation of v per unit area. The left hand side is maximum when the
circle Cr is positioned such that the normal to surface, n( P* ) points in the same direction as
curlv( P* ).
2
359 www.AgriMoon.Com
Stokes’s Theorem
47.2.5 Remark: Stokes’s theorem states that the value of the surface integral is same for any
surface as long as the boundary curve, bounding the projection R on any coordinate plane, is
the same curve C. Hence, in the degenerate case, when S coincides with R, we can take n=k
or j or i depending on whether the projection is taken on the x-y plane or x-z plane or y-z
plane.
Solution: Consider projection of S on the x-y plane. The projection is the circular region
x 2 + y 2 ≤ 16, z =
0 and the bounding curve C is the circle z = 0, x 2 + y 2 = 16.
We have
C C C
=
since z=0. Setting cos θ , y 4sin θ , we obtain
x 4=
2π 2π
3 1
C∫ (3x − y)dx =
∫0 4(3cos θ − sin θ )(−4sin θ )dθ =
−16 ∫ sin 2θ − (1 − cos 2θ ) dθ
0
2 2
1
= 16
= 2π 16π .
2
i j k
∇ × v = ∂ / ∂x ∂ / ∂y ∂ / ∂z = i (−4 yz + 4 yz ) − j (0) + k (1) = k
3 x − y −2 yz 2 −2 y 2 z
Now,
2( xi + yj + zk ) 1 z
=n = ( xi + yj + zk ), (∇ × v=
).n .
2 x2 + y 2 + z 2 4 4
Therefore,
z z dxdy z dxdy
∫∫ (∇ × v=
S
).ndA ∫∫=
4
dA ∫∫ = ∫∫ = ∫∫=
4 n.k
S
4 ( z / 4)
R
dxdy
R R
16π
which is the area of the circular region in the x-y plane. Hence, Stokes’s theorem is proved.
3
360 www.AgriMoon.Com
Stokes’s Theorem
∫ (2 y dx + x dy + zdz
3 3
47.2.7 Example: Evaluate where C is the trace of the cone
C
=z (x + y )
2 2
intersected by the plane x=4 and S is the surface of the cone below z=4.
i j k
curlv =∂ / ∂x ∂ / ∂y ∂ / ∂z =i (0) − j (0) + k (3 x 2 − 6 y 2 ).
2 y3 x3 z
If the outward normal to S is taken, then it points downwards. Then, the orientation of C is
taken in the clockwise direction. Alternatively, if the inward normal to S is taken, then C is
oriented in the counter clockwise direction.
Let f ( x, y, z=
) x 2 + y 2 − z= 0 be taken as the equation of the surface. Then, the normal and
unit normal are given by
xi + yj xi + yj − zk ( xi + yj − zk ) / z xi + yj − zk
=N = −k = and n = except at the
x2 + y 2 z ( x2 + y 2 + z 2 ) / z 2 2z
origin.
(3 x 2 − 6 y 2 ) (3 x 2 − 6 y 2 ) dxdy
We have ∫∫ (∇ × v).ndA =
S
∫∫ −
S 2
dA = − ∫∫ −
S 2 (−1/ 2)
=
since dxdy = (n.k )dA. Therefore, substituting cos θ , y r sin θ , we obtain
x r=
4 0
4 0 4 0
3 3
= ∫ ∫
2 0 2π
[(1 + cos 2θ ) − 2(1 − cos 2θ )]r 3 drd
=θ ∫ ∫
2 0 2π
(3cos 2θ − 1)r 3 drdθ
3 r 4 4 3sin 2θ 0
= −θ 192π .
=
2 4 0 2 2π
=
The bounding curve C is given by x 2 + y 2 = 16, z = 4. Now setting cos θ , y 4sin θ ,
x 4=
4
361 www.AgriMoon.Com
Stokes’s Theorem
∫ 2 y dx + x dy + zdz
= ∫ 2 y dx + x dy
3 3 3 3
C C
0
= ∫π 64[2sin θ (−4sin θ ) + cos3 θ (4 cos θ )]
3
2
We obtain 2π π /2
−256 ∫ [cos θ − 2sin θ ]dθ =
= 4
−1024 ∫ (cos 4 θ − 2sin 4 θ
4
0 0
3 1 π 3 1 π
=
−1024 . . − 2 . . =192π .
4 4 2 4 2 2
Suggested Readings
Courant, R. and John, F. (1989), Introduction to Calculus and Analysis, Vol. II, Springer-
Verlag, New York.
Jain, R.K. and Iyengar, S.R.K. (2002) Advanced Engineering Mathematics, Narosa
Publishing House, New Delhi.
Jordan, D.W. and Smith, P. (2002) Mathematical Techniques, Oxford University Press,
Oxford.
Piskunov, N. (1974) Differentail and Integral Calculus, Vol. II, MIR Publishers, Moscow.
Wylie, C. R. and Barrett, L.C. (2003) Advanced Engineering Mathematics, Tata McGraw-
Hill, New Delhi.
5
362 www.AgriMoon.Com
Module-IV: Vector Calculus
Lesson 48
48.1 Introduction
Let C be a curve in two dimensions which is written in the parametric form r = r ( s ) . Then, the
unit tangent and unit normal vectors to C are given by
dx dy dy dx
T = i+ j, n = i − j.
ds ds ds ds
Then,
dx dy dy dx
fdx + gdy = f + g ds =( gi − fj ). i − j ds =(v.n)ds
ds ds ds ds
where =
v gi − fj. Also
∂g ∂f ∂ ∂
− = i + j .( gi − fj ) =
∇.v
∂x ∂y ∂x ∂y
∫ (v.n)ds= ∫∫ (∇.v)dxdy
C R
The result is a particular case of the Gauss’s divergence theorem. Extension of the Greens’
theorem to three dimensions can be done under the following generalisations.
363 www.AgriMoon.Com
Divergence Theorem of Gauss
Theorem: (Divergence theorem of Gauss) Let D be a closed and bounded region in the three
dimensional space whose boundary is a piecewise smooth surface S that is oriented outward. Let
v( x, y, z ) =v1 ( x, y, z )i + v2 ( x, y, z ) j + v3 ( x, y, z ) k be a vector field for which v1 , v2 and v3 are
continuous first order partial derivatives in some domain containing D. Then,
Remark: The given domain D can be subdivided into finitely many special regions such that each
region can be described in the required manner. In the proof of the divergence theorem, the
special region D has a vertical surface. This type of region is not required in the proof. The
region may have a vertical surface. For example, the region bounded by a sphere or an ellipsoid
has no vertical surface. The divergence theorem holds in all these cases. The divergence theorem
also holds for the region D bounded by two closed surfaces.
or as
Example: Let D be the region bounded by the closed cylinder x 2 + y 2 = 16, z = 0 and z = 4.
2
364 www.AgriMoon.Com
Divergence Theorem of Gauss
4
4 =y 16 − x 2
∫∫∫ (∇.v)dV
=
D
∫∫ ∫
z =0
x=
−4 y =
− 16 − x 2
(6 x + 12 y + 1)dydxdz.
4
4 =y 16 − x 2
∫∫∫ (∇.v)dV
=
D
∫∫ ∫
z =0
x=
−4 y =
− 16 − x 2
(6 x + 12 y + 1)dydxdz.
4 4
=y 16 − x 2
∫∫∫ (∇.v)dV= (4)(2)(2) ∫
D z =0
∫ y =0
dydx= 16 ∫ 16 − x 2 dx
0
1 16 x 4
= 16 x 16 − x 2 + sin −1 = 64π .
2 2 4 0
On S1 =
: z 4,=
n k
∫∫ (v.=
S1
n)dA ∫∫=
zdA 4 ∫∫
S1
= dA
S1
4 (area of circular region with radius 4)=64π.
On S 2 : z = 0, n = −k .
∫∫ (v.n)dA =
S2
=∫∫ − zdA =
0.
S2
2 xi + 2 yj 1
On S3 : x 2 + y 2 = 16, n= = ( xi + yj )
2 x2 + y 2 4
1
∫∫ (v.n)dA =
S3
= ∫∫
4 S3
(3 x3 + 6 y 3 )dA.
Therefore,
4 2π
1
∫∫S (v.n)dA 4=z∫0=θ∫0 [192 cos θ + 348sin θ ]4dθ dz
= 3 3
2π
= 192 ∫ [(cos 3θ + 3cos θ ) + 2(3sin θ − sin 3θ )]
0
3
365 www.AgriMoon.Com
Divergence Theorem of Gauss
Hence, ∫∫ (v.n)=
S
dA ∫∫∫ (∇.v)dV .
D
Divergence theorem can be used to prove some important identities, called Green’s identities
which are of use in solving partial differential equations. Let f and g be scalar functions which
are continuous and have continuous partial derivatives in some region of the three dimensional
space. Let S be a piecewise smooth surface bounding a domain D in this region. Let the
functions f and g be such that v=f grad g Then, we have
∇.( f ∇g ) = f ∇ 2 g + ∇f .∇g
Now, ∇g .n is the directional derivative of g in the direction of the unit normal vector n.
Therefore, it can be denoted by ∂g / ∂n. We have the Green’s first identity as
∂g
∫∫ f (∇g.n)=
dA ∫∫ f = ∫∫∫ ( g∇ f + ∇g .∇f )dV .
2
dA
S S
∂n D
∂f
∫∫ g (∇f .n)=
dA ∫∫ f = dA ∫∫∫ ( g ∇ f + ∇g .∇f )dV .
2
S
∂n
S D
∂g ∂f
∫∫ ( f ∇g − g∇f ).ndA
= ∫∫ f − g dA= ∫∫∫ ( f ∇ g − g∇
2 2
f )dV .
S
S
∂n ∂n D
4
366 www.AgriMoon.Com
Divergence Theorem of Gauss
∂y 2
∫∫ ∇g.ndA =
S
∫∫ S
∂n
dA =∫∫∫
D
∇ gdV .
∂y
∫∫ ∇g.ndA= ∫∫ ∂n dA=
S S
0.
This equation gives a very important property of the solutions of Laplace equation, that is of
harmonic functions. It states that if g ( x, y, z ) is a harmonic function, that is, it is a solution of the
equation
∂2 g ∂2 g ∂2 g
+ + =
0
∂x 2 ∂y 2 ∂z 2
Then, the integral of the normal derivative of g over any piecewise smooth closed orient able
surface is zero.
Suggested Readings
Courant, R. and John, F. (1989), Introduction to Calculus and Analysis, Vol. II, Springer-Verlag,
New York.
Jain, R.K. and Iyengar, S.R.K. (2002) Advanced Engineering Mathematics, Narosa Publishing
House, New Delhi.
Jordan, D.W. and Smith, P. (2002) Mathematical Techniques, Oxford University Press, Oxford.
Piskunov, N. (1974) Differentail and Integral Calculus, Vol. II, MIR Publishers, Moscow.
Wylie, C. R. and Barrett, L.C. (2003) Advanced Engineering Mathematics, Tata McGraw-Hill,
New Delhi.
5
367 www.AgriMoon.Com
****** ☺******
This Book Download From e-course of ICAR
Visit for Other Agriculture books, News, Recruitment,
Information, and Events at
www.agrimoon.com
The information on this website does not warrant or assume any legal liability or
responsibility for the accuracy, completeness or usefulness of the courseware contents.
The contents are provided free for noncommercial purpose such as teaching, training,
research, extension and self learning.