Applications of Differential Calculus
Applications of Differential Calculus
8
Applications of Differential
Calculus
Solution Let f (x) = 2x 3 + 7x + 1, then f (x) = 6x 2 + 7. Example 8.7 Using Taylor series expansion, calculate an
For the evaluation of f (x), a convenient value of x close to approximate value of cos 58◦ . Take π = 3.141593.
2.001 is x = 2.
Solution We write f (x) = cos x and express f (x) as Taylor
We have f (2) = 31. π
series about x = , which is very close to 58◦ , as the values
Also, as x increases from x = 2 to x = 2.001, we have 3 π
d x = 0.001. of cos x and its derivatives are known for x = .
3
Now, f (x) = 6x 2 + 7, so, f (2) = 31. Differentiating f (x) successively, we have
Hence, dy ≈ f (x)d x = 31 × 0.001 = 0.031.
f (x) = − sin x, f (2) (x) = − cos x, f (3) (x) = sin x
Thus, for x = 2.001, an approximate value of f (x) is π
31 + 0.031 = 31.031. For x0 = , we have
3
1 √ √
Example 8.4 Find an approximate value of (8.05) 3 . 1 3 1 3
(2) (3)
1 1 2 f (x0 ) = , f (x0 ) = − , f (x0 ) = − , f (x0 ) =
Solution Let f (x) = x 3 , then f (x) = x − 3 . 2 2 2 2
3 Thus, expanding f (x) in Taylor series about x0 , we can write
Take a = 8 and h = 0.05.
√
π
1
Then, f (a) = f (8) = 2 and f (a) = = 0.083. 1 3
12 f (x) = cos x ≈ − x−
Thus, by the approximation formula 2 2 3
π 2 π 3
1
(8.05) 3 = f (8.05) ≈ f (8) + 0.05 × f (8) x − √ x −
1 3 3 3
= 2 + 0.05 × 0.083 − +
2 2! 2 3!
= 2.00415 58π
Example 8.5 Find an approximate value of tan−1 (0.99). Putting x = 58◦ = , so that
180
Solution Let f (x) = tan−1 x. π 2π 3.141593
x− =− =− = −0.0349
Then 3 180 90
1 we have
f (x) = 1.732 1
1 + x2 f (x) = cos x ≈ 0.5 − (−0.0349) − (−0.0349)2
Take a = 1 and h = −0.01. 2 4
1.732
Then + (−0.0349)3
π 12
f (a) = f (1) = tan−1 (1) = radians and 0.00122
4 ≈ 0.5 + 0.866 × 0.0349 −
1 4
f (1) == 0.5 − 0.1441 × 0.00004
2
By the approximation formula f (a + h) ≈ f (a) + h f (a), we ≈ 0.5 + 0.0302 − 0.000305 − 0.000006 ≈ 0.5299
have Thus, for x = 58◦ , cos x ≈ 0.5299.
tan−1 (0.99) ≈ f (1) − 0.01 × f (1) Example 8.8 Using Maclaurin series expansion, calculate
π √
= − 0.01 × 0.5 ≈ 0.7804 radians. an approximate value of 15.
4
Example 8.6 Find an approximate value of sin 31◦ , assum- √ √ 1
Solution We write 15 = 16 − 1 = 4 1 − .
ing that 1◦ = 0.0175 radians and cos 30◦ = 0.8660. 16
1 1
Solution Let f (x) = sin x, then f (x) = cos x. Consider the function f (x) = (1 − x) 2 , where x = .
16
π Then
Take a = 30◦ = radians and h = 1◦ = 0.0175 radians. 1 1
Then
6 f (x) = − (1 − x)− 2
2
f (a) = f (30◦ ) = sin 30◦ = 0.5 and 1 3
f (2) (x) = − (1 − x)− 2
f (a) = f (30◦ ) = cos 30◦ = 0.8660 4
3 5
Using the approximation formula f (a + h) ≈ f (a) + h f (a), f (3) (x) = − (1 − x)− 2
8
we get For x = 0, we have
sin 31◦ ≈ f (30◦ ) + 0.0175 × f (30◦ ) 1 1 3
f (0) = 1, f (0) = − , f (2) (0) = − , f (3) (0) = −
or sin 31◦ ≈ 0.5000 + 0.0175 × 0.8660 = 0.5152. 2 4 8
Thus, expanding in Maclaurin series, we can write
Approximations using Taylor and Maclaurin Series 1 1 1 x2 3 x3
f (x) = (1 − x) 2 ≈ 1 − x− −
The following examples demonstrate the techniques for find- 2 4 2! 8 3!
ing an approximate values of functions using Taylor and x x2 x3
= 1− − −
Maclaurin series expansion. 2 8 16
170 Engineering Mathematics
1 8.2.2 Errors
Putting x = , we have
16
Let y = f (x) be a differentiable function of x and d x an
15 1 1 1 error in x, then the approximate error in y is dy ≈ f (x) d x.
≈1− − − 4
16 32 256 × 8 16 Approximate relative error in y is given by
Thus, we get dy f (x) d x
=
√ 1 1 1 y f (x)
15 ≈ 4 − − −
8 512 4 × 163 Approximate percentage error in y is given by
≈ 4 − 0.125 − 0.00195 − 0.00006 ≈ 3.87299
dy f (x) d x
100 × = 100 ×
Example 8.9 Using Taylor series expansion, calculate an y f (x)
approximate value of e x cos x for x =42◦ . Assume 1◦ = 0.0175
π π π Example 8.10 The radius of a sphere is measured as 15
radian, e 4 = 1.1934, and e 4 cos = 1.551.
4 cm with an error of 0.01 cm. Find (a) the approximate error,
Solution We write f (x) = e x cos x and express f (x) as Tay- (b) relative error, and (c) percentage error in calculating its
lor series. volume.
Using Result 6.10 of successive differentiation, we have
√ n nπ
Solution If r is the radius of the sphere, then its volume is
f (n) (x) = 2 e x cos x + for n = 1, 2, 3, . . . 4
V = πr 3 .
4 3
π dV
Thus, taking x0 = 45◦ = , we have Therefore, = 4πr 2 .
4 dr
π dV
π Approximate error in volume, d V = dr = 4πr 2 dr .
f (x0 ) = e 4 cos dr
4 4
√ π π When r = 15 cm and dr = 0.01 cm, we have V = π(15)3
3
f (x0 ) = 2e 4 cos =0
2 cubic cm and d V = 4π(15)2 (0.01) cubic cm.
π
π 3π π Thus, the relative error is given by
f (2) (0) = 2e 4 cos = −2e 4 sin
4 4 dV 4π(15)2 (0.01) 3 × 0.01
= = = 0.002
√ π 4π √ π V 4 15
f (3) (0) = 2 2e 4 cos = −2 2e 4 π(15)3
4 3
π Percentage error is given by
Thus, expanding f (x) in Taylor series about x0 = , we can
4 dV
write 100 × = 100 × 0.002 = 0.2 .
V
f (x) = e x cos x
π 2 Example 8.11 The time T taken for a completeoscillation
π
π π
π x − l
≈ e 4 cos − 2e 4 sin 4 of a simple pendulum of length l is given by T = 2π , where
4 4 2! g
π 3
g is a constant. If there is an error of 1% in the value of l, find
π
x− the approximate error and percentage error in the calculated
−2e 4 4
3! value of T .
π π
π π π 2 l π
= e 4 cos − e 4 cos x− Solution We have T (l) = 2π , and hence, T (l) = √ .
4 4 4 g gl
√
2 π π 3 Now, percentage error in l is 100 ×
dl
= 1, and hence,
− e4 x − l
3 4 l
π dl = .
Putting x = 42◦ so that x − = 3 × 0.0175 = 0.0525 radian, 100
4 Therefore, approximate error in T is
we have
π l π l
e x cos x ≈ 1.551 − 1.551 × (0.0525)2 dT = T (l) dl = √ =
1.414 × 1.1934 g l 100 100 g
− (0.0525)3
3 Percentage error in T is
≈ 1.551 − 1.551 × 0.00276 − 0.000057
dT π l 1
≈ 1.551 − 0.00427 − 0.000057 ≈ 1.547 100 × = 100 × = 0.5.
T 100 g l
2π
Thus, for x = 42◦ , e x cos x ≈ 1.547. g
Applications of Differential Calculus 171
the radial vector OP is the limiting value of the angle Example 8.14 For a curve, if f (θ ) = x sin θ + y cos θ and
PQN as Q → P.
f (θ ) = x cos θ − y sin θ , show that its arc-length s is given by
Thus s = f (θ ) + f (θ) + c, where c is a constant.
PN r θ dθ
tan φ = lim = lim =r Solution Differentiating f (θ) = x sin θ + y cos θ w.r.t. θ , we
Q→ P QN r → 0 r dr
Now have
r dθ 2
dx dy
ds f (θ ) = sin θ + x cos θ + cos θ − y sin θ
= 1+ = 1 + tan2 φ = sec φ dθ dθ
dr dr
Thus dx dy
dr dr dθ dθ or sin θ + cos θ = f (θ) − f (θ) = 0 (8.10)
cos φ = and sin φ = cos φ tan φ = r =r dθ dθ
ds ds dr ds
Similarly, differentiating f (θ ) = x cos θ − y sin θ w.r.t. θ, we
Example 8.13 Find the length of the arc in one period have
(t = 0 to t = 2π ) of the cycloid x = a(t − sin t), y = a(1−cos t). dx dy
2π f (3) (θ) = cos θ − x sin θ − sin θ − y cos θ
dθ dθ
Solution Arc-length = (d x)2 + (dy)2
t=0
dx dy
2π 2 2 or cos θ − sin θ = f (3) (θ ) + f (θ ) (8.11)
dx dy dθ dθ
= + dt
0 dt dt Squaring Eqn (8.10) and Eqn (8.11) and then adding, we get
2π 2 2
= a 2 (1 − cos t)2 + a 2 sin2 t dt dx dy 2
0 + = f (θ) + f (3) (θ)
dθ dθ
2π √
=a 2 − 2 cos t dt Thus, the arc-length is given by
0
2π
t 2 2
= 2a sin dt dx dy
0 2 s= + dθ
dθ dθ
π
t
= 4a sin u du putting u =
0 2 = f (θ ) + f (3) (θ ) dθ
π
= 4a (− cos u)|0 = 8a.
= f (θ) + f (θ) + c, where c is a constant.
8.4 Curvature
The shape of a curve depends largely on the rate at which Consider a plane curve ABC as shown in Fig. 8.4. Let
the direction of the tangent changes as the point of contact the tangents at points A and B make angles α and α + α,
moves along the curve. If this rate is high, then the bending respectively, with the positive direction of x-axis.
of the curve is more at the point. This rate of change of Then the average curvature of the arc AB is given by
direction with respect to its arc-length at a point is called
α α
curvature of the curve at that point and is denoted by κ. =
length of arc AB s
For example, a straight line is a curve with curvature equal
to zero and a circle has a fixed curvature. where s = length of arc AB.
Applications of Differential Calculus 175
Illustration 8.14 Consider the parabola y 2 = 4ax. At Illustration 8.15 Consider the curve y = e2x − 1. At the
the origin, the y-axis is tangent to the parabola. origin, neither the x-axis nor the y-axis is tangent to the
Thus, the radius of curvature of the parabola at the ori- curve.
gin is We write f (x) = e2x − 1.
y2 4ax Then, f (0) = 1 − 1 = 0, f (0) = 2e0 = 2, f (0) = 22 e0 =
R = lim = lim = 2a 4, . . . .
x →0 2x x → 0 2x
y→0 y→0 Expanding by Maclaurin’s theorem, we have
Case 3 Neither the x-axis nor the y-axis is tangent to the x x 2 x x2
curve at the origin. y = f (0) + f (0) + f (0) + · · · = 2 + 4 + · · ·
1! 2! 1! 2!
By Maclaurin’s theorem, we have
Here, p = 2 and q = 4.
x x 2 Thus, the radius of curvature at the origin is
y = f (0) + f (0) + f (0) + · · ·
1! 2! 3 √
x 2
x (1 + p 2 ) 2 5 5
= f (0) + f (0) + · · · [since f (0) = 0] R= =
1! 2! q 4
x x2
= p +q + · · · [where p = f (0) and q = f (0)] Polar Form
1! 2! Radius of curvature of a polar curve passing through the
Thus, the radius of curvature at the origin is given by pole can be found as follows.
3 3 3 We know x = r cos θ and y = r sin θ .
1 + {y1 (0)}2 2 1 + { f (0)}2 2 1 + p2 2
R(0, 0) = = =
y2 (0) f (0) q x2 r 2 cos2 θ r
R(0, 0) = lim = lim = lim
x2 x →0 2y r →0 2r sin θ r →0 2θ
where p and q are, respectively, the coefficients of x and y→0 θ →0 θ →0
2
in the Maclaurin’s series for f (x). =
1
lim
dr
2 r →0 d θ
θ →0
a x
2. Find the centre ofcurvature at any point of the catenary curvature.
x
y= e a + e− a . 6. Find the centre of curvature and the equation of the
2 circle curvature of the curve r = 3 cos θ at the point
3. Find the co-ordinates of the centre of curvature of the
parabola y 2 = 4 px (a) at any point on the curve and (3, 0).
(b) at the vertex. Also, find the equation of the circle of 7. For the curve defined by the parametric equations
curvature at the vertex. x = 3t, y = t 2 + 2, find the centre of curvature at
4. Find the centre of curvature and the equation of the cir- any point on it.
cle curvature of the curve r = 1 − sin θ at the point (1, 0).
Applications of Differential Calculus 181
8.8 Envelope
Y
Let f (x, y, α) = 0 be an equation, where α is a parameter.
Corresponding to a particular value of the parameter, there
Envelope
exists a plane curve, i.e., different values of α yield different
B 2 2
curves. Thus, f (x, y, α) = 0 represents a family of curves x 3 ⫹y 3 ⫽1
in the x y-plane, where α is a parameter and therefore, it is
called a one-parameter family of curves. Similarly, an equa- d P
tion of the form f (x, y, α, β) = 0, where α and β are two
parameters, represents a two-parameter family of curves.
O A X
Note If α and β are connected by a relation, i.e.,
one is dependent on the other, the two-parameter family of Fig. 8.9 Envelope of a slipping line
curves reduces to a one-parameter family of curves.
Envelope For a given family of curves, an envelope, if it and
exists, is defined to be a curve which touches every member ∂ x sin θ y cos θ
f (x, y, θ) = − =0
of the family of curves and at each of its points, it is touched ∂θ cos θ
2 sin2 θ
by some member of the family of curves. Solving these equations, we get x = d cos3 θ, y = d sin3 θ .
The following are the steps, without proof, to find the These are the parametric equations of an asteroid.
envelope of a one-parameter family of curves. Eliminating the parameter θ, we can express the equa-
tion of the envelope in Cartesian co-ordinates as
Step 1 Let f (x, y, α) = 0 be the equation, where α is the 2 2
parameter. Construct another equation by differentiating it x3 + y3 = 1
∂
w.r.t. α, i.e., f (x, y, α) = 0. Example 8.18 Find the envelope of co-axial ellipses, the
∂α
Step 2 Eliminate the parameter α from these two equa- sum of whose major and minor axes is a constant, say d.
∂
tions, i.e., f (x, y, α) = 0 and f (x, y, α) = 0. Solution Suppose the major and minor axes of an ellipse
∂α x2 y2
are a and b, respectively, then its equation is 2 + 2 = 1.
Step 3 The resulting equation is the Cartesian equation a b
of the envelope of the family of curves. It is given that a + b = d.
x2 y2
The equation 2 + 2 = 1 represents a one-parameter family
Note If x and y are solved in terms of the parameter a b
α, the parametric equations of the envelope are obtained. of curves, as the two parameters a and b are bounded by the
relation a + b = d. Let t be the parameter.
Illustration 8.18 Consider a segment AB of fixed length Then, we have
d, slipping on x-axis and y-axis as shown in Fig. 8.9. Let
da db db da
us find the envelope of this variable segment. a+b =d ⇒ + =0⇒ =−
dt dt dt dt
The x co-ordinate of A is in [0, d]. When AB makes an
To obtain the parametric equations of the envelope, it is re-
angle θ with AO, the co-ordinates of A are (d cos θ, 0) and
quired to eliminate t (i.e., a and b) from
those of B are (0, d sin θ ).
Thus, the equation of the line AB is x2 y2
f (x, y, t) = + =1 (8.20)
x y a2 b2
f (x, y, θ) = + =d
cos θ sin θ Now
To obtain the parametric equations of the envelope, it is ∂f −2x 2 da −2y 2 db
required to solve for x and y in terms of θ in the following = + =0
∂t a 3 dt b3 dt
system of equations:
x y x2 y2
f (x, y, θ) = + =d or − =0 (8.21)
cos θ sin θ a3 b3
Applications of Differential Calculus 183
dy d2 y
Curve Passing through Origin
Note For a stationary point (where = 0), if
Whether the curve passes through the origin or not can dx dx2
< 0, the curve has a maxima at the point. Similarly, for
be easily found by inspecting the expression f (x, y). If d2 y
f (x, y) does not contain any constant term, then the curve a stationary point, if > 0, the curve has a minima at
dx2
passes through the origin. that point.
If a curve passes though the origin, it is important to find Illustration 8.25 Consider again the curve y = x 3 −12x.
the tangent to the curve at the origin. This can be obtained dy
We have seen that = 3x 2 − 12 = 0 for x = −2 and
by equating the group of lowest degree terms to zero. dx
x = 2, and hence, x = −2 and x = 2 are stationary points
Illustration 8.23 The curve y 2 = 4ax passes through of the curve.
x2 y2 d2 y
the origin, whereas the curve 2 + 2 = 1 does not. Now, = 6x.
a b dx2 2
d y
Sign of First and Second Derivatives For x > 0, is positive and the curve is convex.
dx2
For curve tracing, it is important to know the maxima and d2 y
minima of the curves, as well as the increasing and de- For x < −2, is negative, and hence, the curve is
dx2
creasing properties of it in different sub-intervals. For this concave.
dy d2 y
purpose, it is required to check the first and second deriva- Thus, at x = 2, = 0 and is positive, and hence,
dx dx2
tives of one variable w.r.t. the second one in the following
the curve has a minimum at x = 2. Similarly, the curve has
manner.
a maximum at x = −2.
dy
Sign of
dx
dy Extent of the Curve
From the equation of the curve, find . In an interval
dx For a given curve f (x, y) = 0, solve for y in terms of x
dy
[a, b], if > 0, the curve is increasing in the interval, and find the values of x for which y is defined. This gives
dx
dy the horizontal extent of the curve. Similarly, solve for x in
and if < 0, the curve is decreasing in the interval.
dx terms of y and find the values of y for which x is defined.
dy
If = 0 for x = x0 , then at x = x0 , the curve is sta- This gives the vertical extent of the curve.
dx
tionary.
Multiple Points
Illustration 8.24 Consider the curve y = x 3 − 12x.
dy If multiple branches of a curve pass through a point, it is
We have = 3x 2 − 12. called a multiple point. A point through which r branches
dx
dy of a curve pass is known as a multiple point of order r and
Now, = 0 gives x 2 − 4 = 0, i.e., x = 2 or x = −2.
dx the curve has r tangents at that point.
dy
For x < −2, > 0, and hence, the curve is increasing for A multiple point of order 2 is known as a double point.
dx
x < −2. A double point is further classified as a node, a cusp,
dy or an isolated point according to whether the two tangents
For x > 2, > 0, and hence, the curve is increasing
dx at that point are real and distinct, real and coincident, or
for x > 2.
dy imaginary.
For −2 < x < 2, < 0, and hence, the curve is
dx Multiple points of a curve f (x, y) = 0 can be obtained
decreasing for −2 < x < 2. by solving the equations f (x, y) = 0, f x (x, y) = 0, and
Thus, the curve has maximum at x = −2 and has a mini- f y (x, y) = 0.
mum at x = 2. At a double point, the slopes of the tangents are the
d2 y roots of the following equation in m: f yy (x, y) m 2 +2 f x y (x, y)
Sign of m + f x x (x, y) = 0.
d x2
d2 y Now, the values of m (roots of the equation) are real
In an interval [a, b] for a given curve, if > 0, the curve
dx2 and distinct, real and equal, or imaginary depending on
d2 y 2
is convex in the interval and if 2 < 0, the curve is concave whether f x y (x, y) − f yy (x, y) f x x (x, y) > 0, = 0, or < 0,
dx respectively.
in the interval.
2
d y Thus, a double point (x, y) is, respectively, a node, a
If = 0 for x = x0 , then at x = x0 , the curve changes
dx2 cusp, or an isolated point if
from being a convex one to a concave one, or vice versa. 2
Such a point is known as an inflection point. f x y (x, y) − f yy (x, y) f x x (x, y) > 0, = 0, or < 0.
186 Engineering Mathematics
Illustration 8.26 Consider the equation f (x, y) = x 3 − y 2 (d) Origin does not lie on the curve.
= 0. It describes a plane curve, as shown in Fig. 8.10, known (e) dy
=
d 3
(x − 4x 2 + x + 6) = 3x 2 − 8x + 1 = 0
as the semi-cubical parabola. dx dx √
8 ± 64 − 12
for x = = 2.535, 0.131.
Y 6
We see that the point (1, 0) lies on the curve. Let us y (0.131, 6.06)
check this point of the curve. 6
f x (x, y) = 3x 2 − 3, f y (x, y) = −2y, f x x (x, y) = 6x,
f x y (x, y) = 0, and f yy (x, y) = −2. 5
At the point (1, 0), f x x (x, y) = 6, f x y (x, y) = 0, and
f yy (x, y) = −2. 4
Thus
2 3
f x y (x, y) − f yy (x, y) f x x (x, y) = 02 − (6)(−2) = 12 > 0
Hence, the point (1, 0) is a node on this curve. 2
Example 8.19 Trace the curve y = (x + 1)(x − 2)(x − 3).
1
Solution The curve y = (x + 1)(x − 2)(x − 3) has the fo-
llowing characteristics:
(a) The curve has no symmetry. O 1 2 3 4 x
(b) It has no asymptotes.
(c) The x-intercepts of the curve are −1, 2, 3, and the
y-intercept is 6. Fig. 8.12 Curve for y = (x + 1)(x − 2)(x − 3)
Applications of Differential Calculus 187
2
Example 8.20 Trace the curve x 3 + y 3 = 3ax y. Since f x y − f x x f yy = 9a 2 − 36x y > 0 for x = 0,
y = 0, the point (0, 0) is a node of the curve.
Solution The equation of the curve is f (x, y) ≡ x 3 + y 3 − dy
3ax y = 0. It has the following characteristics: The values of y and for different values of x are not
dx
easy to compute. However, from the above information
(a) As f (y, x) = f (x, y), the curve is symmetric about the
about the curve, the curve can be traced as shown in Fig.
line y = x.
8.13.
(b) (i) The coefficients of highest power of x and also of y
are constants, and hence, the curve does not have any
asymptote parallel to the co-ordinate axes. Y
(ii) Now, let us find oblique asymptotes, if any.
Rewrite the equation of the curve as 3a2 , 3a2
f (x, y) = (x 3 + y 3 ) − 3ax y = 0
y⫽x
If y = mx + c is an oblique asymptote to the curve,
then m is the real value solution of x3⫹y3 ⫽ 3axy
(x 3 + y 3 ) x=1,y=m
= 0, i.e., 1+m 3 = 0, i.e., m = −1
O (0, 0) X
and c is obtained from
d a x⫹y⫽a
(1 + m 3 )c − 3am = 0, i.e., c = = −a .
dm m
Thus, the curve has an oblique asymptote y = −x − a,
i.e., y + x + a = 0.
(c) The curve intersects the co-ordinate axes only at the origin
(0, 0). The points of intersections of the curve with the line
of symmetry y = x are obtained from x 3 + x 3 = 3ax x,
3a
which gives x = 0, . Thus, the points are (0, 0) and
2 Fig. 8.13 Curve for x 3 + y 3 = 3ax y
3a 3a
, .
2 2 1
(d) The origin lies on the curve. Example 8.21 Sketch the curve for the equation y = x + .
x
(e) Equating the coefficient of lowest degree term (i.e., 3ax y)
Solution Rewrite the equation of the curve as f (x, y) =
to zero, we see that at the origin, both the axes x = 0 and
x y − x 2 − 1 = 0.
y = 0 are tangents to the curve.
The curve shows the following characteristics:
Also, differentiating the equation of the curve, we have
dy dy (a) The curve is not symmetric about either the x-axis or
3x 2 + 3y = 3ay + 3ax y-axis. However, it is symmetric about the origin, as
dx dx
f (−x,−y) = (−x)(−y) − (−x)2 − 1 = x y − x 2 − 1 = f (x, y)
dy ay − x 2
or = 2 (b) Equating the coefficient of highest power of y to zero, we
dx y − ax
see that the curve has a vertical asymptote x = 0.
1
3a 3a dy From y = x + , we see that y → −∞ as x → 0 from
At , = −1 x
2 2 dx negative side of x-axis and y → +∞ as x → 0 from
So, the tangent at the point is parallel to the asymptote positive side of x-axis.
y + x + a = 0. Since the coefficient of highest power of x (i.e., coeffi-
dy ay − x 2 x2 x6 x2 cient of x 2 ) is constant, the curve has no horizontal asymp-
(f) = 2 = 0 when y = , i.e., x 3 + 3 = 3ax
dx y − ax a a a tote.
√
6 3 3
i.e., x = 2a x , i.e., x = 3 2a. Now, let us search for oblique asymptotes.
(g) On the curve x 3 + y 3 − 3ax y = 0, x and y cannot be both If y = mx + c is an asymptote, then
negative, because when both x and y are negative, each y 1
m = lim = lim 1 + 2 =1
of the terms x 3 , y 3 , and −3ax y are negative, and hence, x→∞ x x→∞ x
x 3 + y 3 − 3ax y cannot be zero. and
(h) Since f (x, y) = x 3 + y 3 −3ax y = 0, f x (x, y) = 3x 2 −3ay = 1
lim (y − mx) = lim=0
0, and f y = 3y 2 −3ax = 0 has only one solution x = 0, y = x→∞ xx→∞
0, the curve has one multiple point (x = 0, y = 0). Thus, the curve has an oblique asymptote y = x.
188 Engineering Mathematics
(h) Since f (x, y) ≡ y(x − 2) − x 2 − 3x + 1 = 0, f x (x, y) = line passing through the pole and perpendicular to the ini-
y − 2x − 3 = 0, and f y = x − 2 = 0 has no solution, the tial line if
curve has no multiple point. π π
F r, − θ = F r, + θ or F(r, θ) = F(r, π − θ )
dy 2 2
(i) Signs of y and are as follows.
dx A curve F(r, θ) = 0 is symmetric about the pole if F(−r,
θ) = F(r, θ ), i.e., if the equation contains only even powers
Sign
dy
Increasing/ of r .
Interval of x of y Quadrant Sign of decreasing
dx Illustration 8.27 The curve r = 2 cos 3θ is symmetric
−∞ < x < −3.3 −ve Third +ve Increasing about the initial line, whereas the curve r = sin 2θ is sym-
−3.3 < x < −1.0 +ve Second +ve Increasing metric about the line passing through the pole and perpen-
−1.0 < x < 0.0 +ve Second −ve Decreasing dicular to the initial line.
0.0 < x < 0.303 +ve First −ve Decreasing Now consider the curve r 2 = 3 cos2 θ. It is symmetric
0.303 < x < 2 −ve Fourth −ve Decreasing about the pole.
2<x <5 +ve First −ve Decreasing
5<x <∞ +ve First +ve Increasing Pole on the curve
If F(0, θ0 ) = 0 for some θ0 , then the pole lies on the curve.
From all the above information, the curve can be traced as
The tangent to the curve at pole is the line θ = θ0 .
shown in Fig. 8.15.
Illustration 8.28 Consider the polar curve r = 2 cos 3θ.
π
The pole lies on r = 2 cos 3θ, as r = 0 for θ = . Thus,
π 2
50 θ= is a tangent to the curve at the pole.
2
40 Similarly, the pole lies on r = sin 2θ , as r = 0 for θ = 0.
Thus, θ = 0 is tangent to the curve at the pole.
30 Now, consider the curve r 2 = a 2 . The pole does not lie
on it.
20
10 Points of Intersection
Intersection points of the curve with the initial line are
obtained by putting θ = 0 in the equation of the curve and
–4 –2 0 2 4 5 6 8
solving for r .
–10
Intersection points of the curve with any inclined line
–20 θ = θ1 are obtained by putting θ = θ1 in the equation of
the curve and solving for r .
–30
Illustration 8.29 Consider the polar curve r = 2 cos 3θ.
Putting θ = 0, we have r = 2. So, the curve cuts the initial
x 2 + 3x − 1
Fig. 8.15 Curve for y = line at the distance 2.
x −2 π
Consider an inclined line θ = .
π 4
Putting θ = in the equation, we have
8.9.2 Polar Equation 4 π √
3π
The general equation of a curve in polar co-ordinates (r, θ ) r = 2 cos = −2 cos =− 2
4 4 √ π
is F(r, θ) = 0, which may be expressed in explicit form as So, the curve intersects the inclined line at − 2, .
r = f (θ ) or θ = g(r ). Let us define below a few character- 4
istics of the curve as expressed in polar form. Asymptotes
The curve has an asymptote if there exists a θ = θ1 such
Symmetry that lim r = ∞. If the equation of the curve can be ex-
θ→θ1
A curve F(r, θ) = 0 is symmetric about the initial line, 1
i.e., θ = 0 if F(r, −θ ) = F(r, θ). pressed in the explicit form r = , θ1 can be found solv-
f (θ )
The curve is symmetric about a line θ = α passing ing the equation f (θ ) = 0 and the equation of the asymp-
through the pole if F(r, α − θ ) = F(r, α + θ ) or F(r, θ) = tote corresponding to θ = θ1 is given by r sin(θ − θ1 ) =
π
F(r, 2α−θ ). Thus, the curve is symmetric about θ = ,the f (θ1 ).
2
190 Engineering Mathematics
Illustration 8.30 Consider the curve with the equation region bounded by two concentric circles is generally called
3 an annulus.
r= .
1 + 2 cos θ The extent of θ can be found from the condition that
1 + 2 cos θ
Writing f (θ ) = , f (θ ) = 0 implies either θ = for these values of θ, the values of r are real.
3
If r = f (θ ) is a periodic function of θ , to find the extent
2π 2π
= θ0 , say, or θ = − = θ1 , say. of r , it is sufficient to consider one period of the function
3 3
Thus, θ = θ0 and θ = θ1 give the directions of the two only.
3
asymptotes to the curve r = .
1 + 2 cos θ Illustration 8.31 Consider the curve r = 3 + 2 cos θ.
−2 sin θ 2π 1
Now, f (θ ) = . Thus, f (θ0 ) = f = −√ The minimum and maximum values of cos θ are −1 and 1,
3 3 3
2π 1 respectively, and hence, values of r lie between 1 and 5.
and f (θ0 ) = f − =√ . The curve lies within the annulus region bounded by the
3 3
Therefore, the asymptotes are given by circles of radii 1 and 5.
2π 2π 1
r sin θ − = f = −√
3 3 3 Derivative
If
dr
> 0, then r increases as θ increases and if
dr
< 0,
2π 2π 1
and r sin θ + = f − =√ dθ dθ
3 3 3 then r decreases as θ increases.
(c) The curve has no asymptote, as r is finite for every value The values of r w.r.t θ may be noted as follows.
of θ .
(d) Putting θ = 0 in the equation of the curve and solving for Value Value
r , we find that the points of intersection of the curve with of θ cos θ r=a(1+cos θ ) of θ cos θ r=a(1+cos θ)
the initial line are (−a, 0) and (a, 0).
π π π π 7π
(e) cos 2θ ≥ 0 for − ≤ 2θ ≤ , i.e., for − ≤ θ ≤ and 0 1.000 2.000 a
6
–0.866 0.134 a
2 2 4 4
π 3π π 3π π 5π
cos 2θ < 0 for < 2θ < , i.e., for <θ < . 0.866 1.866 a –0.707 0.293 a
2 2 4 4 6 4
The equation of the curve is r 2 = a 2 cos 2θ , and hence, π 4π
π 0.707 1.707 a –0.500 0.500 a
the curve does not exist between the lines θ = and 4 3
4
3π π 3π
θ= . 0.500 1.500 a 0.000 1.000 a
4 3 2
Using the above information about the curve, the curve π 5π
can be drawn as shown in Fig. 8.17. 0.000 1.000 a 0.500 1.500 a
2 3
2π 7π
–0.500 0.500 a 0.707 1.707 a
3 4
3π 11π
–0.707 0.293 a 0.866 1.866 a
⫽Ⲑ2 4 6
5π
⫽Ⲑ4 –0.866 0.134 a 2π 1.000 2.000 a
6
π –1.000 0.000
Using the above information about the curve, the curve can
r 2 ⫽a2 cos2 be drawn as shown in Fig. 8.18.
(⫺a, 0) (a, 0) X
1.50a, 3
O
1.707a, 4
a, 2
1.866a, 6
0.50a, 2
3
⫽⫺Ⲑ4
5
⫽
6
Fig. 8.17 Curve for r 2 = a 2 cos 2θ
O a 2a X
7
⫽
6
Example 8.24 Sketch the curve r = a(1 + cos θ ). [This
4
curve is known as cardioid.] ⫽
3 11
⫽
6
Solution The equation of the curve can be rewritten as 5
⫽
F(r, θ ) = r − a(1 + cos θ) = 0. The following observation can 3
be made about the curve: Fig. 8.18 Curve for r = a(1 + cos θ)
(a) The curve is symmetric about the initial line, as F(r, −θ ) =
F(r, θ ) 8.9.3 Parametric Equations
(b) r = 0 for cos θ = −1, i.e., θ = π . Hence, the pole lies on
the curve. Parametric equations of a plane curve may be written as
(c) The curve has no asymptote, as |r | = |a(1 + cos θ )| = x = f (t), y = g(t), t being a parameter. Sometimes para-
2|a| cos2 θ2 ≤ 2a for every value of θ. metric equations may be reduced to a Cartesian equation
(d) Putting θ = 0 in the equation of the curve and solving for by eliminating the parameter t between x = f (t) and
r , we find that the point of intersection of the curve with y = g(t). In that case, the characteristics of the curve as
the initial line is (2a, 0). discussed in Section 8.9.1 may be used for curve tracing.
(e) Since |r | = |a(1 + cos θ)| ≤ 2a, curve lies within the circle Let us, therefore, consider in this section the cases where
with radius 2a and with centre at the pole. parameter t cannot be eliminated.
192 Engineering Mathematics
Origin O (0, 0) X
Solve the equations f (t) = 0, g(t) = 0 for t. If for any value
of t, both the equations f (t) = 0 and g(t) = 0 are satisfied,
then the origin lies on the curve.
(e) For any finite value of t, the x and y values on the curve (g) For different values of t, values of x and y are as fol-
are finite, hence the curve has no asymptote. lows.
(f) The gradient at any point t on the curve is given by
dy dy/dt −a sin t t t x y t x y
= = = − tan
dx d x/dt a(1 + cos t) 2
t 0 0.00 2.00 a π 3.14 a 0.00
Thus, the tangent at any point is horizontal if tan = 0, π 7π
2 1.02 a 1.87 a 3.17 a 0.13 a
i.e., if t = 0, ±2π, ±4π, . . .. 6 6
t
The tangent is vertical if cot = 0, i.e., if t = ±π, ±3π, π 5π
2 1.49 a 1.71 a 3.22 a 0.29 a
±5π, . . .. 4 4
π 4π
1.91 a 1.50 a 3.32 a 0.50 a
y 3 3
2.0a π 3π
2.57 a 1.00 a 3.71 a 1.00 a
2 2
1.5a 2π 5π
2.96 a 0.50 a 4.37 a 1.50 a
3 3
1.0a 3π 7π
3.06 a 0.29 a 4.79 a 1.71 a
4 4
0.5a 5π 11π
3.12 a 0.13 a 5.26 a 1.87 a
6 6
x π 3.14 a 0.00 2π 6.29 a 2.00 a
0 2a 3.14a 4a 6a
From all the above information, we can draw the curve as in
Fig. 8.20 Curve of x = a(t + sin t), y = a(1 + cos t)
Fig. 8.20.
For sufficiently small h, k, l, . . ., the first-degree terms Solving these two equations, we obtain x = 30, y =
∂f ∂f ∂f −20. (In Illustration 8.40 it is observed that the func-
h +k +l + ··· tion has neither a minimum nor a maximum at the
∂x ∂y ∂z
point x = 30, y = −20.)
govern the sign of the right-hand side, and therefore, of the ∂f ∂f
left-hand side of the above equation. Thus, the conditions = 0 and = 0 are not suffi-
∂x ∂y
It is clear that the sign of the left-hand side will change cient for a maximum or a minimum.
with change of the signs of h, k, l, . . ..
Thus, for a maximum or a minimum of the function 8.10.1 Critical and Saddle Points
f (x, y, z, . . .), we must have
Critical Point A point (a, b) is said to be a critical
∂f ∂f ∂f ∂f
h +k +l + ··· = 0 point of the function f (x, y) if f (a, b) = 0, (a, b) = 0,
∂x ∂y ∂z ∂x
∂f
Since h, k, l, . . . are independent of each other, we must and (a, b) = 0.
∂y
have
∂f ∂f ∂f Figures 8.21 and 8.22 show critical points for a function
= 0, = 0, = 0, . . . (8.26) f (x, y). Fig. 8.21 shows the critical point when a function
∂x ∂y ∂z
If f is a function of n independent variables, then Eqn is minimum at the point P (a, b). Similarly, Fig. 8.22 shows
(8.26) is a system of n simultaneous equations. By solving the critical point when a function is maximum at the point
this system of n simultaneous equations, we shall find the P (a, b).
possible values of x, y, z, . . . for which the function f may
attain one of its maxima or minima.
Thus, Eqn (8.26) forms a set of necessary conditions for
a maximum or a minimum of the function. However, it
should be noted that as these equations are not sufficient, f(x, y)
further tests have to be employed to make any conclusion.
Let us now consider the cases of functions with two vari- P(a, b)
ables. Assume that a, b is a set of values of x, y satisfying
Eqn (8.26).
Illustration 8.37 Consider the functions (a) f (x, y) = Fig. 8.21 Critical point when the function has a mini-
2 − x 2 − y 2 , (b) f (x, y) = 16x + 12y + x 2 + y 2 , and mum at the point P(a, b)
(c) f (x, y) = 20x + 40y − x 2 + 4x y − 2y 2 .
(a) Taking the first partial derivatives, we have
∂f ∂f P(a, b)
= −2x, = −2y
∂x ∂y
∂f ∂f f(x, y)
Putting = 0 and = 0 gives x = 0, y = 0, respec-
∂x ∂y
tively. (Later in Illustration 8.39 it is observed that the
function has a maximum at the point (0, 0).)
(b) Taking the first partial derivatives, we obtain Fig. 8.22 Critical point when the function has a maxi-
∂f ∂f mum at the point P(a, b)
= 16 + 2x and = 12 + 2y
∂x ∂y
∂f ∂f
Solving = 0 and = 0, we get x = −8, y = −6. Saddle Point A critical point (a, b) of a function f (x, y)
∂x ∂y
(Later in Illustration 8.39 it is observed that the func- is said to be a saddle point if f (x, y) has neither a maximum
tion has a minimum at the point (−8, −6).) nor a minimum at (a, b).
(c) Taking the first partial derivatives, we obtain Illustration 8.38 Consider the function f (x, y) = x 2 −y 2
∂f ∂f and the point on the surface corresponding to x = 0, y = 0.
= 20 − 2x + 4y and = 40 + 4x − 4y
∂x ∂y We see that f x (0, 0) = 0 and f y (0, 0) = 0. Hence, (0, 0)
∂f ∂f is a critical point of the function. However, the function
Now = 0 and = 0 gives 20 − 2x + 4y = 0 and has neither a maximum nor a minimum at the point (0, 0);
∂x ∂y
40 + 4x − 4y = 0 respectively. (0, 0) is a saddle point of the function as shown in Fig. 8.23.
Applications of Differential Calculus 195
Lagrange’s multiplier method, another way for optimiza- Illustration 8.43 Consider again the optimization prob-
tion of function with constraint equations, does not require lem of Illustration 8.42 and solve it using Lagrange’s mul-
elimination of dependent variables of the function to be tiplier method.
optimized. Form the Lagrangian function as
L = log x + log y + λ(m − x − y)
8.11.1 Lagrange’s Multiplier Method Taking partial derivatives and equating to zero, we have
Let u = f (x1 , x2 , . . . , xn ) be a function, as defined earlier, ∂L 1
of n variables, m constraint, and only (n − m) independent = −λ=0 (8.27)
∂x x
variables.
∂L 1
The function u has a maximum or a minimum at a point = − −λ=0 (8.28)
∂y y
(x1 , x2 , . . . , xn ) if
and
∂f ∂f ∂f ∂L
du = d x1 + d x2 + . . . + d xn = 0 =m−x−y =0 (8.29)
∂ x1 ∂ x2 ∂ xn ∂λ
Eliminating λ from Eqn (8.27) and Eqn (8.28), we get x =
Also
y.
m
∂ f1 ∂ f1 ∂ f1 Hence, from Eqn (8.29), we get x = y = .
d f1 = d x1 + d x2 + . . . + d xn = 0, 2
∂ x1 ∂ x2 ∂ xn
2
Here, the value of multiplier λ is .
∂ f2 ∂ f2 ∂ f2 m
d f2 = d x1 + d x2 + . . . + d xn = 0, We get the same result by both the methods (Illustra-
∂ x1 ∂ x2 ∂ xn
tions 8.42 and 8.43).
...
∂ fm ∂ fm ∂ fm Example 8.26 Find the optimal value of the function
d fm = d x1 + d x2 + . . . + d xn = 0
∂ x1 ∂ x2 ∂ xn f (x, y) = x 2 −3x y +12x subject to the constraint 2x +3y = 6.
Multiplying these equations by 1, λ1 , λ2 , . . . , λm , respec- Solution Form the Lagrangian function as L = x 2 − 3x y +
tively, and then adding, we get 12x + λ(6 − 2x − 3y). Taking partial derivatives and equating
to zero, we have
P1 d x1 + P2 d x2 + P3 d x3 + . . . + Pn d xn = 0 ∂L
= 2x − 3y + 12 − 2λ = 0 (8.30)
∂x
∂f ∂ f1 ∂ f2 ∂ fm ∂L
where Pr = + λ1 + λ2 + . . . + λm = −3x − 3λ = 0 (8.31)
∂ xr ∂ xr ∂ xr ∂ xr ∂y
and
for r = 1, 2, . . . , n. ∂L
The quantities λ1 , λ2 , . . . , λm are arbitrary, hence, here = 6 − 2x − 3y = 0 (8.32)
∂λ
they are choosen such that Pr = 0 for r = 1, 2, . . . , m. Eliminating y from Eqn (8.30) and Eqn (8.32), we get
The above equation now reduces to
4x − 2λ + 6 = 0 (8.33)
Pm+1 d xm+1 + Pm+2 d xm+2 + · · · + Pn d xn = 0 Solving Eqn (8.31) and Eqn (8.33), we have λ = 1 and x = −1.
Finally substituting these values in Eqn (8.30), we have
8
It does not matter which (n − m) of the n variables are con- y= .
sidered as independent variables; so, assume that x m+1 , xm+2 , 3
Thus, the function f (x, y) = x 2 − 3x y + 12x subject to the
. . . , xn are independent. Since d xm+1 , d xm+2 , . . . , d xn are all
constraint 2x + 3y = 6 has an optimal value when x = −1 and
independent, the coefficients Pm+1 , Pm+2 , . . . , Pn must sep-
8
arately be zero. y= .
3
Thus, we obtain a set of additional (n − m) equations
Pm+1 = Pm+2 = . . . = Pn = 0. Example 8.27 A manufacturer can produce three distinct
From the (m + n) equations f 1 = f 2 = . . . = f m = 0 products in quantities q1 , q2 , and q3 , respectively, and thereby
and P1 = P2 = . . . = Pm = Pm+1 = . . . = Pn = 0, we can derive a profit p(q1 , q2 , q3 ) = 2q1 + 8q2 + 24q3 .
determine the m quantities λ1 , λ2 , . . . , λm and the values of Find the values of q1 , q2 , and q3 that maximize profit
the n variables for which maxima and minima values of the if production is subject to the constraint q12 + 2q22 + 4q32 =
function u = f (x1 , x2 , . . . , xn ) are possible. 4,500,000,000.
198 Engineering Mathematics
Solution The profit p(q1 , q2 , q3 ) = 2q1 + 8q2 + 24q3 has to i.e., 2 + λ 2q1 = 0, 8 + λ 4q2 = 0, and 24 + λ 8q3 = 0 together
be maximized subject to the constraint g(q1 , q2 , q3 ) = q12 + 2q22 with q12 + 2q22 + 4q32 − 4,500,000,000 = 0
+ 4q32 − 4,500,000,000 = 0. Solving q1 , q2 , and q3 in terms of λ, we have
For maximum value of profit p(q1 , q2 , q3 ), q1 , q2 , and q3 1 2 3
must satisfy Lagrange’s conditions, i.e., q1 = − , q2 = − , q3 = −
λ λ λ
∂ p(q1 , q2 , q3 ) ∂g(q1 , q2 , q3 ) Thus, we have
+λ = 0
∂q1 ∂q1 1 4 9
+ 2 2 + 4 2 − 4,500,000,000 = 0
λ2 λ λ
∂ p(q1 , q2 , q3 ) ∂g(q1 , q2 , q3 )
+λ = 0 1 1
∂q2 ∂q2 or λ2 = or λ = − [q1 , q2 , and q3
100,000,000 10,000
∂ p(q1 , q2 , q3 ) ∂g(q1 , q2 , q3 ) being all positive, λ must be negative.]
+λ = 0
∂q3 ∂q3 Hence, for maximum profit, we must have q1 = 10,000,
together with g(q1 , q2 , q3 ) = 0. q2 = 20,000, and q3 = 30,000.
1 Thus, we have
it is passing through the point 1, on this
3 (ω − ω0 )2
track. X (ω) = X (ω0 )+ X (ω0 )(ω−ω0 )+x (ω0 ) +. . .
1 dy 2!
By differentiation of y = x 3 , we have = x 2 and
3 d x When ω is very near to ω0 , the terms with second and
d2 y higher powers of (ω − ω0 ) can be ignored in this expan-
= 2x.
dx2 sion.
Substituting these in Eqn (8.12), we have Thus
2x 1
K = 3 X (ω) = L + 2 (ω − ω0 ) = 2L(ω − ω0 ).
1 + x4 2 ω0 C
1 2 1 This gives an approximate expression of reactance near
At 1, , K = 3 = √ radian per mile. This is the the resonance frequency.
3 22 2
ofchange of direction of the train at the given point
rate 3. Manufacturing problem A CD (compact disc) com-
1 pany manufactures blank CD of radius 5 cm. If the ra-
1, .
3 dius of a manufactured CD is off by more than 0.02
2. Electrical Engineering When a current of frequency cm, it is automatically rejected. The company needs to
ω is applied to a series R LC circuit, the net reactance know the variation in the area of the discs which pass
1
is given by X = ωL − . Find an approximate but inspection.
ωC
simple expression for X at frequencies near the resonant This actually refers to the error in the area of CD
frequency ω0 . which may pass inspection. Writing r for radius and A
1 for area of a CD, we have A = πr 2 .
At resonant frequency, we have ω0 L = .
ω0 C So, dA = 2πr dr .
1
Now, X = ωL − is function of ω. So, we can write Taking r = 5 cm and dr = 0.02 cm, we have d A =
ωC 44
it as X (ω) and expand it as a Taylor series about the × 5 × 0.02 = 0.6286 sq cm.
7
point ω = ω0 .
EXERCISES
Review Questions
1. For a curve given by the Cartesian equation y = f (x), 6. Define circle of curvature of a curve at a point on it.
write the expression for the rate of change of arc- Explain how it may be found for a curve given by a
length w.r.t. the co-ordinate x. Also express the dif- Cartesian equation.
ferential arc-length of the curve. 7. What is an evolute? State the procedure for finding
2. How is the rate of change of arc-length of a polar curve evolute for a curve in Cartesian co-ordinates.
w.r.t. a polar co-ordinate expressed? Express also the 8. Define an envelope. How can the envelope for a given
differential arc-length. family of curves be found?
3. Define curvature of a curve at a point on it. How is it 9. What is a critical point? Define a saddle point. State
expressed in Cartesian and polar co-ordinate systems? the criteria for a critical point to be a maximum point,
4. What do you mean by radius of curvature? How is it a minimum point, and a saddle point.
related to curvature? 10. State Lagrange’s conditions for optimization.
5. Write the expressions for the radii of curvature for 11. Describe the Lagrange’s multiplier method for opti-
parametric curve and also for the curve expressed by mization with constrains.
implicit equations.
1. (b) 2. (d) 3. (d) 4. (d) 5. (a) 6. (c) 7. (b) 8. (a) 9. (c) 10. (d)