M101 – Engineering Mathematics I 2015 WBUT Question Paper — KnowledgeLab
Solutions manual for Kreyszig's book, Advanced Engineering MathematicsFull description
Book on MathematicsFull description
gate psu ies
EngineeringFull description
Mathematics
Mathematics for engineering
Descripción completa
Selected solution to problems from Advanced Engineering Mathematics 4th edition
gate
Mathematics
Free
Engineering Mathematics, volume 1, Algebra, Trigonometry, Geometry, by GellesanaFull description
MathematicsFull description
Engineering Mathematics for GateDescription complète
Full description
ENGINEERING MATHEMATICS -I SECOND EDITION
P.B. Bhaskar Rao M.Sc., Ph.D. Retd. Professor, Former Chairman, Board of Studies, Department of Mathematics Osmania University Hyderabad
S.K.V.S.Sriramachary M.A., M.Phil., B.Ed. Professor & Head (Retd.) Department of Mathematics University College of Engineering (Autonomous) Osmania University Hyderabad
M. Bhujanga Rao M.Sc., Ph.D. Professor, Dept. of Mathematics University College of Engineering (Autonomous) Director of Centre for Distance Education Osmania University Hyderabad
All rights reserved. No part of this book or parts thereof may be reproduced, stored in a retrieval syste'm or transmitted in any language or by any means, electronic, mechanical, photocopying, recording or otherwise without the prior written permission of the publishers,
Published by :
BSP BS Publications =:;;;;= 4-4-309, Giriraj Lane, Syltan Bazar, Hyderabad - 500 095 - A. P. Phone: 040-23445688
CHAPTER -1 Ordinary Differential Equations of First Order and First Degree ..................................................... 1 CHAPTER -2 Linear Differential Equations with Constant Coefficients and Laplace Transforms ...................... 69 CHAPTER-3 Mean Value Theorems and Functions of Several Variables .............................................. 111 CHAPTER-4 Curvature and Curve Tracing ................................................ 213 CHAPTER-5 Application of Integration to Areas, Lengths, Volumes and Surface areas ........................ 313 CHAPTER-6 Sequences of Series .............................................................. 385 _ CHAPTER-7 Vector Differentiation ............................................................. 475 CHAPTER-8 Laplace Transforms ............................................................... 623
"This page is Intentionally Left Blank"
1 Ordinary Differential Equations of First Order and First Degree 1.1
Introduction Differential euqtions play an important role in many applications in the field of science and engineering, such as (i) problems relating to motion of particles (ii) problems involving bending of beams (iii) stability of electric system, etc. For example, Newton's law of cooling states that the rate of change of temperature of a body varies as the excess temperature of the body to that of its surroundings. If 8(t) is the temperature of the body at time 't' and 8 0 is the temperature of the room
de
in which the body is kept, then dt gives the rate of change of temperature with time.
de dt
= K(8 - 8 0) ; K is constant
Similarly Newton's second law of motion for a particle of mass m moving in a straight line can be written as
d 2x m dt 2
=F
Where m is the mass, x is the distance of the particle at time 't' measured from a fixed origin and F the external impressed force.
Engineering Mathematics - I
2
A differential equation is an equation involving an unknown function and its derivatives. Ifthere is only one independent variable and one dependent variable the equation is called (Ill ordinary differential equation. If there are more than one independent variable the equation is called a partial differential equation as this involves partial derivatives. For example:
The first four equations (a), (b), (c) and (d) are ordinary differential equations and the remaining three are partial differential equations.
Order 0/ a differential equation: The order of a differential equation is the order of the highest ordered derivative appearing in the equation.
0/ a differential equation: The degree of a differential equation is the power to which the highest ordered derivative appears in the equation after clearing the radicals if any.
Degree
Ordinary Differential Equations of First Order and First Degree
3
In the above examples:
Example: Example: Example: Example:
1.2
1.1(a) is a differential equation of order 3 and degree I. 1.1(b) is of third order and fourth degree differential equation. 1.1(c) is a second order, first degree differential equation. 1.1(d) is a second order, second degree ditferential equation.
Example Formation of an ordinary D.E : The differential equations ar~ formed by eliminating all the arbitrary constants involved in the functional relationship between the dependent and independent variables.
th~t are
For example:
y
=
cx2 + c 2 where c is an arbitrary constant.
.... (I)
To eliminate 'c': (only one constant) From(l)
dv _0
dx
=
c.2x+ 0
I (~V c= 2x dx Substitution of c in (I) gives
y
=
d ~ ( dx
_1_(dy )2
_I dy x2+ 2x d\: 4x 2
d);
)2 + 2x 3 2d
=
dx
- 4x2y
0
is the required D.E and y = cx2 + c 2 is called the solution of the D.E.
Note: Depending on the number of constants in the given equation differentiate it as many number oftimes successively. Then the elimination of the arbitrary constants from the resulting equations and the given equation gives the required differential equation whose order is equal to the number of constants.
1.3
Example Eliminate the arbitrary constants a, b from xy + x 2 = aeX + be-X and form the differential equation.
Engineering Mathematics - I
4
Solution: The given equation is xy + .x2 = ae-'" + be-x
..... ( I)
The number of arbitrary constants is two. Differentiating (I) w.r., to 'x' two times successively.
dy x dx
_
+ Y + 2x = ae-' - be-x
d2y
ely dy + - + - + 2.1 d-c dx dx
X ---2
=
aex + be-x
.....(2)
.... (3)
From (I), (2) and (3) el imination of a, b gives the D.E. from (I) and (3) we get
d l )' 2dy . . + - - + 2 = xy + x 2 IS the requIred D_E. dxdx
x
1.4
--?
Example Form the differential equation by eliminating the constants a and b from 1
a.x2 + by =
Solution: Differentiating ax2 +
by =
I w.r.t 'x'
dy 2ax+ 2byd-c
=
0
.... (I)
.... (2)
Again differentiating wr.t., 'x'
d 2y dy dy 2a+2by - ? +2b-.- =0 dxdx dx Elimination of a, b from (I), (2) and (3) gives
x-?
i
x
yYl (Y.h + yl2
-I
0 =0 0
Expanding the determinant we get
2 y x d y +x(d )2 2 dx dx
_ y(dY)=o dx
.... (3)
Ordinary Differential Equations of First Order and First Degree
1.5
5
Example Form the differential equation by eliminating the constants from
Example Form the differential equation of all circles passing through the origin and having their centres on the x - axis.
Solution Take any tangent to the circle as y-axis, the centre lies on x-axis. Let 'a' be the radius of the circle. Then centre is (a, 0) :. Equation of the circle is (x - a)2 +
.r
=
a2
.... (1)
y
x
,
x
y' Fig. 1.1 Differentiating (I) w.r. to 'x'
dy dx
2(x - a)
+ 2y -
x- a
dy -ydx
=
=
0
dy dx
a=x+ y -
.... (2)
From (1) and (2)
[x-(x+ .r
y:lJ' y' =[x+ y:J +
(:r +.r
2XY :
+ x2 -
= x2
+.r
(:r +
2XY ( : )
.r = 0 is the required D.E
<
7
Ordinary Differential Equations of First Order and First Degree
1.7
Example Form the differential equation of all central conics whose axes coincide with the axes of coordinates.
Solution The equation of all central conics whose axes coincide with the axes is ax2 +
bl =
.... (1)
1
Differentiating (1) w.r.t., x
dy 2ax+ 2by dx
=
0
.... (2)
Differentiating (2) w.r.t., 'x' again
(d
2 y )2 d y a+by - 2 +b =0 dx dx
.... (3)
Eliminating a, b from (1), (2)(3)
x2
y2 dy ydx
X 2
y~+ ~ d 2 (d dx dx
r
o =0 0
y Xyd2; +x(d )2 _y(dY)=O is the required D.E. dx· dx dx
=>
Exercise - 1(a) 1. Eliminate the arbitrary constants from the following and find the corresponding differential equation : (i) y
= mx + c (m, c arbitrary constants)
(ii) y
=
a e2x + b e-2x
(a, b arbitrary. constants)
d2
[ Ans : ~ - 4y - 0] dx 2
8
Engineering Mathematics - I
(iii)
y
=
ax + bx 2
(a, b arb. constants)
d2 [ADS: x 2
dy ----? 2x -d + 2y dx x
(iv)
(x - h)2 + (y - kf = a2, (h, k a, b arb. constants)
(v)
y
=
y
=
0]
=
OJ
(a + bx)e-X , (a, b arb. constants)
d 2y 2dy [ Ans: - 2 + +Y dx dx (vi)
=
a sin x + b cosx (a, b arb. constants)
d 2y
[ADS: - l 2 (X
+ Y = 0]
2. Find the differential equation of all circles with centre on the line y = x and having radius' I '.
3. Form the differential equation of all the circles with centre on the line y passing through the origin. [ADS. : (xl + y)
(:
-1) =
2(y -
=
-x and
x) (x + y:)]
4. Find the differentiall equation of all the parabolas with vertex at the origin and foci on the x - axis. [ADS:
y
dy -2xy dx
=
0]
5. Find the differential equation of all parabolas with the origin as focus and axis along x-axis.
dy [ADS: 2x dx + Y
?
dy -
(
dx) -Y
=
0]
Ordinary Differential Equations of First Order and First Degree
9
Methods to Solve
1.8
The differential equations of the first order and of the First Degree:
1.8.1 Separation of Variables Sometimes the differential equation /
dy dx
=
q(x, y)
can be written as f(x) dx + g(y)dy = 0
.... (I)
if the variables can be separated. Integration of (1) gives the solution of the equation.
Jf(x)dx+ Jg(y)dy=c
i.e.,
where c is an arbitrary constant.
1.8.2 Example Solve et"tany dx + (1-~sec2ydy
=
0
Solution The given equation
e-'"tany dx + (1-e-'")sec 2ydy =0 can be rearranged as e( sec 2 y - - d x + - - dy=O
I-eX
tany
2
Integrating
eX Jsec y --dx+ - - dy=c
J I_eX
tany
-Iog( e-'" -I ) + log t~ll1Y = c i.e.,
tany -X-I
e -
=
c or tan y = c(e X -I)
Engineering Mathematics - I
10
1.8.3 Example dy Solve - = 1 + xl + dx
Y + xly
Solution The given differential equation can be written as dy dx
=
(1 +xl) (1 + y)
dy x3 -1- 2 = (1+xl)dx => tan-I y = x + - + C +y 3
1.8.4 Example dy Solve dx - 2xy = x, where YCO) = 1
Solution
dy
- =x(1+2y) dx dy I+2y =xdx On integration
x2
I
"2 Iog(1+2y) = 2 +c Giveny= 1 when x = 0 Substituting in (1) 1
Ordinary Differential Equations of First Order and First Degree
11
1.8.5 Example dy Solve - =(4x+y+ If dx
Solution dy = (4x + y + 1)2 dx Substituting 4x + y + I = t in (I)
.... (I)
-
dy 4+dx i.e.,
dy dx
=
dt dx
=-
dt -4 dx
dt - -4 = t 2 dx Integrating
i.e.,
I t -tan-I-=x+c
2
2
tan-I (4X+;+ I) = 2(x + c) The solution can also be written as 4x + y + I = 2tan(2x + c) where C is an arbitrary constant.
Exercise -1(b)
1.9
Solve the Following Differential Equations 1.
:
=~Y
2. (2 - x)dy - (3 + y) dx = 0
[ADS:
~
+ e-Y
[ADS (3 + y) (2 -x)
= =
c] c]
Engineering Mathematics - I
12
[ ADs: yc
dy 6. (x-y)2 dx
= (a + x) (I - ay) ]
[ADs: (x - y) + log ( x- y-a) x-y+a
= a2
7.
dy - =(3x+4y+ 1)2 dx
8.
dy =(2x+y+ dx
9.
dy dx = tan (x + y)
= x + c]
[ADs: 2(3x+4y+ 1)= .J3tan- 1 2.J3x+c]
If
[ ADs:
1 tan-I (2X + .J2 .J2Y + I )
=
x +c]
[ADs: log[sin(x + y) + cos(x + y)] = x - y + c]
10 dy = 2 . dx (x+2y-3)
[ADs: (x + 2y - 3) - 410g (x + 2y + I)
=x + c ]
1.9.1 Homogeneous Equations The differential equation of the form
dy dx
f(x,y) g(x,y)
where f, g are homogeneous functions of same degree in x, y is called a homogenous differential equation. Such a differential equation can be written as
.... (I)
Substitutingy = vx dy dv - =v+xdx dx
Ordinary Differential Equations of First Order and First Degree The D. E (I) becomes
dv \11 (v) v + x dx = ~(v)
~(v)dv
dx
-; + v~(v)-\lf(v)
0
=
Integration yields the solution
dx
J-
X
+
J"'()~(v)dv( ) = c where v = -y v'" v - \11 v
x
\
1.9.2 Example Solve
(xl +
dy I) -dx
=xy
Solution
· . dy dv S U bstltuttngy = vx, -d = v + x x dx
dv v + xdx
xvx =
2
2
x +v x
dv v+ x dx
=--
dv xdx
v - -2- v
=
2
v 1+ v 2
1+ v
dv _v 3 x-=-dx 1+ v 2
dx J I JI-dv=c ~+ J-+ X v3 v
13
14
Engineering Mathematics - I
v-2 logx + + logv
-2
logxv- -
1
=
e
e
=
2V2
(y)
X2
logx - ~ x 2y i.e.,
2ylogy
=
=
e ~ 2ylogy - xl
2ey + x 2
1.9.3 Example Solve
dy y2 x-+-=y dx x
Solution
dy xl-+y=xy dx i.e.,
Substitutingy = vx, dy dv -=v+xdx dx
dv v+x - =v-vl dx dv x - =-vl dx dx
dv
-+=0 X v2 Integration yields V-2+1
logx+ - - =e -2+1
1
~logx--
v
=c
=
2ey
Ordinary Differential Equations of First Order and First Degree
x logx - y
ylogx -x
=
=
c cy
or
ylogx
=
x + cy
1.9.4 Example Solve [x+ YSin(Yx)]dx = xSin(Yx)dy Solution The given differential equation can be written as dy
x+ ysin(X)
= __
--,---.0.-,.-----.:-
dx
Substitutingy = vx
xSin(j~)
~
dv v+x dx
dy dy - = v + x. dx dx x+ vxsinv xsinsin(v)
= ----
dv l+vsin v+x- = - - dx sinv dv l+vsinv 1 -v=-.x-= dx sin v SI11 v dx sin v dv = x Integrating -cosv = logx + c logx + cos(Yx) = c
15
Engineering Mathematics - I
16
Exercise - 1 (c)
1.10
l.
Solve the Following Equations dy y2 dx - xy-x2
2. (2 - 2xy)dx
=
[ADS : y = ceix ]
(x 2 - 2xy)dy
[ ADS : xy(y - x) = c ]
3. 2xy + ~ - r) dy = 0 dx
[ADS:
r + y = cy ]
[ADS: logx = 2tan- i (Yx) + c]
5. xdy - ydx = ~ x 2 + y2 dx y [ ADS: cos - = logcx ] x
7. xcos(Yx) (ydx + xdy) = ysin (Yx) (xdy - ydx)
8. (ry - x 3 )dy - ~
9. (r + y) dx
10.
1.10.1
X
+ xy) dx
=
0
[ ADS:
= 2xydy
[ADS: sec (Yx) = yxc]
y~ x + y = cx.e 2
2
( ADS: (r -
dy dx = y[logy - logx + 1]
tan
-I(YI) Ix
]
y) = xc ]
( ADS : y = xecx ]
Non-homogeneous Differential Equations The D.E of the form
dy dx
=
ax+by+c Ax+By+C
.... (J)
where a, b, c A, B, C are constants, is called a non-homogeneous dif.ferential equation.
Ordinary Differential Equations of First Order and First Degree
17
Case (i) If
a
b
A
B
--::;:.-
Substituting x = X + h, Y = y + K where h, k are (constants) to be chosen so as to satisfy. ah + bk + c = 0, Ah + Bk + C = 0
Solving these equations, values of h, k are obtained. The given D.E then reduces to a homogeneous D.E. dY
aX +bY
dX
AX +BY
which is then solved taking Y = VX and then substitute X
=
x - hand Y
=
y - k in the solution.
Case (ii) If
a
b
A
B
Then the differential equation will be of the fonn dy _ (ax + by)+c dx - m(ax+by)+C
since Ax + By will be constant m times ax + by. Now substitute ax + by = t,
Differentiation gives a+bdy=dt dx dx
i.e.,
dt --a dy dx -=-dx b
D.E (2) then reduces to dt
--a
t+c
b
mt+C
~=---
.... (2)
Engineering Mathematics - I
18
Then the solution is obtained by using the method of separation of variables.
1.10.2
Example Solve
dy dx
x+2y-3 2x+ y-3
Solution dy x+2y-3 = dx 2x+ y-3
... - (I)
Substituting x = X + h, y = Y + k where h, k are chosen to satisfY h + 2k - 3 = 0 and 2h + k - 3 = 0 solving we get h=I,k=I
.... (2)
i.e., we take x = X + I, Y= Y + I The D.E (I) reduces to dY
dX
X-2Y 2X+Y
=---
dY
dV
Substituting Y = VX, dX = V + X dX dV X+2VX V+X-=--dX 2X+VX V+XdV =1+2V dX 2+V
~
XdV =1-V
dX
2+v dv=dX I-v X dX +(1 +_3_) dv = 0 X v-I
10gX + v + 3 log(v-I) = loge
2
2+V
Ordinary Differential Equations of First Order and First Degree
v + log( v-I )3 + logX = loge v-log(v-I)3 X = loge
7. (4x-6y-l)dx+(3y-2x-2)dy=O 3 [ADS: x - y + 41og(8x - 12y - 5) = c ]
8.
dy dx
x-y+3 2y + 5
= 2x -
[ ADS: (x - y) + log[2 + x - y ]
x = -
2
+ C]
1.12 Linear Differential Equations A differential equation of the from :
+ py = q where p, q are functions of 'x',
alone or constants is said to be a linear differential equation offirst order. Multiplying both sides of the equation by el pta [called the integrating factor (I.F)] we get,
el ptlr dy + el ptlr py = q .e Iptlr dx The left hand side is the differential coefficient of y. e (i) can be written as
d(y.e
1ptlr
)
=
qe 1ptlr
Integrating
which gives the desired solution.
.... (i)
l pcb:
Engineering Mathematics - I
24
Note: In some cases a differential equation can be reduced to the linear form by taking 'y' as independent variable and x as the dependent variable. The D.E is written as
p" q I are functions of y or constants Now the I.F
=
eJ
Pldy
Solution is
x.e
1.12.1
JPldy
=
fql·e
Jpl"Y
•
+c
Example Solve
(I + x 2 )
dy + 2yx - 6x2 = 0 dx
-
Solution Rearranging the given differential equation to the form
dy dx + py = q We have
dy 2x 6x 2 -+--y=-dx I +X2 1+X2 Here
6x 2
2x p=
I+X2
,q=
I+X2
~dx
I.F= e
Jpdx =e J I+x2
log(I+X2 )
I.F = e
=
(1 + x 2 )
Solution is given by
y(1. F) =
fp(I.F)dx + c
Ordinary Differential Equations of First Order and First Degree
1.12.2
Example Solve
xlogx
dy - +y=2Iogx dx
Solution dy y 2 -+--=dx xlogx x
Here
I 2 p= - - andq=xlogx x
Solution is
ylogx =
J~ logxdx + c
ylogx=2
1.12.3
(IogX)2 2
+c
Example Solve
dy Y (xsmx+cosx . ) =I xcosx -+ dx
Solution xsin+ cosx dy -+y.---dx xcosx
p=
xsin x + cosx
xcosx
xcosx
I ,q=--
xcosx
25
Engineering Mathematics - I
26
J
~Slll X~.COS~ dx
I.F
=
e
XCDSX
= e(log(xsecx)
I. F
= e1og(x sec x)
xsecx
=
Solution is
y(xsec x)
xysec x xysecx
1.12.5
f_lxcosx
=
= =
x
x sec x dx + c
2
Jsec xdx E c tanx+c
Example Solve
dy + 2ytanx = sinx given that y ~
=
0 where x
=
~, 3
Solution dy + 2ytanx = sinx dx p
=
2tanx
q
=
sinx ef
2 tan xd!:
IF
=
~
= e210gsecx
ylF=
JqxIF.dx+c 2
ysec 2x = Jsecx.sec xdx + c
Ordinary Differential Equations of First Order and First Degree
ysec 2x ')
Given that y
Jsec x tan xlir + c
=
ysec-x
secx + ('
=
.... (I)
;7j
0 when x =
=
1C
o=
sec - + c => 3
C =
Substituting c = -2 in (I)
ysec 2x
=
secx - 2
is the required solution
1.12.6
Example dy Solve (x + 2.v) -,
=
(X
y
Solution
+
2.v
dx dy
x y
X
=
--- =
y yd dx
2.v
is in the form of
P
-1
I
= -
Y
y
27
'
q
1
=
2,1 Y
-2
Engineering Mathematics - I
28
x-1
f2y
=
y
2
1 x-dy+c y
x
- =y+C Y
1.13.7 Example dy dx
(x + y+ I) -
Solve
=
1
Solution dy dx- X =y+l PI=-lq\,=y+1
IF =
e
f-1dy
e- Y
=
Solution is given by x(IF)
xe-Y
= fql(/F)dy+c
=
f(y + I)e-Y dy + c
or
i.e.,
x +y + 2
=
ce Y
Exercise 1(e)
1.13 Solve the Following Differential Equations I. (I + y)dx
= (t~n-Iy
- x)dx [ADS: xetan- I y = tan-Iyetan-Iy --etan- I y + c]
dy 2. cos2x - + Y
dx
= tanx [ADS: Y = ce-tanx - tanx -1]
29
Ordinary Differential Equations of First Order and First Degree
dy 3. x - +2y-x2 1og=0 dx (ADS: y
4.
dy dx + ycot x = 4x cosec.x, if y
=
c =
-?
x-
I .r2 +- x210gx --I 4
16
n/ 0, when x = ~2 n2 (ADS: ysinx = 2x2 - 2
5. yeYdx
=
(y
+ 2x&)dy
6. (x + 31) dy dx
7. (xy -I )
=
dy dx+ y
3
=c-
( ADS : x
21 + cy ]
+ y3
3
=
= 0
I ADS: x dy dx
e~Y
y
I
8.
I
( ADS : xy~2
= x 3 - 2xy if y = 2 when x =
.
= ce Y
I
+ - + II Y
I
[ADS: 2y - x 2 + 1 = 4 el~~ ]
dy x+ ycosx 10. dx = I +sinx [ADS: y(l+sinx)
=
x2
c - -] 2
30
Engineering Mathematics - I
1.13.1 Non-linear Differential Equation of First Order Ber noulli j. equation: The differential equation of the form ely dx + PY
=
llyn
.... (I)
where p, q are functions ofx alone is said to be a Bernoulli's differential equation. Dividing (I) throughout by yn
dy =q dx Substitutingynt-l(coefficient ofp) "" v y-n_ + py-n+l
.... (2)
1 dy dv (1-n) - - = y" dx dx (2) reduces to I dv - - - +pv=q
(I-n)dl'
dv - + (1-n)pv = (I-n)q dx which is linear in v. The avove D.E can be solved by using the method given in 1.12.1 example.
1.13.2 Example Solve dy - ytanx = ysecx dx Solution
dy
dx - ytanx = ysecx
dy 1 y-2 -d - - tan x = sec x x Y Substituting
-y =v,
.... (I)
Ordinary Differential Equations of First Order and First Degree
1 dy
dv
dx
dx
+--=/
(1) reduces to dv dx + vtan x
=
sec x
is linear in v. Here
p
tanx, q
=
IF
=
e
=
e
secx
=
fpd< flanXd<
= elogsec x =
secx
SolutIon of th [) E (I) IS therefore
v(lF) = Jq(IF}dx + c
v. secx
=
Jsecx.secxdx + C 2
vsecx.= J sec xdx + c Substituting
1
v=-Y 1
- - sec x y i.e.,
=
tan x + c
y(tan x + c) + sec x
1.13.3 Example Solve (3xy + .0) dx - 3i2dy = 0
Solution dy 3x2 - - 3xy dx
=
.0
=
0
31
32
Engineering Mathematics - I
is in Bernoulli's fonn 1 dy _ dx - xy - 3x 2
"7
Substituting
-1 =v
y
D.E reduces to
dv I 1 ---v=dx xy 3x 2 which is linear in v Here
1
1 q= 3x 2
p=~,
IF = e =
e
Jpdx
J!d\ x
=x Solution of the D.E (1) is
v(lF) = fq(lF)dx+c -..!..x= f_l_. xdx + c y 3x 2
x 1 --=-Iogx+c Y 3 i.e.,
y(110g x + c) + x = 0
1.13.4 Example dy Solve dx + (2xtan- 1y -.x3) (1 + y) = 0 Solution
dy 1+ y2 dx
dv
= dx
.... (1)
Ordinary Differential Equations of First Order and First Degree
33
The given differential equation becomes dv
-
d-c
+ 2xv=x3
.... (1)
which is linear in v Here
p = 2x, q = x 3
IF
e
=
J
2xdx
=e
x2
Solution of the D.E (1) is
v(lF) = Writing
Jq(IF}dx+ c
x2 = t xdx
1
= 2dt 2
r
I
dt
v.e x = {e 2+c
1.14.5 Example dy Solve tany dx + tanx
= cosy cos3x
Solution dy tany d'C + tanx
= cosy cos 3x
Dividing by cosy throughout
dy secy tany dx + secy tanx = cos 3x Substituting secy = v, we get
Ordinary Differential Equations of First Order and First Degree
is linear in v
p
Here
=
~
IF
COSX,
= e
q = SIl1 X cos X
fpdr = e feo') '(tl\
=e
~1Il X
Solution of the D.E (I) is v.e SlIlt
=
J. xcosx.e m"'. I + c SII1
(X
sinx = t=> cosx dx = £II in the RHS
write
siny
e;lIlt
sinyeSlI1\"
= [te t =
el] + c
eSItlX [sinx-l] + c
Exercise - 4(f)
1.14 Solve the Following Differential Equations I. (ylogx -I )ydx =
x(~y
1 (Ans: y 2.
dy . - cosx + YSlnx = dx
r.::::::::
( Ans : dy tany ---£Ix l+x
=
1 + logx + ex )
...; ysecx
(Ans :2y y
4.
=
(I+x)
~ y-
= -
12
-Jsecx
=
tanx + 2c )
sin 2x -sinx -~ + ce 2SIIl\"
2
)
~secy
[Ans: siny = (1 + x) (eX + c) )
36
Engineering Mathematics - I
dy 5. x - + ylogy = xyeX dx
I
+c
J
x3
2
dy
ADs : xlogy = (x - I)e"
6. 3-+--y=-) dx x+ I y
I ADs: 7.
dy + ytanx dx
=
(x +
Iii =
dy dx
X
2
6
2x 5
+-
X4
+- +C 5 4
isecx
I ADs: 8.
x6 -
cos2x
=
y(c + 2sin x) J
1
Y +xy
1.14.1 Exact Differential Equations Let us consider the differential equation Mdx + Ndy = 0 where M, N are functions ofx,y. If this equation is to be exact, then it must have been derived by directly differentiating some function ofx,y. Hence
Mdx + Ndy = du, say
.... (I)
But from differential calculus
au
au
du = -dr:+-dy ax ay From (I) and (2) we get
au
N=
M= ax'
Now
aM
a2u
ay
ayax
-=--
au
ay
and -ax- - -ax-ay-
.... (2)
Ordinary Differential Equations of First Order and First Degree
aM
-
oy
oN.IS tIle con d·· ItlOn ..lor exactness.
= -
AX
:. The differential equation Mdx +
N~v =
0
.IS exact 'f aM aN I -=-
oy
ax
Then the solution is expressed in the form
+ (integrate w.r.t y those terms that are independent of x)
(treatingyas constant integrate w.r.t x)
Note: IfN has no term independent ofx then the solution is fMdx
=
e
1.14.2 Example Solve (x + 2y - 3) dy - (2x - Y + \ )dx = 0
Solution (x+2y - 3)dy - (2x - Y + I) dx = 0 M = -(2x - Y+ \) N = (x + 2y - 3)
aM
-=\
oy
aN
-=\
ax
The given differential equation is exact The solution is - f(2x - y + \ )dx +
fe x + 2y - 3)
_2X2 2i - -2- + yx - x + 2
=>
I - x 2 + xy -
x - 3y
=
e
=
- 3y =
(ry = e
e
37
38
Engineering Mathematics - I
1.14.3 Example Solve (x 2+ y.)dr: + 2xy dy = 0
Solution
(x 2 + .V) dx + 2xy dy
=
0
M=x2+y.
N
aM av
-
=
2xy
aM ax
-=2y
=2y
The given differential equation is exact 2
Solution is f(x + y2
}h + f2xydy = c
1.14.4 Example Solve (I + e
-,:~, )dx + /Y
[I.'... ;
1
dy = 0
Solution (I + e'x Y ) dx + e x'Y
(
J -;;, ) dy = 0
x'y M = J + e:t Y , N = e,i
aM
-x
~
i
-=-
[
x)
J- Y
,y' e'Y
oM oN oy
ox
The given differential equation is exact
39
Ordinary Differential Equations of First Order and First Degree
Solution is
X
+ e
x
Y
(y)
=
c
Exercise - 4(g) 1.15 Solve the Following Differential Equations 1. (el' + 1)cosxdx +
1.15.1 Integrating factors If the differential equation Mdx + Ndy = 0 is not exact, it can be made exact by multiplying it with some function of x, y. Such a function is called an integrating
jactor. Rule!t' for fillt/illg the illtegrtltillg factors : 1. Integrating factors found by inspection:
Example Solve x dy- ydx
=
0
Solution
xdy- ydx
=
0
40
Engineering Mathematics - I
Dividing by x2
xdy- ydx
---'-,=-=--
2
x
=
0
On integration
Yx
=c
First method /0 find an integrating jac/or : If the differential equation Mdx + Ndy is not exact, but is homogenous and Mx + Ny
7:-
0, then the integrating factor is
1 . Multiply the differential Mx+Ny
equation by IF. The DE becomes exact.
1.15.2 Example Solve
(x 2y - 2xy2)dx - (x 3
-
3x2y)dy
=0
(x 2y - 2xy)dx - (x 3
-
3x2y)dy
=0
Solution .... (I)
The differential equation (I) is homogeneous
M
= x 2y- 2xy
aM
-
oy
N
= -
aN
=x2-4xy
-
ax
(x 3
-
3x2y)
=-3x2 + 6xy
The DE is not exact and
Mx + Ny = x 2y IF
=
7:-
0
) 1 Mx+ Ny - x 2 y2
Multiplying the DE by the integrating factor ~ • x y
(
y2 X2Y -2X 2
X
Y
2
)d _(X2 -3X3Y )dY X
2
X
Y
2
=
0
.... (2)
Ordinary Differential Equations of First Order and First Degree
1
2
Y
X
-x 3 N =-+I y2 Y
write
M
then
--=-=--
= --I
aMI
aNI
ry
ax
y2
DE (2) is exact Solution is
x - - 210gx + 3logy + c y i.e., x
or
1.15.3 Second Method to Find the Integrating Factor If the differential equation Mdx + Ndy = 0 is not exact and is of the form
MUltiplying (1) by IF (0/,+2 + 30+2 /,+I)dx + (2x h+3/,
oM
-
ox oN ox
-
0+ 1/,+I)dy
=
0
.... (2)
3(k + 2)0/,+1 + 2(k + 1) 0+2yK
=
= 2(h + 3) x h+2/, ---(h + 1)0/,+ 1
The DE(1) is to be exact
oM
ON
ox
ox
2(k + 1) = 2(h + 3) => h = -
%k = -/i
Comparing the coefficients
k + 2 = - (h + J) and
2(k+ 1) = 2(h + 3) => h =
-,% k= -/i
Substituting h, k values in (2) -5/
(
3/
X /2 y/2
II) ( II -II + 2x -1/ 12 yl2 dx + 2X 72 Y 12 I
(3) is an exact DE Solutionis
2
--x 3
-2/
3/
13 yl2
II
II
+ 4X 72 y72 = c
-3/
II) dy = 0
_ X /2 y72
.... (3)
51
Ordinary Differential Equations of First Order and First Degree
Exercise - 4(j) I
Solve the following differential equations I. x(3ydx + 2xdy) + 8y4(vdr + 3xdy = 0
I Ans:
4
3x
3.
(y2 + 2x2y)dx + (2x 3 -
xy)dy
=
10 +1
2 3+l
4
3
y
3
X 1
---y 10 - - +I 3
+1
7
3
=c]
0 I
( Ans : 4x 2 y
I
2
3
2 __ X
3
2Y 2 =
ex ]
3
4. (2.~y - 3y4)dx + (3x 3 + 2xy3)dy
=
0
(Ans: 5x
- 36, 13
24
y 13_12x
_10 ' I]y
15 I]
=c]
1.17.4 Applications to geometry, law of natural growth and Newton's law of cooling: When some action is applied on a quantity the changes and the action affects all the parts equally. The rate of change depends on the original quantity. For example the total population 'p' of a country increases with respect to time 't', say, then its rate of change with time is
~ . Under ideal condrtions the rate of
change of the population will be proportional to the total population at any given time and is called the law of natural growth. The growth of the population satisfies the differential equation
dp =kp dt
-
where k > 0 is a constant Solving the differential equation we have
p
= cekl
dx The rate at which a quantity x is decreasing is given by dt and this is negative.
52
Engineering Mathematics - I This rate of decrease or decay is found to be proportional to x itself
dx =-kx dt
-
(where k> 0 is a constant)
is the law of decomposition
= ce-kl
Solving the DE x
Newton s law of cooling: Newton's law of cooling states that the rate of decrease of the temperature of a body is proportional to the difference of the temperature 9fthe body and that of its surrounding medium. Let 0 be the temperature of the body at time 't' and 00 the temperature of its surrounding medium. :. Difference between the temperatures = 0 - 8 0 (0 > 00) and the rate of decrease
dO
of the temperature of the body is - dt' since the body is cooling oc. Using Newton's law ofcooling
or
-dO
dt dO
-
dt
=
= k(O -
de - oc(O - 0 ) dt 0
00 ) where k> 0 is a constant of proportionality
-k (0 - 0 ) 0
f~=-k fdt+c 0-0 0
10g(0 - 00) = - kt + c
or
0=0 o + Ce-kl
1.17.4 Example The mass of crystalline deposit increases at a rate which"js proportional to its mass at that time. The deposit has started around a crystal seed of 5 grams. Find an expression of its mass at time 't'. Ifin 30 minutes the mass of the deposit increases by'l' gram, what will be the mass of the deposit after 10 hours.
Solution Let m be the mass of crystalline at time 't' then by law of growth
dm=k' dt '
dm -=dt
km
Ordinary Differential Equations of First Order and First Degree
53
on integration log m = kt + c
.... (I)
initially when t = 0, m = 5 (1) =>
log5 = 0 - c c=-log5
SubstitutilllJ; c value in (1) log m = kt-log5
I.e.,
kt= 10g(
~)
.... (2)
When t = 30 minutes mass deposit increases by , I ' gram m = 5 + 1 = 6 grams Substituting t=
"21 (hours), m =6
Substituting in (2) log(1Y t= 109; Now to find the mass after 10 hours (i.e t = 10) from (3) we get
IOg(H IO~ log 7~ log (7) ~ IOg(%)" x
(
6)20
m = 5"5
grams
.... (3)
54
Engineering Mathematics - I
1.17.5 Example The rate at whi~h a certain substance decomposes in a certain solution at any instant is proportional to the amount of it present in the solution at that instant. Initially, there are 27 grams and three hours later, it is found that 8 grams are lett. How much substance will be left after one more hour.
Solution If m grams is the amount of the substance left in the solution at time 't', then the rate at which it decomposes is dm , which is proportional to m. dl By law of decay
dm
dt
= -
f~
km (k> 0)
=-k fdt+c
logm = -kt+ c
.... (I)
Initially when t = 0, m = 27 From (1) we have log27 = - k.O +c
=>
c= log27
Substitution of'c' in (I) gives log m = - kt + log27
10g(;; )= -
kt
It is given that m = 8 when t = 3 .. From (2) 8 log (2 ) = - k, 3 7
8 -k= log ( 27 2 -k=log3
)X
.... (2)
Ordinary Differential Equations of First Order and First Degree
55
Then (2) becomes
10g~ =IOg(%}
... (3)
when t = 4
III
log 27
2
4
= log ( "3 )
m=27 x (%r grams
m=
16
3
grams.
1.17.6 Example The number x of bacteria in a culture grow at a rate proportional to x. The value ofx was initially 50 and increased to 150 in one hour what will be the value ofx after
1
12" hour. Solution
dx -=/0: dt
dx =kdt x
-
logx = kt + c c is the constant of integration when
t= O,x= 50
..
log 50 = k.O + c
or
c= log50
(1) =>
logx = kt + log50
.... (I)
Engineering Mathematics - I
56
x log- =kl 50 x = 150, when 1= 1 ..
150 log 50
or
k= log3
= k.I
(2) then gives log (
x
5O) = flog3
we want to find x when I
3
=
X
50
"2
3
=
(3)2
3
X
= 50 (3)2 grams
1.17.7 Example The rate of cooling of a body is proportional to the difference between the temperature of the body and the surrounding air. If the air temperature is 20°C and the body cools for 20 minutes from 140°C to 80°C, find when the temperature will be 35°C. Solution If 8 is the temperature of the body at time '1' then from Newton's law of cooling
-d8 -a(8-20) dl = f~ 8-20
~
de
- =-k(8-20) dl
k rldl + c
J'
log(8 -20) = - kt + c Initially when t = 0, 8 = 140 log(8 - 20) = 0 + c,
.... (I)
Ordinary Differential Equations of First Order and First Degree
57
c = log( 120), (I) reduces to log(140 - 20) = - kt + log 120 or
+kt = log( 120) - 10g(S - 20)
.... (2)
It is given that q = 80 when t = 20 minutes k. 20= logI20-log(80-20) I (120) k= 20 log 60 I k= -log2 20
Substituting in (2) 1 (2 0 IOg2)t= logI20-log(S-20) It is required to find t when
In a certain reaction, the rate of conversion of a substance at time "t' is proportional to the quantity of the substance still untransformed at that instant. At the end of one hour 60 grams while at the end offour hours 21 grams remain. How many grams of the first substance was there initially? ( ADs: 85 grams approximately I
2.
The rate of growth of a bacteria is proportional to the number present. If initially there were 100 bacteria and the amount doubles in '1' hour, how many bacteria will
1 be there after 2"2 hours. (ADs: 564
I
58
Engineering Mathematics - I
3. Under certain conditions cane sugar in water is converted into dextrose at a rate which is proportional to the amount unconverted at any time. 1f75 grams was there at time t = 0.0 and 8 grams are converted during the first 30 minutes find the amount I converted in 12 hour. [ADS: 21.5 gms]
4. The rate of cooling of a body is proportional to the difference between the temperature of the body and the surrounding air. If the surrounding air is kept at 30°c and the body cools from 80°c to 60°c in 20 minutes. Find the temperature ofthe body after 40 minutes. [ADS: 48°c J
5. If the air is maintained at 30°c and the temperature of the body cools from 80°c to 60°c in 20 minutes. Find the temperature of the body after 40 minutes. [ ADS: 48°c]
6. The rate at which a heated body cools in air is proportional to the difference between the temperature of the body and that of the surrounding air. A body originally at 80° cools down to 60°c in 20 minutes the temperature of the air being 40°c what will be the temperature of the body after 40 minutes from the original temperature. [ADS: 50°c J
7. The rate at which bacteria multiply is proportional to the instantaneous number present. If the original number doubles in 2 hours in how many hours will it triple. [ ADS: 210g3 J
log2 8. Water at temperature 100°C cools in 10 minutes to 80°C in a room of temperature 25°C. Find the temperature of water after 20 minutes. [ ADS: 65.5°c]
9. A cup of coffee at temperature 100°C is placed in a room whose temperature is 15°C and it cools to 60°C in 5 minutes, find its temperature after a further interval of 5 minutes. [ADS: 38.8°c ]
Ordinary Differential Equations of First Order and First Degree
59
Orthogonal Trajectories 1.18.0 Definition (1) A trajectory of a family of curves is a curve cutting all th.e members of the system according to some law. For example a curve cutting a family of curves at a constant angle is a trajectory.
Definition (2) If a curve cuts every member of a given family of curves at right angles, it is called an orthogonal Trajectory. The orthogonal trajectories of a given family of curves _ themselves form a family of Curves. If the two families of curves are such that each member of either family cuts each member of the other family at right angles then the members of one family are known as the orthogonal trajectories of the other. In two dimensional problems in the flow of heat, the curves along which the heat flow takes place and the isothermal curves or loci of points at the same temperature are orthogonal trajectories. In hydrodynamics, the flow of water from a lake into narrow channel produces a family of streamlines which are orthogonal trajectories to the curves of equal Velocity Potential. In the flow of electricity in thin conducting sheets, the paths along which the current flows are the orthogonal trajectories of the equipotential curves and vice - versa.
be a family of curves, where c is a parameter. We can form a first order differential equation by eliminating 'c' from (1) \
F(X,y,:)=o
i.e.,
..... (2)
is the differential equation whose general solution is (1 ) If the two curves are orthogonal (curves intersecting at right angles) the product of the slopes of the tangents at their point of intersection must be equal to -I. Suppose
(x, y)
trajectory .
is the point of intersection of the curve (I) and its orthogonal
60
Engineering Mathematics - I
At this point slope of the tangent to the curve (I) is dy and -dx is the slope of
dx
dy.
the tangent to the orthogonal trajectory. Therefore on replacing dy by - dx in (2)
dx
dy
, the equation thus obtained is the differential equation of the family of orthogonal trajectories of the family(l) .
. . ¢(
xJI,
-
:)=0
...
(3)
is the differential equation of the system of orthogonal trajectories, and its solution is the fami Iy of orthogonal trajectories of (I) .
is the family of curves where c is the arbitrary constant. We can form a differential equation
F(r,o, :~ )= 0
..... (2)
of the family (I), after elimination of the constant 'c'. Let ¢ be the angle between the radius vector and the tangent at any point
(r, 0) .
on a member of the family of curves. Then
dO tan¢=rdr
.... (3)
Let ¢I be the angle between the radius vector and the tangent at any point (lj, ~ ) on the trajectory. Then
tan;/, Y'I
dBI = r,·1 _ dr,
..... (4)
I
At a point of intersection of the given curve and the orthogonal trajectory
lj=r,
01 =0
Hence tan¢1
From (3) and (4), lj
dOl dr l
1 dr
= --; dO
1 dr i.e., -; dO
dOl
= -lj dlj
= -cot¢
Ordinary Differential Equations of First Order and First Degree
61
Hence the differential equation of the orthogonal trajectory is obtained by
·. dB fi 1 dr. J dB fi dr Sll bstltutmg - r - or - - , I.e., -r- ordr
r dB
dr
Therefore the differential equation F
dB
(r, B, _r2 ~~) = 0
..... (5)
gives the differential equation of orthogonal trajectory of the family of cllrves(l) Solution of (5) is the orthogonal trajectory of the family of curves.
1.18.3 Example Find the orthogonal trajectory of the family of curves a/
= x 3 , where a is variable
parameter.
Solution: Given family of curves is
a/ = x
3
Differentiating (I) w.r.t. 'x'; 2ay dy
dx
..... ( I)
= 3x 2
..... (2)
Eliminating 'a' from (1) and (2)
. (X3) .y dy _ dx -3x
.. 2 /
2x
2
dy =3y dx
is the differential equation of the family (I). Now replacing :
..... (3)
by - : in(3),
gives the differential equation of the orthogonal trajectory to (I) as -2x dx = 3y
dy
Integrating both sides 2X2
23 2 -2~ = L + c . Therefore 2 2
+ 3/ = c is the equation of orthogonal trajectory of (I)
Engineering Mathematics - I
62
1.18.4 Example Find the orthogonal trajectory of the family of parabolas y2 parameter.
= 4ax
where 'a' is the
Solution:
y2 =4ax
..... (1)
differentiating (I) w.r.t. 'x';
y dy dx
Eliminating 'a' from (I) and (2).
.. y=2x
= 2a
..... (2)
y2 = 2X(Y:)
dy dx
..... (3)
is the differential equation of the family (1) Replacing
dy by _ dx in(3) dx dy
Y=2X( -: J
..... (4)
is the differential equation of the orthogonal trajectory to (I). Integrating (4)
fydy = -2 fxdx + c; 2X2 + y2 = c is the orthogonal trajectory of (I) 1.18.5 Example 2
Find the orthogonal trajectories of ~2 + (
a
2
a +A
= 1 whereA is the parameter.
Solution:
X2
y2
-+ =1 a2 a2 +A Differentiating (I) w.r.t. 'x';
..... (1)
2x2 + 22y dy a a +A dx
=0 ..... (2)
Ordinary Differential Equations of First Order and First Degree
63
)
xxy -2- - - = 1 a a 2 -dy dx
Eliminating A from (I) and (2)
:. ( x
2
-
a 2 ) : = xy
..... (3)
is the differential equation of the family (I) ·· -dx - fOor -dy.111 (3) S ub stltutmg
dy
y
dx
dY~-( x' x a
2
;
Y (2 X -a 2)( -d- ) = xy dx
. I.e,
)dx
is the differential equation of the orthogonal trajectory to (1) 2
Integrating we get
x +l 2
= 2a
2
2
Y 2
X
2
-=a logx---
2+c
log x + C is the orthogonal trajectory of the family of curves( 1)
1.18.6 Example Find
the
orthogonal
trajectories
of
the
family
of
coaxial
circles
x + l + 2gx + C = 0 where g is the parameter. 2
Solution: Given
x 2 + y2 + 2gx + C = 0
..... (1 )
Differentiating (I) w.r.t. 'x' ;
dy 2x+2y-+2g+0=O dx
..... (2)
Eliminating 'g' from (I) and (2)
x + y2 - 2x ( x + y : ) + c = 0
2
y2 _ x 2 _ 2xy dy + c = 0 dx
..... (3)
is the differential equation of the family of (1 ) Substituting· - dx for dy in (3) we get,
dy
dx dx y2 _ x 2 + 2xy - + C = 0 dy
..... (4)
---
Engineering Mathematics - I
64
Which is the differential equation of orthogonal traject{)ry. SimplifYing
dx 1 2 c+ / 2x---x = - - dy y y
..... (5)
is a linear equation of the Bernoulli's form (non - linear differential equation of first order and first degree) Substituting x 2 = t,
2x dx = dt in (5) dy dy
dt t ---= dy
y
(c+ / ) y
..... (6)
is linear in t. I· --dy 1 I :.I.F=e Y =-
Y General solution of(6) is
t.! = y
f +yy2 .!dy+k y C
t
C
y
y
-=-y+-+k.
x 2 + Y - ky -
C
=0
( since t = x 2 )
is the equation of the orthogonal trajectory of (I)
1.18.7 Example Find the equation of the system of orthogonal trajectories of the parabolas
Ordinary Differential Equations of First Order and First Degree
1 dr ---tanB/2 = 0 r dB 2
65
..... (2)
£10 elr. for m (2) we get dr dB
Substituting -r -
1(
- -r-?dB) -tanB/ 2 =0 r dr
dr
-
=
r
-cotB/2dB
..... (3)
is the differential equation of the 0I1hogonai trajectory . Integrating (3)
log r
= -2 log sin B/2 + log c
log(~) = logsin c
.
2
/
(
2
B/2
1- cos B)
-=sm B 2=--'--------'r 2
:. 2c = r (1 - cos D) is the equation of orthogonal trajectory of the fami Iy of ( I ) 1.18.8 Example Find the orthogonal trajectory of
rill
= alii
cosmB where 'a' is the parameter.
Solution: Given
rill = a cosmB ln
..... (1)
mlogr = mloga+ log cos mB Differentiating w.r.to 'B' ;
~ dr r
dB
=
m -dr= -
r dB
1
cosmB
( -msmmB . )
-tan(mB)
is the differential equation of the family of curves (I) Replacing
dr
-
dB
2
dB . m(2) d,.
by -r -
..... (2)
66
Engineering Mathematics - I
~.(_r2 da) = -tan (rna) r
dr
dr = cot (rna)da
..... (3)
r
is the differential equation of the orthogonal trajectory. Integrating (3) 'log r = J.-Iogsin (rna) + loge
rn log rln = log e sin (rna) ; rill lll
= e lll sin (rna) is the orthogonal trajectory of the
family (I) Exercise 4 (s) 1.
Find the orthogonal trajectories of the following family of curves, where a
IS
a
variable parameter.
(ii) (iii)
x 2 -xy+ l =a 2
(iv)
2X2 + y2
(v) (vi)
2.
2
y=ax 2 xy=a 2
(i)
[ADS: x +21 =e] 2
?
[ADS: x - y-
= kx
x 2 _ y2 = a2 X2/3 + y2/3 = a2/3
=e ]
[ADs: (x-y)=e(x+y)2] [ ADs: x
2
= - y2 log ( ;
)]
[ADs: xy = e] [ADs: y4/3 _ X4/3 = e 4/3}
Show that the system of confocal and coaxial Parabolas
l
=
4a (x + a)
is self
orthogonal.
3.
3 Find the orthogonal trajectories of the curves 3xy = x 3 - a , a being the parameter.
4.
Prove that orthogonal trajectories of a system of circles x 2 + y2 - ay = 0 is
x 2 + y2 -bx =0
Ordinary Differential Equations of First Order and First Degree 5.
67
Find the orthogonal trajectories of the following family of curves. -
(i)
r=aB
(ii)
r =e aO
(iii)
rn sinnB = a"
2 (iv) reosB = asin B (v)
r = a(l +eosB)
r = aeosB 2 (vii) r2 = a eos2B (vi)
0' --
= ce 2 I (log)2 + B2 = c I
[ADS: r [ ADS:
[ADS: rll eosnB = e"
I
[ ADS: r = e (3 + eos 2B) ] [ADS: r = c(l-eosB) I [ADS: r [ADS: r ?
= esinO I
= C- Sin 20 ) ?
"This page is Intentionally Left Blank"
2 Linear Differential Equations with Constant Coefficients and Laplace Transforms 2.1.1 Let (i) the differential operator" ~" be donoted by 'D'
dx
(ii) PI' P2' P3' ............. Pn be either functions of x or constants (iii) R be a function of x then the general form of a linear differential equation (L.D.E) of order n is given by D nY + PI Dn- I Y + P2 Dn- 2Y + ................. PrJ! -- R
.... (I)
or simply (Dn + ppn-I + ............ +Pn)Y = 0, If PI' P2' ..................., Pn are constants then (I) is called a L D E with constant coefficients. Denoting the differential operator (Dn + PI Dn - I + .............. + Pn) by f (D), (I ) can be written as f(D)y
Suppose p is a constant = -m, say, then solution is y (i.e.) the solution of (D - m)y = 0 is y = Ce rnx
2.1.3
..... (3)
= Cernx .... (4)
In this chapter the attention is mostly confined to L.D.E.S with constant coefficients. If R = 0 then equation (2) becomes f(O)y
=
0
We shall take, here afterwards, f(m)
= 0 as the auxiliary equation, (A.E).
2.1.4 Consider a second order L.O.E (D2 + a\ D + ~)y
A.E. is m 2 + aim + a2
=0
..... (6)
= o.
Let m\, m 2 be the roots of this equation. Now four cases arise. (i) m\, m2 are real and distinct
(ii) m\, m 2 are real and equal (iii) m\, m2 are complex and distinct (iv) m\, m2 are complex and equal
2.1.5 Case i : (6) can be written as [0 2 - (m\ + m 2 ) D + m\m 2 ] y
Call (0 - m l ) y
=
= 0 or (D - m2 ) (0 - m\) y = 0
.... (7)
y
(7) now becomes (0 - m 2) Y = 0
From (4) it follows that Y = Ce D12X From (7) (D - m\)y = Ce D12X I.F.
=
e-rn\X
Solution is y . e-m\x = C fe(m,-m')x dx + C\
y.e-rn\x=
C m 2 -m!
e(rnrrn\)x dx + CI' Call
.... (8)
C m 2 -m!
as C 2 .
-
Linear Differential Equations with Constant Coefficients ...
(i.e.,) y
=
71
C2 e-1Il 2x + C, elll,x where C, and C 2 are arbitary constants
Similarly ifm" m2, .......... I11n are real and dinstinct roots off(l11) = 0 then y = C, e'n,~ -I- C 2 el11 2x -I- ....•....•..•..•.. + C n elllnx
.... (I)
(I) is the solution of f(O)y = 0 where C" C 2, ........... , C n are arbitrary constants.
2.1.6 (Case ii) : m,
=
m2,
From (8) it follows that the solution of(6) is given by
:. y
= (Cx
+ C,) elll,x
Similarly if m, is repeated say 'r' times then solution of(D - m,Y y
=
0 is
y = (C, + C 2x + C 3x2 + .............................. + c;-r-') elll,x
..... (II)
(II) is the solution corresponding to a root m, repeated r times
2.1.7 Case (iii) : Ill" m2 are complex, say (a ± if3) Then solution of(6) is given by
y
=
A eta ± I~)X + B e(a-
(i.e.,) y =
=
IfJJX
where A and B are arbitrary constants
Ae ax . e Vlx + Beax.e-1jJx = eax [A(cosjJx + isinjJx) + B (cosjJx' - isinjJx)]
eax [C, cos/lr + C 2 sin,lk]
..... (III)
where A + B = C, and i(A - B) = C 2 (III) gives the solution corresponding to two complex conjugate roots (a± ifJ)
2.1.8 Case (iv) : If (a ± ifJ) are repeated say's' times then the corresponding solution (follows from case (ii) and (iii) is given by y = ~[(C, + C 2x + C 3x2 + ........................ + CsXS- ' ) ] cosf3x
+ (d, + d 2x + d 3x2 + ....................... + dsxS- ' ) + sinf3x] Example 2.1.9 Solve
d2 y
dy
dx
dx
- 2 +5-+6y=O
Sol. A.E is m2 + 5m + 6 => (m + 2) (m + 3) = 0 :.m=-2,m=-3. Solution is y
Linear Differential Equations with Constant Coefficients ...
Exercise 2 (c) I.
2.
Ans: y
d 2y dy Solve22 +4- +3y=O dx dx
= Acos
J7x
+ Bsin
J7x
x x Ans : y = e-x (Acos.fi + Bsin .fi )
3.
4.
d4 y Solve dx4 -64y = 0 Ans : y
=
C\e2x + C 2e-2x + e-x (C 3cosfix+c 4 sin fix) + ~(C3cos.J3x+C6sinfix)
2.2.1 Consider the O.E.f{O) y
=
R
.... (I)
Operating both sides with f(IO) (called inverse operator)
1 we have f(O) [f(O)]y
=
I f(O) R
I
or y = f(O) R This solution is called the particular Integral of(l) while the solution off(O) y = 0 is called the complementary function (C.F.) Suppose y\ is C.F and Y2 is P.I for (I) Thenf{O) y\ andf{O) Y2
= 0 from (2)
= R from (1)
Hencef(O) [y\ + Y2]
=
f(0) y\ +f(O) Y2
=O+R=R
..... (2)
76
Engineering Mathematics - I
Which shows thatYI + Y2 is also a solution of(l)'Yl + Y2 (i.e.) C.F. + P.I is called the most general solution of(l).
2.2.2 Calculation of
I D-I11,
R
D _ nl, R = y, say, Operating both sides with (0-m 1), we get (0-m 1) Y = R
dy
(i.e.)
dx - m1y= R •
f-mjtU
I.F.lse =e-ml x :. Solution is given by
y.e-mf = fRe-lIl/ dx+c (i.e) y = Ce-ml X + ell/IX fRe-1Il IX dx = C.r. + P.I P.I.=e mI X fH£lIllx dx
m
If
I
=O~ '0
R= fRdx
I Thus the operator 0 stands for integration
I Find the value of - D x
Example 2.2.3
~
Sol. D 5
-5
Fe" fe".xdx ~ <' 0" [ -
+('
m} e'+ xe;" - ~;'l
x I ----5
25
Example 2.2.4
1
Sol.
I eX 02 _50+6 = (D-2){D-3) .(eX)
f
I I (eX) I 3x -3x x dx = (0-2) . (0-3) = (0-2) e e.e
=
I 3x (0-2) . e
2
f
I -I - (eX) = - -I e 2x Je r -2x .exdx = -e I 2 x e-xdx=-e I 2 x.e- t 2 (0-2) 2 2 2
• 2
Linear Differential Equations with Constant Coefficients ...
77
Example 2.2.5
Solve (02+a 2) y
=
tan ax
Sol. A.E is m2 + a2 = 0 => m = ±ai C.F. is y = C I cos ax + C 2sin ax
I
I - I- - -I tan ax = (0+ ai)(O- ai) 2ai [ O-ai D+ ai
P.I =
2ai [e lax J e
r
-a/\
tan axdx - e -IllX J eat" tan axdxJ 2
alX
Je - tanaxdx
=
J. . sinax (l-cos ax) (cos ax -Ismax) - - dx = JSin ax- i dx cos ax cos ax cosax
= - -- -
a
lilly
II
J tanax
=
Je
an
tanax
dx
i isin ax -\og(sec at' + tan ax) + - a a
cosax i isinax - - - + - \og(sec ax + tan ax) - - a a a
I eh PI e2x sin2x = sin2x .. - OZ-40+1 (D+2)2-4(0+2)+1
e 2x 0 -3 Complete solution is
= - sin2x = 2
e 2x sin2x -4-3
1"
= --
7
eX sin2x
Exercise 2(g) I. 3x
Ans : y = (C 1 + C 2x) ~ + 2.
8e
(al -
4x + 3)
Solve (03 - 70 - 6) Y + e2x + x 2e3x 1 2 +-x+-) 5 169 Ans:y= C 1e-x +C 2e- 2x+C 3e 3x_e 2X_(x 12 6 72
3.
Solve (0 2 - 20 + 4)y = eXcos 2x Ans: y =
~(C, cos3Fx + C
2
sin3Fx)+le
x
-~sin 2x
92
Engineering Mathematics - I
4. Ans. y = Clcosx + C2 sinx -~ ~(2cosx - sinx) 5 5.
13 C3sm-x ' 13) - - 1 e 2x.(11cosx-3slnx) . Ans.y= C le- x +e7i (C 2COS-X+ 2 2 130 2.7.1 To find
I~O) [xv] where v is a function ofx
O(xv l) = xOv l + vI 02(xv l) = O(xDvl + vI) = x02v I + 20v I = x0 2v I + (0 2) vI Similarly 03(xv l ) = x03v I + (03) vI in general on((xv l ) = xOnv l + (on) vI (here 1 indicates differentiation w.r.t '0') ..... (1)
Thus, f(O)(xv l ) =o.x(O)vl + fl(O)vl now let f(O)vl = v
1 vI = I(O)'v From(1) j(0)
1
1
[x'/(D) v]=xv + fl(O)[/(O) v]
Operating both sides with I/O) _1_ _ _ 1_
__1_
I
_I_
x /(0) v - 1(0) (xv) + I(O)·f (0) [/(0) .v]
Linear Differential Equations with Constant Coefficients"""
93
2.7.2. Note: For finding [xlll.(sinax or cosax)] when m > 1 the above method becomes very lengthy. Hence we can take it as [Xlll. (R.P or I.P of e/a,"] and apply the method (2.6.1 ). Solved Examples Example 2.7.3 Solve (D 2 + 9) y = xsinx
Sol. A.E. is m2 + 9 = 0 =>m
=
± 3;
1 .) 1 I . -2-(XSIllX = {x--,-.2D}.-,-.sll1x D +9 D-+9 D-+9 [using the result of2.7.1. v = sitlX] P. I.
2.8.1 Cauchy's Homogenious Linear equation Most general form of the homogeneous equation is
(xnon + K1x ll- 1 on-l + ..... + Kn- 1 xO +
=X
.... (I)
Where kl' k2' ........ kn are constants and 'X' is either a constant or a function of x To solve such equations it is convenient to transform them into linear equations with constant coefficients with the substitution x = eZ i.e. z = logx.
dy d\:
= dy
dz dz' dr
=
= ~ dy
.... (2)
x dz
I 'dy 1 d 2 y -l----t----] x x dz X dz2
dy ( .: dz
= x)
__ , d2y dy - 2(d 2 - d )
z
x
..... (3 )
z
..... (4)
d
By taking '0' to stand for dz' dy (2), (3), (4) ~ x dx
=
dy dz
= Oy
..... (5)
Xll
dn dx"
-Z = 0(0-1) (D-2) ......... (O-n+l)y
By substituting (5) in the given equation (') it transforms into a linear differential equation with constant coefficients.
97
Linear Differential Equations with Constant Coefficients ...
Solved Examples 2
2.8.2 Solve x 2
d y
m"2
dy -x d~" + 4y = 0
Sol. Writing x = el.
(i) reduces to
d
[D(D-I) -D + 4tv = 0, where D = -
dz
(D2_2D+4)y=00
A.E. is m2 - 2m + 4
=
:. C.F. is given by y
=
substituting
0 :::) m = 1 ± fi i
(C 1cos,,'3z+C
2 sinfiz)e
I
z = logt,
Y= Clcos (J3logx) + C~sin (J310gx) e logx y = [C ICOS (fi log x) + C 2 sin (fi logx )lx 2.8.3 Solve
d 2y
dv
x2_~ -
3x" - 4y dx
dx-
Sol. Writing x = e Z
2x2
=
..... (i)
(i) reduces to
[D(D-I) + 3D (D 2 - D
+ 3D -
4lY =
2e 2/ ,
4)y = 2e 21.
A.E. ism 2 +2m-4=0 ~m=-1 C.F. is y
=
(C1e
(I.JS)L
2
±J5
+ C 2e
")_
(h'5)L)
2
P.1. = -D-2-+-2-D - 4 (e--) = 4 + 4 _ 4 e Most general solution is
y
= (
C
e-(Il /"S)Z
1
+C
2
e (I JS)7
+ e~z
1
2_ - =
2e 2z e 2z -4- = 2
98
Engineering Mathematics - I
2
2.8.4 Solve.xl Sol.
d --?2y = x d\'
1
2
+-
X
Taking x = eZ equation (i) reduces to [D(0-1)-21v
=
A.E. is m2-m-2
:. C.F. is y
e 2x + e-Z =
(D-2)(D+l)
e2z
1 ~_ ---e-- 3(D-2) =e2z
=
3
I (D-2)(D+l)
e--Z
=
1 (D-2)(2+1)
__ e-
3(D+I)
I (I) e-3 (I) 3(D+2-2) 3(D-I+l)
3
1
+
I
.!. e2z . .!. (I) - .!.
= -
= 2,-1
C l e2x + C2e-z
=
1
P.1. =
0 ~m
D
_
e-z
3
.!. (1 ) D
I
ze2x - - ze-z 3
Complete solution C.F + P.1. I.e.
Y
=
1
C e2z + C e-= + - z( e2x I 2 3
-- e-Z )
I
Y = C1 e210gx 2 + C e-Iogx + - (e2logx
3
Y
= C
I
-
e-1ogX)logx
1 1 1 x 2 + C - + - (x 2 - - )Iogx 2 x 3 x
Sol. Taking x = ez the equation reduces to [D(O-I) (0-2) + 2D(0-1) +21v = 10ez + IOe-z
(D 3-D 2+2)y
=
IO(ez+e-Z )
e2z
+
1 (-1-2)(D+I)
e- Z
Linear Differential Equations with Constant Coefficients ...
A.E. is m3-m 2+2
=
99
0
=> (m+ I) (m 2 -2m+2) = 0 (i.e.) m =-1, m = :. C.F. is y
10 e -L = 5e -L + 2e -L . - -I - - I 5( D + I) (D - I) + 1
=
= C.F. + P.1.
C1e-z + e=(C 2cosz + C sinz) + 5e-z + 2e-z . z logx 3
C1
= -
x
x
•
+ x(C 2cos(logx) + C3sm(logx)) + 5x + 2
d 2y
dy
.
2.8.6 Solve x 2 dx 2 + x dx + Y = logx sm(logx) Sol. Writing z = logx the equation reduces to [D(D-I) + D+ IlY = zsinz (D 2 + l)y = zsinz A.E. is m2 + 1 = 0 => m
C.F. is y
P.I.
= ±i
= C1cosz + C2 sinz 1
= -2-
o +1
•
zsmz 1
ze 'Z
=
I.P. of
=
I.P. of e'z ---:-,----z
-2-
o
+1
=
J.P. of e'z
1 0- -1+20i+l
1 (0+i)2 +1
z
100
Engineering Mathematics - I
=
I.P. of e=
I
Z
D2 +2Di
0)'-1
=
1 ( I P of e1z 1+ . 2Di 2i
=
I.P. of e1z
_I_(z_~) 2Di
2i
_e .i(z2 IZ
=IPof-- - - -Iz .
2
2
= Imaginary part of Z2
= - - cosz
2
Z
J
2i
[2 .(cosz + -I
I.
-- ZSInZ = 4
z. 4
- - (SInZ
. .
[22 + 2.J] Z
ISInZ)
U
+ 2zcosz)
Hence complete solution is
=
C,cos(logx) + C2sin(logx) -
lo;X
[sin(logx) +2(logx) cos(logx)]
Example 2.8.7 Solve(x2 02-xO+4)y = cos(logx) + xsin(logx)
Sol. Taking x = eZ or Z = logx the equation reduces to [0(0-1) -0 + 4lY = cosz + eZ sinz.
A.E. is m2 - 2m + 4 = 0 => m =1±J3i C.F.
=ez (C,cosJ]z+C 2sinJ]z)
P.I. =
I 2
D -2D+4
COSZ
+
2
I . Z
D -2D+4
e smz
101
Linear Differential Equations with Constant Coefficients ...
1 3-20
--.cosz+
=
e
•
(D+I)--2(0+1)+4
smz
(3+20)
eZsinz
3cosz-2sinz
9+4
3
13
---cosz+ - - = The complete solution is y
=
L
1
=
eZsinz
+-3
C.F. + P.1.
~
.
x[C1cos J3 (Iogx) + C2Slll v(logx) ] +
1
.
13 [3cos(logx) -2sll1(1ogx)] +
xsin{logx) 3
Example 2.8.8 Legendre's Linear equation: An equation of the form d ll
dll-1y
(ax+b)11 dx~ + K1(ax+b)n-l dx n - I + ...... + koY
=
X
Where K's are constants and X is either a constant or a function of x is called Legendre's linear equation. To solve such an equation we transform it into linear equation with constant coeficients by the substitution. ax + b = eZ or z = log (ax +b) Now
dy dx
=
= _a_
dy dz dz'dx dy dx
(ax + b)2 -
I
Similarly (ax + b)3
dy ax+b'dz dy dz
= a-
3
d ---f dx
=
= aOy
d
where 0 = dz
a3 D(O-l) (0-2)y
by substituting these values the given equation reduces to linear equation with constant coefficients.
Engineering Mathematics - I
102
Solved Examples Example 2.8.9 Solve
d 2y dy (J+x)2 dx 2 + (I +x) dx + Y
.
= sm2[log(l +x)]
.... (I)
I +x = e => log( I +x) = Z
Sol. Write
dy
(I+x) dx
=
Dy, (I+x)2
d2y -?
dx-
= D(O-l)y where D =
Substituting these values (1) reduces to (D2_D+D+ I)y = sin2z; (i.e) (D2+1)y = sin2z A.E. is m2 + I C.F.
=
0 => m =
±i
= C ICOSZ + C2sinz 1 0 +1
P.1. = -2 - sin2z =
sin 2z
-~4+1
sin 2z
= -~3
· ·IS = C .F + PIC sin32z C omp Iete so IutlOn . = lCOSZ + C' 2Sll1Z - sin 2(1og(l + x))
.
= C1cos[log(l+x)] + C2slI1[log(l+x)]-
"
.J
Example 2.8.10 d 2y
dy
Solve (2x+3)2 dx 2 +6(2x+3) dx + 6y= log(2x+3)
Sol. Writing 2x+3 = eZ or z= log(2x+3) the given equation reduces to [2 2.D(0-1) + 6.2D + 6]y = z i.e.
(4D 2+8D+6)y = z
where D =
d dz
-
A.E. is 4m2 + 8m + 6 = 0 => 2m2 + 4m + 3 = 0 =>m= C.F. is y
~2±J2i
2
=
1
=-l+J2i
e-Z [Acos
~ z + Bsin ~ z]
-
d
dz
Linear Differential Equations with Constant Coefficients
Pol. =
I 1 z = 40- +80+6 6 1
1 6
= -
r1--40 ]z= 3
1 6
1
I
103
000
40
:m 2
40 2D2'Z = -6 [1+- +-3-rl z 3 1+-- + ---3 3 z 6
4 3
2 9
- [z- - ] = - --
Complete solution is y = C.F. + P.1. 1 B' 1 ] + z -2 =e-Z[ Acos r;::z+ SIl1-Z ..;2 ,fi 6 9 I
1
1
e x+..;2..;2
= -2 3 [Acos( r;:: log(2x+3»+ Bsin(
1
')
log(2x+3))] + -6 log (2x+3)- ::.. 9
Exercise 2(h) Solve the following differential equations I. (x 202 - 3xO + 4) Y = 2x2
2.
Ans. (C I + C2logx~ + x 2 (Iogx)2 (x 2D2 - 3xD + 5)y = sin(log(x»
3.
Ans. x 2(C Icoslogx + C2sinlogx) + (coslogx + sinlogx) (x 2D2 - xD +2)y = xlogx
4.
Ans. x(Clcoslogx + C 2sinlogx) + xlogx (x 2D2-3xD+5)y = x 2 sinlogx Ans.
C I coslogx + C2sinlogx + (Iogx? -x[2coslogx - sinlogx] d 2y
-2 +
dx
dy (2x-l) - +2y = 0
dx
C
y= C 1(2x-l)+ ~ 2x+ I
d2 dy (5+2x)2 dx; - 6 (5+2x) dx + 8y = 6x
Ans.
y
=
C (5+2x)2+.fi + C (5+2xi t fi _ 3x _ 45 I
2
2
8
2
11.
(2x+I)2 d y _ 6 (2x+l) dy + 16y = 8 (l+2x)2
dx 2
Ans. 12.
y
dx
= C 1 + C2 log(2x+ 1) + [log(2x+ 1)2](2x+ 1)2
d2y dy (2x+3) - 2 -(2x+3) --12y=6x dx dx
Ans.
C 1(2x+3)2 +
Ci2X+3f~
- 2x +
%
2.9.1 Linear Differential equations of second order d 2y dy The general form is dx 2 + P dx + Qy = R
where P, Q, R are functions of x.
Method of variation of Parameters Let the C.F. of the equation 2
d Y
dx 2
be
+ P dy + Qy = R dx
y=Au+ Bv
where A and B are constants and u, v are two independent solutions-,of d 2y
dy
dx2 + P dx +Qy = R
.... (I)
Linear Differential Equations with Constant Coefficients ...
y
be
Au + Bv
=
105
.... (2)
where A and B are constants and u, v are two independent solutions of d 2y dy - 2 +P-+QllJ= 0 dx
dx
.... (3)
.J'
+ PU I + Qu = 0 and v2 + PV I + Qv = 0
such that
1I2
where
tl2 =
d 2u dx 2
ul
'
=
du dx etc.
In order to obtain the solution of the equation (I) the arbitrary constants A and Bare treated as arbitrary functions of x and are chosen in such a way that
y = A(x)u + B(x)v satisfies (1)
.... (4)
Differentiation of(2) gives = All,
+ lIAI + BVI + vB, (suffixes indicating the order of the derivative) .... (5)
Now we chose A and B such that Alu + BI V = 0 So, Again differentiating d2y dx 2 = (Au 2 + lIIAI) + (Bv2 + vIB,)
.... (6)
substituting (4), (5) and (6) in (I) we get [AlI 2 + BV2 + Alu, + Blv,] + P[Au, + Bvd + Q(Au + Bv] = R
(or)
A(u2 + Pu, + Qlt) + B(v2 + PV I + Qv) + Allt, + BI vI = R
Alu, + Biv i = R (u and v are solutions of(3)
.... (7)
Solving (5) and (7) we get dA
-
dx
=A = I
-vR
uv, --u, v
dB
and -
dx
uR
= B = --I
1/V,
-II, v
integrating we obtain A(x)=
J
-vRdx +CI;B(x)= (uv, -1I,v)
J
uR +C 2 lIV, -u,v
.... (8)
where C I and C 2 are arbitrary constants Substituting (8) in (4) we get the complete general solution of the equation (I)
Engineering Mathematics - I
106
Working Rule: First find two independent solutions 1I and v of(3). Then C.F. is given by y = Au+Bv where A and B are arbitrary constants. Treating A and B as functions of x, we have the solution of (1) as y = AI(x) 1I + BI(x) u where A(x) and B(x) are given by (8)
Solved Examples Example 2.9.2 d2
Solve by the method of variation of parameters ~- y = dx 2
2 --x
1+c
Sol. A.E. is m2 - I = 0 => m = ± I C.F. is y
= A~
+ Bc-x
Assuming A and B as functions of x such that the given equation is satisfied by
y
=
dy dx
Ac-'" + Bc-X we have =
dA
dB
dx-
dx
AcX - Be-x + ~-+e-x-
= A~
dA
- Be--x
.... (I)
dB
Choosing A and B so that ~~ + e--x dx = 0
.... (2)
.... (3)
Substituting (I) to (3) in the given equation, we get dA
~-
dx
dB
+e-x -
2
=--
.... (4)
l+c x
dx
dA
Solving (2) and (4) we get -dx dB
and -
dx
C
-x
=
_c_ eX + I
.... (5)
X
=---
eX + I
= log (~+ I) -e-x-x Integrating (5) and (6) A = _e- x + log( 1+e-') - x + C" B = -Iog( 1+~) + C 2
.... (6)
Linear Differential Equations with Constant Coefficients ...
107
Hence the solution is
y=[-e-x+log(ex + y=C,ex
(or)
Remarks: If(i) I
e-x
+C 2
1)-x+C,1~+[C2-log(1 +~)]e-X
-1 +elllog(e'+I)-x]-e-'\'log(el+I)J
+P+Q=Otheny=~(ii)
I-P+Q=Otheny=e-'x
and (iii) P + Qx = 0 then y = x are solutions of D2y + POy + Qv = 0
Example 2.9.3 Solve by the method of variat ion of parameters d 3y (ix"
- \ + (I 3
Sol:
d
cOlt") ~ dx
cotx = sin 2x
.... (I)
d y dv dx 3 + (I - cotx)
dx - cotx = 0
C.F.P
=
I - cott
.... (2)
Q = -cou comparing with original eqn.
:. I - P + Q = I - I + colx - cotx = 0 showing that v = e-X is a solution or (2) To find another independent solution of(2) let
11
= eO-x . w
By this substitution, the equation (2) reduces to d 2 w dw 2 -x -+-(l-cotx--e )=0 dx 2 dr c- x
or
d2w dw - 2 = (l+cotx)dx dx
~
i(~) dw
= I
+ cou
dX dw or - =eTsinx dx
dw
log - = x + logsinx dx
:. w= JeXsinxdx=- e; (cosx-sinx) eX
.
I
.
u = e-x [--(cosx-smx)] = - - (cosx - Slnx)
2
2
The second independent solution can be taken as cosx - sinx
Engineering Mathematics - I
108
Thus solution of(2) is
y
=
A (COS X
-
silu) + Be-x
For finding the solution of (I) we treat A and B as functions of x
i.e. y = A(x) (cosx - sinx) + B(x).e-x is the solution of (I) where A(x) and B(x) are obtained by solving (cosx - sinx) A I + e-x . BI
0
=
.... (4)
and
.... (5)
dA solving (4) and (5) ~
I
-2 sinx
=
sin2x -dB = ~--ex
and
dx
4
(I-cos2x)
--'-----'4
integrating we get
f'slflxdx+C 1 =--+C cosx 1
A = - -I
2
.... (6)
2
.... (7)
and Hence the complete solution is
y
=
C 1(cosx - sinx) + C2e-x -~ (sin2x - 2cos2x) 10
Example 2.9.4 Solve by the method of variation of parameters (I -
x)
2
d y + x dy _ Y dx 2 dx
= (I
_ x)2
The given equation can be written as d 2y dx 2
x
+ I-x
dy I dx - l-x Y = I-x
x
p= _ . Q= I-x'
I
I-x
.... (I)
andX= I-x
Clearly P + Qx = 0 Hence y
.
.
d2y
= x IS a-solutIon of - 2 + dx
X
-I-
-x
dy I dx - --y = 0 I-x
.... (2)
Linear Differential Equations with Constant Coefficients ...
109
To obtain the second independent solution of(2) take y = vx d 2v
x
2
I-x
x
dv
Then (2) reduces to -'} + dx (-+-.1) dx-
=
0
d . dv dv x 2 - [-]+- [-+-] = 0 dx dx dx I-x x d dv ~~
x
2
I-x-I
2
1
2
- - =- - - - =+( ) -- =+1---d,,1- x x I- x x I- x x dx
lo.g(:) =x+ log(x-I)-logx2 dv
-
dx
x-I
1
(-I)
X
x-
= -.e-t= e-t [-+-, ] 2 x
:. v
1
= e"'" .x
The second independent solution is xv = e"'" solution of equation (2) is Y = Ae"'" + Bx
.... (3)
To find the solution of (I) treat A and B as functions of x such that
dA dB e"'" - + x -
=
0
.... (4)
dA dB e"'" - + x -
=
I -x
.... (5)
dt"
and,
dx
dx
dx
dA dB Solving (4) and (5) we get d; = -xe-X and dx = I
Hence the complete solution is y
=
A e"'" + B.x
= [C1+e-x(I+x)]e"'" + (x+C 2 ) x = C1e"'" + Cr + x 2 + (I+x)
Engineering Mathematics - I
110
Exercise (i) Solve the following by the method of variation of parameters I.
d 2y
-
dx 2
+ a2y
=
sec ax
· x . I Ans. y = C,cos a:r + C'2sm ax + -SIl1 ax + -cos ax Iog(cos a,·) a a
3 Mean Value Theorems and Functions of Several Variables 3.1.0 This chapter deals with (i) Rolle's theorem (ii) Lagrange's mean value theorem also called as first mean value theorem. (iii) Cauchy's Mean Value theorem (iv) Higher Mean Value theorems. (v) Curvature (vi) Centre of curvature. (vii) Evolutes (viii) Envelopes 3.1.1
Rolle's Theorem Ifj{x) is (i) continuous in [a, b], (ii) differentiahle in (a, b) and (iii)j{a) = j{b), then there exists a'c'E(a,b) Proof: Suppose (i)
(ii)
3f'(C)=0
f (x) is a constant function throughout the interval [ a, b ], then f' (x) = 0 VX E ( (I, h) Hence theorem is proved ... f (x) is not a constant function in [ (I ,h ]. As f(x) is continuous in [a,b], there exists a maximum value say, at 'c' sayat'{f (a~d~h) for
f(x)
(a ~ C ~ b)
in(a,b).
and a minimum value,
Engineering Mathematics - I
112
f (c) 7:- f (d) Suppose
f( c) 7:- f( a) then
and Further
and at least one of them is different from
f(x)
and
'c + h' be a point in the neighborhood of 'c' ,.
/(c+II)- /(c) $;
h
/(c+I1)- /(c)
.... ( 1)
.... (2)
is differentiable in (a, h).
/(c+h)- /(c) h-~O h LI
1'( c) $; 0
0 when II> O.
;::: 0 when h < O.
Iz
Lt /(c + h) - /(c) ... / ,(c) == h~O h
i.e.,
f ((l) == f (b)
and
$;
0 and
As h -t 0 we get from (I) and (2) that -
/(c+h)- /(c) h~O h Lt
.
;::: 0 respectively.
1'( c) ~ 0 simultaneously => 1'( c) == 0
Similarly the theorem can be proved when
f (d) 7:- f (a)
3.1.2 Geometrical interpretation of Rolle's theorem Let P and Q be two points on the curve y=f(x). AP==BQ ordinates f(a)=f(b» and the curve is continuous from P to Q. It can be shown that there is at least one point on the curve y = f(x) between x = a and x = b at which the tangent to the curve is parallel to x-axis. . y
x=a
y=a
P
Q
.i(b)
f(a)
x'
0
A
y'
fix) is a constant function
B
x
Mean Value Theorems and Functions of Several Variables
113
x=c y
Q f(a)
o
f(b)
f(d)
x
B
A
y'
j(c) andj(d) are both different fromj(a)
=
j(b)
x=c y
x=a
x=b
f(a)
o
f(c)
f(d)
B
A
x
y'
j(c)
-:F j(a)
andj(d)
=
j(a)
=
j(b)
x=b
x=a y
x=d f(b)
f(a) f(d)
x'
A
0
B
y'
j(d)
* j(a) and j(c) = j(a) = j(b)
x
Eng~neering
114
3.1.3
Mathematics - I
(x) = log { :;a:a:) }in ( a,b )
Verify Rolle's theorem for / Solution: Consider the function
J (x) = 10g{ x(2 + ah)} in the interval ( a, b ). x a+b
J(x+h)- J(x) h~O± h
Lt
2
=
I[ (x+h)2 +llh I x +abj - log - og--+ h~O± h ( a + b) ( x + h) ( a + b) x
Lt
~rlOg (x +_~h) + 2xll+h2 -log x+ hj 2
=
Lt
+h~O±
hl
(x 2 +ah)
x
1
= Lt -1 [ log {'1+ 2xh + h
}
7
r"~O± h
x- + ab
-log {h}] 1+x
Lt [2X
=
+h~O± x 2 + ab
1 ] --+O(h) x
where O(h) indicates terms of order h and higher powers of h.
I
I(
X)
which indicates that
=
2x - ~ x +ab x 2
I (x)
is differentiable in ( a,b ) and hence continuous also.
2
I
Further
(a +ab) (a ) = log ( ) a a+b
= log 1 = 0
2
I (b) = log
(b +ab) (
b a+b
) = 0.
Thus
I (a ) = 0 = I (b)
All the conditions of Rolle's theorem are satisfied. Hence 3c( a < c < b) such that
I
I(
C)
=0 II () C =0
Clearly c
= +ab
. gIves
2
2c
c +ab
1 c
--=O~c
is the G.M of a and b and so
The theorem is thus verified.
E
2
=ab or c=±ab
(a, 1)
Mean Value Theorems and Functions of Several Variables
3 .1.4
115
j"() S1l1X.10 x = ---
. Ven'fy Rolle's t I1eorem '"lor the function
eX
(0 ,Tl )
Solution:
f(x) =~lIl
(O,Tl)
f(x+h)- f(x) h
Lt h-... O±
= Lt -1 [Sin(x+h) -Sinx] h--.O±
h
e( ail)
e'
e< [Sin ( x+ h) - e" sin x] = Ltx iI--.O±
h
ex+1I .e
2
Sin(x+h)-Sinx{l+ h + h + .. ..1 1 I! 2!
=
f
-----------X+-il~x--------~
LI h
e
}HO±
.f!
~
~
1 2 cos ( x + ). sin ( ) - h {sin x + 0 ( h )}
=h--.O± Lt h
----~--~--~-----------ext-II
. (h)
S1l1
=II--.O± Lt !h 2COS(X+!!'). ( h ) -{sinx+O(h)} 2 2
f' ( x ) = cos x ~ S1l1 x e Thus
f (x)
Further
f
is differentiable in
(0) = 0 = f
and hence continuous there.
(Tl) .AIl the conditions of Rolle's theorem are satisfied.
3.1.5 Example Verify Rolle's theorem for the function.f{x) = Ixl in (-1, I). Solution Here
.f{x) = - x for -I < x < 0
= 0 for x = 0
=x for 0
.f{-I)=.f{I)
/'{x}
=-1
/'{x} =
1
for
-1.:sx.:sO
for
:. f'{x} does not exist at x = 0 and hence j(x) is not differentiable in (- 1, I)
.. Rolle's theorem is not applicable to the functionj(x) = Ixl in (-1, I). y y = f(x) = Ixl
------------------~~------------------x
y'
Fig. 3.1 (The curve is not smooth at x = 0).
117
Mean Value Theorems and Functions of Several Variables
Exercise - 3(A) I. Verify Rolle's theorem for the following functions : 1. lex)
= x (x + 3)e-xI2
2. j(x) = sin x
3. j(x)
m
(-3.0)
m
(O.n)
(.'s;+6:
= (x- a)m (x- b)n in (a,
6. fix)
= xl -
3.2.1
(2,3)
m
5. fix)
[Ans : c = n/2]
(:.?~)
= e" (sin x - cos x m
4. j(x) = log
[Ans:c=-2]
b), m > 0,
[Ans : c = n]
[Ans : c =
11 > 0
5x + 7 in (2, 3)
[Ans : c =
J6]
{mb+ na} {m+ 11}
]
[Ans : c = 5/2]
Lagrange's Mean Value Theorem (First mean Value Theorem): Ifj(x) is (i) continuous in [a, b] and (ii) differentiable in (a, b) then ::3 at least one "alue 'c' in (a, b)
3
f'{c)
= f{bi- f{a) -a
Proof: Define a new function «j>x where A is a constant
3
fix) + A.\
..... ( I)
«j>(a) = «j>(b)
i.e.,
fia) + A (a)
i.e.,
A
= _
=
=
fib) + Ab
[f{b)- }'(a)] b-a
..... (2)
j(x) and Ax are continuous in [a, b]
Hence «j>(x) is continuous in [a, b]
..... (3)
j(x) and A x are differential in (a, b),
Hence «j>(x) is differential in (a, b) andj(a)
=
«j>(b)
..... (4) ..... (5)
Engineering Mathematics - I
118
(3), (4), (5) show that ~(x) satisfies all the conditions of Rolle's theorem.
3 atleast one value c in (a, b) ;) ~'
(x)
~' (c)
A
~' =
0
f' (x) + A
=
= f' (c) + A = 0 f'(c)
= -
..... (6)
From (2) and (6) it follows that
f'{c) = j'(b)- f{a) b-a 3.2.2 Geometrical Interpretation of Langrange's Mean Value Theorem P and Q are two points on the continuous curve y = j(x) corresponding to x = a and x = b respectively. .. P[a, j(a)],
Q [b,j(b)] are two points on the curve.
· ... h . PdQ· SI ope on t he Ime JOll1lng t e pomts an IS
f{b)f{a) ' R·IS a pomt . h b_a on t e
curve between P and Q corresponding to x = c, so that f'{c) is the slope of the tangent line at R [c,j(c)].
f'{c) =
f{b)- f{a) b_a
means that the tangent at R is parallel to the chord P Q.
y Q
p
f(a)
x'
0
x=a
f(c)
x=c
f(b)
x=b
x
y'
Fig. 3.2 Hence this theorem tells that there is at least one point R on the curve PQ where the tangent to the curve is parallel to the chord PQ.
119
Mean Value Theorems and Functions of Several Variables
3.2.3 Example Verity Lagrange's theorem for the functionf(x) = (x - I) (x- 2) (x- 3) in (0, 4).
Solution j(x)
=
(x - I)(x - 2)(x - 3) and a
=
0, b
=
4
j(x) = x 3 - 6 Xl + II x - 6 is an algebric polynomial and (0, 4) is a finite interval. Hencej(x) is differentiable in (0, 4) and is continuous in [0,4] showing that the conditions of Lagrange's Mean Value theorem are satisfied.
:3 atleast one value 'c' in (0,4), such that f'{c)
f{b)- f{a) b-a j(0) = - 6, f(4) = 6 =
..... ( I)
f'{x) = 3 x2 - 12 x + II f'{c) = 3c2 - 12 c + 11
3c2 _ 12 c + II
From (I)
=
3c2 - 12 c + 8 = c=
°
6±2fj 3
c = 6 ± 2fj
The point
6 - (- 6) 4-0
3
W h·ICh CIear Iy
°)
I·Ie .111 ( ,4 .
3.2.4 Example Verify Langrange's Mean Value theorem for the functionj(x)
= ~
in (0, I) :
Solution j(x) = ..
~
is differentiable in (0, I) and continuous in [0, I] :3 atleast one value 'c' in (0, 1) such that 3
3
l'{c)
'{ ) f c
=
f'{x)
=
..... ( I)
f(l) - f(O) 1-0
=
j(0) = eO
f{b)- f{a) b-a
=
..... ( 1)
I, f(l) = e
~ gives
f'{c)
=
e
C
Engineering Mathematics - I
120
eC =
From (I)
e-I ' 1-0
-
c=log(e-I), thisc EO, I)
3.2.5 Example Verify Langrange's Mean Value theorem for the functionj(x) in (3,4).
= 5x2 + 7x + 6
Solution j(x) is an algebraic polynomial and the interval (3,4) is finite, j(x) is differentiable in (3,4) and continuous in [3,4]. . . 3 atleast one value c r.(3, 4) such that j'(e) = f(ll) - f(3) 4-3 3
f'{c}
..... (1 )
f{h}- f{a} b-a
=
j(3) = 72,j(4) = 114, f'{x} = lOx + 7 114-72
IOc+7=
c = 3, 5
4-3 (3, 4)
E,
Exercise - 3(8) I. Verify Langrange's Mean Value theorem for the following functions : I. j{x)
=
x(x - I) (x - 2) in (0,
X)
2. j(x)=logxin(1,e) x 2 _ 3x _
I
in ( -~
1, I;)
3. j(x)
=
4. j(x)
= a2 - 7x + lOin (2, 5)
3.3.1
Cauchy's Mean Value Theorem If two functionsj(x) and g(x) are (i) continuous in [a, b] (ii) differential in (a, b) and (iii) g'(x) =F- 0 in (a, b) then 3 atleast one value 'c' in (a, b) 3
f'{c} g'{c}
f{b}- f{a}
= g{b}-g{a}
121
Mean Value Theorems and Functions of Several Variables
Proof: Define a new functionj(x)
= j(x) + A g (x)
.....( I)
Where A is a constant such that 4>(a) + 4>(b)
..... (2)
j(a) + A g (a) = j(h) + A g (h)
[j{b}- j{a}]
..... (3)
A=- [g{b}-g{a}] j{x), A g (x) are continuous in [a, b)
Hence 4>(x) is continuous in [a, b) j(x), A g (x) are differentiable in (a, h)
..... (4)
Hence 4>(x) is differentiable in (a, b)
..... (5)
and 4>(a) = 4>(b)
..... (6)
(4), (5), (6) show that 4>(x) satisfies all the conditions of Rolle's theorem 3 atleast one value c in (a, b) ;) 4>' (c)
3.3.2 Example Verify Cauchy's Mean Value theorem for j(x)
=
1 and g(x) x-
-?
Solution j(x), g(x) are differentiable in (a, b) and continuous in [a, b)
1 in (0, h) x
= -
122
Engineering Mathematics - I
::3 atleast one value 'c' in (a, b) :)
f'{c)
g'{c) Here
f{b)- f{a) g{b)- f{a)
=
f'{x) = ~ X3 -I
g'{x)
=7
-2 C 3
=--IT
=
1
//c 2
a+b ab
2
c c
a
b
2ab a+
= --b which is the Harmonic Mean of , a' and 'b'
c E(a, b)
3.3.3 Example Verify Cauchy's
Mea~
Valve Theorem for fix)
=
e, g(x) =
Solution fix), g(x) are differentiable in (3, 7) and continuous in [3, 7]
::3 atleast one value 'c' in (a, b) :)
f'{c) f{b)- f{a) g'{c) = g{b)-g{a) Here
f'{x)
=
e
g'{x)
=
e-x
c = 5 E (3, 7)
e-X in (3, 7)
Mean Value Theorems and Functions of Several Variables
123
Exercise - 3(C) I. Considering the functionsj{x):::: x 2 , g (x):::: (x) in Cauchy's Mean value theorem for (a, b) prove that 'c' is the arithmetic mean between a and b. 2. Verify Cauchy's Mean Value theorem for j{x):::: sin x, g(x):::: cos x in (a, b) 3. Verify Cauchy's Mean Value theorem for j{x)::::
3.4.1
I
fx, &>(x):::: fx
in (a, b)
Higher Mean Value Theorem with Lagrange's form of remainder (Taylor's theorem with Lagrange's form of remainder): If a functionj{x) is such that
(i) j(x); f'{x}, r{x} ..... fn - I (x) are continuous in [a, a + h] (ii)
fW{x} exists in (a, a + h), then :3 at least one number '8' between '0' and' I' h hn - ' h" 3j{a + h):::: j{a) + ,f'{a} + ..... + - ( \"f(n-I)(a) + -F (ll + 9h) 1. n -I,. n!
Proof: Define a new function. j{x):::: j{x) + (a+h-x) f'{x} + (a+h-xY r(x) + ... J! 2!
+
(a + h - X ),,-1 /"-1 (a + h - x)" (-I) (x) + .A n . n.,
where A is a constant
3
..... (1 )
cp (a) :::: cp (a + h) 2
h j" (a) + 2! h fW (a)+ cp(a) :::: j{a) + 1! hn - I
hn
..... + - ( I",,r-I(a) + n- J!
n!
A
..... (2)
and
cp (a + h) :::: j{a + h)
..... (3)
but
cp (a) :::: cp (a + h)
..... (4)
Thus
j(a + h) :::: cp (a + h) :::: cp (a)
Engineering Mathematics - I
124
Hence using (2) we get
h2 I" (a) + ..... . 2!
h J!
j(a + h) = lea) + - j(a) + -
..... (5) j(x), I'(x) , I" (x ), ... .fn-I(x) and (a + h - x) (a + h - x)2 etc continllolls in [a, a + 17] and differentiable in (a, a + h)
Hence
~(x)
is continuous in
ra, a + 17] and differentiable in (a, a + 17)
..... (6)
(4) and (6) show that ~(x) satisfies all the conditions of Rolle's theorem.
:.
:3 atleast one value 'c' (a < c < 0 + 11) ;)
~' (c) =
0
Write c = a + e 17 where 0 < e < I
:.
:3
e E (0,
I) ;)
~'
(0 +
e 17) = 0
..... (7)
Differentiating (I) with respect to 'x'
.,.
+ [(a+h-x)"J
(n-I)
I" (x)- (n-IXa+h-xt-
2
(n-I)
In-J/(x)j_ n(a+h-xtF;'n!
J
A ..... (8)
From (7) we get ~' (0
+ e 17) =
(h - eh),,-J
(n-I)
[rea + eh]
=0 ..... (9)
From (5) and (9) it follows that 2
j(a + h)
=
hi' (a) + -, h I j(a) + -1' .
2.
"
(a) + ...
(0 < e < 1)
Mean Value Theorems and Functions of Several Variables
The last term
3.5.1
125
~ .!' (a + 0 h) is called Lagrange'sform of remainder. n!
Higher Mean Value Theorem with Cauchy's form of remainder (Taylor's theorem with Cauchy's form of remainder): If a functionf{x) is such that (i)f(x), f'{x) , r{x) ..... fll(X) arc continuous in
[a, a + 11] and (ii),f'1(x) exists in (a, a + h) then 3 atleast one number '0' between '0' and' 1' such that 2
f(a + h)
=
I!
h '" «(I) + ... .f{a) + h f' () a + j!f
Proof: Define a new fUllction ~(x) = .f{x) +
(a+h-x) () (a+h-x)2 () f' X + f" x + .... I! 2! '
(a + h - x)" where A is constant
From (1)
~
(a + h)
I
{n-l}!
+
..... (1)
fll-l(x)+(a+h-x)A
(a + h) = ~ (a)
3 ~
..... (2)
= f(a + 11) h .
~(a) = f(a) + -1'.1
'I
{a} +
h2
f 2!
"
h"-' , {a} + ... + -(-1)'.1 11-1 (a) + II
n- .
A
..... (4)
From (2), (3), (4) f{a + h)
=
h " j
.f{a) + -"
.
,J
n- f"
(0) + -,
2.
(a) + .......
.. ... (5) f(x), f'{x) , r(x) ..... jil- I (X) and (a + h - x), (a + h - x)2 ...... are continuous in [a, a + h] and differentiable in (a, a + 1/)
126
Engineering Mathematics - I Hence (x) is continuolls in [a, a + h]
..... (6)
(x) is differentiable in (a, a + h)
..... (7)
(2), (6) (7) show that (x) satisfies all the conditions of Rolle's theorem. :. 3 atleast one number '0' in between '0' and' I' Differentiating (I)
fll (a + 0 h) is called Cauchy's form of remainder.
Mean Value Theorems and Functions of Several Variables
3.5.2
127
Alternate form of Lagrange's Mean Value theorem If f(x) is i) continuous in the closed interval [a. a + h] and ii) derivable in the open interval ( a . a + h ) then there exists at least one number B, 0
< B<1,
such that f(a+h)=f(a)+f'(a+Bh).
Proof: In ( 2.1 ) put b exists 'c' ,
=a + h
f'(c)
such that
,then from the proof of L . M . V theorem there
= f(a+h)- f(a)
a+h-a writing c = a + Bh we have 0 < B < 1 and (I) becomes
....... ( I)
f'(a + Bh) = f(a+ h)- f( a) => f( a + h) = f(a)+ hf'(a+Bh) a+h-a 3.5.3 Example: If f(x+h)=f(x)+hf'(x+Bh), O
Solution: Let
f (x) = ax 2 + bx + c , then f(x+h)=a(x+h)2 +h(x+h)+c
Alternate form of Cauchy's Mean value theorem If f (x) and g (x) are (i) continuQ,Us in the closed interval [ a, a+h] (ii) derivable in the open interval ( a . a+ h ) and (iii) g' ( x) *- 0 at any point in the open interval
( a, a+ h ) , then there exists at least one number B , 0 < B < 1, such that
f'(a+Bh)
f(a+h)- f(a)
g'(a+Bh) - g(a+h)-g(a) Note: Lagrange's Mean Value theorem can be deduced from C . M . V theorem replacing g (x) with x in the interval ( a, b )
by
Engineering Mathematics - I
128
= b , g ( a ) = a and g X ) = 1 f'(C) f(b)-l(a) -- = . which is L. M. V. theorem. 1 b-a
We get g ( b)
I(
Example: Using Cauchy's Mean Value theorem, show that sin b - sin a < b - a given that 0 < a < h <
Jr 2
Solutio,,: Taking
f (x)
3e
=
sin x and g (x) = x and applying C.M. V theorem we get that,
such that
. Furt her cose < ] 111
sinb-sina
b-a
(0 Jr)2 ;-
:.
= cos
Jr
e where 0 < a < c < b < -
sin b - sin a <
b-a
3.5.5 Example: I f f '( x) is continuous in [ a, b] and
l' .
2
' (I.e.,) Sin b -SIna < (b -a)
fll (x) exists in ( a, b ) then show that I
f(b)-f(a)=(b-a)f '(a)+(b-a)2 fll(a) 2!
where
a
Solution:
g( x) = h( x) _(b - X)2 h( a) where b-a h(x) = f(b)- f(a) -(b -X)f'(X) ........(1) g(a)=h(a)-h(a)=O and g(b)=h(b)-O=O from (I) and L.M.V
Choose a function
Now
g( x)
such that
theorem.
g(a)=O=g(b) . It is given that f'ex) is continuous in [ a, b] and differentiable inC £I, b) .
h(x) and g (x) . :. g (x) satisfies all the conditions of Rolle'stheorem.Hence 3e (a
g
3.6.1 Students are quite familiar with functions of a single independent variable. Functions which depend on more than one independent variable are called Functions of Several Variables. (i) Volume V = nr2h ofa cylinder, r = radius, h = height, curved surface area (ii) Area of a triangle A =
A = 2prh.
I
"2 xy, x = base, y = altitude
(iii) Volume of a rectangular parallelopiped v = xyz, where x = length, y = bredth z = height are all examples offunctions with more than one variable. (i) and (ii) are examples offunctions oftwo variables and (iii) is a function of three variables. Let z =j(x, y) is a function of two independent variables. Here x, yare independent variables and z is the dependent variable and let the function f(x, y) be defined in a region R. Suppose (x, y) be a movingpoint and (a, b) a fixed point in the region R ofthe xyplane. The point (x,y) may approach (a, b) along different paths (see figure) PI' P2' P3 etc.
3.6.2 Concept of a Limit Let
and
(i)
f(x, y) be defined in a region R
(ii)
(x, y) tend to (a, b) along any path
(iii)
E
< 0 be given.
Mean Value Theorems and Functions of Several Variables
If ::3 0 > 0 such that Ij(x, y) -II <
133
V Ix - al < 0, Iy - bl < 0, thenj(x, y) is said to
E,
tend to I.
j(x, y)
Lt
We write
=
'I
x~a
y~b
I(x, y)
Lt
or as
=
I
(x,y)~(a,b)
3.6.3 Example
x2
x2 + y2
Lt
x-~2
Consider
---'--, it can be seen that if 4xy
y->!
+ y2
5 - -8 <
4 xy
E
whenever Ix - 21 < 0.02,
E =
lv -
0.01 (say) then
II < 0.02 thus 0
=
0.02
satisfies )
)
x- +,v-
4xy
5
-"8
5
< 0.0 I. Hence the desired limit is
"8 .
Note: 1 Let then
Lt
«,y)-->(a,b) .
(i)
f(x, y)
=
Lt
[f ± g]
Lt
[fg]
(x,y)->(a,b)
(ii)
I and
Lt
=
=
(x,y)->(a,b)
and
(iii)
Lt (x,y)->(a,h)
g("(, v)
=
(X,Y)-'(a,b)'
m
I±m
/m
1
I
g
III
Note: 2 (x,y) LI ~ (a, b) exists iff
LI
[
LI
(x-~a) y~b
1
I(x,b) ==
Lt
[
Lt
(y~b) x~a
f(x,h)
1
Engineering Mathematics - I
134
3.6.4 Example ~
f(x, y)
=
x- + V . 2 ' find Lt as (x, v) 2x+ y .
~
(2, I)
Solution
LI [Lt
y~ x~2
1 [4
2
y]
x +y LI + 5 2X+;2 . = y~ 4+ y2 =5=1 ..... (i)
Lt [Lt
x~2 y~l
2
x
Lt
+y ]
2x+ y2
=
x~2
[X2 +
I] 5
2x+1 =5=1
The two limits (i) and (ii) are equal. 'J
Lt
Henee
x' + Y ---=-;;-=1
(x,y) ~ (2,1) 2x + y2
3.6.5 Example J
f(x, v) .
=
)
x- - Y' ? J find whether the limit exits as (x, y) x- + y-
~
(0, 0)
Solution
I,I
Thus
X
L/
~ 0 [Y ~ 0
2 2]
x - Y J
')
x- + y-
lIenee the limit does not exist.
3.6.6
Example
Lt Find
xy ~
(x,y)
~
')
(0,0) y- - x-
LI
7:-
y
I[""'
I,/
x] - Y
~ 0 ly ~ 0 x- + y ?
2] )
..... (ii)
135
Mean Value Theorems and Functions of Several Variables
Solution Puty =
11/X
Lt
Lt
xy
(X,y)~(O,O)
/
__ X2
mx
2
X~O m 2 x 2 __ X2
It is clear that limits will be different for different values of m i.e., the limit depends upon the slope of the path along which (x. y) approaches (0,0). Hence the limit does not exist.
Exercise - 3(0) I. Examine whether the following limits exist. Find them if they exist. x 2 ~ y2 + 4
Lt (i)
(x.y)
~
Lt (iii)
(x.y)
3.7.1
yx~2y
J
(x, J~
-~
(0,0) y2
(iv)
X
4
~ 2X2
Y4 (x,y) ~ (0,0 ) x-J ~ y-J Lt
~(2,2) xy~2x
Lt (v)
(ii)
3xy2
(2,1)
x 2 + 4y2
Lt
~
J
x-y-
(x,y) ~ (0,0) x 2 + y2
Concept of Continuity Suppose
Lt
(i)
j(x. y) exists
(~,y) ~(a,b)
and
Lt
(ii)
f(.~, y)
=
j(a, b)
(x,y)~(a,b)
thenj(x. y) is said to be continuous at (a. b).
Note: 1 Ifj(x. y) is said to be continuous at every point of a region R, then it is said to be continuous in R.
Note: 2 Letj(x. y) and g(x, y) be continuous at (a, b) then f ± g,fg, and are all continuous at (~ b).
f
g
(g -:t 0)
Engineering Mathematics - I
136
3.7.2 Example Consider the function/ex, y) Lt
(x,y)
~(I,I)
= x 2 + Y - 2x when (x,
f(x,y) =
Lt [Lt x~
~
1 y
I
(x 2
y) 7= (0, 0) and.l( I, I)
=0
+ y-J -- 2x) ]
=/(I, I) Hence the function is continuolls at (1, 1)
3.7.3 Example 2
2
Consider the function f(x, y) = ~ y J when (x, y) 7= (0,0) and.l(O, 0) = 0 . . x + yLet (x, y)
~
(0, 0) along the path y
=
Lt (.Y,y)
~
mx Lt
(0,0)
2
2
~y 2 (x, y) ~ (0,0) x- + y
f(x,y) =
x 2 m 2 y2
Lt
(x, y) ~ (0,0) x 2 + 1112 x 2
m2x 2 (.r, y) ~ (0,0) I + m Lt
-2 =
°
=.1(0 0) ,
lienee .l(x, y) is continuous at (0, 0)
3.7.4 Example Consider the function.l(x, y) Lt
=
xy x2 _ y2 (x, y 7= 0, 0), .1(0, 0)
(x, y) ~ (0,0) x 2 _ y2 along the path y
•
Lt
xy
mx
=
°
2
(x, y) ~ (0,0) x 2 (I - m 2 )
= niX III
= - - 7=
1-JIl
2
.1(0, 0) except when
Hence the given function is discontinuous at (0, 0)
III
=
°
Mean Value Theorems and Functions of Several Variables
137
3.7.5 Example X
Consider the functionj(x, y) = =
2
Y
2
(x, y) cf. (0, 0)
~X2 _ y2
°
at (x, y)
=
(0, 0)
In this case it is convenient to introduce (polar co-ordinates). Substitutions x = rcosB, y
=
rsinO.
x 2 y2
r4cos20sin2B
l"(2 + y2 - Jr2 (cos
2
0 + sin 20)
r'
-(sin 2 20) 4
4
{x2 + y2
=
y2
4 Now
. I (x-, +y-, )3 2 ~ -E provided /x/ < E· 3
thus when
4
2
/x - 0/ <
E
I· i3
and
lY - 0/ <
I
E3
..... (i)
Hence
x 2 y2
Lt
---;===== =
(x,y)~(0,0)~X2 + y2
°
frol11 (i)
=j(O, 0) 2 2
j(x, y)
=
Jxx Y+ y 2
2-
is continuous at (0, 0)
Engineering Mathematics - I
138
Exercise - 3(E) I. Investigate the continuity of the following: (i)
j(x, y)
=
2x + ),2
(x, y) t= (2, 3)
}
=0
(x, y) = (2, 3)
at (2, 3)
Ans: Not continllous (ii)
xy j(x.y) = - - , ~Xl + y2
(x, y) t= (0, 0) at (0,0)
=0
(x, y)
=
(0, 0)
Ans : Not continuolls. x-y (iii) j(x,y) = - - , 2y+x
(x, y) t= (0, 0)
=0
3.8.1
(x, y) = (0,0)
Ans : Not continuolls
Partial Differentiation Consider z = j(x, y) where x and yare independent variables and z is the dependent variable. Keeping one of the two variables x and y as constant and allowing the other to vary we get a partial derivative of 'z' with respect to the variable that is varied. Keep 'y' constant and allow x to vary, then partial derivative of z w.r. to 'x' is obtained and is denoted by az or f (x, y) or Dj{x, y) or 8f ax x ax az So Similarly
ax az
By
= =
Lt
j(x+Jx,y)- j(x,y)
ax~o
bx
Lt
f(x,y+Jy)- j(.x:,y) by
ay~o
Note: 1
z = j(x, y) represents a surface in the cartesian co-ordinates (x, y, z) system. The section of the surface z = jex, y) with the plane x = k (parallel to yz plane) is a curve. (Similarly the sections with planes parallel to xy - plane and zx - planes also will be curves).
Mean Value Theorems and Functions of-Several Variables
139
Note: 2
az) ( Ox [
gives the slope of the tangent at (x, y, z) to that curve x=k
~z 1 ~
gives the slope of the tangent at (x. y. z) to the curve obtained as
y=k
the section of the surface z =f(x. y) with the plane y
3.8.2 Example z = x3
(i) To find
+ x2y + I
az , we keep 'y' as constant ay
partial derivative w.r. to x of x 3 = 3x2 partial derivative w.r. to x of x 2y partial derivative w.r. to x of I
az
.. (ii) To find
ax
=
=
=
2xy
0
3x2 + 2xy
az , we keep 'x' as constant
ax
partial derivative w.r., to y ofxJ = 0 partial derivative w.r., to y of x 2y
=
x2
partial derivative w.r., to y ofy2 = 2y
az
-
ay
=x2 + 2y
3.8.3 Example f= x 2 + 1- 2xy
af
- =2x-2y ax
af
-
ay
=2y-2x
3.8.4 Example f= (x - y) (x + 2y)
=
k.
Engineering Mathematics - I
140
This is a product, and the usual product rule applies
at' = (x - y) (1 + 0) + (x + 2y)( I - 0) -' ax
= x - y + x + 2y = 2x + y aj
-
ay
=
(x - y) (0 + 2) + (x + 2y) (0 - 1)
=
2x - 2y - x - 2y = x - 4y
3.8.5 Example y Iff= x-2 , Find af and aj x+y ax ay Solution
Applying the Quotient rule
aj
(x+ y)(I-O)-(x-2y)(I+O)
ax
(x + y)2
aj
(x+ y)(-2)-(x-2y)(O+I) (x + y)2
ay
3y (x+ y)2 -3x
= (x + y)2
3.8.6 Example Iff
=
sin(ax + by), sino aj and aj ax ay
aj
-
ax
a
= cos(ax + by) x -
ax
a
ay
=
cos(ax + by)
x
a
ay
(ax + by)
= acos(ax + by)
(ax + by) = bcos(ax + by)
3.8.7 Example Consider f Then
=
ax2 + hxy - by aj
ax
at
= 2ax + hy - ' = 2by + hx ' ay
The expressionix = 2ax + hy is itselfa function ofx andy. We could therefore find its partial derivative with respect to x or y.
Mean Value Theorems and Functions of Several Variables
141
If we differentiate it partially w.r. to x, we get
(i)
2
I]
~ [Of] ax ax
and this is
a a [- , = - (2ax + hy) = 2a ax-
ax
This is called the second partial derivative of 'f \V.r. to 'x'. If we differentiate Of partially w.r. to 'y' we get
ax
!L (a~) ay
ail ayax
written as - -
Thus
al ax
-
= 2ay
+ hy,
-8f ay
=
hx - 2by
Similar steps can be carried out with the expression for
a~ 2
ay
=
2
aI axay
=
~ (8f) ~ ay al ax
=
(8f) ~
-2b
=
h
8fay
ax
and this is -
Engineering Mathematics - I
142
3.8.8 Example 2
2
2
2
al 8f a 1 a f a 1 a 1 Iff= ax3 + hx2y + bl, find - ' - ' - 2 ' - 2 , - - , - ax ay ax ay ayox axay
au ) a2u -=c--at ax 2 = .I(x + Ky) + g(x - Ky) then show that a2u = K2 a2u ay2
ax 2
I
Engineering Mathematics - I
150
V
II. If
1I
=~
z
z + - prove that x
x au + y ay + z
ax . uy
12. I f
II = log(x3 t-
y3
-f :;3 ~3xyz),
au =: 0
oz
prove that
all au
UII
3
ax
az
x+y
-+-t----=:~--
13. I f
1I
=.I(r)
14. [1' 11 =
W Ilere
x
~-
y+z
r
~)I
t- :;
I
I,,)
.jx" + y + z- prove t HIt
=
y
z
z+x
y+z
+ _.- -+ - - then prove that
lADs: -I)
16. If u = xyf (;
l'
prove that
au
all
x-- + y-- = 211 Dx . ~v 17. If x = rcosO, y
18. If 1-1 = logr, r
=
=
rsin0, then prove that
0 20 0 1 0 - , +-) =0 ox- ay 3 x + );J - x 2y - xy2, prove that 4
19. lf 11
=
x 2y + J,2z + z2x , prove that
au+ -au+ -au= (x+ V+z-)~ ax
~
az
.
Mean Value Theorems and Functions of Several Variables
20. If u = .x3 +
151
i - 3ary then prove that a2u
a211
21. If u = e ax + by j(ax - by), prove that
au
all
ax
0;
h-+a-=2abu X
22. If
2
- ?-
a- +u
y-)
z-)
+ - 2- + - 2 - - = I , prove that b +u c +u
-2( allax
+ (au)2 + (au)2 (-au)2 ax 0; az 3.9.1
au au) x-+ y-+zBy
az
Composite Functions If u is a function of two variables x and yare themselves functions of an independent variable I. then
au au ax au av at ax at ay al
-=-.-+-.-'
..... (i)
Suppose now that 'u' is a function of the variables x and y, and x and yare themselves functions of two other variables rand s. then
au au ax au ay ar ax ar ay ar all au ax au 0; -=-.-+-.'as ax as 0; as
-=-.-+-.-
..... (A)
..... (8)
II Ily if u is a function of rand s, where rand s are themselves functions of x and ythen
au au or ax as ax ar ax as ax au au or au as -=-.-+-.0; ar ay as ay -=-.-+-.-
..... (C)
The second and higher-order partial derivatives of'u' can be obtained by repeated application of the above formulae. Also, the formulae can be extended to functions of three and more variabl~s.
Engineering Mathematics - I
152
3.9.2 Example If z = fix, y), x = ell + e- v, y = e-l/ eV then show that
az az az az ---=x--yau av ax ay Solution
az GZaX azOy az u az au ax .au Oy .ax ax . av .
-=--+--=-e --e
az az ax az Oy
az
-=--+--=--e av fu·av Oy·av fu·
v
u
az Oy.
--e
v
By substraction
az - az)_( - - e (au av
v)aZ - - (" e -e v)aZ ax ay az y-az =x--ax ay
3.9.3
11
+e
Example If u
=
j(x, y), x
=
rcos(J, y
=
sinO prove that
Solution we know that
au = au. ax + all .ay = au cosO + au sin 0 ar ax ar Oy ar ax Oy au au ax au Oy au . au - = - . - + - . - = -(-rsll1 0) + -(rcosO) as ax as ay ae ax ay au lau au +2cosOsmO-.-+ . auau ()2 =cos-O.()2 (-Or )2+2 rae ax axOy 1
Mean Value Theorems and Functions of Several Variables
.
~
au 2 .
SIIl-O(- ) +SIIl 2 c3y
au -~ 0(-) ax
153
. au au ax ay
-2cosOsIIlO-.-+ COS 2
3.9.4 Example If u
=
j{x - y, y - z, z - x) prove that
au + au + al~ = 0 ax c3y az Solution Let X = x - y, y = y - z, Z = z -x then u
=
f(x, Y, Z)
all all ax au ay au az au au ax ax ax ay ax az ax ax az au au all au au az ay = ay - ax' az az ay
-=-.-+-.-+-.-=---
lilly
Adding all these gives
au + au + au =0 ax c3y az 3.9.5 Example If
z = j{u,
u x2 -y and v = 2xy az y-=2(x· az x-- y)az )ax ay au
v) Where
=<
prove that
?
(a)
(b) Solution
u = x2 - y, v = 2xy au au =-2y, av ax =2x - 'c3y ax =2y, av c3y =2x
0(-au )2 c3y
Engineering Mathematics - I
154
..... (i)
..... (ii)
From these, we get
at at 2 2 aZ x - - y - = 2(x + y ) ax ay au From (i) and (ii) we obtain
(-J
az -
ax
7
az +( By
az az az az)2 = {[ 4 x-+y- J? + [-y-+x- J? } au av all av
3.9.6 Example If x = reos Bandy = rsill B, prove the following:
ax 1. -=eose· as ' 2.
3.
ar ax
ax . e. -=-rsll1
as
x
~X2 + y2
,
(:r +(: r
=1
ay = sine· By = rsinE> ar 'as
'
ar By
y
~X2 + y2
,
as
-
ax
=-
y x- + y-? ?
as
'
By
x x2 + y2
Mean Value Theorems and Functions of Several Variables
155
Solution I. We have x we get
=
reos O. Differentiating this relation partially w.r.t., rand w.r.t., '0,
ox = cosO· -o=r . or ' ao -rsl110
-
..... (i)
Similarly, the relation y = rsinO yields
oy = sin {} oy = reosO or 'ao 2. From the relations x = reosO, y
.... (ii)
rsinO, we get
=
x 2 +1=,2
..... (iii)
y = tall 0 x
.... (iv)
or
Dr Differentiating (iii) paliially w.r.t., 'x' and w.1:l. ), 'we get -0 ;-0 x y This gives
or y or x x ox r Jx 2 + y" '
y ..... (v)
Differentiating again (iv) partially w.r.t., 'x' and 'y' we get
y( __1 ) = sec" 0 ao , x2
ox
~ = sec" 0 ao ~v
x
which gives
y = cos " 0 = --, y =, -ao = - -,
ox
x·
00
cos" 0
x
x
,.2
r-
-=--=-
oy
y , x- + y
..... (vi)
x
.... (vii)
3. From (v) we get
2 or)2 +(or)" =X +y2 =1 ( ox oy r- rr2
=< =<
(::l\(:r~1
..... (viii)
Engineering Mathematics - I
156
4. From (v) we also get ..... (ix)
..... (x) Adding (ix) and (x)
a2 ,. a2r ax ay2
I r
-+ 2 -=-
..... (xi)
5. From (ix) and (x), we get ..... (xii)
Differentiating
ar = Y
ay
partially w.r.t., 'x' we get
r
So that
..... (xiii)
From (xii) and (xiii) we get that ..... (xiv) 6. From (vi) we get
Similarly from (vii), we get ..... (xvi)
157
Mean Value Theorems and Functions of Several Variables
Adding (xv) and (xvi) we get 2
~
a 0 a-A 2xy 2xy 0 ----:? + - 2 = -4- - - . 1 = ax ~v r r
Note: I. The variables (I', 0) related to (x, y) through the relations x = rcosO, y = rsinA are the polar co-ordinates of the point whose cal1esian co-ordinates are (x. y). 2. In obtaining the resllits in (i) rand 0 are regarded as independent variables and x and yare regarded as dependent variables. 3. In obtaining the results (ii) - (vi), x andy are regarded as independent variable~ and rand 0 are taken as dependent variables. 4. After partial differentiation is carried out, the final expression for the partial derivatives are to be expressed in terms of the independent variables.
.
or ax ax ao ax ar as ax
.
ar is not the ax
5. From the expressIOn for - , - ; - , - , It may be observed that -
. I reclproca
0
fax ax.IS not the reclproca . I 0f or an d ae
ao ox'
ar = ax and ax = 1'2 as . Similarly ar = 0' ay = 1'2 as ax ar ao ax 0' ar ' as 0' 6.
I n Clact we note IJere
.
The results (i - vi) arc useful in calculations involving polar co-ordinates.
Exercise - 3(G) I. If u = x + y, v = xy and f = feu, v), show that
aj Of Of Of x-+ y-=u-+2vax 0' au av 2. If z = 1(r, y) and x = e" + e- X and y = e- II -ev show that
az az az az ---=x--yau av ax 0'
158
Engineering Mathematics - I
3. If z = ttu, v) where
II
'--=
OZ x - --
x2 - y2 and v = 2xy prove that UZ
2
()z
2
v- == 2(x + y ) -
ox . oy
4. If II
=
(!"siny. v
5. If Z
=
/og(u 2 + v).
OU
e'cosy and z =.f(1I, v) prove that
=
U =
eX' +/.
V =
!!..~ _ 2.\"(21/ 2 + I)
ox -
/1
6. If w =.f(u. v) and u = x 2
x2
-f
y, show that 2
OZ _ 4Y1l + 1 • ~v 1/ + v
-
2+V
-I, v = 2.\:v prove that
r 1 l--+-ov-
J o2 W U-W 4(.2 2) ) ) .\ +Y --+--=
olr
7. Ifu
x
~v-
J ;-,J o-W u-W 2)
OU
V and w =zandf=}(u, v, z
= -;V=-'--
z
ll~
show that
Of Of Of Of x-+y-+z-=w-
ox
8. If z = }(x, y), x
9. If z = }(u, v),
= 11
It =
oy
oz
cosh v, y =
Ix +
1I
ow
sinh v prove that
II~Y, l' =- ~v
-
I11X,
show that
a-J + z 0 2z (/ 2 0 Z+02Z - = +111 2)(2 2 ox ~2 0,,2 ov 2 10. If z =f(u, v),
U =
x2
-I,
v
=
1- x 2 show lhat
OZ oz x-+ y-=O ~v
ox
1
Mean Value Theorems and Functions of Several Variables
II. If x
ell /e/1/V.)' C~ ell sec)'
=
159
fino the vallie of
vv + }'-~l ( x~l!-ax + .}' Dy(~)' 1x (x ax'
lADS: 0
~v
(
12. If w
I Iog(y-? -. x-) ) . . II == 2 Ilmt an d v::= Sill =
I(
yx
I
II
y, )~ .. z2, 22 -- x 2) prove that
.f(x2 -
-.!. ow + _~_ ow + _~ all' == 0 x ax )' ~y z Dz 13. If w
=
j(x, y),
11
0=
e\ v
D2 w
eY show that
3 2 \1'
--:=
3x~v
14. If x::= -!vw,y::=
~
---.(111')
ouvv
j;;,z = -J;;;,
Prove that
o
ax . ~Y
oZ
all
ov
all'
where
3.10.1 Homogeneous Function
Definition A functionj(x, y,
Z, II .. .)
written in the form xl1
Ex.
I. z = sin-I
is said to be homogeneous in x, y,
g(Yx '!-,!~ ..... ) and is said to be of degree x X
(~ ) + tan 1(~)
Z=xO[sin-I(~J+tan 2•
Z
= x
2
+
y
x-y
Z, 11 .....
2
. IS
if it can be
11.
is homogeneolls and is of degree '0'. Since
I:]
IlOmogeneous an d'IS 0 f'd egree Isll1ce '
Engineering Mathematics - I .
160
3.
~ = x 3sin(~) + A)!ZIOg(;1+ XZ2 tan I(~ 1
~ x'ISin(~)<~lot H~r tan'(; llence
~
lJ
is of degree' 3' and homogeneous in x, y, z.
3.10.2 Euler's Theorem The following theorem, known as the Euler's theorem gives a very lIseful formula for a specific combination of partial derivatives ofhomogeneolls functions.
Theorem :
ffz is a homogeneous function ofx andy with degree '11' then
oz
oz
x-+ y-=nz ox
.
..... (i)
~y
Proof : Since Z is a homogeneous function of degree '11' we can write 'z' in the form.
z=
Now
xllf(~)
~: = Il - ~) + x l /(:)( ~;) I1X
=
1f
(
nxlllf(~)_xll 2y/(~)
and
xoxoz + yoz = f( xy) -x" 'y/ (y) + x" 'y/ (:Ix) x x ~y
I1X"
= nxllf(:Ix) =nz
[.:Z=Xllf(;)]
Mean Value Theorems and Functions of Several Variables
uz
oz
oX .
~V
x-+ V-=IIZ
Thus
·3.10.3 Ir z is a homogeneous functions of x and y with degree 'n' then 2 ]OZ
J
o-z uxoy
)
X - , + 2xy - - + y
ax -
20-Z
-i = uy
(
11 II -
)
1z
..... (ii)
Pr(}(l: Diflerentiating (i) pUl1iaily w.r., to x, we obtain
02 Z OZ 02 Z OZ x--+-+ y--=ll-2 ox ox oxoy ax 02 Z 02 z UZ x-) + yx--=(n-I)-
I.e.,
ox'
oxoy
ax
.... (iii)
Differentiating (i) par1ially w.r., to 'y', we obtain
02 Z 02 Z oz y-+x--=(n-1)oyJ oyox oy
..... (iv)
Multiplying (iii) by 'x' and (iv) by 'y' and adding we get
02 Z 02 Z ) 02 Z (OZ Oz) x--? +2xy--+ y-? =(n-I) x-+y- =n(n-l)z ox' oxuy oy' ox ry ?
3.10.3 General Formula The Euler's theorem may be generalized to homogeneous functions of more. than two variables. For a homogeneous function j(x, y, z) of degree 'n' in three variables, the theorem can be read as :
x oj + y OJ + Z oj ax ry oz
= nj
..... (v)
In general, for a homogeneous functionj(xl ....x,,) of degree 'n' in 'k' independent variables x\ .... x k the theorem reads
xl -
OJ
aXl
+x2 -
OJ
aX 2
+ ..... +xk -
oj
ax k =l1j
The proofs of (v) (vi) are left to the student
..... (vi)
Engineering Mathematics - I
162
3.10.4 Example
.
1/ 1I :~. cosec
-[;-;+JY] au __ ~ if;: + if;; ,thcn prove that x au ax + y ay - 6 tanu I
Solution We note that from the given cxpression for' u' we can write cosecll
;-;+JY . ;1.: f (say)
= "il
'Y.\"
+vY
=:
..... (i)
Then
:. fis a homogeneous function of degree
~,
Hence by Euler's theorem
W WI. I x-+ y-=- j =-cosecu ax
I.e.,
a
Oy
6
6
a
I x-(cosecu) + y-(cosecl~ =: -cosecli ax ay 6
..... (ii)
163
Mean Value Theorems and Functions of Several Variables
Carrying Ollt partial differentiation, we get
I
--cosecu 1 au a 6 x-+ y-=-=----=--tanu II
~v
ax
- cosec II. cot u
6
au + yau- = - -1 tanlt
x-
ry
ax
6
3.10.5 Example
x y z Ifu= j o( -,-,yzx
J prove that x-+y-+z-=O au au au ry
ax
az
Solution 11=
o(x
z)
y j -,-,y z x
all au au 0'( -----+---+x z x y y z x-+y-+z-=j ax ry az y x y z z x =0 3.10.6 Example If u = sin
1; + tan -1 ~
au
Solution
u = f(x. y) is of degree '0' By Euler's theorem
all
prove that x ax + Y ry = 0
J
Engineering Mathematics - I
164
au
ali
ax
0'
x-+ y-=o.u=o
.
,x - = v, y
Sill
Alter:
,y
tan - -
x
=w
u=v+w
au avow -::::-+au av Ow ax ax ax' ay 0' 0' x . 1 av -:::: Sill V, - = cosvy y ax y - y ? aw - = tan w, -} = sec- w.x xax -x av I ) aw --:::: cosv- - = sec"w.-
-=-+-
0"
y2
x
0'
=0
3.10.7 Example If
ZI =
X3 + 1'3)
tan' - __ ( x+y 0
(i)
prove that
au au . ax oy 2 a-u? au. . 2 a-u X - 2 +2xv--+ y-} =sm4u-sm2u ax . axay ay-
x - + y - = s1112u J
(ii)
-
?
Solution
Proof: For the given expression for 'u' we have
x3 + y3
tanu = ----'-- = f, say x+y
Mean Value Theorems and Functions of Several Variables
Ilere
~f'
165
is a homogeneous functions of degree'2'
By Euler's thcorem
at" of
.\"-' Ox +Oy I.e.,
.1' 7f' = 1l.J =~ --:
X~(lanl/) I~ y~~~(tanll)= 2tanll ax
~v
Carrying oul partial derivatives, we get
:"\ vII
? xsec-II-;-~
vx
all
?
+ ysecII
au
all
aX
(~y
x - + y~~~
2 tan 11
= --,,-
sec ~
- 2tanll
~v
.
ceo
S1l1211
..... (i)
u
Next dirterentialing bL)th sides of the resull (i) partially \Y.r., to 'x' and 'J" wc get alII
all
ax-
D.\·
Ou
alII
x ~-, + ~- + v .~-- =-- 2 cos 211 -
..... (ii)
a"11 D) II UII all x -~ + )' - , + -~ = 2 cos 211 ~~
..... (i i i)
(~VaX
. axD)'
~v'
Dx
~v
(Iy
Multiplying the result (ii) by x, the result (iii), by 'y' and adding the results we get
,a"II +
a~1I + y' x', ~-,(1.\'-
~-:,~
ay
'1
= (2cos211 =
a"II
~xy-- =
OXG)'
(2cos211 ~-l (all all) x - + y--Dx
-I ) sin211,
sin411 - sin211
3.10.8 Example
Solution Let
and
(from i)
~v
Engineering Mathematics - I
166
Then v and ware hOl1logeneous./llllclioll.\' (!ldegree '2 '. Applying 3.10.3 to 'v' and 'w' we get ..... (i)
!Illy
..... (ii)
(i) - (ii) gives
Note: Letj(z) be a homogeneous function ofx and y of degree Il ; then from Euler's theorem
a" =I/f + y_:'l ax Dy
('JI"
x_:'l
i.e.,
.' f)z , az x..f . ax + xi ay =
.
oz
OZ
Ilf
1?/ => X ax + y ay = 7
DifTerentiating (i) w.r., to y'partially
j
--JJ a'z aZ }I--=Il a'z [j") 'fj'" aZ \'-+-+ ' 2 . ax Dx . oxoy /2 -ax J
J
Differentiating 0) w.r., to 'y' partially 2
2
[/2/2-jt' 1
a z a z az x--+y-+-=n . axay ay2 ry
-az ay
..... (i)
..... (ii)
..... (iii)
(ii) x + (iii) Y gives
=
Il[j '2
~ f( ][x.az + y_~~]
j2
OX~}I
..... (iv)
Mean Value Theorems and Functions of Several Variables
167
Using (i) in (iv)
..... (v)
llere
3.10.9 Example If z = log(,3 + .vI __ x 2y __xy2), prove that (i)
a=
oz
x-+ v- =3
ax . ry
Note: (i) and (iv) can be llsed as formula in such cases
"a z a z )a z then show that x- - 2 + 2xy-- + y - " + 1 == () ax
j then show that
x5 _ y5
')
2
a-z
2
a z )a z 1 2 X - - 2 +2xy--+ y - , =--tanz(tan z-19) ax ax~ ay 400 2
~ x 2+ y2] 15. If z = log (
x+y
;\2;\ 2 ; \2UZ 2 2uZ uZ
find the value of x
GX2
+ 2xy axay + y ay2
3.11.1 Jacobian Ifu and v are functions of two independent variables x and y then the determinant
au au ax Ov ax
~ i.e.,
Ov
ux
uy
v,
vy is called the Jacobian of u, v with respect to x and y
au
and is denoted by
a(u, v) a(x,y) If u,
l~
or
.J(~) x,y
ware functions of three independent variables x, y, z then the determinant
Mean Value Theorems and Functions of Several Variables
lI(
UI'
1I,
V(
vI'
Vz
W(
Wy
lV,
is called the Jacobian of 1I.
V,
w with respect to x.
a{/I, v, w) a{x,y,=)
)~
173
z and is denoted by
)(U,\" w)
or
x,y,z
Similarly aUI
uU I
a{U p 1l2·····uJ _ aX I aX 2 a{X I'X 2.....x,,) - all" OU" oX I aX2
(J/I~
ax" au" ax"
SOME PROPERTIES OF JACOBIAN 3.11.2 Property
a(u,v) a{x'Y)_1 a(x,y)' a{u, v) where
J = a{lI,
vI and
a{x,y)
_ a{x,y) .J - 0(11, v}
Proof: Let u = u(x. y) and v = lI(x. y) so that u and v are functions of x and y. Then differentiating these partially each with respect to II and v, we get I :=
au . ax + all . 0' ax all
ay au
Engineering Mathematics - I
174
av ax av ay ax au ay au
O= - . - + - . -
1=
Now
av .ax + av .ry ax av ay av au ax av ax
a(lI, v) a(x,y) a(x,y)' a(u, v)
all oy av ay
au -au ax ay av av ax ay
ax ax all av ry ay au av
x
-
ax au ax av
-
x
ry
all ry
av
(Interchanging rows and columns in 2 nd determinant)
ali ax au oy ax au oy au av ax av ry -.-+-.-au au ry au
au ax au ay ax' av ay' au av ax av ay --+--' ax' av ~v' av
-.-+-.- --+--
=I~ ~1=1
Proved
3.11.3 Property
a{u,v) _ a{u,v) a{r,s) a{x,y) - a{r,s)' a(x,y) where u, v are functions of r, sand r, s are functions of x, y.
Proof By differentiation offunction ofa functions, partially; we get
au au ar au as ax ar ax as ax
-=-.-+-.-;
Mean Value Theorems and Functions of Several Vanables
175
au -+--'---_. or oy OS
ou
-=-
or .~ as· oy ,
~
..... (i)
oval' - + avos -_. ox or .ox os· ox ' Ov 01' or 01' as -=- -+--_. ~y or .oy os· ~y , 01'
-=~-
and
all au
o(u, v) a{r,s) (1(r,s)· a(x,y)
(lr
=?----=
uv i"lr
os elv us
x
(lr (Ix us
(Is (Iy us
ilyl
(!y.
or as or os x ax ax ov -oV (II' as or as oy oy
Oll
all
(Interchanging rows and columns of the second determinant)
or ou os or ox os ox Ov or OV os -.-+-.-or ox os ox
or Oll as or oy os oy ov or ov as -.-+-.or oy as ~
all
ou
-.-+-.- -.-+-.-
ou -ou ox ~ _ o(lI,v) ov Ov - o(x,y) ox ~
-
using (i)
-
3.11.4 Property If u, v are functions of two independent variables x and y thell u, v, are independent if
.
Otherwise they are dependent. oo((U,V))oFO. x,y
.
.
. o(u, v, w)
(functions of x, v, z) are II1dependent If and only If
(
a x,y,Z
)
oF
Similarly if
0
II,
v,
W
Englneenng Mathematics - I
176
Theorem If the functions 11 1, 1I::...... 1I11 of the variables xI' x 2 ....• xn be defined by the relations Il, ==.ll'",), 112 =flx ,. x 2), 113 = .li x ,. x 2• x 3 ) ..... 1/11 .I;,
O(lI p 1/2 .....U,,}
Then
OU, UU] OU 3
o(x"x2.....x,,) -
ox, . oX1 'a~~
OU II OXII
We know that Similarly lor
Oil,
Oil 2
ox,
ox,
Oil,
0111
0(U"l/2 .. ·..U,J _ aX 2 o(x, 'X 2 ..... X II )
Then
oX1
all,
~1I2
aU
OXII
OXII
OXII
OU,
Oll]
OU"
ox ,
ox ,
ox ,
Olll.
Ollil oX1
0
oX 2
o
OU, Olll ox , . oX 1
o
3.11.5 Example
If x
= 112
-vl . Y
_ o(x,y)
=
2uv ,. tllld --( --) 0 1l,V
Solution We have
x
=
u2 - v2
~
-ox
all
=
2lI.
ox
--=-2v OV
Ii
OU II OXII
177
Mean Value Theorems and Functions of Several Variables
cYv
0)1
v = 2uv => --- = 2 v -"-- = 211 V I I ' 01'
-
ox AX V{X,y) _ or ()I~ ~ 1211 - 21'1 O(II.~) - VI' ~F 21' 211 au VV
:=
4(u-» +) 1'-
)
3.11.6 Example
V(x, y) V r,O
a( t,
°)
If x ~ reose. y ~ rsintJ find - (- ) lllld-(--)
,
a x,y
Solution
x = reosq, (It
=>
~ =
y
rsinO
=
cos()o
(1,.
-ax = -rsiIlO,
ov = rc()sf)
---
0(1
of)
ax ax
Then
o(x,y) _ or
ao
12 = x 2
- sinO
ao
or Again
-rsinol . r Sill 0
_Icoso
~v
0(1',0) - ~v
+ y2 ; 0 =
tal1- 1 ~
x or
2r-- = 2x
ox
-- y
(l()
ax
0
r+y or
-
ax Dr ay
-y 0 r-
0
y
,.
+x
= -
I'
(1()
ilx
x2
-y
-y
+)'2
r
)
=
r(cos 2 0 + sin 20)
= r
Engineering Mathematics - I
178
x
x
lar
ar
x
a(,., 0) ~ ax ~}! aCy,y) ~ ao ao --
Y y
~
vy
lay
3.11.7 Example If x
=
uv, y
=
11 -
then show that .1./ ::.: I
v
Solution ax
We have
-=V,
all
al'
~F
ay
v
VII
But
')
lr =
aI'
=
11
--)
v·
VU )I, ax .
DII 2u-=x
2u- =
xy
,
and
at
-=11
v·
x
=-, 2v
J =
y
av =!...
aI'
ax y'
21'-= ay
y
x
211
II"
lly
2u
Vr
v)'
I
~x
2vy
21:J!
X
X
4uvy
411vy
ay
X
2uv· Y
I? =---v211v
V
2u
x --?
y.
Mean Value Theorems and Functions of Several Variables
179
- 3.11.8 Example Prove that JJ
1 for x
=
e V secu, y
=
e V tanu
=
Solution
ax = e secu tanu, all and
oy
au
=
=
- = ('
"
I
I'
)
e sec1I =
ye-
.
tanu
v
e secutanu
2v =
e secu
=
e" tan 11
II
1
Y
eVtanu
--e" secll
x
y
SlIlli = X
U
(Y) an d v '21 Iog(x
·-1 = Sill -;
Du
y
X~X2
i1x
av ax
=
x- -
X~X2 -y 2
(Jy
- /
Dv y
? '
- Y-))
all
,
x 1
2
y ?
?
x- - y-
i1y
Y J
lIx
lIy
Vt
VI'
X~X2
_ y2
X
x2 _ y2
~X2
_ y2
Y
x-1 -
1
y-
2"
-e secu = -xe
1
y, also
e2v = x 2 -
v
av
sec 2u - lan 2
x 2 e- 2v -
e V secli
0;
eV sec 2 u;
x" Yv Now
ax av
-=
I'
Engineering Mathematics - I
180
, \' I JJ =xe ,-=1 \'
xe
3.11.9 Example VIV Wlt lIV I I I O(X, y, z) If x=-,y=--,z=-- t Jen SlOW tlat - ( - '-)=4
v
II
W
01l,V,\I'
Solution
ax Oll
o(x,y,z) _ oy 0(1/, V,-H')
Oil OZ 011
ox ov oy ov oz OV
oX
VI\'
all'
2
u
lI'
v
It
1/
11'
Wit
II
V
2
oz
V V
v
ow
IV
IV
oy ow
- VlI'
Wli
VlI'
- Wli
VlI'
lIW
II
uv
ltV ltV '-2 '-2 ' 2 II v W --l/V
-I -I -I =-1(1-1)-1(-1-1)+ 1(1+1)=4
Exercise - 3(1) I. If y + x + Z = 1I, Z + Y
=
1Il: Z = lIVW, find o(x,y,z) o(u, v, w)
[ADS:
(1I
2
v) J
lSI
Mean Value Theorems and Functions of Several Variables
2.
If
II
2/ - Z2, V = x2yz, W = 2X2 -xy find
= x+
o(u, v, w)
at the point (1,-1,0)
o(x,y,z)
Ans:19
3.
If x==rcosB,}J==rsinB, find .
o(x,y)
o(r,O)
o(r,O)
o(x,y)
--'--.
~(x,Y)
and
!JS!~~~ 8(x,y)
and show that
= 1 (Hint :relcr 3.6.6) X X,
XOXI
8(r,0)
XIX)
o(YI,VnY,)
4.
I ,y, = - - - ,show that ----'----'--== 4 If YI = - --- ,Y) = --XI x2 x, 0(X I ,X2,X3)
5.
If
6.
If x=a(coshu)cosv,y=a(sinhll)sinv, then show that
YI ==1-XI'Y2
o(x, v)
a
=xl (l-x2 ),y] ==xlxl(l-x,),provcthat
O(YI' Y), y, ) 0(X I ,X2 ,X,)
= (1)3 - XI 2 x"
2
--,-'- = - (cosh 211 - cosh 2v) 2
0(11, v)
7.
If II = /(x l ), v = ¢(x l , x 2 ), W = If/(x!, X 2' Xl) then show that
Oell, v, w) 0(X I ,X2 ,X3 ) 8.
If
X
. 1')
= r S1l1 r.
Solution:
011 Ov
ox
do • 0 . do 0 I oC X, y, z) ). 0 cos 'I', Y = r Slll S1l1 or, Z = r cos - S lOW that - - - = ,.- S1I1
o(x,y, z) o(r,O,¢)
=
uw
=--' oX • ox} 2 i
o(r,a,¢)
=
ox
ox
uy
or
00
o¢
ay
oy
uy
--
or
UO tJ¢
OZ
OZ oz --
or
oe
sinOcos¢ reosOcos¢
= sin Osin ¢ cosO
r cos Osin ¢ -rsine
- r sin Osin ¢ rsinO eos¢
o
o¢
sin ecos¢[ 0 +,.2 sin 2 Oeos¢ ] + r coseeos¢.r sineeosecos¢ +
+( -r sin esin ¢)[ -r sin 2 Osin¢ - r cos 2 esin ¢ ] == r2 sin 3 Oeos 2 ¢ + r2 sin Ocos] Oeos 2 ¢ +,.2 sin 3 Bsin 2 ¢ + r2 cos 2 Osinesin 2 ¢ == r2 sin 3 0 + r1 sin () cos 2 () = ,.2 sin O.
Engineering Mathematics - I
182
3.12.0 JACOBIAN OF COMPOSITE FUNCTIONS: Let (u,v) be functions of x,y where x,y themselves are functions of r,e .Then 1I,V arc
. fi .. . e 8(1l, v) composite unctlOlIls of r, and
= 8(u, v) . 8(x,y)
.
(from 3.6.3)
8(r,B) 8(x,y) 8(r,B) Note 1: Ifu,v are functions ofx,y we may regard x,y as functions ofu,v. Taking r = ll, = v in the above result, we have
e
au
8u
8u 8v
8v
1 0
av
0
811
8v
-
8(1l, v) 8(x,y)
---
8(x,y) . 8(1l, v)
..
=
8(u,v) 8(1I, v)
8(1l, v) 8(x,y)
=
8(x,y) . 8(1l, v)
=1
1
1 => .1.1 1 = 1
I where J = 8(1l,:i,J = 8(x,y)
8(x,y)
Note 2:
8(u, v, w) 8(r,e,¢)
are functions of
8(1l, v)
=
8(u, v, w) 8(x,y,z)
.
8(x,y,z) 8(r,e,¢)
where u,V,W are functions of x,y,z and x,y,z
r, B, ¢ .
Also
J./' = 1; where J= 8(1l, v, w) , .II 8(x,y,z)
3.12.1 Example: If
II =
= 8(x,y,z)
8(u,v,w)
xyz, V = x 2 + y2 + Z2, W = x+ y+ z find .I = 8(x,y,z) . 8(1l, v, w)
Solution: we first evaluate .I
I
=
8(1l, V, w) 8(x,y,z)
,
since u,v,w
yz zx xy functions of x,y,z.
.II
= 2x 1
2y 2x
1
Hence from the relation that .1.1
1
= -2(x -
y)(y - z)(z - x)
1
= 1,
we have
J = 8(x,y,z) = -] 8(u, v, w) 2(x- y)(y-z)(z-x)
are explicitly given as
183
Mean Value Theorems and Functions of Several Variables
Exercise - 3(G) I.
)
1
C(ll, v)
•
o(r,O)
2.
If x =
;;;;,y = Jwu,z =-r;;;;
where II = r sin OcosrjJ, v =,. sin Osin rjJ, W = rcosO, show that
3.
3
Ifu=2xy,v=x--y, whcrex=rcosO,y=rsmO show that ---=-4,..
·
If x =
ltV )' ,
v(x,y,z) a(r,O,rjJ)
ll+V
_
= -1I-V'
1 4
= --,.
2
•
Slll
0
v(u,\')
(1I-V)2) (A liS : --'------'-
fll1d - - o(x,y)
411v
·
yz zx xy 8(x,y,z) V = HI = prove that = x ' y , z 8( ll, V, lV) 4
·
x +Y
4.
II II = -
5.
1I1I=--,O=tan-l(x)-tan-- (y) show that =0. 1- .xy v( x, y)
6.
If lI=x--2y,v=x+y+z,w=x-2y+3z show that ----=1Ox+4
•
D( ll, 0)
I
v(u,V,lI') 8(x,y,z)
1
3.12.2 Jacobian of implicit functions Let
1I1' 1I2' 1I3,
be implicit functions of x,,
jJu" Then
then
1I2' 113'
x 2' x3
x" x2, x3 ) = 0,
so that r= 1,2,3
oj; + a.t; .aU I + a.t; .aU 2 + a.t; .aU 3 aX I aUI oX I oU 2 aX I aU 3 oX I
= 0 and so on
a(.t; ,j~ '/1) 0(11" u), uJ -a(u"U 2 ,1I))' a(X I ,X2 ,X3 ) I
aj~
all,. =
ax!
I af2 . OUr Ollr
I
ax!
a.t:l }u au,.
(by the rule of prodllct of determinants)
oj~ . OUr I a.f! . Ollr au, oX2 all, oX3 I of2 . all,. I a.f~ . all r au,. aX 3 OUr aX 2 I af3 . au,. I 'fJj~ . OUr all, aX3 au,. oX2
x y z If u = r:--:;' v = r:--:;' W = r:--:; prove that v'1-r2 v'1-r2 v'1-r2
a(u, v, w) _ 1 a(x,y,z) - (l-r2 where
-? =
x2 +
(1-/ 2 )
=x2
r+
y;
z2
Solution we have
II
W, x, y,
=>
II(u, v,
Similarly
./; == v2 (1 - x 2 ~
13
z)
== w2 (1 - x 2 ~
= u2 (I _x2 -
r - z2) - r
r-
r
_z2) -
=
0
z2) - z2 =
0
x2 = 0
.... (i)
(see 3.13.2)
Mean Value Theorems and Functions of Several Variables
aCt; ,.I~ ,.I;) _
J
- 2u 2 x- 2x
- 211-Y
- 2112 z
-2v 2 x
-2v 2 y-2y
- 2v 2 z
8(x,y,z)
-2w 2 x
- 211'- y
)
(-2x) (-2y)(-2z)
=
- 2w 2 z -2z
)
u2 + 1 v)
w-
185
?
)
u-
11-
2
v +1
V
2
w2 +1
)
11'-
u-?
0
)
v-
-8xyz -I
w2 +1
-I
0
by C 1 -C 2
C2 -C 3 =
-8xyz(u2 + v2 + w2 + I)
=
-8xyz
[
X2
y2
?
')
')
x- + y + z- + 1-
=-8xyz
;also.
2u(l-x2 - / 8(ui' v2 , w3 )
x
o
_Z2)
o o
=
8l1vw(l- Xl
8
,.,
1'-
1-r-
=--2
8U;,.I;, .I;)
]
J
-8xyz I-r
Z2
+--- 2 +--? + 1 I-r- I-r l-r-
--J
~
.
o
2u(l_x2 -
0
y2 _Z2)
o
2u( 1- Xl
_ yl _ Z2 )
_ y2 _ Z2 )3
Y
JI-?
.
z
(l_r2)3=8~yz(l_r2)112
.J1-r 1
_. 8(11, V, w) 1)3 -8xyz ~ Hence trom (I), we have - - - - = (- - , 8 (I 2)1/2 8(x,y,z) I-rxyz -r
=
I
r:-? ',"1_1'2
X Y Z Cor: Given U=- V=- W=k' k' k
where k
= VII -
X
2
-
' f iIn d •,(x,y,z) Y, - Z-, (u, v, w)
in terms ofk. lAns: k
S
]
Engineering Mathematics - I
186
Exercise - 3(H) I.
Ifx+y+z = u ,y+z = uv, z = uvw, then prove that
8(x,y,z)
----'----=--~
=- II
2
V
8(u, v, w) ~)
8(u,v) x--y2. Ifx + l + u _v = O,uv + xy = 0, prove that 8(x,y) u 2 + v 2 then 3. If u1 =x1 +X2 +x3 + x,pU 1U 2 =x2 +X3 +X4 ,1lIU2113 =x4 8(X p X2 ,X3,X4 ) 3 2 =1I1 .il2 ·113 8(u l , 112 , u3' 114) 4. If u 3 +V3 +W3 =X+ y+z,u 2 +V2 +W2 =X3 + / +Z3, 8(1I,v,w) (y-z)(z-x)(x-y) t Ilen SlOW t Ilat = ---"'----'---'..-----'-----'-----'----"-I 8(x,y,z) (v- w)(w-u)(u - v) 2
2
2
prove
that
3.7.4 Use of Jacobians in determining functional dependence and independence of functions. Let u and v be two functions of x and y connected by the relation v = f(u), then we say that u and v are functionally dependent. We shall prove that the condition for functional dependence is
o(u, v) = 0 o(x,y) Consider
w = v - feu) = 0
Now w is a function of u and v where u, v are functions of x and y. w is a composite function of x and y
ow ow au Ow Ov -=-.-+-.-
ax ax ax ov ax ow ow au ow Ov -=-.-+-.ax au ax ov ax But w considered as a function of x, y is identically zero. i.e.,
Ow -=0
ax
'
..... (i)
..... (ii)
Mean Value Theorems and Functions of Several Variables
187
From (i) and (ii) we get
all ow ov au ax av ax
Ow
-.-+---.--=0
Dw
all ow av
au
~v
--.--+-.-=
and
Eliminating
D1'
..... (iii)
0
ay
..... (iv)
Ow aw -a-·-o from (iii) and (iv) we have 11 I'
all aI' ax ax all -av Oy
=0
~v
au au ax ay ou au =0 ax Oy
--
a(u,11=o
=>
a(x,y)
The concept of functional dependence can be extended to any number of variables. Thus if u, 1\ 11' are functions of x. y, z the 1I, v, w will be functionally dependent (i.e., there will exist a relation between ll. 1\ 11') if
Q~!_,V, w) = 0 a(x,y,z) 3.14.5 Example Are x -I- y - z, x - y + z, x 2 + y2 + z2 - 2yz functionally dependent? I f so, find a relation between thelll.
Solution Let
lI=x+y-z
..... (i)
1'=x-y-l-z
..... (ii)
11' =
Xl
-I-
I
-I-
i2 -
2yz
..... (iii)
Engineering Mathematics - I
188
o(U, v, w) = o(x,y,z)
Then
\ \
2x
-\
-\
2y-2z
2z-2y
-\
-\
2y-2z
2y-2z
= 2x
=0
(since 2nd and 3 rd columns are identical)
:. u, v, ware functionally dependent To find the relation between u, v, w we notice that, from (i) and (ii) 11
+ v = 2x,
v = 2(y-z)
U -
From (iii)
u2 + v2 = 2w, the required relation 3.14.6 Example
i
11,
If 11 = x + Y + z, v = x 2 + y2 -I- z2, W = x 3 + + z3 - 3xyz then prove that l~ ware not independent and find the relation between them.
We shall now find the relation between them We have
x+y x-y
u=--
x-y x+ y
u
By componendo and dividendo
x
1+u
Y
I-li
-=-
A
191
Mean Value Theorems and Functions of Several Variables
Now
Y
xy
o X
v = (x - y)' = (;
r
V
1
x
1+11
y
l-lt
'V=C~:~2 J u2(4v+ 1) = I The required relation of functional depedence between
3.14.8 Example Examine for functional dependence of u = sin-I x + sin-1y
v=x~l-y2 +y,JI-x 2 and find the relation between them, if it exists.
Solution We have
u =-===
,
~1_x2
lilly
a{u, v)
1I (
It )'
a{x,y)
Vx
Vy
lI,
v.
by (A)
192
:. ll, V
Engineering Mathematics - I
are functionally dependent we shall now determine the relation between u and v Let A=sin-1x,B=sin-1y; =>x=sinA y=sinB The given function can then be written as u = A +B 2
2
v = sin A)1-sin B +sin BJl-sin A = sin(A+B)= sinu Hence the required relation of functional dependence between u and v is v = sinu
3.14.8(a) If u
= eX sin y, v = eX cos y,
show that u,v are functionally independent
au ax
au ay
av
av
ax
ay
-
a(u, v) Solution: The Jacobian a(x,y)
X
•
X
e smy e cosy eX
cos y
_ex
sin y
=
_ex
* o.
:. x, v are functionally independent.
.
x-y x+z
x+z y+z
3.14.9 VerifY whether u = - - , v = - - are functionally dependent, and if so find the relation.
Solution: Treating u,v as functions of x,y regarding z as an absolute constant, the condition for functional dependent in
a(u, v) a(x,y)
=0
au ax
Now - =
,
y + z au 1 av 1 av = - = - - - and -=~-, (x+z/' ay x+z ax y+z ay
x +z (Y+Z)2
Mean Value Theorems and Functions of Several Variables
B(u,v) B(x,y)
=
y+z (x + Z)2
193
---
X+Z
=0
x+z (y + Z)2
1
y+z
u, v are functionally dependent.To lind the relation between u and v, we eliminate 'x' between the given relation viz: u
x-y
= ------------- (A),
x+z
v = ---------------(B)
X+Z
y+z
+Y 1-11
liZ
from(A) ux+vy = x-yo =>(I-u)x=uz+y => x = - liZ
+Y
Substituting for 'x' in 'B', v = 1-11
+z
y+z
=
(y+Z) I-l/' y+z
I-u
Hence the required functional relation between 'u' and 'v' is v( I -u) = I.
Exercise - 3(1) I.
Verify whether the following are functionally dependent and if so, lind the relation between them.
x- y x+ y x+y x x-y ( b) u = - - , v = tan -\ x + tan -\ y I-xy (a) II = - - , V = - -
(c)u=x-y,v= xy x+y (X+y)2 2.
(Ans: u( I+v) = 2) (Ans: v = tan (Ans: 4v + 1I
2
-\
II)
=1)
Prove that the following functions are not independent.
= X + 2y + Z, v = x- 2y + 3z, W = 2xy - 2z + 4yz - 2Z2 2 (b) II = x+ y- z, v = x- y+ z, w= x + l + Z2 - 2yz 3 (c) u = x + y + z, v = xy + yz + zx, w = x + l + Z3 - 3xyz
(a) II
3.
If u
= x + y "v = y + z Z
functional relation.
x
w
= y(x + y + z) , xz
show that u,v are connected by a
Engineering Mathematics - I
194
4.
x-y l-xy
Examine for functional dependence between u = - - , v = tan-I x - tan -I y.1f ( Ans: u=tan v)
dependent, find the relation. 5.
Prove that functions
II
=Y + z, v = x + 2z 2 , W = x -
dependent and find the relation between them.
4yz - 21 are functionally
(Ans: V=W+2l?)
3.15.1 Maxima and Minima The method of finding the maxima and minima of functions of two and three independent variables is discussed in this topic. The maximum or minimum values are also called Stationary or extreme or turning values. For finding these values we use Taylor's theorem for functions of two variables.
3.15.2 Taylor's theorem for a function of two variables: We know that Taylor's theorem for a function f(x) of a single variable 'x' is
f(x+h) = f(x) + hj\x) +
~~ /I(X)+ ...... .
....( I) Now, consider f(x,y), a function of two independent variables. Keeping 'y' constant and following (I) we get
. a h a j(x+h,y+k) = f(x,y+k)+h-f(x,y+k)+---J j (x,y+k)+ .........(2) 2
ax
2
.
2! ax-
Keeping 'x' constant and by applying (I), we get
a. ea f(x+h,y+k) = f(x,y)+k-j(x'Y)+---2 f(x,y)+ ...... . 2
ay
2! ay
.... (3) Substituting (3) in(2) we get
. a. ea j(x,y)+k a;f(x'Y)+2T ay2 f(x,y)+ ...... . 2
f(x+h,y+k) = [
]
2
a [ f(x,y)+k-j(x,y)+---) a. e a f(x,y)+ ...... .] + hax ay 2! ay 2 2 112 a [ aj. a ] + 2T al f(x,y)+k iY (x,y)+2T al f(x,y)+ ....... +..... .
e
Mean Value Theorems and Functions of Several Variables
195
Substituting (iii) in (ii) we get 2
2
f(x + hy + k)
=
k a f(x,y)+ ....] a [f(x,y)+k-f(x,y)---2 ay 2! ay
2 a [ f(x,y)+k-f(x'Y)+---2 a a f(x,y)+ .....] +hax ay 2! ay
e
2 2 2 h a [ f(x,y)+k-f(X'Y)+---2 a a f(x,y)+ .....] +---2 <2 ax ay < 2 ay
e
+ .....
[
aj
[aj
2 2
2
k a f + .....] + h-+hk--+ a f .....] = f(x,y)+k-+-2 ay
< 2 ay
+[~ a <2 ax
+. . .]
2
{
f(x + h. y + k)
Thus
=
ax
f(x, y) + (h! + k
axay
~ )f
(a
a)2 f+·····
1 +~hax+kay
..... (iv)
Expression (iv) if known as Taylor's expansion for functions of two variables. Writing x
=
a, y
=
b expansion about (a, b) is given by
a
a)2
[a
a
If(a+h,b+k)=f(a,b)+ ( hax +kay
ira,b)
1 h-+k- ]2 i((jb)+ ..... +1! ax ay , Now if we call h = x - a, y - b = k we get
f(x,y) =f(a, b) + ((x-a)!
+(y-b)~)i((j'b)
Engineering Mathematics - I
196
1 ( (x-a)-+{y-b)a a )2 fiab) + ..... . +-
ax
2!
ay'
..... (v) a = 0 = b, we have
Writing
f(x, y)
j(0, 0) + (x
=
~+y;
)./(0,0) + ......
.. ... (vi)
This is Maclaurin's series for two varialbes. 3.15.3
Example
Solution
Ix (x, y)
=
fyy (x, y)
=> f{I,I) = e 2
f(x, y) =- ~2 + y2
We have
=
2xe x
e
2+ 2 Y
X2 + y2
=> Ix (I, I)
2
+ 4Y eX
2
7
=
2e
2
=> fyy (I, I)
+Y-
=
6e2
Hence putting these values in second form of Taylor's theorem, we get e,2+y2
' /x Y')/ = () => j'»)')'(I , -2) = 0 j .Y.ly\" Substituting a = -1, b = 2 and above values in (i) we have
x 2y + 3y - 2
10 -4(x -I) + 4(y + 2) - 2(x-l;2
= -
+ 2(x- -1) (y + 2) + (x -
1)2
Cv + 2)
3.15.6 Example Evaluate log[(l.03) I] + (0.98) '4 -1] ,approximately.
Solution
h)~ + (y + k)~ -
Then
f(x + h. y + k)
Assuming
x = 1, y
1, II
=
0.03, k
We get
F(x+h, y + k)
=
IOg[(l.03)~ +(0.98)~ -IJ
=
=
of ox
I -x 3
2
3
IOg( (x +
=
[
oy
x 3 + y4 _I
=
0 .02
oF
OF)
.
F(x, y) + h-+ k - approximately ..... (A) ox oy
1 --3/4 -y
of
1
= -
I)
4 I -
-
I
x 3 + y4 -1
F(x, y) + h
. ox + k [OF) oy approximately (OF)
Mean Value Theorems and Functions of Several Variables
199
~ IOg[ (1.03)~ + (0.98)~ -I] =
F(I, I) + 0.03
(aF)
ax
+
(_O.02)[aF)
(1,1)
ay
(1,1)
~O.03 ll+I-IJ r j to.oJ ~ 1~O.005 ll+I-IJ
approximately
(-: F(l,I) =
0)
Exercise - 3(J) I.
Obtain the expansionllsing Taylor's theorem of the following:
I. xy2 +
cos(xy) at
(I'~)
2. xY at (I, I) lIpto second term
I ADS: 3. sinx. siny in powers of
xY = 1
+(x -I ) +(x -I) x (y -I ) + -21 (x_I)2 I
(x -~) and (y -~) ~+~(x_~)+~(V_~)_~(X_~)2
( ADS. ·2242'
4.
444
eax sinby at (0, 0) ( ADS: hy + xyah + '" )
200
Engineering Mathematics - I
.
5. Prove that sm(x +y)
( x+
=x +Y-
y)
<3
3
+ .....
2
6. Show that eY 10g(1 + x) = x + xy - ~ approximately 2
7. Expandf(x, y)
=
sinxy in powers of (x -1) and (y [ ADs • I - -1t
•
~) upto second degree term.
2 (x - 1)2 - -1t(x -1 \
8
2
y - -1t) - -1 ( y - -1t)2 + ..... I 2 2 2
3.16.7 Maxima and Minima The function F(x, y) is maximum at (x, y) if for all small positive or negative values of hand k; we have
f(x + h, y + k) - F(x, y) < 0
.... (i)
Similarly, F(x, y) is minimum at (x, y) iffor all small positive or negative values of hand k; we have
F(x + h, y + k) - F(x, y) > 0 Thus ifF(x, y) has a minimum or maximum value then [F(x + h, Y + k) -F(x, y)] keeps the same sign for all small positive or negative values of hand k. Now using Taylor's expansion, we get
F(x + h, y + k) - F(x, y)
=
+
-f
(
aF) haF -+k ay
ox
2 2 2 -2-(h2 a F +2hk a F +k2 a F) 21 ax 2 axay 0'2 + .....
..... (iii)
Assuming h, kto be sufficiently small; the sign of the expression on the left can
aF aF be made to depend on that of h ax + k 0' Hence the necessary condition for f(x, y) to be a maximum or minimum is that
Mean Value Theorems and Functions of Several Variables
of
of
ox
oy
201
- = 0 and - = 0
of p
I.e.,
=
0 and q
then (iii) ~ f(x + h y .
-I-
0
=
where p = ox'
k) - j(x y)
= -
I[
'2!
2 2 2 ,0 F 0 F , 0 F] 11- --) + 2hk-- + K- --) + .....
or
I 2!
2
=r- [ h-,
2!
r
=-
2!
z= 0
ox8y
oy-
r rh 2 + 2shk + tk2 ]
= -
02F s=-ox8y'
of q = oy
~l
~y-
s t ,] + 2hk-+-K-
r
r
s - s ) k, )+ (2 -,-k- +1-) [( h+-k r ) rr-
II
=~[(h+~k)2 +(rt-.~.2)k:l 2! r rIt has the same sign for all h, k if and only if (rt - s2) is positive. When
(i)
(rt - s2) positive,j(x, y) is maximum for negative 'r' and minimum for positive 'r'.
(ii)
(rt - s2) is negative, we have neither a maximum nor minimum and such stationary points are called saddle points.
Solvingp = 0, q
=
0, we get the extreme points
Extreme value :
.I(a, b) is said to be an extreme value ofJ(x, y) if it is either maximum or minimum. 3.16.8 Example Find the maximum and minimum values of x 3 + 3xl- 31 + 4.
Solution We have
J(x, y)
=
x 3 + 3xl- 3x2
fx = 3x2 + 31- 6x,
fy
-
31 + 4
= 6xy - 6y,
f)y = 6x- 6,
1;)' = 6y
Engineering Mathematics - I
202
We now solve
Ix = o,f;, =
0 simultaneously
I.C.,
x2+y~2x=0
and
xy
~Y =
..... (i)
0
..... (ii)
xy ~ y = 0 => y = 0 or x = I From
x 2 + y ~ 2x
when
y = 0, .~ ~ 2x = 0
0 we get
=
x = 2 or x = 0
when
x= l,y~1
=
0 i.e.,y= I or~1
Thus the points are (0, 0), (2, 0), (I, I), (I, ~ I) (a) For
x = 0 , y = 0 , r =j'xx = ~6 ' s = 0 ') t = ~f) rt ~
s2 =
36 ~ 0 > 0
j(x, y) is stationary at x = 0, y = 0 But
r=/xx =
~6
<0
fix, y) is maximum at x = O,y = 0 and the maximum valuej{O, 0)
=4 (b) Por Again
x = 2 , .y
= 0,
r =/xx = 6 ' s =/.xy = 0 , t =j'n' = 6
rt ~s2 > 0
j(x, y) is stationary at x But
=
2, Y
=
0
r =/xx = 6 > 0 j{x, y) is minimum at x = 2, y = 0 and minimum value j(2, 0) = 0
(c) At
x
=
I, Y = I or at x = I, Y = ~ I
rt ~ s2 > 0 and we can reject these points as they are not stationary points.
3.16.9 Example Discuss the maximum and minimum of f('(,y)
=
x2 ~
Y + 6x -- 12
P = 2x + 6, q = ~2y, r = 2, s = 0, t = ~2 p = 0 gives x = ~3, q = 0
~
y =0
Mean Value Theorems and Functions of Several Variables
203
Stationary point is (-3, 0); and here r = 2 >0; and rl-s~ == 4-0 == 4 > 0 ) ) - 1,... x- -.V + 6x - .!. is minimum at (-3,0) thc minimum valuc is
j'(-3,O)==-21. Exercise -3(k)
....
Find the maximum and minimum valucs of: x' + 3x/ - 3/ + 4 jAns : (0,0) max. its valuc 4,(2,0) min. its value 'W] 3xy ) - 3.c)-3Y) + 7 x'+ [Ans :max . value 7 at (0,0) min.value '3' at (2,0))
3.
x~ + /
I. I. ')
+6x+ 12
[Ans: mill. value '3' at (-3,O)j
3
4.
x + x/ + 21x -12x2 - 2./
5.
X4
6.
X4
+x/ + / + / _ x2 _
l Ans : max.value
10 at ( 1,0) min. value -98 at (7,0) 1
jAns: minimum value '0' at (0,0)1 y2
+1
lAns: max.value I at (0,0) min.value Y2 at t(Hlr pOints( ±
7. 8.
9.
)
3
0
2
~
~ ± ~)]
1 (1 1)
3
jAns: max.valuc --at -,- ] 324 2 2 ) ) ) 3 Discuss the stationary values of II == x-y-3x- - 2y -4y + [Ans: u is a maximum at(O,-I). The maximumvalue is 5, u is neither maximulll nor minimum at(±4,3) Find the maximulll and minimum values ofthc following functions. (a) x'/(12-3x-4y) jAns:maximulllvalucat(2,1)] X
Y -x Y -x Y
(b)
Xl
+ / - 3axy( a > 0)
[Ans: minimum value is _(/' at (a,a)1 1
[Ans:maximum value is ~ at (a/3,a/3)]
(c) xy(a-x-y) (a>O) 10.
27
Determine the values of x,y for which the following functions are maximum or mllllmum .. (a) x ' y 2(1-x- y) l Ans: maximum for x = 112 , Y = I/3 j (b) (x 2
+ y2)"
-2a\x~
- /)
[Ans: minimum for x == ± a,y = OJ [Ans : maximum for x = y = 7[ /3]
(c) sinx siny sin(x+ y) (d)a[sinx+siny+ sin(x+y)] [Ans: maximum for x = y =7[ /31 (e) x/(3x+6y-2)
[Ans: minimum at x = y = 1/6J
204
Engineering Mathematics - I
Lagrange's method of undetermined multipliers: 3.17.0Theorem: Find the maximum and minimum values of
f(x"X 2'X3' ......xlI )
where x" x2 ' x 3 ' •••••• XII are connected by the following m equutions.
f(xl'x2, .........x,J .Since n variables are connected by the m given
relations, only n-m of the variables are independent. The maximum and minimum values ofu can be found by the method of Lagrange's method of undetermined Multipliers .For u to be max or min of u , du i.e.
=
0
au au au -dx, +-dx2 + ........ +-dxn =0 ax, aX2 aXil also from, we have
..... ( 1) multiplying the above equations by 1, ~, ~, ....Am respectively and adding, we get
The values ~,~, .... Am are at our choice .We can, therefore choose them so as to satisty 'm 'linear equations .
.. .. (3) I(is immaterial which n - m of the n variables are regarded as independent. Let these bexl/I+"x"H2' .......x".Then since the n - m quantities dxm+"dxm+2, .......dx" are all independent, their coefficient must be separately zero. This gives additional equations.
205
Mean Value Theorems and Functions of Several Variables
..... (4)
au axil
+
A, atA aXil
+
~ a¢2
+ ........ + Alii
- axil
a¢m = 0 aXil
Now we have m+n equations in all and given relations ¢,
= 0, ¢2 = 0, .......¢III = 0
The equations (3) and n-m equations (4) are sufficient to determine the m + n quantities A,,~, .... AII/ and
x"x2 ' ••••••••• x
lI
for which maximum and minimum values of u
are possible .The m multipliers are called Lagrange's Multipliers.
3.17.0 (a)Theorem: If u=
discuss the
maximum and minimum values ofu
Proof: For a maximum or minimum, we have du I.e.
=0
af" af" af" -'-' d" + -'-' dy + _0_' dz = 0 ax ay oz
Also from
j; = 0 andJ; = 0
FaJ; dx+ aJ; dy+ a.t; dz=O ax iY az
'" .. ( 1)
aJ; dx+ aJ; dy+ aj; dz = 0 ax iY az Multiplying these equations by 1, A" ~ respectively and adding
a¢ + A, aJ; ( ax ax ...
+
~ aJ;) dx + ( a¢ + A, a.t; + ~ aj ;) dy + l( a¢ + A, a.t; + ~ a.t;) dz = 0 ax
liY
iy
iy
fu
~ ... (2)
The lagrange's multipliers
A"
~ are choosen such that
fu
fu
Engineering Mathematics - I
206
o¢ +:1/11 ()I; + /l, ~ 01; oy oy - ~v
=
0
.... (4)
then equation (2) becomes
o¢ +:1/11 ~/~ + /L1) ~/; ()z oz - oz
=0
•.••
These equations (3) ,(4), (5) and given
(5)
1;<'t,y,z)=0, .I;(x,y,z)=O are
sufficient to find ~, ~ ,x,y,z which give maximum and minimum values of u ..
3.17.0(b) Theorem: If u
=
¢(x,y, z) and f(x,y,z) = 0 find the maximum or
minimum values ofu.
Proof: For a maximum or a minimum I.e.
o¢ ox
o¢ oy
o¢ oz
~.dx+~.dy+~.dz
II
,we have du
=
0
=0
.... ( I)
and
f(x,y,::) = o. of ol ~l - .dx + -.Jy + -.Jz = 0 ox oy OZ
from given
.... (2)
Multiplying (1) by 1 and (2) by A and adding
o¢ + A of)dx+( o¢ + A ol)dY+(o¢'+ A ol)dZ = 0 ( ox ox oY oy OZ oz
.... (3)
Lagrange's multiplier A is choosen such that
o¢ +A ~l =0 ox ox .... (4) and
o¢ +A ~l Oy oy
=0
.... (5)
The equation (3) gives
o¢ +A ~l =0 oz oz The equations (4),(5),(6) and
.... (6)
1 (x,y, z) = 0
which u is a maximum or a minimum.
enable us to tind the values of
A "x,y,z for
207
Mean Value Theorems and Functions of Several Vanables
Exercise - 3(1) 3.17.1 Example
Find the point upon the plane ax + by + cz
=
p at which the function
$ = x 2 +1+ z2 has a minimum value and find this minimum $. Solution
ax+by+:::=p
p-ax-hy e
~:::= ~---~
..... (i) )
Again
$=
')
. .2
J
+ .V"" + z-
x~
= r
')
J
+ ), +
(
17 - axc- ·bV)-
Then for min./max., we have
ox = 2.x -
0$
2a (p _ ax - by) = 0 e2
~$ = 2y- 2~ e-
uy
From (8)
x
(p _ ax - hy) = 0
p-ax-hy
a
) c-
y
p-lLr-hy
h
c2
~-=
~-=
x a
.
bx y=. (l Now, substituting this value ofy in (ii) we obtain y
~-=-'--~
b
2
X) =0 x- 2a ( P - l I Xb- e 1I c 2x - (lP
Then,
and
+ a2x + b2x
=
0
+
x
ap bx bp ) ) and Y = =) ) ) a +h- +c a a- +b- +c-
=
2
b2
( 2) =
x( a 2
+
(lP
..... (ii) ..... (iii)
208
Engineering Mathematics - I
Hence 3 a min, when
p -x =-y =----:---=-::----c2 2 2 Also
a
b
Z =
p-ax-by c
Z
I.e.,
a +b +c
=~[p_ c
2
a p
a2 + b2 + c2
P
=
Hence 3 a min. at
x
y Z abc
-==-
4> . =.xl + y2
+ z2 =
P 2 a +b 2 +c 2
a 2 p2 + b 2 + p2 +C2p2
(a 2 + b2 + c2 r
----'---,---...:...--.---:'--
mm
3.17.2 Example Divide 24 into the three parts such that the continued product of the first, square of the second and cube of third is maximum.
Solution Let 24 be divided into parts x, y, z then x + y + z
x 2 y2z
Given
4>(x, y, z)
then
4>(x, y) =.xl;? (24 .- x - y) from (i)
or
4>(x, y) = 24x3y2 - x4y2 - x 3y3
=
Differentiating, we get
p
= ~! = n.xly2 -
a 4> r =-?
4~y2 -
3x2y3
2
ax-
=
144xy2 - 12x2y2 - 6xy3
=
24
..... (i)
209
Mean Value Theorems and Functions of Several Variables
Since rl - s2 > 0 and r < 0, therefore ~ (x, y) is maximum at (12, 8). Putting x
==
12, and y
==
8 in (i), we get z == 4
The values of x, y, z are 12, 8, 4 respectively. This is the division of 24 for maximum ~(x, y, z).
3.17.3 Example A rectangular box, open at the top, is to have a volume of 32 c.c. Find the dimensions of the box requiring least material for its construction.
Engineering Mathematics - I
210
Solution Let I, band h be the length, breadth and height of the box respectively. Then, wc have 2(1 + h) h + Ih
v
=
Ihh
=
32 ; surface
S
=
2(1
I-
b)h + Ih and b
=
=
32 =
-
III
32) h + I (32) 32 s = 2 ( 1+= 2111 + -64 +III liz I h
~~
Now
=
2/-(:;)
For maximum and minimum ofS, we get
as =0= as al
ah
as = 0 ~ 2h _ 642 = 0 h = 32 0/
'
12
as = 0 ~ 2/- 32 = 0
and
and
1
ah
82~ = +2
8h-
h2
so 'S' is minimum for
1= 4,b = 4,h = 2.
S (say) .....(i)
Mean Value Theorems and Functions of Several Variables
I
211
1
1
3.17.4 Example: If u=a'x"+h'/+c3z1 where -+-+-=1 show that the x y z stationary point of u is given hy x
a+h+c
a+h+c
a
h
=- - - , y =
a+h+c
,,"" - - - c
Solution: For a stationary value of 'u', du =0 I.e.
1
a ' 2xdx+h'-2ydy+c 2:::dz
Also we have the given relation
=0
..... ( I )
~~ + l + l = 1 x
y
:::
1 1 1 -) dx+-) dy+-) dz = 0
y
x~
..
~
..... (2)
Equating the coefficients of dx , dy ,dz from (1)+ A(2) separately to zero
A
I
A y
I
1
A --
We get a x + -) = 0, h x + -) = 0, c x + -, = 0 . x-
1
or
:. a x ' =h ' /=c ' z'-=(-A) ax = by = cz = (-At] = k(say) x = k I a; y = k I h; z = k I c ..... (3)
· . j'or x, y, Z SU hstltutll1g
.
1 + -1 + ~~1 = I Y z
In X
abc k k k :. k = a+h+c
-+-+-= 1
..... (4)
Substituting in (3) we get X=
a+b+c II
,y=
a+b+c a+b+c , z = - - - for which u is stationary. c b Exercise 3(M)
I.
Find the maximum and minimum distances from the ongll1 to the curve 5x" + 6xy + 5 y2 - 8 = O. [Ans : max. 4, min .16,6,3]
2.
Fin~the dimensions of the rectangular box, with out a top of maximum capacity rAns:6",6",3" J whose surface is 108 Sq. inches.
Engineering Mathematics - I
212
3.
If
f = u 1 +j'1 +W1
where
11
+ v + w = 1, show that
f
IS stationary when
1
U=V=W=~
3 4.
Use Lagrange's method of multiplies to determine the minimum distance from the origin to the plane x + 2y + 3z = 14.
Jl4"1
[Ans:
5.
Find the maximum value of / = xyz when xy + yz + zx = k. [Ans: (k / 3)3/?]
6.
If the distance from the origin of any point
3x+2y+z-12 = 7.
°
Lagrange's method). Find the maximum
is P
=~X2 + y2 +Z1
volume
of a
p(x,y, z) on the plane
,find the minimum value of p(use the
parallelepiped
inscribed
b-
h
8a c] [Ans: 3 (:;3
c-
-V j
Given x+y+z=k,findthemaximumvalueof XII,yll,ZY [Ans: If
r2
=
x
2
+
l
+
Z2
12.
1
aa./3II.rY.k
and x+y+z =30 find the values of x,y,z for which
r
IS a
[Ans: x = y = z = 10]
mll111nUm 11.
]
Find the volume of the largest rectangular parallelepiped that can be inscribed in
a
10.
2
3J3
, , ) · ·d x- Y z- 1 t I1e e II IpSo! - 2 +-) +-) = . 9.
a sphere
8r [Ans:
' x-, + y-) +z 2 = r.
8.
III
Find the dimensions of a rectangular box without a top of maximum capacity whose surface area is 108 square inches. [Ans: 6,6,3 inches] Find the dimensions of a rectangular box with open top, so that the total surface area S of the box is a minimum, given that the volume 'V' of the box is constant. fAns: x = y = 2z = (2\,)1/3]
13. 14.
Show that the rectangular solid of maximum volume of that can be inscribed in a sphere is a cube I f the total surface area of a closed rectangular box is 108 sq.cm. Find the dimensions of the box having maximum capacity
[Ans:
.Ji8,.Ji8,.Ji8]
15.
Perimeter of a triangle is constant, prove that the area of this triangle is maximum when the triangle is equilateral.
16.
If u=¢(x,y)where/(x,y)=O, find the maximum and minimul.1 values of 11 using Lagrange's multipliers Method
4 Curvature and Curve Tracing
4.1.1
Curvature The curvature of a curve (bending of a curve) varies from point to point on the curve. Let P be a point on the curve and Q be
c.l
point nearer to P, are OQ = 0 s
Let tangents at P and Q make angles, \.11, and \.11 + 0\.11 with the X- axis.
LT R T' = d\jJ is the angle through which the tangent at P turns as a point moves along the curve form P to Q through a distance os along the curve and hence 0\.11 depends on the arc length
os.
0\.11 is defined as the total bending or total curvature of the arc P Q. 0\.11
8s
is defined as the average curvature of the arc P Q.
Now os
~
Lt
~\.11
&s~O uS
0 as Q ~ P =
dd\.l1 is defined as curvature of the curve at the point P and is denoted S
by K.
Thus
d\.l1 K=ds
Engineering Mathematics - I
214
x'
0
T
Fig. 4.1
Radius of Curvature The reciprocal of the curvature of a curve at any point' P' is called the radius (?f curvature at that point and is denoted by 'p'.
ds
Thus p = -d is called the radius of curvature of the curve at the point P. \1' .
4.1.2 Theorem If 'r' is the radius of a circle then the radius of curvature of the circle is same as its radius.
Proof: Let 'c' be the centre and 'r' be the radius ofa circle. P be any point on it and Q be a neighbouring point, arc PQ = ds. The tangents at P and Q make angles \If and \If + 8~, respectively. From the figure
LT R T' = 8 \If LP R Q'
= 180 - 8\1f
LP C Q'
= 8\jf,
8s
~ = u
~,
from the sector PCQ.
r and as Q ~ P, 8s
~
0
Curvature and Curve Tracing
215
Lt 8 \1' = ~ 8s r
01->0
1
I.e.,
P
r
p=r
y
o
x'
x
y'
Fig. 4.2 4.1.3 Cartesian form of the Radius of Curvature Let y = j{x) be the Cartesian ~quation of a curve and tan (\lJ) be the slope of the tangent to the curve at a point P.
dy tan \lJ = dx Differentiating w.r to 'x'
We know that
ds dx
Engineering Mathematics - I
216
... ,..
,
1+(ddxy )2
[I+(~~)' I+(~)' 1 p=
3/
. _I[ + y( ) '2 I.e.,pwhereYI y~
dy d 2Y = -d 'Y2=-2 X dx
is the radius of curvature of the curve (Cartesian form).
4.1.4 Parametric form of Radius of Curvature Let x = fil), Y = get) be the Parametric form of the equation of a curve.
-d'( = f\t() -dy = g ,(t ) dl
i.e.,
x'
' dl
= /'(t),
y'
= gl(t),
dy dy dl -=dx dx dt
'., dy dx
y'
- = -I X
.
-,
217
Curvature and Curve Tracing
Differentiating w.r. to 'x' on both sides
" .- y·x '" •,{2y X 'Y --=
{X')3
3/
dy
d'y
Substituting -d ' - I 2 values in p = X
ex
[l+(~)T d2
.
Y dx 2
[l+(Hr p=
-
'----'..,..
x'y" - y'x"
{X')3 [(XI)2 + {y')2 y~ P=
where
x'y" - y'x"
I dx x =-
dt '
I
dy dt
y=-
is the parametric form of the radius of curvature,
4.1.5 Polar form of Radius of Curvature Let P (r, 8) be a point on the curve r =f(O). If '0' is the pole, 0 X is the initial line.
LXOP =0 The tangent at P makes an angle \I' with the initial line.
$ is the angle between the radius vector OP and tangent PT.
Engineering Mathematics - I
218
From figure
\If=e+cj> d\lf de d -=-+ds ds ds
de d de +-.ds de ds
1 p
-=-
!
p
[1
= de + d] ds
..... (1 )
de
x
Fig. 4.3 We know
ds de
and
tan
i.e.,
tan cj>-=dr de
Differentiating w.r. to
r2+
(dr)2 de
de
= r-
dr
r
'e' on both the sides
..... (2)
Curvature and Curve Tracing
(1 + tan
219
2
d~
de
<1»
=
d~
de
d de
The value of
I+
~~
r-~~
t~
I + r' +
de
1+
d,h
-'I'
r2
r
(dr)2
+ de
r2 +2(_dr)2 _rtl_2_r 2
de de = --'-----'------:---
de ds Substituting the value of de'
r2
+(-:r ,
I+
d~
de
from (2), (3) in (I)
..... (3)
Engineering Mathematics - I
220
is the polar form of radius of curvature.
4.1.6 Pedal form of Radius of Curvature Let P be any point on the curve. the tangent P T and P makes an angle \jJ with the initial line OA, LP 0 X = e (see the figure in 4.1.5) \jJ=e+~
..... ( I)
Differentiating w.r. to s
d\jJ de d~ -=-+d\' ds ds -
I
p
I p
-
We have
=
I de d~ dr -.r- + - r ds dr cis d~
I r
= - (sin~)
de
+ dr I p
sin~=r-, ~
p
I [ sm~+rcos~. d~] ,cos~=dr
= -
~
r
I d(r sin ~) r dr
_ = ~dp p
cos~,
r dr
:. p
~
p = r sin ~
=
dr rdp
where p is the perpendicular distance from pole to the tangent.
4.1.7 Examples 2
Prove that the radius of curvature of the curve y = a cosh (x/a) is
Solution y
=
a cosh (x/a)
L
a
Curvature and Curve Tracing
dy
-
=
dx
221
I a sinh(x/a)a
d 2y -2 =
dx
I cosh (x/a)a
dy d 2 y Substituting the values of - , - - ? in p = dx dx-
( dy)2l·~2
[ I + ~j; d2
Y
dx 2
]3 2
?
P =1[ + sinh - x / a cosh{x / a ).~ p
= aeos h2(x/a)
p
~ a( ~:),
as y
~ cosh(x/a)
p~ ( : ) 4.1.8 Example Prove that for the rectangular hyperbola xy = c2 the radius of curvature at any point r3
is given by p =
-2
2c
where 'r' is the distance of the point from the origin.
Solution Differentiating xy = e 2 w.r. to 'x'
dy
x - +y=O dx
222
Engineering Mathematics - I
..
dy y .
[
( ~b:dy)2l~~
d 1+ Substltutmg for -d ,--? m P = /2 X dx( Y 2
dx 2
p=
as ,.2 = x 2 + rand xy = c 2
4.1.9 Example Prove that the square of the radius of curvature at any point on the parabola
y = 4 ax varies as the cube of the focal distance of the point. Solution From
y=4ax
dy dx d 2y dx 2
-
=
~ f~ I -3/ lax /2 2 -va.
= -
..... (1 )
223
Curvature and Curve Tracing
p=
Hence
y
x'
o
x
y'
Fig. 4.4
p=
-1 r -~~ x 2
-~"a
p
=
1a(x+a)~
..... (2)
The focal distance of the point p (x, y) sp=
~(x-aY +(y-oy
= ~X2 +a 2 -2ax+ y2 =
sp=
~ x 2 + a 2 - 2ax + 4ax
using(l)
~(x+aY
SP = focal distance = (x + a)
..... (3)
224
Engineering Mathematics - I
From (2) and (3)
-2
P = ~ [Focal distance]3/2
4
p2 = - [Focal distance]3
a
p2
00
= cube of the focal distance.
4.1.10 Example
Show that the radius of curvature at an end of the major axis of the ellipse Xl
a2
y2
+---;;: =
1 is equal to the semi
latu~-rectum
of the ellipse.
Solution
i ..e.,
..... (1 )
Differentiating (1) w.r. to 'x'
dy dx
2
2fy-xdyl dx
d y b dx 2 = - --;;
/
y2
Curvature and Curve Tracing
225
I
y
x'
o
A'(-a,O)
(0,0)
A(a,O)
y'
Fig. 4.5
d 2y dx
1
b2 2
( using ax'2 + b2 y' = I
-2 ----
a
i
[1+(:)']" p=
Hence
d 2y dx 2
[ l+---~ h" 4 a i p= _b 4 2 a
r;
(a 4i+b 4x2) a 4b 4
i
The coordinates of the ends of the major axis are A (a, 0) and A' (-a, 0) Taking the end A (a, 0) we have from (2)
pl(a,o)
b2 = -;; =
Semi latus-rectum of the ellipse.
..... (2)
226
Engineering Mathematics - I
4.1.11 Example J
J'
5<
Prove that p at any point of the curve x;; 3 + y/3 =a 3 is three times the length of the perpendicular from origin to the tangent at that point.
Solution The parametric equations of the curve are
x
=
acos 3e,
y = asin3e
dy
3a sin 2 e cose dy - 3a cos 2 e sin e dx
de
dy dx
-------,,---- - - = -
dx
=
tan e
de and
d 2y dx
de -sec2edx = - sec2e / (-3acos 2 e sin
-= 2
e)
3a cos 4 e sin e
gives
Hence
p=
(I + tan 2e
f2
I 4
3a cos e sin e p
=
..... ( I)
3a cos e sin e
The perpendicular distance from the origin to the tangent
dy y-xp
= ---;===dx==
I+(:J
P = a sin e cos e
From (1)&(2) p=3p
2
a sin 2e - cos e(- tan e)
2 JI + tan e
..... (2)
227
Curvature and Curve Tracing
4.1.12 Example Find p at any point for the cycloid
x = a (8 - sin 8) y = a(1- cos 8) and also find p at 8 = 90°
Solution dy dy _ de _ a(sin8) dx - dx - a(l-cos8)' de 2
d y
and
-
dx
2
=
dx
2 8 1 d8
cosec - - 2 2 dx
1 28 -2cosec 2
d 2y 2
=
-1
a(1-cos8) ,
[1+(:)'r
..
p=
d 2y dx
p
8 dy = cotdx 2
-
= - 4a sin
( at 1+cot
p=
2
2
2
-1 4a sin4 ~ 2
8
2
at 8 = 90°, P = - 4a sin 45° i.e., p = -
212 a
4.1.13 Example Find p at (r, 8) for the curve r = a (1 + cos 8) and also find at 8 =
Solution r
= a (1 + cos 8)
dr de = r) = a(- sin 8) d 2r - 2 =r =acos8 de 2
~.
228
Engineering Mathematics - I
[a 2(1 + cos e)2 + 2a 2 sin 2 e + a 2cos e(1 + cos e)] [2a 2 (1 +cose)7i] 2
2
2
a + 2a + 3a cose
2
[2a .2cos
2
~~]7i
3a .2cos ~~ 2
2
=
.iacos e/ 3 /2
e=~
at
4 P = -acos 45° 3
=
2 r;:; -,,2a 3
4.1.14 Example Find p for the curve r'" = if! cos me Solution r'"
= if! cos me
mlog r = m log a + log (cos m e)
1 dr m -. r de r)
= -r
=
1 cos me
(-
m sin m e)
tan me
..... ( I)
ar r2 = -r m sec2 me - tan m e. de r2 = - rm sec2 me + r tan 2 me From (1)
=
3 {r 2+ r2 tan 2me)2 r2 + 2r2 tan 2me + mr2 sec 2me - r2 tan 2me
..... (2)
229
Curvature and Curve Tracing
r
{m + J)cos me am
r m
{m+J( m a
For the curve ,-In = d" cos m
P=
e
(m + J)r m- 1
4.1.15 Example Find the radius of curvature of the curve p2 = ar
Solution
p2
=
ar
..... ( 1)
Differentiating w.r. to 'r'
2p dp dr dp dr p
=
a
=
a
2p
r dr
= r.
dp
2p a
From(l) 2r
p=-j;;; a
4.1.16 Example 1 Find the radius ofcurvature at the point (p, r) on the ellipse
-2
P
J
J
= -2 + -b2 a
r2 ~b2
a
Engineering Mathematics - I
230
Solution
Differentiating 2 dp p3 dr
-2r a 2b 2 2
=>
dr a 2b r-=-dp
p3
4.1.17 Example
Find the radius of curvature p of the curve ,2 cos28
=
if
Solution
,2 = cos28 = a2
..... ( I)
Taking logarithms on both sides we get
2 log r + log cos 2 8 = 2 log a Differentiating w.r.to '8' 1 dr 2 -;: d8 - 2 tan 28 = 0
d8
r-=cot28 dr
d8
tan <\> = r dr tan <\> = tan
=
cot 28
(~ - 28 )
1t
<\>=2- 20 We have
P = rsin<\> => p = r sin P = r cos 28
Substituting
cos28 =
a2 -2
r
from (1)
(~-28)
Curvature and Curve Tracing
231
dp = dr
~ ,.2
dp ,. p=r-=---) dr -(r
4.1.18 Curvature at Origin (i)
When y can be expanded in powers ofx by some algebraical or trigonometrical method from the equation of the curve as qxl
y = px + -
2!
dx
y
x=O
y=o
d2 q = --) dx- x=o
y=o
3
..
p (at ongll1) =
(ii)
dyl
+.... then p = -
(I + p).i
--'----~
q
Newtonian Method: If a curve passes through (0, 0) and the tangent at (0,0) IS
X - axis, then we have x = 0, y =
°and d,Y = °at origin. {X
The expansion of y by Maclaurin's theorem reduces to Xl
Y
= (yo)
+ (YI)O +
° °+ °. +
2
(Y2)O
Taking upto terms of (x 2 ) Y=
x
x2 q - +... 2!
Dividing each term by x 2 and taking limits 2y q= L t x->o x 2
232
Engineering Mathel\latics - I
.,
We have
P = (at ongln)
as
p = 0 at (0,0)
and
Lt 2 q = x~Ox
and
P (at origin) =
3
(I + p)2
= -'-------=--<-q
2y I
q
p (at origin) = Lt
x~o
(~) 2y
(iii) Similarly, if a curve passes through (0, 0) and the tangent at (0, 0) is Y-axis 2
then
p (at origin) = Lt (Y2 x~o
X
)
(iv) Curvature at pole: If the initial line is the tangent to the curve at the pole then
2 X2 ) r2 cos 8 (.: x = r c.os 8) p= Lt = Lt x~o ( 2y x~o 2rsin8 y=rsm8 r
8
2
P = Lt -·--·cos 8 e~o 28 sin8
-
Lt - r )
p - e~o ( 28
8
(as Lt - .- ~ 1.cos28 ~ I) e~osm8
4.1.19 Example Find p at the origin of the curve y4 + xl + a (xl + y) - a2y method.
= 0 by Newtonian
Solution Equating to zero the terms of the lowest degree in the equation of the curve
- a2y =?
=
0
y= 0
i.e.,
x - axis is the only tangent at origin
Then
'p' at (0,0) when X - axis is the tangent is given by
x2 p= Lte~o2y
..... (1)
233
Curvature and Curve Tracing
Dividing the given equation by y x-J
I
+ x· -
x2
+ a. - + ay = y y
a2
x~O
Let So that
y~O
and
-
From (1)
o + 0.2p + a2p
x
2
=2p
y
+ a.O = a2 => 2 a p = a2 => p =
a
"2
4.1.20 Example Find p at the origin for the curve x = a (8 + sin 8), y = a (I - cos 8) by Newtonian method.
Solution
x = a (8 + sin 8),
y = a( 1 - cos8) is a cycloid
X-axis is the tangent to the curve at (0, 0)
p at (0, 0) = Lt
(~l
x---+o 2y
2 Lt a (8 + sin 8)2 0---+0 2a(1- cos 8)
a (8 + sin 8)2 L t - -'-----'-0---+02 l-cos8 ' =
(Of the form % )
Lt a (8 + sin 8)2(\ + cos 8) sin 2 8
0---+02
=
Lt
!!.(~+ 1)2 (1 + cos 8)
0---+02 sm8
a
8
= -2 (1 + 1)2 (I + 1) using Lt = I 0---+0 sin8
a [8] 2 P = 4a
p=
-
Engineering Mathematics - I
234
4.1.21 Example Find pat (0, 0) for the curve 2x4 + method.
31 + 4x2y + x Y - y2 + 2x = 0 by Newtonian
Solution Equating to zero the terms of the lowest degree in the given equation of the curve 2 x = 0 => x = 0 y-axis is the tangent to the curve at (0, 0) 2
P at (0, 0) = Lt L x->o2x ..... ( I)
So that Dividing the given equation by 'x' 2x3 + 3y2. L
2
x
Taking
y2
+ 4x y + Y - -
x
+2
=
0
x~O,y~O
2.0 + 3.0· (2p) + 4.0 + 0 - 2 P + 2 = 0; p = 1
Exercise - 4(A) I. Find p for the following curves : 1.
y
=
4 sin x - sin 2x at x
n
=
2" (Ans:
515 ..
-I 4
3J2a
[Ans: - - ] 16
235
Curvature and Curve Tracing
5.
E +JY =J;
at
%,~ a
lADs: J21
= 1t
6.
x = a (0 - sin 0), y = a (I - cos 0) at 0
7.
= a sin 2 t (I + cos 2t) Y = a cos 2 t (I - cos 2 t) at any point 't'
8.
x = a (cos 0 + 0 sin 0) y = a (sin 0 - 0 sin 0)
9.
? = a 2 cos20
lADs: - 4 al x
lADs: 4 a cos 3 tl
lADs: a 01
a2
[ADS: );]
10.
r =
~
(\
+ cosO) at
(~,~) J2 [ADS·. -3a ]
an
lADs:
4.2
{n + 1)rn-1 1
Centre of Curvature, Circle of Curvature and Evolute 1. Definition: The centre of curvature at a point 'P' of a curves is the point 'C' which lies on the positive direction of the normal at 'P' and is at a distance 'p' (in magnitude) from it. 2. Definition: The circle of curvature at a point' P' of a curve is the circle whose centre is at the centre of curvature 'C' and whose radius is 'p' in magnitude.
236
Engineering Mathematics - I
y
o
x'
y'
Fig. 4.6
4.2.1
Centre of Curvature Let P (x, y) be any point on the curve and C (X, Y) be the centre of curvature at 'P' which is on the normal PC at P. Then PC =: p = radius of curvature at 'P'. The tangent P T at P to the curve makes an angle \jJ with X - axis. P B, C A are perpendiculars to X - axis, and P R is perpendicular to C A. L.PTX=L.RPT=L.PCR=\jJ I cos \jJ = - sec \jJ
I
(I + tan 2 \jJ2 )2" ..... ( I)
dy dx
..... (2)
237
Curvature and Curve Tracing
..... (3)
and we have
(i.e.,)
x = 0 A = 0 B - A B = OB x = x - P C sin \jJ X
=
x - p sin
PR (From
\jl
From (2) and (3)
X=
X=x-
d~[1 +(tfl)2] dx
dx
d 2y dx 2
Y=CA=RA+RC=Y+PCcoo\jJ From (I) and (3)
~
P R C, PR = PC sin \jJ)
238
Engineering Mathematics - ,
y
o
x'
x y'
.. The coordinates of the centre of curvature of a curve
p (x, y) is C
at
2
2
2
_ dy _ d Y [ _ Yl (I + Yl ) 1+ Yl ] , . Denoting Yl - d 'Y2 - - - 2 we can WrIte C x ,Y + x dx Y2 Y2 C (X, Y) is the centre and 'p' is radius ofthe circle of curvature .. The equation of the circle of curvature is (x -
xi + (y - Yi = p2
4.2.2 Example Find the centre of curvature of the curve Y
= x 3 - 6x2 + 3x + 1 at (1, -1)
Solution
dy = 3x2 - 12x + 3 dx
-
Curvature and Curve Tracing
239
dy
_
dX(I,I)
--6
=-6 dX2 (1,_1)
X=
YI (I +
X -
2
yn
Y X= I-
(-6XI+36)
-6
I+ YI2
and Y = Y + - Y2 1+36 and Y = I + - -
-6
-43 X = - 36 and Y = 6 -43) The coordinates of the centre of curvature are ( - 36, -6-
4.2.3 Example Find the centre of curvature for the cycloid x
=
ace - sine),
Y = a(1 - cos e)
Solution
dy = a sin e _ cot e/ dx a(l-cose) 12 d 2y 2 e 1 de 1 1 -=cosec -.--=x---,------,dx 2 2 dx 2· 2 e a(l cos e) Sin 2
. 4 4aSIn
e
--
2
240
Engineering Mathematics - I
• 4
4a SIO
X=
ace -
e
-
2
. 4e e 4 aslO -cot-sine) + _ _-'2"'--------=-2
. e
SIO-
2
ace - sine) + 2a sin e X = ace - sine) X=
2
1+ YI
Y=y+-Y2
(l+cot2E!) Y=a(l-cose)12 . 4 - 4 a SIO
e
-
2
Y = a(l- cose)- 4a sin 2 ~ 2 Y = a(l - cose) - 2a(1 - cose) Y = -a(l - cose) .. The coordinates of the centre of curvature are
[ace + sine),
-
a(1 -
cose)]
4.2.4 Example Find the centre ofthe curvature at any point (x, y) on the eIlipse.
Solution ..... (1 )
241
Curvature and Curve Tracing
Diff. (I) w.r. to 'x'
2b 2x + 2Q2y dy dx
=
0
..... (2)
..... (3)
Using (I) (2) and (3)
X= x _ YI (1 + yO Y2
[~: (a' - h'), -;.' (a' - h' 4.2.5
ij
is the centre of curvalure,
Example Find the circle of curvature of the curve
fx + fY
=
fa
at the point
(~, ~)
Solution ..... (1)
Engineering Mathematics - I
242
Diff (1) w.r. to 'x' 1 1 -1 -I 1 --I dy -x 2 +_y2 -=0
2
2
dx
dy
JY
..... (2)
dx--.r; and
:10/
=
ai
14' /4
[
d 2y dx 2
~1
, I ~-I ----..jy-x dy C 1 ~-I] ..;x-y 2 dx 2
=
X
d 2y
dx
..... (3)
4
..... (4)
a
0/ (// /4' )14
y~ ¥= [1 + 1]% = ~
p = [1 + Y2
i
.fi
a a a is the radius curvature of the curve at ("4'"4 )
X= x _ YI (1 + YI2 ) Y2
a 1(1+1) 3a X=-+--=4 4 4 a
Y=y+
(1 + y2) 1
Y2 a
(I + 1)
Y="4 + -4- =
3a
4
a 3a 3a)
Coordinates of the centre of curvature ( 4' 4
..... (5)
243
Curvature and Curve Tracing
Equation of the circle of curvature of the curve (x - X)2 + (v _ y)2
3a)2 ( 3a)2 _ ( X -4- + y -4- -
=
p2
(~)2 Ji
Exercise - 48 I.
Find the coordinates of the centre of curvature of the curve y = x J at the point
7 31 IAns: - - ) 64' 38 2
2. Show that the centre of curvature for the curve y
=
x x+ 9 at
(3,
6) is
(3, I;)
3. Show that the centre of curvature at the point 't' on the ellipse x = a cos t, .
4.
a- - b 2) . 2 , )cos-I,- () b Slll-'
a - ba
.
y=bSllltlS
[(
1
Find the centre of curvature at (
1
1
3;1, 3;1) of the folium x
3
+I
=
3a xy
21(1 21a IAns: - - --I 16 ' 16
5. Find the centre of curvature at the point '0' of the curve
+ a sin 0
x
=
(I - aO) cos 0
y
=
(I - aO) sin 0 + a cos 0 IAns : a sin 0, - a cos 01
6.
Find the coordinates of the centre of curvature of the cycloid x = a (0 + sinO), y
=
a( I - cos 0) (Ans: a (0 -sin 0), a (I + cos 0) 2
7.
For the curve y x2
+I
= 111X
= a 2 (I
+ ~ show that the circle of curvature at the origin is a
+ 1112) (v -
I1lx).
Engineering Mathematics - I
244
4.3 4.3.1
Evolutes Definition .' Corresponding to each point on a curve we can find the curvature of the curve at that point. Drawing the normal at these points we can find centre of curvature corresponding to each of these points. Since the curvature varies from point to point, centres of curvature also differ. The totality of all such centres of curvature of a given curve will define another curve and this curve is called the evolute of the curve. Thus the locus of centres of curvature of a given curve is called the evolute of that curve.
Notes: (i)
Elimination of x, y from the coordinates of the centre of curvature (X, Y) gives an equation in terms of X, Y which is called evo/ute of the given
curve y = fix). (ii)
When the equation of the curve is given in the parametric form (say t) elimination of '1' from the coordinates of the centre of curvature (X, Y) gives an equation in terms of X, Y which is called evo/ute of the curve x = fit), y = get).
4.3.2 Example Find the evolute of the parabolay = 4 a x
Solution Differentiating
y
=
4ax
2ydy dx dy dx
=
4a
Fa Fx -Fa
=
2
d y dx2
w.r. to 'x'
..... (1)
=-3-
..... (2)
2x2
X= X _
YI
(1 + yl2 ) Y2
X=x
-~(1 +~) Fx -Fa 3
2x 2
x
Curvature and Curve Tracing
245
X=3x+2a
..... (3)
I-Fa
I + Cl ) 1 + y2 X Y = Y+ _ _I = Y+ ~----=c'y)
,
2x 2
..... (4\
Squaring y2
4 -.x3
=
a
X-2a
From (3) substituting x = --3- in (5)
4(X-2a)'
y2= _ a 27 ay2 27 a
r
3 =
4(X - 2 a)3
=
4(x - 2 0)3 is the evolute of the curve.
4.3.3 Example
Solution ..... ( I)
Differentiating (I) with respect to 'x' dy 2b 2x + 2a2y dx
dy dx
=
=
0
2
-b x 02y
..... (2)
Engineering Mathematics - I
246
d 2y dx2 X
=
=
_b 2 a 2b2 a 2 y2· a 2;.
=
_b 4 a 2y 3
..... (3)
(I
x _ YI + yl2 ) Y2
X
Substituting c?y2
=
=
x
x - -4 (a4y2 2
a b
+ b4;x2)
a 2b2 - b2x 2 from( I)
..... (4)
b4 x 2 1+4""2
=y+
ay
_b 4
a
2
i
247
Curvature and Curve Tracing
Y
Substituting b2x 2
=
y - ~ (0 4; a 2b 4
= 02b 2 -
Y
=
Y
=
y
a 2;
from (I)
-+ a b4
a 2 _b 2
+ b4x 2 )
[a 4;
+
b
2 2
a (a 2
-;)]
y3
b4
..... (5)
Eliminating x, y from (4) (5)
(aX)32+ (bY)}2= (a 2 - b 2 )2:I is the evolute of the given curve.
4.3.4 Example Show that the evolute of the curve 2
2
2
X3 +y3 +a 3
2
2
~
= is(X+Y)3 +(x-Y)l =2a 3
Solution The parametric equation of the curve is x = acos 3S, y = a sin 3S dy dy dx
= de dx
3 a sin 2 ScosS - 3a cos 2Ssin S
de dy - =-tan S dx
d2y de - 2 =-sec 2 S dx dx
248
Engineering Mathematics - I
30 cos 4 esin e ..... ( I)
(I + tan
2
e)
1
..... (2)
From (1) and (2) (X + Y) = a(cose + sine)3 and X - Y = a(cose - sine)3 ')
1 2 -.t 1 -dy =acost -dx =a -smt+--.sec dy t t 2 2 ' dt an2
j [.
r·
-dx = a -smt+ =a -smt+-.IJ dt 2 sm-. t1cos-t sm t
2
acos 2 t dy . and -d SIn t t
dx dt
2
= a cost
dy dy
.J!L dx dt
dx
--
acost
= tan t
acos 2 t sin t
X = a( cos t + log tan
t
t
"2) - a cos t = a log tan "2
..... ( 1)
1+ yf Y2 y
=
a sin t+
a
y=sin t From (1)
t tan 2
From (2)
SIn t=
.
= eX/a a
Y
1+ tan 2 t sin t
..... (2)
Curvature and Curve Tracing
251
2
(i.e.,)
2 tan t 2
a
1+ tan 2 t . 2
y
x
2e
a
2X
I+e
y y
=
a cos II
=
([
Y
a
~( /
a
+ e - X,,)
= ((
cosh
(%)
(%) is the evolllte of the given curve. Exercise - 4(C)
I. Find the evolllte of the curve x
2.
=
c t and y
c = -
t
Prove that the equation to the evolute of the parabola x 2 27 a x 2
3. Show that the evolute of the hyperbola
x2
4ay is 4(y - 2(/)3
=
I is (ax)2/3 - (by)213
=
=
(a 2 + b2 )2/3
4. Show that the evolute of the cycloid x equal cycloid.
=
a (q - sin q) y
=
a (I - cos q) is another
4.4 4.4.1
Envelopes Family of curves: Let us consider the equation of a straight line x cosa + y sina = p. Where a is a parameter. For different values of a, x cosa + y sina = p gives different straight lines but all of them are at a constant distance 'p' from the origin. The set of all of these straight lines is said to form a family of straight lines and 'a' is called the parameter of the family. For a given p, different values of'a' give different straight lines which touch the circle x2 + .Y = p2. The circle x2 + 1 = p2 is called the
envelope of the family of straight lines.
252
Engineering Mathematics - I
A curve, which touches each member ofa given family of curves and each point is the point of contact of some member of the family is called the envelope of the family of curves.
Letj(x, y, a)
=
0 be the given family of curves and a is the parameter.
. . WIt . h respect to " . II y Now d ·f'" I Jerentlatmg a partla
a/(x,y,a) aa
=
0
· · 0 f' a' f romJ'"(x, y, a ) = 0 an daf(x,y,a) = o· . .III EI· Imlllatlon gIven an equatIOn terms of x, y which is called the envelope
aa 0/ fhe family of curves.
If the given equation of the family ofcurvesj(x,y, a) = 0 is the quadratic in 'a', say of the form Aa2 + Ba + C = 0 where 'a' is the parameter and A, B, Care functions of x, y. Then the envelope of the family of curves is B2 - 4AC = o.
4.4.2
Example Find the envelope of the family of straight lines y
= mx + ~ m
where 'm' is the
parameter.
Solution a
y=mx+ -
m
III
..... ( I)
is the parameter
Diff(l) w.r. to 'm' a
o=x--2 m
111=
~
Substituting for' m' in (I)
y
=
2.[;;;
Squaring we gety = 4
ax is the envelope of the family of curves
..... (2)
Curvature and Curve Tracing
Aliter:
253
y
a
=
mx + -
III
m2x -
111 Y
+a= 0
is quadratic in parameter '/11' (A = x, B = - y, C = a)
8 2 - 4 AC = 0 then gives the required envelop. .. Envelope of the family of curves is
i
i - 4(a)(x) = 0
= 4 ax is envelope.
4.4.3 Example Find the envelope of the family of curves I J J J y=mx+ va-nr+bw here'm'ist heparameter.
Solution I
J
J
J
y=mx+ vlrm-+b(y - I1lx)2 = ([2m 2 + b 2
m 2(x 2 - a 2) + m(- 2 y x) + (V
- b2) =
0
is quadratic in parameter '111' The envelope of the family of curves (-2 x y)2 - 4(x 2 - a 2)(1- b2) = 0 b 2x 2 + a 2 = a 2b 2
i
2
2
.;. + -=;- = I is the envelope of the family of curves. a b 4.4.4 Example Find the envelope of the family of curves x cos 38 + y sin 3 8 = c when 8 is the parameter.
Solution Given x cos 3 8 + y sin 3 8 = c
(8 is parameter)
Diff. (I) partially w.r. to '8' - 3 x cos 2 8 sin 8 + 3y sin 2 8 cos 8 = 0 Hence
tan8= y
x
..... ( I)
Engineering Mathematics - I
254
Substituting these values in (]) X
[ ~ X2 Y+ y2
x 21 4.4.5
= C(X2
3 ]
[
+y
13
j=c
x ~ X2 + y2
+ I) is the envelope of the family curves.
Example 2
2
4
Find the envelope of the ellipse -;- + a b
=
I where
a, b are parameters connected by the relation a2 + b 2
=
c2
Solution Given
x2
y2
a
b
-2+ 2
=]
..... (1 )
a2 +b2 =c2
..... (2)
Assume that' a' and' b' are functions of' t' Diff. (]) and (2) w.r. to '1'
2) da + I a3 dt
2 (-
x
x 2 da
y2
(- 32)db
db
b
-3 - + -3 - =0 a dt b dt Diff.
(i.e.,)
a2 + b2 = c2 w.r.
=
0
dt
..... (3)
to '1'
da db 2a- +2b- =0 dt dt da db a- + b- = 0 dt dt
Eliminating da, db from(3),(4) dt dt
..... (4)
255
Curvature and Curve Tracing
Using (I) and (2) we get
x
2
/
£=~=-2 b
a
(i.e.,)
[/4
(14
c
c2
=x2c2
(12= XC
/
-=-
Substituting (12, h2 , values in a 2 + h2
x 2c2 + Ic 2
xl + J,z = I 4.4.6
=
= ("2
c2
is the envelope of the family of curves.
Example Find the envelope of the family of curves
~ + Eh
= I where (I,
b are parameters
{/
connected by the relation ah
= ('2.
Solution Given
~+E a b ab
=
I
= ('2
..... ( I)
..... (2)
Assume that a, b are functions of 't' Diff. (1) and (2) w.r. to 'I' - I-) + da y.-( -I -) = db O x( a 2 dt h 2 dt
x da
y db
- - - + 2- - =0 (/2 £II h dl
Diff.
ab
= ('2
..... (3)
w.r. to 'I'
db da a-+b-=O dt dt
..... (4)
256
Engineering Mathematics - I
da db Eliminating dt ' dt from (3), (4)
x
y
x
y
-
-
ab
ab
x
y
-+!L.=lL~-.!L=JL=~ 2
2
b
a
2ab
Using (I) and (2) we get y
x
£=~=-2 b
2c
a
x
a = 2c
2
ab
Using (2) c2 ~= 2c 2 X
..... (5)
a=2x
a y
ab 2c 2
-=--
b
Using (2)
b=2y
..... (6)
Substituting a, b values from (5), (6) in (2)
(2 x)(2 y)
=
c2
4 xy = c2 is the envelope of the family of curves.
257
Curvature and Curve Tracing
Exercise - 4 (0) 1. Find the envelope of the family of curves: (a) y
(b) ;
= mx + a ~I + m 2 when =
111
is parameter.
2a(x - (I) where a is the parameter.
2; =
IAns : x 2 (c) x cos
(d)
~ +~ a
b
=
=
p where
I when a + b = e
(a, b parameters).
(Ans : (e)
~ + ~ = 1 when a2 + b2 = e2 a
b
ff H +
=
I when a + b = e
fx
+
.JY = -!c ]
fx
+
.JY = -!cl
(a, b parameters).
IAns:
(f)
01
(a, b parameters). 2
2
(Ans : x 3 + y3 In
2
=c 3 ]
m
~ + L = I when dnb n = elll + n ani bill
(g)
Find the envelope of the family of curves
(h)
Find the envelope of the family of straight lines
~+~
=
1 where a & bare
related by the equation d l + bn = en, e being a constant. n
(Ans :
xn+1
II
+ yll+1
/J
= C"+ I ]
Engineering Mathematics - I
258
4.5
Curve Tracing Generally a curve is drawn by plotting a number of points and joining them by a smooth line. If an approximate shape of the curve is sufficient for a given purpose then it is enough to study certain important characteristics. This purpose is served by curve tracing methods. The points to be observed for tracing of plane algebraic curves are given below.
1. Whether the curve is passing through the origin, if so the equations of the tangents to the curve at the origin. Suppose F(x,y) = F(O,O) =
° =>
°
is the algebraic form of the equation of the curve.
the curve is passing through origin (i.e.,) If there
IS
no
constant term in F(x,y) then origin lies on the curve. The equations of the tangents to the curve are obtained by equating the lowest degree terms in F(x,y) to zero. If at 0(0,0) the tangents are: (i) real and coincident then '0' is called a cusp (ii) real and different then '0' is called a node (iii)imaginary then '0' is called a conjugate point 2.
Whether the curve is symmetric about an axis or about other any line. If (i) F(x, y) = F(x-y)
=> curve is symmetric about x-axis
(ii) F(-x, y) =F(x, y)
=> curve is symmetric about y-axis
(iii) F(-x,-y) =F(x, y)
=> curve is symmetric in opposite quadrants.
(iv) F(y, x) = F(x, y) curve is symmetric about y = x (v) F(-y,-x) =F(x, y) curve is symmetric about y =-x
Curvature and Curve Tracing
259 y
y
----+-----------~~x
--------~~~----------.. x
y' =4ax
(i)
(ii)
y
y
------~~~----~x
----~~~----------------~x
x
(iii)
y
3
+y'
~~V
(iv)
-----..;:::.~"'-------~x
(v)
3. Weather the curve intersects the axes, if so the tangents to the curve at these points. 4.
Find the region in which the curve exists (i.e.) the curve is defined. The values of x for which y is defined gives the extent parallel to x-axis and the values of y for which x is defined gives the extent in a direction parallel to y-axis. If y is imaginary for values of x in a certain region then the curve does not exist in that region.
260
Engineering Mathematics - I
Similarly with respect 10 values ofy . 5.
Finding the asymptotes. An asymptote is a line that is at a finite distance from (0, 0) and is tangential to the curve at infinity (i.e.) the curve approaches the line at infinity. (i) Sum of the coefficients of the highest degree terms in x equated to zero gives the equations of the asymptotes parallel to x - axis. (ii) Sum of the coefficients of the highest degree terms in y equated to zero gives the equations of the asymptotes parallel to y-axis. (iii) To find the asymptotes that are neither parallel to x-axis nor parallel to yaxis (i.e.) oblique asymptotes, the following method is suggested. Substitute y =mx + c equation in 'x'as
$
n
(m)x
n +$
F(x,y) = 0 and rewrite the equation as a polynomial
111
n-
l(m)x
n-I + .... +$ =0 n
The slopes of the asymptotes are given by $I/(m) = O. Let the slopes be
The values of'c' can be obtained trom $I/-I(m)
= 0,$1/_2(m) = 0
(if necessary).
Let the corresponding values of c be c1' c2 , •••••
Note: if $1/ (m) is a constant then there are no oblique asymptotes to the curve. 6. Obtain dy from F(x,y) = 0 by differentiation.
dx
If in the interval for
and :
x, dy > 0 then the curve is increasing dx
< 0, then the curve is decreasing in that region.
. · ·· . by -dy 7. TIle po liltS 0f maxima an dminima are given
dx
=
0.
Curvature and Curve Tracing
261
y
~o+-------------------~X
dy dx
-
=
o.gIves
XI,X)
-
If in an interval for x 2
(i)
~
el
~-
> 0 then the curve is concave upwards
(X), X 4 ' .... in -
the figure) and has
a minima in that interval. 2
(ii)
el
~
llx-
< 0 the curve is concave downward (XI' xw .... in the figure) and has a
-
maxima in that interval. 2
(iii)
If at a point
X)
-
(say) in the interval d
~
dx-
=
0 then that point is called a point
of inflection and at that point the curve changes it concavity to the opposite direction.
4.5.1 Example Trace the curve y
X2
+2x
.... ( 1)
=
x+ 1 Solution: Substitution of
X
=0
in (1) gives y
F(x,-y):f:- F(x,-y):f:- F( -x, y) etc 1.
=
:. The curve is not symmetric about any line
can be written as y (x + 1) - x (ie) x
2
-
0:. origin lies on the curve.
2
-
2x = 0
.... (2)
xy - (y - 2x) = 0 Lowest degree term is y - 2x = 0
Equation of the tangent at (0,0) is Y Substitutingy
=
0 in (1) we have x(x+ 2)
= 2x .
= o,x = O,x =-2
Engineering Mathematics - I
262 :. The curve crosses x -axis at (0,0) and (-2,0)
Since x =0 in (I) gives y = 0 only, there is no intercept on y-axis
x :. The curve is defined in the region R -
The
curve
is
not
defined
-I.
at
Hence
it
IS
discontinuous.
{-I}
Coefficient of highest degree term in x «i.e.) of Xl) Hence there is no asymptote parallel to x -axis.
IS
I which is a constant.
Coefficient of highest degree term in y is x + I. :. x + 1=0 is the asymptote parallel to y- axis. From
(2)
dy dx
which
is >
0 always except
at
x
Hence the curve is increasing in the regions. 00
< x < -2, -2 < x < -1, -1 < x < 0 and 0 < x < 00 \ \
\~
\..-
\j.:, \ \ \ \ \ I \
/
\ /
1
2
Fig. 4.7
4.5.2 Example 2
Trace the curve
y-
x +1
- Xl
-1
Solution: ....... ( I)
= -1
Curvature and Curve Tracing
263
°
x= in (I) gives y y intercept is -I.
r (-x,y) = r (x,y) y=
°
gives x
2
=
-I :. Hence curve does not pass through (0,0), and
curve is symmetric about y-axis
°
+ 1= =>
= ± i curve does not cross
x
Coefficient of highest degree term in x is y-I
x -axis.
(from( I»
llence y = I is the asymptote parallel to x -axis Coefticient of highest degree term in y is to y-axis
dy =
Hence y
dx
I
=0
x2 -]
.
gives x
=0
x = ±] are asymptotes parallel II
and Yx=o =
-34 < 0 and =I:- 0
:. The maximum point on the curve is (0,-1) and there are no points of inflection.
y
~
~
l~Y.,
I'J
Ix
I I I
I I
--------1---_ I
----r--------I
I I I
I
I I
I
I I
0
I --------1---
X
I I __ ..JI ________ _ (0,-1) Y=-1
I I I
I
I I I I
I I I
I
Fig. 4.8 The curve is decreasing in the interval
.. an d ·· IS IIlcreaslllg 111
-00
°
<
Ixl < ] and I< x <
00
as dy is -ve there
dx
. t I" i < X < 1 as -dy.IS +ve III l1S mterva dx
4.5.3 Example: Trace the Folium x
3
+ / == 3axy
.... ( I)
Engineering Mathematics - I
264
Solution: Clearly (0,0) lies on the curve .No intercepts on y-axis.
F (y,x)=/+x'-3ayx=0=F (x,y)
:. curve is symmetric about x-axis
Coefficient of highest degree term in y is I. a constant, and hence no asymptote parallel to y-axis. Lowest degree term equated to zero gives xy = 0, x = 0, y = 0 (i.e) the axes are the tangents to the curve at the origin.
y = x in (I) gives 2x' - 3ax 2 = 0 x = 0, 3;1 (ie) the line y = x meets the curve at (O,O)and F (-x, - y)
(3~, 3~)
= x' - /
+ 3axy
= F (x, y) :. curve is not symmetric about (0,0)
To obtain oblique asymptotes: Substitute y 3
= mx + c in (I) and
rewriting in terms
2
2
of 'x' we get (I+m )x'+(3m c-3am)x + ..... =0 (ie) ¢3(m)=I+m'and
¢2 (111) = 3m 2c - 3am, ..... ¢, (m)
= 0 gives m= -1, there is only one oblique asymptote with slope-I
Equation of the tangent to the curve at this point is y
+ x = 3£1 which is parallel to (2)
4.5.4 Example: Trace the curve a 2y 2- x 2( a 2- x 2) = 0 ........ ( I )
Solution: x
=0
F
(-x, y)
in (I) gives y = 0,0 (ie) double point on the curve =
F
(x, -y) = F (x, y) curve is symmetric about both the axes.
Lowest degree terms, in (I), equated to zero gives a are the tangents to the curve at (0.0). y
=
°
in (I) gives x
The curve intersects the x-axis at
(a, 0)
and
Ixl > a
l
is negative (ie) the curve does not exist for
Tangents at the points
ddx2y ) at x ( 2
=
°
gives
x
2
)
= 0 ~ y = ±x
(-(1,0)
From (I) it can be seen that
Ixl > a
= 0, ± .i2' points of maxima and III inima.
a x = ± J2 to the curve are parallel to x-axis .Further
a . . = J2 IS negative.
in 0< x <.i2
x
(/ -
= ±a
for
YI
2
:. x
while decreasing in
a .IS a maxImum . = J2
. d" . POlllt an Y IS IIlcreaslIlg
.i2 < x < a .Thus after tracing the curve in
the first quadrant, the symmetry about both the areas can be utilized for tracing the complete curve. Substituting y
= mx + c
in (I )gives
X4
+ a 2 (mx + c
r -a
2 2
x
=0
266
Engineering Mathematics - I
¢4 (m) = 1 which is a constant. Hence there are no oblique asymptotes to the curve. y
~(_-a~,O~)~--------~----------+-~~X
, Fig. 4.10
4.5.5 Example: Trace the curve y2 ( x 2 - a 2) + a 2x 2 = 0
..... (1)
Solution:
x = Oin (1) gives y = 0 :. (0,0) lies on the curve. F (x,-y) =F (x,y) =F (x,-y) :. curve is symmetric about both the areas. Highest degree terms of y equated to zero gives (X2 - a 2) l
= 0 => x = ±a
:. x = ±a are asymptotes parallel to y-axis. Highest degree terms of x equated to zero gives x 2(y2 + a 2) = 0 => x = ±ia ... No asymptotes parallel to x-axis. Equating lowest degree terms to zero we get a 2(y2 - x 2) = 0 => y ± x are tangents at (0, 0) writing y
= mx + c
in (1) and rearranging in terms of I x I
m2x 4+2mcx3+(c 2 _a 2m2+a 2)x2
-2a 2mcx-c 2a 2 =0
=> c is indeterminate since m = 0 ¢2 (m ) = 0 gives c 2+ a 2 = 0 => c = ±ai oblique asymptotes do not exist. Further there are no intercepts on either of the axes.
267
Curvature and Curve Tracing 2
( 1) can be written as
l
=
a 2x x _a
-2 - which indicates that curve does not exist 2
when Ixl > a :. The curve exists only in the region -a < x < a. y
x
Fig. 4.11
dy a4 x = and y> 0 in first quadrant. Hence curve is increasing in 0 < x < a . dn y 4.5.6 Example: Trace the curve
l (2 - x) - x2 (1 + x) = 0 ....... ( 1)
Solution: x=Oin (1) gives
y=O
:.(0,0) lies on the curve
F(x,-y)
=
F(x,y) :.curve is
symmetric about x-axis . Lowest degree terms equated to zero gives
2l- x 2 = 0 => y = ±
12
tangents at (0,0) to the curve y
=
0 in (\) gives -1,0
curve passes through (-1,0) in addition to (0,0) Coefficient of highest degree terms of y equated to zero gives 2 - x = 0 (ie) x = 2 is an asymptote parallel to y-axis 2
No asymptote parallel to x-axis as coefficient of x is 1,a constant
Engineering Mathematics - I
268 Substituting Y
x
3
(
= mx + c
in (I) and rearranging
2 2 2 2 2 m + 1) + x ( 1+ 2mc - 2m ) + ( c - 4mc ) x - 2c
¢3 (m) = m 2 + 1 => ¢3 (m) = 0 gives m = ±i :.
=0
There are no oblique asymptotes.
--~----~~--------~----~x ...
~...,
I (2,0) I '" J2 I "'",. I I I -x
y=-
I I
Fig. 4.12 From (I) Y1 = y2
2
+ (2x + 3x h' h h O· 0 < X < 2 were h ) . T IS sows t at YI > 111 t he curve
2y 2-x
is increasing and x = 2 is an asymptote. YI
= 0 => x = 0, x =
5±J5j 4
.
.. In the secon d quadrant y IS maxImum at x
5-£1 =-4
where the tangent IS
parallel to x-axis It can be seen that the curve forms a loop in
[-1,0]
due to symmetry about x-axis.
4.5.7 Example Trace the curve /
-(x-2)(x-4)2 =0
Solution: y = O. (I) gives x = 2, 4 Curve intersects the x-axis at (2,0) and (4,0)
.. , .. (1)
Curvature and Curve Tracing
269
F(x, -y) = F (x,y) :. curve is symmetric about x-axis For x < 2 , /
is -ve :. therefore curve does not exist for x < 2
From (I) it can be seen that x = 2 is a tangent at (2,0) Coefficient of x 3 = l,a constant :. No asymptote oarallel to x-axis Differentiating (1) w.r.t 'x'
dy _ (x-4)(3x-8)
dx At
2y
(4.0) dy
= Lt ± (x-4)(3x-8)
dx
Lt± 3x-8 = ±J2 2.Jx-2
2.J x - 2 ( x - 4)
H4
H4
Equation of tangents at (4,0) are y =
±J2(x -
4)
Further dy --}- 00 as x --}- 2 :. x = 2is a tangent at (2,0)
dx
y
--~--~----~~~~~x
,,
(2,0)
Fig. 4.13
dy = dx
°
gives x = 8/3. Thus the tangents at x = 8/3 are parallel to x-axis. For the
region. 2 < x < 8/3, dy >
dx
°
above x-axis and dy <
dx
° °
extreme values. For the region 8/3
dx
x -axis. :. the curve forms a loop in 2 ~ x ~ 4
below x- axis. At x = 8/3, y takes above x -axis and dy >
dx
°
below
Engineering Mathematics - I
270
4.5.8 Example 2
Trace the curve y ( x + 4a
2
) -
8a 3
=0
.... (1)
Solution:
y
* 0 when x = O.
:. origin does not lie on the curve.
Curve is symmetric about x-axis ,since F
(-x,y)=F (x,y)
y is +ve for all real x :. Curve completely lies above x-axis. Substituting y = mx + c in (I)
m3 + cx 2 + 4a 2 mx + 4a 2 c - 8a 3
¢3 ( m ) = 0 => m = 0 and
=0 ¢2 ( m ) = 0
gives c = 0
dy
(ie) y = 0 is an asymptote to the curve From (I) -
dx In first quadrant dy < 0 and hence y
IS
dx
maximum at
16a 3 x
t
= - ----
(X2 + 4a2
decreasing. dy
dx
o at (0, 2a). y s
(0, 2a) .
----------~--------------~x
Fig. 4.14
4.5.9 Example: Trace the curve y- c cosh x / c = 0
Solution:
x= 0 in (1) gives y = (c/2)[e O +e-D] = r
.... (1)
Curvature and Curve Tracing
271
:. Curve does not pass through (0,0) and intersects y-axis at (0, c). F
(-x,y)= y -(c/2)( e- x / c +e
curve is symmetric about y-axis
It can be seen from (I) that y>o for all real x .
dy= Sill . h -
dx
( x /) C = -1 [x/c e -x -Xlc] 2
which is >0 for x >0 and <
°
for x <0. y
Y=c
(D,c) ------~~------------~x o
Fig. 4.15 Thus the curve is increasing in
°
(0,00), and decreasing in (-00,0).
When x = we have dy = 0. Further y = c for x = 0.
dx
(i.e.) The tangent to the curve at (0, c) is parallel to x-axis.
4.5.10
Example
Trace the curve
(x / a )2/3 + (y /a )2/3 -1 = 0
Solution:
y
:I;
0 when
x =
° :.
It can be observed that F
(0,0) does not lie on the curve.
(-x,y)= F (x,y)= F (x,-y)
The curve is symmetric about both the axes.
.... ( 1)
Engineering Mathematics-I
272
--~---------+----------~----~x
(0, - b)
-
Fig. 4.16
x = 0 gives / = b 2 => Y = ±b y = 0 gives x = ±a. Thus the curve meets x axis. at (-a,o),( a,o) and y axis at (o,b ),( o,-b). For
Ixl> a or \YI >b
(x/a) 2/3 +(y/b)2/3 -1> 0
Equation is not satisfied.
Ixl> a or IYI >b. :. Curve completely lies in the region Ixl ~ a and Iy ~ bl :. Curve does not exist for
y -!!..( a )I/3. In the first quadrant y is decreasing as a bx
From (I) dy
dx
dy < 0 for dx
O 0
Solution: Clearly y F
(x, - y)
= 0, =
0 when x = O. Origin is a double point on the curve.
F ( x, y). Curve is symmetric about x-axis.
.... (1)
Curvature and Curve Tracing
273
dy
Differentiating (I) w.r.t, 'x' we get -
Fx (3a - 2x )
= -
dx
3jJ
2 (a-x) -
dy dx
1(0,0) =
o.
So x axis is a tangent to the curve at origin.
x = a is an asymptote parallel to y-axis. 3
A" coefficient of x is a constant, there is no asymptote parallel to x -axis. For x < 0 and x> a curve does not exist (as y2 is negative) Substituting y
x
3
(
= mx + c
in (1) and rearranging
m 2 + 1) + x 2 ( 2mc - am 2 ) + x ( c 2 - 2ame ) - ac 2
=0
fA (m) = m 2 + 1=> m = ±i . I I x=a
PI
I I I --~~--------.---------~x
(a,O)
Fig. 4.17 Hence there no oblique asymptotes.
3x2 + y2
From (I) YI
= (
)
2 a-x y
For 0 < x < a, YI is positive. Hence y is increasing in 0 < x < a .
. 4.5.12 Example Trace the curve x
2
(a
- x) = l
(b + x)
.. .,(\ )
Engineering Mathematics - I
274
Solution: Clearly curve passes through origin F ( x, - y ) = F
(x, y) .'. Curve is symmetric about x-axis.
Equating lowest degree terms to zero
by2 _ax 2 = 0 => y
= ±~.x
There are two real, distinct tangents to the curve at (0,0)
y
=0
in (1) gives
The curve intersects x-axis at
x2 ( a - x) = 0 => x = 0, a
(a, 0) 3
There is no asymptote parallel to x-axis as coefficient of x is a con y
I I I I
--~I~--------~~--~~~~----~x I I I l X=-b I
Fig. 4.18 Coefficient of highest degree term in y is b + x
:. x + b = 0 is the asymptotes parallel to y-axis Substituting y