Engineering Mathematics
Engineering Mathematics
Calculus
CHAPTER HIGHLIGHTS
Limit oF a Function Now, let us see what happens when x is greater than 2.
Let y = f(x) be a function of x and let ‘a’ be any real number. When x = 2.1, f (x) = 9.261
We must first understand what a ‘limit’ is. A limit is When x = 2.01, f (x) = 8.12
the value, function approaches, as the variable within that When x = 2.001, f (x) = 8.01
function (usually ‘x’) gets nearer and nearer to a particu- When x = 2.0001, f (x) = 8.001
lar value. In other words, when x is very close to a certain As x decreases and approaches 2, f (x) still approaches 8.
number, what is f(x) very close to? This is called ‘right-hand limit’ and is written as lim+
x→ 2
→ ←
Meaning of ‘x → a’ x 2 2 x
Let x be a variable and ‘a’ be a constant. If x assumes values We get the same answer while finding both, left and right
hand limits. Hence we write that lim x = 8.
3
nearer and nearer to ‘a’, then we say that ‘x tends to a’ or ‘x x →2
approaches a’ and is written as ‘x → a’. By x → a, we mean
that x ≠ a and x may approach ‘a’ from left or right, which Meaning of the Symbol: xlim
→a
f(x) = l
is explained in the example given below. Let f (x) be a function of x where x takes values closer and
Let us look at an example of a limit: What is the limit of closer to ‘a’ (≠ a), then f (x) will assume values nearer and
the function f (x) = x3 as x approaches 2? The expression ‘the nearer to l. Hence we say, f (x) tends to the limit ‘l’ as x tends
limit as x approaches to 2’ is written as: lim Let us check to a.
x→2
out some values of lim as x increases and gets closer to 2, The following are some of the simple algebraic rules of
x→2 limits.
without even exactly getting there.
1. xlim kf ( x ) = k lim f ( x )
When x = 1.9, f (x) = 6.859 →a x →a
When x = 1.99, f (x) = 7.88
2. lim[ f ( x ) ± g ( x )] = lim f ( x ) ± lim g ( x )
When x = 1.999, f (x) = 7.988 x →a x →a x →a
When x = 1.9999, f (x) = 7.9988 3. xlim[ f ( x ) ⋅ g ( x )] = lim f ( x ) ⋅ lim g ( x )
→a x →a x →a
As x increases and approaches 2, f (x) gets closer and closer
to 8 and since x tends to 2 from left this is called ‘left-hand f ( x) x→a
4. lim = ( lim g ( x ) ≠ 0)
limit’ and is written as lim− . x →a g ( x ) lim g ( x ) x → a
x→ 2 x →a
Continuous Functions
Let f: A → B be any given function and let c ∈ A. We say f is
continuous at c, if given ∈ > 0, there exists d > 0 such that
|f (x) - f (c)|< ∈ whenever |x - c| < d
(d, m )
In words, this means that, if x is very close to c in domain,
x
then f (x) is very close to f (c) in range. a d c c1 b
Equivalently f is continuous at c. If lim f ( x ) = f (c)
x →c NOTE
We observe
1. c ∈ A, i.e., f (c) must exist 1; 0 < x ≤ 1
The converse may not be true as f ( x ) = is
2. lim f ( x ) exists −1; 1 < x ≤ 2
x →c
bounded on [1, 2] but it is not continuous at x = 1.
3. f (c) and lim f ( x ) are equal.
x →c
If any of these three conditions fail, then f is discontinuous Intermediate-value Theorem
at x = c. If f is continuous on [a, b] and f (a) ≠ f (b) then f takes every
value between f (a) and f (b).
Equivalently, if f is continuous on [a, b] and f (a) < k < f
Algebra of Continuous Functions (b) or f (b) < k < f (a), then there exists c ∈ (a, b) such that
If f, g be two continuous functions at c, then f + g, f - g, fg f(c) = k.
are also continuous at x = c. Equivalently, If f (a) and f (b) are of opposite signs then
To solve a problem of continuous functions at a point a, there exists c ∈ (a, b) such that f(c) = 0.
you can take the following approach. y
NOTES d 1
(c) ( x ) = ; x≠0
dx 2 x
1. dy is the rate of change of y with respect to x. d
dx 2. [axn + b] = an · xn-1
dx
2. If the function y can be represented as a general curve, d
and a tangent is drawn at any point where the tangent 3. [ax + b]n = n a (ax + b) n-1
dx
makes an angle θ with the horizontal (as shown in the d
4. [eax] = a · eax
dy dx
figure), then = tan θ , In other words, derivative of
dx d 1
5. [log x] = ; x > 0
a function at a given point is the slope of the curve at dx x
that point, i.e., tans of the angle, the tangent drawn d x
6. [a ] = ax log a; a > 0
to the curve at that point, makes with the horizontal. dx
d
7.
(a) [sin x] = cos x
Y y = f (x ) dx
d
(b) [cos x] = –sin x
dx
θ
d
(c) [tan x] = sec2 x
dx
O X
d
(d) [cot x] = –cosec2 x
dx
d
Standard Results (e) [sec x] = sec x · tan x
dx
If f (x) and g(x) are two functions of x and k is a constant, d
(f) [cosec x] = –cosec x · cot x
then dx
1.
d
(c) = 0 (c is a constant) Inverse Rule
dx If y = f (x) and its inverse x = f–1(y) is also defined, then
d d dy 1
2. k ⋅ f ( x) = k f ( x ) (k is a constant) = .
dx dx dx dx
dy
d
3. ( f ( x ) ± g ( x ))
dx Second Derivative
d d If y = f (x), then the derivative of derivative of y is called as
= f ( x) ± g( x)
dx dx d2 y
second derivative of y and is represented by .
dx 2
Product Rule d2 y d dy dy
= f ′′( x ) = where
4. d { f ( x ) ⋅ g ( x )} = f ′( x ) ⋅ g ( x ) + f ( x ) ⋅ g ′( x ) dx 2 dx dx dx is the first derivative of y.
dx
d 1
8. (a) sin −1 x =
Quotient Rule dx 1 − x2
d f ( x ) g ( x ) ⋅ f ′( x ) − f ( x ) ⋅ g ′( x ) d −1
5. = (b) cos ec −1 x =
dx g ( x ) ( g ( x )) 2 dx | x | x2 −1
d −1
Chain Rule (c) cos −1 x =
dx 1− x2
6. If y = f (u) and u = g(x) be two functions, then
dy dy du d 1
= × (d) sec −1 x =
dx du dx dx | x | x2 −1
d −1 1 dn nπ
(e) tan x = (f) n sin (ax + b) = an sin + ax + b
dx 1+ x2 dx 2
(f) d cot −1 x = −1 dn nπ
dx 1+ x2 (g) n cos (ax + b) = an cos + ax + b
dx 2
9. (a) d sinh x = cosh x
dx dn
(h) n (eax sin bx)
d dx
(b) cosh x = sinh x b
dx = ( a 2 + b 2 ) n / 2 eax sin bx + n tan-1
a
(c) d tanh x = sech 2 x d n ax
dx (i) (e cos bx)
dx n
d b
(d) coth x = −cosech 2 x = (a2 + b2)n/2 eax cos bx + n tan-1
dx a
d d 1 ( −1) n n+1
n n
(e) sech x = − sech x tanh x
dx (j) n 2 = n + 2 sin θ sin (n + 1)θ
dx x + a 2 a
d
cosechx = − cosechx coth x
where θ = tan-1
(f) x
dx
a
d 1 dn
10. (a) sinh −1 x = (k) n (tan-1x) = (-1)n-1 (n - 1)! sinnθ ⋅ sin nθ
dx dx
1 + x2
where θ = cot-1x.
(b) d cosh −1 x = 1
dx x2 − 1 Application of Derivatives
(c) d tanh −1 x = 1 Errors in Measurement
dx 1 − x2 Problems relating to errors in measurement can be solved
using the concept of derivatives. For example, if we know
(d) d coth −1 x = −1
dx x2 −1 the error in measurement of the radius of a sphere, we can
find out the consequent error in the measurement of the
d −1 volume of the sphere. Without going into further details
(e) sech −1 x =
dx x 1 − x2 of theory, we can say dx = error in measurement of x and
d −1 dy = consequent error in measurement of y, Where y =
(f) cosech −1 x = dy
dx x x2 + 1 f (x). Hence, we can rewrite = f ′( x ) as dy = f ′(x) · dx.
dx
Successive Differentiation Thus, if we know the function y = f (x) and dx, error in
measurement of x, we can find out dy, the error in meas-
If f is differentiable function of x and the derivative f ′ is also
urement of y.
a differentiable function of x, then f ″ is called the second
derivative of f. Similarly 3rd, 4th ... nth derivative of f may NOTES
be defined and are denoted by f″′, f ″′′, ... f n or y3, y4 … yn. 1. An error is taken to be positive when the measured
11. The nth derivatives of some special functions: value is greater than the actual value and negative
when it is less.
dn
(a) n x n = n ! dy
dx 2. Percentage error in y is given by ×100.
y
dn m!
(b) n x m = x m − n s(m being a positive
dx ( m − n)
integer more than n) Rate of Change
d n ax While defining the derivative, we have seen that derivative
(c) n e = an eax is the ‘rate of change’. This can be applied to motion of bod-
dx
ies to determine their velocity and acceleration.
dn 1 ( −1) n n !
(d) n = ; x ≠ −a
dx x − a ( x + a) n +1 Velocity If we have s, the distance covered by a body ex-
pressed as a function of t, i.e., s = f (t), then rate of change
dn ( −1) n −1 ( n − 1)! ds
(e) n log( x + a) = ; (x + a) > 0 of s is called velocity (v). v = = f ′(t ).
dx ( x + a) dt
Acceleration Rate of change of velocity is defined as Another Form If f is defined on [a, a + h] such that
acceleration. Since v = f ′(t) itself is a function of t, we can 1.
f is continuous on [a, a + h].
write v = f ′(t).
2.
f is differentiable on (a, a + h) then there exists atleast
dv d 2 s i.e., acceleration is the second derivative
=a = , one q ∈ (0, 1) such that f (a + h) = f (a) + hf ′(a + qh).
dt dt 2
of the function s = f (t). Meaning of the sign of the derivative
h2
Mean Value Theorems f ( a + h) = f ( a) + hf ′( a) +
2!
f ′′( a) +
Rolle’s Theorem Let f be a function defined on [a, b] such
that hn −1 n −1
+ f ( a) + Rn .
1. f is continuous on [a, b]; ( n − 1)
2. f is differentiable on (a, b) and
Where
3. f (a) = f (b), then there exists c ∈ (a, b) such that f ′(c)
=0 hn f n ( a + θ h)
Rn =
Lagrange’s Mean Value Theorem Let f be a function de- n!
fined on [a, b] such that (Lagranges’ form of remainder)
1. f is continuous on [a, b], hn (1 − θ ) n −1 f n ( a + θ h)
2. f is differentiable on (a, b) then there exists c ∈ (a, b) Rn =
( n − 1)!
f ( b) − f ( a)
such that f ′(c) = .
b−a (Cauchy’s form of remainder)
xn n Example 7
+ f (0) + is known as
n! If f (x) = 2x2 + 3x + 4, then find the value of q in the mean
value theorem.
Maclaurin’s infinite series.
Series expansions of some standard functions Solution
f (a) = 2a2 + 3a + 4
1. x 2 x3 xn
ex = 1+ x + + + + + f (a + h) = 2(a2 + 2ah + h2) + 3a + 3h + 4
2 ! 3! n!
f (a + h) − f (a) = 4ah + 2h2 + 3h = 2(2ah + h2) + 3h
x3 x5 ( −1) n x 2 n +1
sin x = x −
2. + − + + f ( a + h) − f ( a )
3! 5 ! ( 2n + 1)! = 2( 2a + h) + 3
h
x2 x4 ( −1) n x 2 n
3. cos x = 1 − + − + + h
2! 4 ! ( 2n)! = 4 a + + 3 (1)
2
x3 x5 x 2 n +1
sinh x = x +
4. + + + + Now f ′ (x) = 4x + 3, f 1 (a + qh)
3! 5! ( 2n + 1)! = 4a + 4hq + 3 (2)
Comparing Eqs. (1) and (2) we have 4 a + + 3
x 2
x 4
x 2n h
5. cosh x = 1 + + + + +
2! 4 ! ( 2n)! 2
x 2 x3 x 4 ( −1) n −1 x n h
6. log(1 + x ) = x − + − − + = 4 a + 4 hθ + 3 ⇒ a + hθ = a +
2 3 4 n 2
(1 + x)-1 = 1 - x + x2 - x3 + …
7. ⇒ θ=
1
(1 - x)-1 = 1 + x + x2 + x3 + …
8. 2
(1 + x)-2 = 1 - 2x + 3x2 - 4x3 + …
9.
Partial Differentiation
−
1
x 1⋅ 3 2 1⋅ 3 ⋅ 5 3 Let u be a function of two variables x and y. Let us assume
10. (1 − x ) 2
= 1+ + x + ⋅ x +
2 2⋅3 2⋅ 4 ⋅6 the functional relation as u = f (x, y). Here x alone or y alone
or both x and y simultaneously may be varied and in each Similarly, f (x, y) is said to have a local minimum at a
case a change in the value of u will result. Generally the point (a, b), if f (x, y) has least value at (a, b) in a neighbour-
change in the value of u will be different in each of these hood of (a, b).
three cases. Since x and y are independent, x may be sup-
posed to vary when y remains constant or the reverse. Procedure to Obtain Maxima and Minima
The derivative of u wrt x when x varies and y remains Let f (x, y) be a function of two variables for which we need
constant is called the partial derivative of u wrt x and is to find maxima and minima.
∂u
denoted by ∂f ∂f
∂x Find f x =
1. and f y =
∂x ∂y
Take fx = 0 and fy = 0 and solve them as simultaneous
2.
∂ 2 u ∂ ∂u ∂ 2 u ∂ ∂u
= , = . equations to get pairs of values for x and y, which are
∂x 2
∂x ∂x ∂x∂y ∂x ∂y called stationary points.
∂2 f ∂2 f
Total Differential Co-efficient Find r = f xx =
3. , s = f xy = and
∂x 2 ∂x∂y
If u be a continuous function of x and y and if x and y receive
small increments Δx and Δy, u will receive, in turn, a small ∂2 f
t = f yy =
and find rt – s2.
increment Δu. This Δu is called total increment of u. ∂y 2
Δu = f (x + Δx, y + Δy) - f (x, y)
4.
At a stationary point, say (a, b)
In the differential form, this can be written as
(a) If rt – s2 > 0, then (a, b) is called an extreme point
∂u ∂u of f (x, y) at which f (x, y) has either maximum or
du = dx + dy.
∂x ∂y minimum which can be found as follows.
Case 1: If r < 0, then f (x, y) has a local maximum at
du is called the total differential of u. If u = f (x, y, z) then a, b)
Case 2: If r > 0, then f (x, y) has a local minimum
du ∂u dx ∂u dy ∂u dz
= ⋅ + ⋅ +⋅ ⋅ at (a, b).
dt ∂x dt ∂y dt ∂z dt
(b) If rt – s2 < 0, then (a, b) is called as saddle point
of f (x, y) where f (x, y) has neither maximum nor
Implicit Function minimum at (a, b).
If the relation between x and y be given in the form f (x, y) =
c where c is a constant, then the total differential co-efficient Example 8
wrt x is zero.
Find the stationary points of the function f (x, y) = x2y + 3xy
– 7 and classify them into extreme and saddle points.
Homogeneous Functions
Let us consider the function f (x, y) = a0xn + a1xn-1y + Solution
a2xn-2y2 + … + anyn. In this expression the sum of the indi-
Given f (x, y) = x2y + 3xy – 7
ces of the variable x and y in each term is n. Such an expres-
sion is called a homogeneous function of degree n. ∂f ∂f
∴ fx = = 2 xy + 3 y and f y = = x 2 + 3x
∂x ∂y
Euler′s Theorem
Now fx = 0 ⇒ 2xy + 3y = 0 and fy = 0
If f (x, y) is a homogeneous function of degree n, then
⇒ x2 + 3x = 0
∂f ∂f
x +y = nf . −3
∂x ∂y ⇒ y = 0 and x = ; x( x + 3) x = 0 and x = −3
This is known as Euler’s theorem on homogeneous 2
function. 3
But for x = , fy ≠ 0
2
Maxima and Minima for Function \ The stationary points of f (x, y) are (0, 0) and (-3, 0)
of Two Variables Now r = fxx = 2y; s = fxy = 2x + 3 and t = fyy = 0
A function f (x, y) is said to have a local maximum at a point And rt – s2 = 2y × 0 – (2x + 3)2 = -(2x + 3)2
(a, b), if f (a + h, b + k) ≤ f (a, b) for all small values of h and \ rt – s2 < 0 at (0, 0) as well as (-3, 0)
k, i.e., f (x, y) has a local maximum at (a, b), if f (a, b) has a Hence the two stationary points (0, 0) and (-3, 0) are saddle
highest value in a neighbourhood of (a, b). points where f (x, y) has neither maximum nor minimum.
Example 9 ax
5. ∫ a dx = +c
x
1
Indefinite Integrals 18.
∫ 1+ x 2
dx = tan −1 x + c or − cot −1 x + c
If f (x) and g(x) are two functions of x such that g′(x) = f (x),
then the integral of f (x) is g(x). Further, g(x) is called the 1
19. ∫ dx = sec −1 x + c or − cos ec −1 x + c
antiderivative of f (x). x x −1 2
dx x 3.
If f (x) is continuous function of x over [a, b], then
30. ∫ = sin −1 +c
a −x
2 2 a a b
∫b
f ( x )dx = − ∫ f ( x )dx.
a
dx x
31. ∫ = sin h−1 + c or 4.
If f (x) is continuous in some neighbourhood of a,
a +x
2 2 a a
then ∫ f ( x )dx = 0.
a
log | x + a 2 + x 2 | + c
dx x 5.
If f (x) and g(x) are continuous in [a, b], then
32. ∫ = cos h−1 + c or b b b
x −a
2 2 a ∫a
[ f ( x ) + g ( x )]dx = ∫ f ( x )dx + ∫ g ( x )dx.
a a
log | x + x 2 + a 2 | + c b b b
6. f ( x )dx = f ( z )dz = f (t )dt
1 1 x
∫ a ∫ ∫ a a
33. ∫ 2 dx = tan −1 + c a a
x + a2 a a ∫ ∫
7. f ( x )dx = f ( a − x )dx
0 0
1 1 x−a
34. ∫ 2 dx = log +c a
8. f ( x ) = 0, if f ( x ) is odd
x −a 2
2a x+a ∫ −a
1 1 a+ x a a
35. ∫
a −x
2 2
dx =
2a
log
a−x
+c 9. ∫− a f ( x)dx = 2∫0 f ( x)dx if f ( x) is even
2a a
x a2 + x 2 a2 x 10. ∫0 f ( x )dx = 2 ∫0 f ( x )dx, if f ( 2a − x ) = f ( x )
36.
∫ a − x dx = + sin −1 + c
2 2
2 2 a
= 0 if f (2a - x) = - f (x)
x a2 − x 2 a2 x
37.
∫ a + x dx = + sin h−1 + c
2 2 na a
2 2 a 11. ∫0 f ( x )dx = n ∫0 f ( x )dx, if f ( a + x ) = f ( x )
38. x x 2 − a2 a2 x
∫ x 2 − a 2 dx =
2
− cos h−1 + c
2 a Applications of Integration
x Area as a Definite Integral
39. ∫ log x dx = x(log x − 1) = x log + c 1. The area enclosed by a curve y = f (x), the lines x = a
e and x = b and the x-axis is given by:
40. ∫ e [ f ( x ) + f ′( x )] dx = e f ( x ) + c
x x
b f ( x )dx, if f ( x ) ≥ 0, a ≤ x ≤ b
b ∫a
Definite Integrals A = ∫ | f ( x ) | dx = b
a
− ∫ f ( x )dx, if f ( x ) ≤ 0, a ≤ x ≤ b
The difference in the values of an integral of a function f (x) a
for two assigned values say a, b of the independent variable
x, is called the Definite Integral of f (x) over the interval [a, y
b y = f (x )
b] and is denoted by ∫a f ( x )dx.
The number ‘a’ is called the lower limit and the number
‘b’ is the upper limit of integration.
A
Fundamental Theorem of Integral
o x=a x=b
Calculus x
2. If f(x) is continuous function of x over [a, b], then the area enclosed by the curve y = f (x), the lines x = a
b b
∫ Kf ( x )dx = K ∫ f ( x )dx. c b
a a
and x = b and the x-axis is A = ∫ f ( x ) dx − ∫ f ( x ) dx
a c
y x = f (x ) 3
= 2
23
A
(c, 0) x=b 3
=
o x=a x
3
4
Example 11
Find the area enclosed by the curve y = x2 and line y = 4?
4. The area enclosed by the curves y = f (x) and y = g(x)
and the lines x = a and x = b is given by,
Solution
y f (x )
The area enclosed by the curve y = x2 and the line y = 4 is
A
the region OAB.
\ The region OAB is bounded by line y = 4 and the curve y
g(x ) = x2 from x = -2 to x = 2 and 4 ≥ x2 for all x ∈ [-2, 2]
o x=a x=b x
y
y g(x ) B (2, 4)
( −2, 4 ) A y=4
A A
y = x2
f (x ) o x
o x=a x=b x
2
b ( f ( x ) − g ( x ))dx, if f ( x ) ≥ g ( x ), \ The required area = ∫ ( 4 − x 2 ) dx
∫a
x =−2
a ≤ x ≤ b 2
= 2 ∫ ( 4 − x 2 )dx (∵ 4 − x 2 as even )
b
A = ∫ | f ( x ) − g ( x )| dx = b
a
∫ ( g ( x ) − f ( x ))dx, if f ( x ) ≤ g ( x ); 0
a 2
a ≤ x ≤ b x3 32
= 2 4 x − =
3 0 3
Example 10
Find the area enclosed by the curve y = x3, the line y = 2 and Rectification
the y-axis in first quadrant? The process of determining the length of arcs of plane
curves is called Rectification. The length of the arc can be
Solution
calculated by any one of the methods given below.
The area bounded by y = x3, y = 2 and the y-axis is the area
OAB as shown in the figure. Cartesian Equations Let y = f(x) be a function of x. The
1 length of arc between the points with x-coordinates ‘a’ and
So, the region OAB is bounded by the curve x = y , the 3
2
b dy dy
‘b’ is given by S =
∫ 1 + dx, provided is con-
1
y=2 B NOTE
A
If the equation of the curve is given in the form x =
o
f (y), then the length of the arc between the points with
x y-coordinates ‘c’ and ‘d’ is given by
y =x3
2
\ The required area d dx dx
S=∫ 1 + dy provided is continuous on [c, d]
2 c
dy dy
2 1 3 4
= ∫ y 3 dy = y 3
y =0 4 0
Parametric Equations Let x = f (t) and y = g(t) be paramet-
3 4 ric functions of ‘t’.The length of the arc between the points
= × 23
4 {f (t1), g (t1)} and {f (t2), g (t2)} is given by
dx 2 dy 2 dx dy π
t2
g ( x ) ≥ 0 in 0,
∫t1 dt + dt dt provided dt and dt are both 2
\ By first mean value theorem,
continuous on [t1, t2].
π
π
Polar Equations Let r = f (θ) be a function of θ, the length ∫0
2
x cos dx = µ ∫ cos x dx = µ
0
of the arc between the points {f(θ1), θ1} and {f(θ2),θ2} is
π π
2
\There exists is µ ∈ 0, such that
given by S = ∫
θ2 dr
r2 + dθ provided
dr
is continu- 2 ∫
0
2
cos x dx = µ
θ1
dθ dθ
ous along the arc. Example 13
If the equation of the curve is given in the form θ = f(r), Verify second mean value theorem for f(x) = x2 and g(x) =
then the length of the arc between the points (r1, f (r1)), (r2, x2 on [-1, 1].
f (r2)) is given by
Solution
r2 dθ
2
dθ Given f(x) = x2 and g(x) = x2 on [-1, 1] both f and g are
S=∫ 1+ r2 dr provided is continuous along
r1
dr dr continuous and integrable on [-1, 1] but g is a decreasing
the arc. function on [-1, 0] and increasing function on [0, 1] \ g is
not monotonic.
Theorems on Integration 1 1
\ ∫−1 f ( x ) g ( x )dx = ∫−1 x ⋅ x dx
2 2
1. If f is a continuous function on [a, b] then there exists
b
c ∈ (a, b) such that ∫a f ( x ) dx = f (c)(b − a)
1
x5 1 1 2
2. If f, g ∈ R [a, b] and g keeps the same sign on [a, b] = = + = (1)
5 −1 5 5 5
then there exists m ∈ R lying between the infimum and
b b But by second mean value theorem,
the suprimum of f such that ∫ f ( x ) g ( x ) = µ ∫ g ( x )dx b µ b
∫ f ( x ) g ( x ) dx = g ( a) ∫ f ( x )dx + g (b) ∫ f ( x ) dx
a a
NOTE a a µ
1 µ 1
This is called the first mean value theorem. \ ∫−1 x dx = g ( −1) ∫−1 x dx + g (1) ∫µ x dx
4 2 2
(Singularity at upper limit)
Solution
b b
Take f(x) = x and g(x) = cos x 2. ∫−∞
f ( x )dx = lim
a →− ∞ a ∫ f ( x )dx.
π and g is integrable on π ∞ lim
\ f is continuous on b
also
0, 2
0, 2
3. ∫−∞
f ( x )dx = b → ∞ ∫ f ( x )dx. Or
a →− ∞ a
∞ 0 b Example 14
4. ∫− ∞ f ( x )dx = lim ∫ f ( x )dx + lim ∫ f ( x )dx.
a →− ∞ a b →∞ 0
∞ dx
Or = lim ∫ f ( x )dx.
r Examine ∫
1 xp
for convergence/divergence.
r →∞ − r
If these limits exist then it is convergent otherwise it is Case 2: If the limits of x are function of y, say x1 = g1(y) and
divergent. x2 = g2(y) and the limits of y are constants, say y1 = c and y2
Improper Integral of Third Kind If the limits of the inte- = d, then integrate wrt x first treating y as constant and then
gral are infinite or f (x) may be discontinuous or both then integrate wrt y.
the improper integral is known as third kind.
x2 = g 2 ( y )
NOTES y2 = d
That is, ∫ ∫ f ( x, y )dR = ∫ ∫ f ( x, y )dx dy
R y1 = c
∞1
1. ∫ p dx is convergent when p > 1 and it is divergent x = g ( y)
1 x
1 1
when p ≤ 1. This result is used in comparison test for Case 3: If both the variables x and y have constant limits,
testing the convergence or divergence of the integral then one can follow any order of integration.
of first kind.
c 1 Change of Order of Integration Evaluation of some of the
2. ∫ dx is convergent for p < 1 and is divergent double integrals can be made simple by changing the order
a ( x − c) p
of integration. In change of order of integration, we take the
for p ≤ 1. This is used for convergence or divergence limits of the variables for the given region of integration in
of an improper integral of second kind. such a way that the order of integration reverses.
y x −x
2
1
∴ I = ∫ ∫ e y dx dy (2) (b) The volume of the solid of revolution obtained by
y =0 x = y y
revolving the area A about y-axis is
x2 2x Volume = V = ∫∫ 2π x dx dy
Put =t ⇒ dx = dt
y y A
x =2 ∫
6
y =4
2π ydxdy = ( ∫ dx )( ∫
5
x =2
6
y =4
2π y dy )
\ In D, y varies from y = − 4 − x 2 to y = 4 − x 2 and x ( )( )
= x ]5x = 2 π y 2 ]6y = 4 = 3 × 20π = 60π
varies from x = -2 to x = +2.
\ The volume under the surface x + 2y + z = 4 and above
the circle x2 + y2 = 4 in xy-plane is Change of Variables Evaluation of some of the double
(or) triple integrals can be made simple by changing the
2 4 − x2
V =∫ ∫ zdxdy = ∫ ∫ ( 4 − x − 2 y )dxdy (2) variables.
D x =−2 y =−2 4 − x 2
1. In a double integral: Let a double integral ∫Rxy∫∫
Evaluation of this double integral can be made simple by f ( x, y )dxdy in x and y is to be converted into the
changing it into polar coordinates.
In polar coordinates, x = r cos θ, y = r sin θ and variables u and J where x = f(u, J) and y = Ψ (u, J).
Then
∂x ∂x
∂( x, y ) ∂r ∂θ cos θ − r sin θ
∫ ∫ f ( x, y)dxdy = ∫ ∫
Rxy R1uϑ
f (φ (u, ϑ ), ψ (u, ϑ )) | J | dudϑ
J= = =
∂( r , θ ) ∂y ∂y sin θ r cos θ ∂x ∂x
∂r ∂θ ∂( x, y ) ∂u ∂ϑ
Where J= =
\ J = r. Also, in the circle x2 + y2 = 4, r varies from r = 0 to ∂(u, ϑ ) ∂y ∂y
r = 2 and θ varies from θ = 0 to θ = 2p ∂u ∂ϑ
\ From (2), Is the Jacobian of x and y wrt u and J and Ru′ϑ is the
region of integration in u, J-plane corresponding to
V = ∫ ∫ ( 4 − x − 2 y )dx dy Rxy in xy-plane.
D
2π
=∫
2
2. In a triple integral: Let a triple integral ∫Rxyz∫∫
∫
θ =0 r =0
( 4 − r cos θ − 2r sin θ ) | J | drdθ
2π 2
f ( x, y, z ) dxdydz in x, y and z is to be converted into
=∫ ∫ ( 4 − r cos θ − 2r sin θ )rdr dθ the variables u, J and w, where x = f(u, J, ω), y = Ψ
θ =0 r =0
(u, J, ω) and z = h (u, J, ω)
=∫
2π
θ =0 (∫ r =0
2
( 4 r − r 2 cos θ − 2r 2 sin θ )dr dθ ) Then ∫ ∫ ∫ f ( x, y, z )dxdydz = ∫ ∫ ∫ f (φ (u, ϑ , ω ),
Rxyz R1
uϑω
2
2π r 2r 3 3
ψ (u, ϑ , ω ), h(u, ϑ , ω )) | J | dudϑ dω
= ∫ 2r 2 − cos θ − sin θ dθ
θ = 0
3 3 r =0 ∂x ∂x ∂x
2π 8 16 ∂u ∂ϑ ∂ω
=∫ 8 − 3 cos θ − 3 sin θ dθ ∂( x, y, z ) ∂y ∂y ∂y
θ =0
where J =
= is the
∂(u, ϑ , ω ) ∂u ∂ϑ ∂ω
2π
8 16 ∂z ∂z ∂z
= 8θ − sin θ + cos θ = 16π
3 3 θ =0 ∂u ∂ϑ ∂ω
Multiplication of Vectors | a ×b |
4. The angle between two vectors: sin θ =
Scalar or Dot Product If a and b are two non-zero vectors | a || b |
and q is the angle between them (0 ≤ q ≤ p), then their dot or 5. A unit vector perpendicular to the plane of a and b
scalar product is given by a ⋅ b =| a || b | cos θ ⋅ a ⋅ b is a scalar. a ×b
is given by n where n =
| a ×b |
NOTES 6. Area of parallelogram whose adjacent sides are
1. If one or both of a , b , are 0 , then a ⋅ b = 0 a and b is given by | a × b |
7. When the diagonals are given, the vector area of par-
2. a ⋅ b = | a | ⋅ (scalar component of b along a )
1
allelogram ABCD is ( AB × AC )
= |b | (scalar component of a along b ), 2
3. a ⋅ b = b ⋅ a 1
8. The vector area of the triangle ABC = ( AB × AC )
4. If a , b , c are any three vectors, then a ⋅ (b + c ) 2
= a ⋅b + a ⋅c 9. If a = a1 i + a2 j + a3 k and b = b1 i + b2 j + b3 k ,
5. Two non-zero vectors a and b are perpendicular if
i j k
a ⋅b = 0
Then a × b = a1 a2 a3
6. i ⋅ j = j ⋅ i = j ⋅ k = k ⋅ j = j ⋅ k = k ⋅ i = 0
b1 b2 b3
7. a ⋅ b is positive, negative or zero according as 0 ≤ q
< 90°, 90° < q ≤ 180° or q = 90° 10. Vector product is distributive with respect to vector
addition a × (b + c ) = a × b + a × c
If a = a1 i + a2 j + a3 k , b = b1 i + b2 j + b3 k ,
4. d dG dF
(F ⋅ G) = F ⋅ + ⋅ G,
dt
dt dt
a1 a2 a3
c = c1 i + c2 j + c3 k , then abc = b1 b2 b3 d dG dF
( F × G) = F × + × G.
c1 c2 c3 dt
dt dt
5.
The volume of a tetrahedron with co-terminus edges 6.
Partial derivatives: If F is vector function dependent
1 on x, y and z, say F = F ( x, y, z ), then partial
a , b , c is [abc ] cubic units. derivative of F with respect to x is defined as
6
∂F F ( x + ∆x, y, z ) − F ( x, y, z )
[abc ] = ( a × b ) ⋅ c = a ⋅ (b × c )
6. = lim .
∂x ∆x → 0 ∆x
Vector Triple Product If a , b , c are three vectors, then the ∂F ∂F
Likewise, one can also define and .
triple product a × (b × c ) is called the vector triple product. ∂y ∂z
If a , b , c are any three vectors, then a × (b × c ) It is also possible to define higher order partial
derivatives as:
= ( a ⋅ c )b − ( a ⋅ b )c
∂2 F ∂ ∂F ∂ 2 F ∂ ∂F
= , = .
Vector Variable ∂x 2
∂x ∂x ∂y 2
∂y ∂y
A variable of the form r = xi + y j + zk is called a vector ∂2 F ∂ ∂F
variable and x, y, z are scalar variables. = , etc
∂x∂z ∂x ∂z
Scalar Function If t is a scalar variable on a range a ≤ t ≤ b
and a function f defined as f = f (t) for t ∈ [a, b] is called a Differential Vectors
scalar function of t.
If G = G ( x, y, z ) then
1.
Example: f (t) = 9t3 + 4t2 + 7,
∂G ∂G ∂G
f (t) = sint + 5cost + et, etc. dG =
dx + dy + dz
∂x ∂y ∂z
Vector Function If t is a scalar variable defined on a do-
main [a, b], and a function F (t ) = x(t )i + y(t ) j + z (t )k is If F = F1 i + F2 j + F3 k , then
2.
called a vector function of the scalar variable t.
dF = dF1 i + dF2 j + dF3 k
NOTE 3. d ( F ⋅ G ) = F ⋅ dG + dF ⋅ G
t is generally taken as ‘time’.
4. d ( F × G ) = F × dG + dF × G
Differentiation If F (t ) is a continuous single valued vec-
tor function of the variable t, then the derivative of F (t ) is Vector Differential Operators ∇ is to be read as del or nabla
dF F (t + ∆t ) − F (t ) ∂ ∂ ∂
defined as = lim where Dt is a small
dt ∆t →0 ∆t ∇ = i + j + k
∂x ∂y ∂z
increment in t.
One can also look at second and higher order derivatives ∂2 ∂2 ∂2
∇2 = + + is called Laplacian.
in a similar way. ∂x 2 ∂y 2 ∂z 2
NOTE
C2
If div ( F ) or ∇ ⋅ F = 0, then F is called ‘solenoidal’ C3
x
0
Curl of a Vector C1
Let F ( x, y, z ) is a vector field defined for all (x, y, z) in ∫c F ⋅ dr = ∫C ( xydx + y 2 dy) + ∫C ( xydx + y 2 dy)
a certain region of space and is differentiable, i.e., F is a 1 2
i j C1 C2 C3
k
y=0 y=1-x x=0
∂ ∂ ∂
curl F = ; F = F1 i + F2 j + F3 k 0<x<1 1<x<0 dx = 0
∂x ∂y ∂z
dy = 0 dy = - dx 1<y<0
F1 F2 F3
1 0 0
NOTE =∫ [ x(0)dx + 0 + ∫ x(1 − x )dx + (1 − x ) 2 ( −dx ) + ∫ y 2 dy
x =0 x =1 1
Standard Results 0 1
= ∫ ( −2 x 2 + 3 x − 1)dx − ∫ y 2 dy
1. div (f F) = f div F + F ⋅grad f or ∇ . f F = f ∇ ⋅ F + 1 0
F ⋅ ∇ f 2 3 1 −1
= − + 1 − =
2. curl (f F) = ∇f × F + f curl F 3 2 3 6
Cartesian Form Let F ( r ) = F1 i + F2 j + F3 k , where, F1, F2, B = − xi + x 2 j + x 3 k , then find the values of
F3, are continuous and differentiable functions of x, y, z. d d
(i) ( A ⋅ B ) and (ii) ( A × B ).
If cosα, cosb and cosg be the direction cosines of the unit dx dx
normal N, then
Solution
N = i cos α + j cos β + k cos γ . d d d
(i) ( A ⋅ B) = A ⋅ ( B) + B ⋅ ( A)
dx dx dx
∴ ∫S F ⋅ N ds = ∫S ( F1 cos α + F2 cos β + F3 cos γ )ds
d
= ( x 3 i + x 2 j + xk ) ⋅ ( − xi + x 2 j + x 3 k )
But then ds cosa, ds cosb and ds cosg are the projections dx
of ds on yz, zx and xy planes. If dx, dy, dz are the differentials d
along the areas then + ( − xi + x 2 j + x 3 k ) ⋅ ( x 3 i + x 2 j + xk )
dx
ds cosa = dy dz; ds cosb = dz dx; ds cosg = dx dy.
∴ = ( x i + x j + xk ) ⋅ ( − i + 2 x j + 3 x 2 k )
3 2
∫S
F ⋅ N ds = ∫ ∫ ( F1dy dz + F2 dz dx + F3 dx dy )
S
+ ( − xi + x 2 j + x 3 k ) ⋅ (3 x 2 i + 2 x j + k )
NOTE = - x3 + 2x3 + 3x3 - 3x3 + 2x3 + x3 = 4x3.
If R1 is the projection of S on xy-plane, then d dB dA
dxdy (ii) ( A × B) = A × + ×B
dx dx dx
∫S F ⋅ Nds = ∫R1 ∫ F ⋅ N cos γ
d
= ( x 3 i + x 2 j + xk ) × ( − xi + x 2 j + x 3 k )
dxdy
= ∫∫ F ⋅N ⋅ (| N ⋅ k | = cos γ ) dx
S
|N ⋅k| d
+ ( x 3 i + x 2 j + xk ) × ( − xi + x 2 j + x 3 k )
Equivalently, dx
dydz dzdx = ( x 3 i + x 2 j + xk ) × ( − i + 2 x j + 3 x 2 k )
∫∫ F ⋅ Nds = ∫ ∫ F ⋅ N
| N ⋅ i |
= ∫∫ F ⋅N
| N ⋅ j | + (3 x 2 i + 2 x j + k ) × ( − xi + x 2 j + x 3 k )
S R2 R3
i j k i j k
Volume Integral
x2 y2 z2 = x3 x 2 x + 3x 2 2x 1
∫x1 ∫y1 ∫z1 f ( x, y, z )dz dy dx
−1 2 x 3 x 2 −x x2 x3
x2 y
2 z2 = i(5 x 4 − 3 x 2 ) − j (6 x 5 + 2 x ) + k (5 x 4 + 3 x 2 )
= ∫ ∫ ∫ f (x, y, z) dz dy dx
x1 y1 z1 Example 20
If f = x3 - 6xy² - 9xyz is a scalar function, then find
Gauss’ Divergence Theorem ∂2 f ∂2 f
, .
∂x 2 ∂x∂y
If F is continuously differentiable vector function in the
region bounded by a surface S, then ∫ ∫ F ⋅ N ds =∫∫∫ div Fdv Solution
S V
where N is the unit normal to the surface. f = x3 - 6xy² - 9xyz
∂f
Green’s Theorem If P and Q are scalar point functions, ∴ = 3 x 2 − 6 y 2 − 9 yz
possessing continuous derivatives of the first order, in a re- ∂x
gion S of the xy plane bounded by a closed curve C then ∂2 f
∴ = 6x − 0 = 6x
∂Q ∂P ∂x 2
∫C Pdx + Qdy = ∫ ∫S ∂x − ∂y dxdy. ∂f
= −12 xy − 9 xz
Stoke’s Theorem If S is an open surface bounded by a ∂y
closed curve C and F is a continuously differentiable vec- ∂2 f ∂ ∂f ∂
= = ( −12 xy − 9 xz )
tor point function, then ∫ F ⋅ d r = ∫ curl F ⋅ Nds, where N ∂x∂y ∂x ∂y ∂x
C S
is unit outward drawn normal at any point on the surface. = - 12y - 9z.
So, n1 ⋅ n2 = | n1 || n2 | cos θ ⇒ 6 − 16 + 4 2
2 −Z 3 −8
∴ ∫ F ⋅ dr = ∫ − z 2 dz = =
= ( 1 + 16 + 4 )( 36 + 16 + 4 ) cos θ PQ 0
3 0 3
−6 −3 3 (3) Along QR:
∴ cos θ = = =
21 56 7 6 7 6 y = 2, dy = 0 and x changes from 3 to 0
0
x2
0
−27
\ θ = cos 3
−1 ∴ ∫ F ⋅ dr = ∫ 3 x dx = 3 =
7 6
QR
3 2 3 2
Solution y
3 3x 2 27
3
( x 3 − x 4 )dx
∫ ⋅ = ∫0 = =
x =0 x =0
F dr 3 x dx
2 0 2
x 4 x5 1 1 1
1
= − = − =
Z 2 4 5
(0, 2) 0 4 5 20
R Q(3, 2)
Example 29
O P X By applying Gauss theorem, evaluate ∫ ∫ (x 3 dy dz + x 2 ydz
(3, 0) S
Solution π
5a 4 b sin 4θ 2
We have = θ−
F1 = x3; F2 = x²y; F3 = x²z 2 4 0
2 ∴ ∫ (∇ × F ) ⋅ ( N ⋅ k )ds
Lower limit: x = 0 ⇒ a sinq = 0 ⇒ q = 0]
π
= 20b ∫ 2 a 2 sin 2 θ a 2 (1 − sin 2 θ ) a cos θ dθ
= ∫ ( − 4 ky ) ⋅ N ds = ∫ − 4 y ( N ⋅ k ) ds ∫ ∫ − 4 y dx dy
R
0
π Since N ⋅ k ds = dx dy
= 20 a 4 b ∫ 2 sin 2 θ cos 2 θ dθ And R is the region bounded by the rectangle.
0
π b
1 2 − 4 y2
= 20 a b ∫
a b a
sin 2θ dθ =∫
x =− a ∫ y = 0 ∫− a 2 dx
4 2
( − 4 y ) dy dx =
0 4
0
π
1 − cos 4θ
= 5a 4 b ∫ 2 dθ
a
= −2 ∫ (b 2 − 0)dx = −2b 2 [ x ]−a a = − 4 ab 2 .
0
2 −a
Exercises
3. Evaluate lim (x − [x]), where [x] is the greatest inte-
1. lim {3 x − 9 x 2 − x } = ______. x→2.7
x →∞ ger less than equal to x.
1 (A) −0.3 (B) 0.7
(A) (B) 3 (C) 4.7 (D) 2
6
1
(C) 6 (D) None of these 4. Evaluate lim 189 .
x →0 x
24 cos x − 24 + 12 x 2 − x 4 (A) 0 (B) ∞
2. lim = (C) −∞ (D) None of these
x →0 24 x 6
1/ x
1 −1 2 x + 3x
(A) (B) 5. lim =
720 120 x →0 2
(A) 1 (B) 3
1 −1
(C) (D) 2
(C) 6 (D)
120 720
23. Let f (x) = eax and g(x) = e−ax be two functions defined 31. ∫ sec3 x dx = ______.
in [ p, q], If the functions satisfies Cauchy mean value sec x tan x
theorem then the value of ‘c’ is ______. (A) + log (sec x + an x )
3
p+q
p + q (B)
(A) sec 2 x tan x 1 π
2 (B) + log tan + x
3 3 4
(C) 2( p + q) (D) None of these
sec x tan x 1 π x
∂z (C) + log tan +
24. If x = cos (z + y2), then = 2 2 4 2
∂y
(D) None of these
(A) 1 (B) y
π /2
(C) 2y (D) -2y
6
32. ∫ sin 4 x cos6 xdx = ______.
0
4 x + 4 y ∂u ∂u
2 5. If u = 6 , then x + y = 3π 2π
x + y
6 ∂x ∂x (A) (B) ⋅ ⋅ c
128 425
u 4
(A) (B) 3π 3π
2 u (C) (D) ⋅ ⋅
2560 512
(C) 4u (D) 6u
33. Area bounded by the curve y = −3x2, x = 2 and the two
26. The stationary points of the function f (x, y) = x3 + y4 − coordinate axes is ______ sq units
27x + 32y + 100 is/are (A) 2 (B) 3
(A) (3, 2), (3, −2) (C) 6 (D) 8
(B) (−3, 2), (−3, −2) 34. The volume of the solid obtained by revolving the area
(C) (3, 2), (−3, −2) bounded by the parabola y2 = x - 4, x-axis and the lines
(D) (3, −2), (−3, −2) x = 4 and x = 7, about x-axis is ______ cubic unit
27. For the function f (x, y) = 2x2 + 4y2 + 4xy + 2x + 9 11
(A) π (B) π ⋅ ⋅
10y + 7. 2 2
(A) Local maximum exists, but no local minimum. 13 15
(B) Local minimum exists, but no local maximum. (C) π (D) π ⋅ ⋅
2 2
(C) Neither local minimum nor local maximum exists.
35. The length of arc of the curve y = ln (cos x) from x = 0
(D) Both local minimum and local maximum exists.
π
28. For the function xyz, if x + y + z = 3, then the local to x = is ____.
⋅
4
maximum occurs for xyz at the point ______.
(A) ln (1 + 2 ) (B) ln ( 2 − 1)
1 1
(A) 4, 2 , 2
(C) ln ( 2 + 3 ) (D) ln ( 2 − 3 )
(B) (5, -1, -1)
(C) (1, 1, 1) π /4 π /4
(D) (7, -3, -1) 36. Evaluate ∫ ∫ (3 cos θ + 4 sin θ )dθ dφ ______.
0 0
29. The ratio of the dimensions of a rectangular box of vol-
ume 64 cubic units and open at the top that requires 2 −1 ( 4 2 − 1)π
(A) π (B)
least material for its construction is 2 4 2
(A) 2 : 2 : 1 (B) 2 : 4 : 5
(C) 2 : 3 : 4 (D) 1 : 2 : 3 ( 4 2 − 1)π ( 4 2 − 1)π
(C) (D)
2 4 2
∞∞ ∞ 1
along the curve x = 3 sin 4t, y = 3 cos 4t, and z = 5t
(t > 0). The magnitude of its velocity and acceleration
∫ ∫ f ( x, y)dxdy (D)
(C) ∫ ∫ f ( x, y)dxdy respectively are
0 0 1 y+2
(A) 13, 45 (B) 12, 48
39. By changing the variables in the double integral (C) 13, 48 (D) 12, 45
dxdy
∫ ∫ xy , where x = eu+v and y = uv, it changes to
⋅
df
45. f (t ) be a vector function and f × = 0 implies
R
dt
∫ ∫ φ (uv)dudν then f(u, v) is (A) f is a vector function with constant magnitude.
R (B) f is a vector function both in direction and
R eu + v magnitude.
(A) u + v ) (B)
(e (uv ) uv (C) f is a vector function of constant direction.
(D) Either A or C.
1 1 1 1
(C) + (D) − 46. The directional derivative of f = x3 y + y3 z + z3 x in the
v u v u
direction of i + 2 j + 2 k at (0, 1, -1) is
40. By changing the variables from x, y to u, v where x = u
5 4
+ 2v and y = 4u + 3v, the given integral ∫ ∫ f ( x, y )dxdy ⋅ (A) (B)
R
3 3
changes to ∫∫ f (u + 2v, 4u + 3v )ψ (u, v )dud v then Y −4 −5
(C) (D)
R 3 3
(u, v) is ______.
(A) 5 (B) -5 47. If r = x iˆ + y ĵ + z k̂ and | r | = r, then ∇ rn =
1 1
(C) (D) − (A) n(n - 1)rn-1 r (B)
n(n - 2)rn-2 r
5 5
n ⋅ rn-2 × r (D)
(C) n(n - 1) r
41. The area bounded by the circle x2 + y2 = 6 and the
parabola y = x2 is given by:
2 x 2 −6
Direction for questions 48 and 49:
(A) ∫ ∫ dydx Two equations f = x y2 z - 2y + z2 and g = x2 + yz - x - 2
x =−2 y = x represents two surfaces
2 6− x2 48. Find normal vector to ‘g’ at (1, -1, 2)
(B) ∫ ∫ dydx (A) i + 2j + 2k (B) i + 2j - k
x =− 2 y = x
(C) 2 i - j - k (D) i - j - 2k
2 x 2 −6
49. The acute angle between the surfaces f and g at (1, -1,
(C) ∫ ∫ ( x 2 + y 2 )dydx 2) is
x =−2 y = x
15 15
2 6− x2 (A) cos–1 (B) cos–1
(D) ∫ ∫ ( y − x 2 )dxdy 390 390
x =− 2 y = x (C) 60° (D) 30°
∇ × ( r × a ) = −2a
(C) C
C is the triangle with vertices (0, 0), (2, 0) and (2, 1).
(D) All of these
11 11
5 5. Compute the value of div(∇f × ∇f ). (A) (B)
24 12
(A) ∇f curl (∇f) (B) ∇f curl (∇f)
−11 11
(C) curl (∇f × ∇f ) (D) 0 (C) (D)
6 4
56. For what value of p the vector f = (2x + 3y)i + (z + 2y) j
+ (x - pz) k is solenoidal? 64. Find the area of the region in the first quadrant bounded
(A) 4 (B) -4 1 x
by the curves y = 4x, y = and y = using green’s
(C) 2 (D) 0 x 4
theorem.
57. For what values of p, q and r the vector f = (x + ry - z) 1
(A) log 2 (B) log 2
i + (3x - y + qz)j + (px + y - z) k is irrotational? 2
(A) p =1, q = -1, r = 3 (C) log 4 (D) log 16
p = -1, q = 1, r = 3
(B)
(C) p = -1, q = 1, r = -3
⋅
∫∫
65. Evaluate F nds where F = 2xziˆ - yz ĵ + yxk where S
s
(D) p =1, q =1, r = -3 is the cube bounded by x = 0, x = 3, y = 0, y = 3 and z =
0, z = 3.
58. If ∇f = yziˆ + zx ĵ + xy k̂, then f(x, y, z) =
(A) xyz + f (y, z); f ≠ constant 27 81
(A) (B)
(B) xyz + g(x, z); g ≠ constant 2 4
(C) xyz + h(x, y); h ≠ constant 27 81
(C) (D)
(D) xyz + k; k is a constant 4 2
59. If F = (5xy - 6x2)iˆ + (2y - 4x) ĵ, compute the line inte- 66. For the force field F = x2 i + xyj in the square region in
gral ∫ F d r where C ≡ y = x3 in the xy-plane joining
⋅
the xy-plane bounded by the lines x = 0, y = 0, x = 2, y
C = 2. Using stokes theorem, find the value of ∫ F , dr . ⋅
2
3
π π
(B) (D)
10
⋅ ⋅
2x
Q 1
16. lim 1 + is equal to [GATE, 2015]
x →∞ x
(A) e-2 (B) e
R
(C) 1 (D) e2
17. While minimizing the function f (x), necessary and
P sufficient conditions for a point, x0 to be a minima
are: [GATE, 2015]
(A) f ′(x0) > 0 and f ″(x0) = 0
O (B) f ′(x0) < 0 and f ″(x0) = 0
(C) f ′(x0) = 0 and f ″(x0) < 0
[GATE, 2012]
(D) f ′(x0) = 0 and f ″(x0) > 0
ad - bc
(A)
18. The directional derivative of the field u(x, y, z) = x2 -
ac + bd
(B)
3yz in the direction for the vector ( iˆ + ˆj − 2kˆ ) at point
ad + bc
(C)
ab - cd
(D) (2, -1, 4) is _______. [GATE, 2015]
19. The optimum value of the function f(x) = x2 - 4x + 2
11. There is no value of x that can simultaneously sat-
is [GATE, 2016]
isfy both the given equations. Therefore, find the least
(A) 2 (maximum) (B) 2 (minimum)
square error solution to the two equations, i.e., find
(C) -2 (maximum) (D) -2 (minimum)
the value of x that minimizes the sum of squares of the
errors in the two equations 20. The quadratic approximation of f(x) = x3 - 3x2 - 5 at
2x = 3 the point x = 0 is [GATE, 2016]
4x = 1 [GATE, 2013] (A) 3x2 - 6x - 5 (B) -3x2 - 5
π /6 (C) -3x2 + 6x - 5 (D) 3x2 - 5
12. The solution for ∫ cos 4 3θ sin 3 6θ dθ is[GATE, 2013] xy
0 21. What is the value of lim ? [GATE, 2016]
x →0 x2 + y2
y →0
1
(A) 0 (B) (A) 1 (B) -1
15
(C) 0 (D) Limit does not exist
8
(C) 1 (D) 22. The area between the parabola x2 = 8y and the straight
3
line y = 8 is _____. [GATE, 2016]
x + sin x 23. The area of the region bounded by the parabola y = x2
13. lim equals is [GATE, 2014]
x →∞ x + 1 and the straight line x + y = 3 is [GATE, 2016]
(A) -∞ (B) 0 59 9
(A) (B)
(C) 1 (D) ∞ 6 2
10 7
xα − 1 (C) (D)
14. The expression lim is equal to [GATE, 2014] 3 6
x →0 α
(A) lnx (B) 0 24. The angle of intersection of the curves x2 = 4y and y2
(C) xlnx (D) ∞ = 4x at point (0, 0) is [GATE, 2016]
(A) 0° (B) 30°
15. With reference to the conventional cartesian (x, y)
coordinate system, the vertices of a triangles have the (C) 45° (D) 90°
following coordinates: (x1, y1) = (1, 0): (x2, y2) = (2, ∞ 1 ∞ sin x
2): and (x3, y3), = (4, 3). The area of the triangle is 25. The value of ∫0 1 + x 2 dx + ∫0 x
dx is
equal to [GATE, 2014]
[GATE, 2016]
3 3 π
(A) (B) p
(A) (B) ⋅
2 4 2
4 5 3π
(C) (D) (C) ⋅ (D) 1
5 2 2
Answer Keys
Exercises
1. A 2. D 3. B 4. D 5. C 6. D 7. D 8. C 9. A 10. A
11. D 12. C 13. B 14. A 15. D 16. D 17. A 18. C 19. B 20. D
21. A 22. D 23. B 24. D 25. A 26. D 27. B 28. C 29. A 30. A
31. C 32. D 33. D 34. A 35. A 36. B 37. C 38. B 39. D 40. A
41. B 42. A 43. D 44. C 45. C 46. D 47. C 48. B 49. A 50. C
51. D 52. D 53. D 54. D 55. D 56. A 57. B 58. D 59. A 60. A
61. C 62. D 63. C 64. C 65. D 66. A 67. C 68. B 69. C 70. D
CHAPTER HIGHLIGHTS
IntroDuction d3 y d 2 y dy
6. + 8 + + 9 y = 16 x 2
Familiarity with various methods used in evaluating indefi- dx 3 dx 2 dx
nite integrals or finding anti-derivatives of functions [or, in ∂u ∂u
other words, evaluating ∫ f(x) dx] is a pre-requisite. 7. x + y = 8u
∂x ∂y
Differential EQuations 8.
∂2u ∂2u
+ = 10
An equation involving derivatives of a dependent variable ∂y 2 ∂x 2
with respect to one or more independent variables is called ∂2u ∂2u
a differential equation. The equation may also contain the 9. = 25 2
∂y 2 ∂x
variables and/or their functions and constants. If there is
only one independent variable, the corresponding equation ∂4u ∂2u ∂4u
is called an ordinary differential equation. If the number of 10. + 6 2 2 + 4 = e3 xy
∂x 4 ∂x ∂y ∂y
independent variables is more than one, the corresponding
equation is called a partial differential equation. We note that in the given examples, Eqs. (1) to (6) are
ordinary differential equations while Eqs. (7) to (10) are
Examples: partial differential equations. We refer to these examples
dy later on in next chapter.
1. = x 4 + e− x + y
dx
2
Certain Geometrical Results may also be
d2 y dy Expressed as Differential Equations
2. x2 + 3 + 3 y 4 x = sin x + 6
dx 2 dx
Illustration 1 Consider a family of parallel lines. All these
dy lines have the same slope. If k represents the slope, we may
3. + 5 y = x 3 − tan x
dx interpret the family of parallel lines as curves having the
d2 y dy
4. + 4y = 0 same slope. As represents the slope of the tangent to
dx 2 dx
2 4 a curve at any point (x, y), we may say that the differential
d3 y dy dy
5. 3 + 5 + e 2 xy = 6 equation = k represents a family of parallel lines.
dx dx dx
dy dn y d n −1 y d n−2 y
Illustration 2 The differential equation y = k (a constant) P0 + P1 n −1 + P2 n − 2 +
dx dx n
dx dx
may be said to represent the family of curves having the
dy
length of subnormal equal k at every point (x, y) on the + Pn −1 + Pn y = Q
curve. (We may note that the family of curves is the family dx
of parabolas). Our study is confined to ordinary differential where P0, P1, P2, . . ., Pn, Q are functions of x or constants.
equations. In what follows, differential equation means If an equation is not linear, it is called a non-linear differ-
ordinary differential equations. ential equation. In examples, 1, 3, 4, 6 are linear differential
equations, while examples 2 and 5 are non-linear differential
equations.
Order of a Differential Equation
It is defined as the order of the highest derivative present in Solution of a Differential Equation
the equation. Examples (1), (3) are of first order; (2), (4) are A function y = f(x) or F(x, y) = 0 is called a solution of
of second order and (5), (6) are of third. a given differential equation if it is defined and differenti-
able (as many times as the order of the given differential
Degree of a Differential Equation equation) throughout the interval where the equation is
The degree of a differential equation is defined as the degree valid, and is such that the equation becomes an identity
of the highest order derivative present in the equation. (It dy d 2 y
when y, , 2 ,… are replaced by f(x), f ′(x), f ″(x), …
is assumed that the various order differential co-efficients dx dx
or derivatives present in the equation are made free from respectively.
dy d 2 y
fractional powers). [In the case of F (x, y) = 0 one has to get , , … by
dx dx 2
Examples (1), (2), (3), (4), (6) are of first degree while successive differentiation of F(x, y) = 0 with respect to x].
Example (5) is of second degree.
Examples:
Consider the differential equation,
dy
1. y = e7x is a solution of = 7 y, since on substitution
dy 2 5/ 2 dx
d3 y
1 + =4 . of y = e7x, both left and right sides of the differ-
dx dx 3 ential equation become identical. We find that
−1
Taking the square on both sides (to free it from fractional y = e 7 x , 3e 7 x , e 7 x or, in general, y = Ce7x, where
powers), the differential equation is 2
C is an arbitrary constant represents solutions of
5
dy 2 d3 y
2 dy
= 7 y.
1 + = 16 3 . dx
dx dx
2. y – x = 4 is a solution of the differential equation
2 2
y = e7x, y = 3e7x, … are its particular solutions. The general solu- SOLVED EXAMPLES
tion represents a family of exponential curves.
In Example (2) y2 – x2 = C represents the general solu- Example 1
dy x
tion of the differential equation = and the solutions Form the differential equation representing the one-
dx y parameter family of curves
y2 – x2 = 4, y2 – x2 = 5, … are its particular solutions. The
general solution in this case represents a family of rectan- x3 – Ay = 0.
gular hyperbolas.
Solution
3. y = 2e–3x + 5e6x is a solution of the second order differ-
d2 y dy Given, x3 – Ay = 0 (1)
ential equation − 3 − 18 y = 0. Ay = x3
dx 2 dx
(which can be verified by actual substitution). Also, dy 3x 2
y = 4e–3x – 10e6x, e–3x + e6x, … or, in general, y = Ae–3x A = 3x 2 ⇒ A= (2)
dx dy
+ Be6x where A and B are arbitrary constants represents dx
d2 y dy
solution of 2
− 3 − 18 y = 0. Substituting A in the Eq. (1),
dx dx we have
4. y = 2 cos 4x + 3 sin 4x or, in general, y = A cos 4x + B
sin 4x where A and B are arbitrary constants represents 3x 2 dy
x3 − ⋅y=0 ⇒ x − 3 y = 0.
d2 y dy dx
solutions of + 16 y = 0.
dx 2 dx
In Example (3), the general solution is y = Ae–3x + Be6x and
in Example (4), the general solution is y = A cos 4x + B sin 4x. Example 2
By assigning particular values to the arbitrary constants Obtain the differential equation of all the circles in the first
one can generate particular solutions. quadrant, which touch the co-ordinate axes.
From Examples (3) and (4), we infer that the general
solution of a second order differential equation is a relation Solution
between x and y involving two arbitrary constants such that The equation of any circle in the first quadrant, which
the differential equation is satisfied by this relation or the touches the co-ordinate axes may be represented as (x – h)2
general solution of a second order differential equation is + (y – h)2 = h2.
a two-parameter family of curves where the parameters are Differentiating with respect to x,
the arbitrary constants.
dy
To sum up, the general solution of an nth order differ- 2( x − h) + 2( y − h) =0
ential equation is a relation between x and y involving n dx
dy
arbitrary constants, such that the differential equation is sat- x+ y
isfied by this relation or the general solution of an nth order or h = dx
differential equation is an n-parameter family of curves dy
1 + dx
where the parameters are the arbitrary constants. For the
first and second order differential equations, we have Substituting the above expression for h in the equation
First Order Equation of the circle
One parameter family of curves: dy
2
dy
2
dy
2
Example 3
1
3∫
− dt = ∫ x 2 dx
dy 1 + y2
Solve: = .
dx 1 + x2 −1 x3
t= +c
Solution 3 3
dy 1 + y2 1 3 x3
= − e− y = + c.
dx 1 + x2 3 3
1 1 Given: When x = 1, y = 0;
dy = dx
1 + y2 1 + x2 1 1
− e° = + c
3 3
Integrating on both sides,
2
1 1 c=−
∫ 1+ y2
dy = ∫
1 + x2
dx. 3
1 3 x3 2
∴ The solution is − e − y = − .
sinh y = sinh x + c.
–1 –1
3 3 3
Example 4 x3 + e–y3 – 2 = 0.
dy
Solve: ( x − xy 2 ) + ( y + x 2 y ) = 0. Homogeneous Differential Equations
dx
Homogeneous differential equation will be of the form f (x,
Solution y)dy = g(x, y)dx, where f (x, y) and g(x, y) are homogeneous
dy functions in x and y of the same degree.
( x − xy 2 ) + ( y + x 2 y) = 0
dx
(x – xy2) dy + (y + x2y) dx = 0 Definition
A function F(x, y) in x and y is a homogeneous function in
x(1 – y2) dy + y (1 + x2) dx = 0
x and y of degree n(n, a rational number), if F(x, y) can be
1 − y2 1 + x2 y x
dy + dx = 0 expressed as x nφ or y nψ .
y x x y
4 y y3 Solution
1. x 3 + 4 x 2 y − y 3 = x 3 1 + − is a homogeneous
x x3 dy y
Given: x = y + x sin
function in x and y of degree 3. dx x
y dy y y
2. x 3 tan is a homogeneous function in x and y of = + sin (1)
x
dx x x
degree 3. dy dv
Put y = vx, =v+x .
x+ y dx dx
3. is a homogeneous function in x and y of
2x − 3y Substituting in (1) we get,
degree 0. We change the dependent variable y to v by
dv
dy dv v+x = v + sin v
the substitution y = vx. Then, =v+x . dx
dx dx
dy xdv 1 1
On substitution y and in the given homogeneous ⇒ = sin v ⇒ dv = dx
dx dx sin v x
equation, it reduces to the variables separable form. 1
= ∫ cosec v dv = ∫ dx
x
Example 6 ⇒ log (cosec v – cot v) = log x + log c
dy ⇒ cosec v – cot v = cx
Solve: x 2 = x 2 + 7 xy + 9 y 2 .
dx
y y
cosec − cot = cx.
Solution x x
dy
x2 = x 2 + 7 xy + 9 y 2
dx Example 8
2
dy 7y y Solve 3y2 dx + (2xy + 3x2) dy = 0.
= 1+ + 9
dx x x
Solution
dy dv
Put y = xv ⇒ =v+ 3y2 dx + (2xy + 3x2) dy = 0.
dx dx
dy −3 y 2
dv =
v+ = 1 + 7v + 9v 2 dx 2 xy + 3 x 2
dx
dy dv
dv Put y = vx ⇒ =v+x
x = 9v 2 + 6 v + 1 dx dx
dx
dv −3v 2
v+x =
1 1 dx 2v + 3
dv = dx
9v + 6 v + 1
2 x dv −3v 2
x = −v
Integrating on both sides, dx 2v + 3
1 1 dv −3v 2 − 2v 2 − 3v
x =
∫ 9v 2 + 6v + 1 dv = ∫ x dx dx 2v + 3
1 1 1 2v + 3 1
dv = dx
∫ (3v + 1)2 dv = ∫ x dx − 3(3v + 1) = log x + log c −5v − 3v
2 x
2v + 3 1
1 −x ⇒ dv + dx = 0
=− = loge cx = = loge cx v(5v + 3) x
3y 9 y + 3x
3 + 1 Integrating on both sides,
x
1 3 1
where C is an arbitrary constant. ⇒ ∫ v − 5v + 3 dv + ∫ x dx = 0
Example 7 3
⇒ log v − log (5v + 3) + log x = log c.
dy y 5
Solve x = y + x sin
dx x ⇒ 5 log v –3 log (5v + 3) + 5 log x = 5 log c.
v5 Example 10
⇒ log x 5 = log c5 Find the solution of (ey + 1) cot x dx + ey log(sin x) dy = 0.
(5v + 3)3
y5 Solution
⇒ 3
= c1 , where c1 = c5
y Given (e y + 1)cot x dx + e y log (sin x)dy = 0
5 x + 3
Let M = (e y + 1)cot x and N = e y log (sin x)
y5 x3 ∂M ∂N
⇒ = c1 ⇒ x 3 y 5 = c1 (5 y + 3 x )3 = e y cot x and = e y cot x
(5 y + 3 x )3 ∂y ∂x
x + 2y y − 2x 1 ∂M ∂N 1
2 dx + 2 dy = 0 − = [−2 y − 2 y ]
x +y
2
x +y
2
N ∂y ∂x 2 xy
−2
The solution is U + ∫ φ ( y ) dy = C = = f ( x)
x
x Integrating factor (IF)
x + 2y
U = ∫ M1dx, where M1 =
x2 + y2 = e ∫ f ( x ) dx
−2 1
1
= e∫
dx log
= e∫ −
2 log x
x
x + 2y x =e x2 =
=∫ 2 dx x2
x + y2 1
∴ Multiplying the given equation with 2 , we get
x 1 x
=∫ dx + 2 y ∫ 2 dx x2 − y2 2 xy
x2 + y2 x + y2 2 dx + dy = 0
x x2
1 1 x x2 − y2 y
= log( x 2 + y 2 ) + 2 y tan −1 2 dx + 2 dy = 0 (2)
2 y y x x
x2 − y2 2y
1 x M1 = and N1 =
= log ( x 2 + y 2 ) + 2 tan −1 x2 x
2 y ∂M1 −2 y ∂N1 −2 y
y − 2x = 2 , and = 2
Since in N1 = 2 there is no term independent of x, ∂y x ∂x x
x + y2
the solution is ∂M1 ∂N1
=
∂y ∂x
1 x
log ( x 2 + y 2 ) + 2 tan −1 = C ∴ Eq. (2) is an exact equation and its solution is
2 y x
Method 2
∫ M1dx + ∫ (the terms of N1 not containing x) dy = C
x 2
x − y2
If the differential equation Mdx + Ndy = 0 is of the form ∫
x2
dx + ∫ 0 dy = C
1 x
y2 y2
y, f(xy)dx + x g(xy)dy = 0 and Mx – Ny ≠ 0, then is ⇒ ∫ 1 − 2 dx = C ⇒ x + = C.
Mx − Ny x x
an integrating factor of Mdx + Ndy = 0.
Example 13
Method 3 Find the solution of xy 2dx + ( y + y2)dy = 0.
In the equation Mdx + Ndy = 0, Solution
1 ∂M ∂N Given xy2dx + ( y + y2)dy = 0 (1)
= f ( x ), then e ∫
f ( x ) dx
if − is an integrating
⋅
Mdx + Ndy = 0
N ∂y ∂x
factor of the given equation. M = xy2; N = y + y 2
∂M ∂N
1 ∂N ∂M = 2 xy and =0
Similarly if − = g ( y ) then e ∫g ( y ) dy is an ∂y ∂x
M ∂x ∂y
∂M ∂N
integrating factor of the given equation. ≠
∂y ∂x
Example 12 1 ∂N ∂M 1
− = [−2 xy ]
Find the solution (x2 – y2)dx + 2xy dy = 0. M ∂x ∂y xy 2
Solution −2
= = g( y)
Given (x2 – y2)dx + 2xy dy = 0 (1) y
M = x – y and N = 2xy
2 2 Integrating factor is e∫g(y)dy
−2 1
∂M ∂N ∫ dy log 1
= −2 y and = 2y =e y = e −2 log y dy = e y2 =
∂y ∂x y2
∂M ∂N 1 xy 2 dx y + y 2
≠ Multiplying Eq. (1) by , we get + 2
dy = 0
∂y ∂x y2 y2 y
d
Now, ( ye ∫ Pdx ) = 0 General solution
dx
dy Pdx dy y ⋅ IF = ∫Q ⋅ IF dx + c.
That is, e ∫ Pdx + ye ∫ × P = 0 or e ∫ + Py = 0.
Pdx
dx dx sin 3 x
y(1 + x 4 ) = ∫ (1 + x 4 ) dx + c
This means that if we multiply both sides of Eq. (2) by 1 + x4
e∫Pdx, the product 3 sin x − sin 3 x
= ∫ sin 3 x dx + C = ∫ dx + c
dy d
e ∫ Pdx + Py is {ye ∫Pdx}. The factor e ∫Pdx is called 4
dx dx cos 3 x 3
an integrating factor of Eq. (2). y(1 + x 4 ) = − cos x + c
12 4
Suppose we multiply both sides of Eq. (1) by e ∫Pdx, it is 12y(1 + x ) = cos3 x – 9 cos x + c
4
2 x2
x− 1
IF = = e ∫ Pdx e x = 3 + Ce 2
y
General solution is y ⋅ IF = ∫Q ⋅ IFdx + c 1
2 2 y=
x− 8 x− x2
ye x = ∫ 4 + 2 e x dx + c 3 + Ce 2
x
2
2 x− Example 18
= 4 ∫ 1 + 2 e x dx + c
x dy y y
x−
2 Solve + log y = 3 (log y ) 2 .
(Put e x =t dx x x
2
x− 2 Solution
⇒ 1 + 2 dx = dt )
e x
x dy y log y y(log y ) 2
= 4 ∫ dt + c = 4t + c Given + =
dx x x3
The general solution is
1 dy 1 1 1
x−
2
x−
2 ⇒ 2
+ ⋅ = 3 (1)
ye x = 4e x +c y(log y ) dx x (log y ) x
1
Let = u,
Bernoulli’s Linear Equations log y
dy
An equation of the form + Py = Qy n is called Bernoulli’s Differenting wrt x
dx −1 1 dy du
linear equation, where P, Q are continuous functions in x. 2
⋅ =
(log y ) y dx dx
Example 17
∴ Eq. (1) becomes
dy
Solve + xy = − (3 xy 2 ). −du 1 1 du 1 −1
dx ⇒ + u= 3 ⇒ − u= 3
dx x x dx x x
Solution
It is a linear equation in u.
dy
Given + xy = − (3 xy 2 ) −1 −1
dx Here P = and Q = 3
Throughout the equation dividing with y2 we get x x
1
dy − ∫ dx 1
y −2 + xy −1 = −3 x (1) IF = e ∫ Pdx = e x = e − log x =
dx x
dy du ∴ Solution is u . IF = ∫QIFdx + c
Let y −1 = u ⇒ − y −2 =
dx dx 1 −1 1
u = ∫ 3 ⋅ dx + c
−du x x x
The Eq. (1) becomes + xu = −3 x
dx 1
u = − ∫ x − 4 dx + c
du x
− xu = 3 x
dx 1 1
= 3 +c
The above equation is a linear differential equation in u. (log y ) x 3 x
− x2
∴ IF = e ∫ Pdx = e − ∫ xdx = e 2
Second Order Linear Differential Equations
with Constant Co-efficients
∴ Solution is u ⋅ IF = ∫ QIF dx
The standard form of a second order linear differential equa-
− x2 − x2 tion with constant co-efficients is
u⋅e 2 =∫ 3 xe 2 dx.
d2 y dy
a02
+ a1 + a2 y = F ( x )(1)
−x2 dx dx
= − ∫ 3e −t dt when t = where a0, a1, a2 are real constants and F(x) is a function of
2
only x. The second order equation,
−3e −t
= = 3e −t d2 y dy
−1 a0 + a1 + a2 y = 0(2)
− x2 − x2 dx 2 dx
u⋅e 2 = 3e 2 +C represents the corresponding homogeneous equation.
Let y = u (x) represent the general solution of Eq. (2) Since the roots of the auxiliary equation are equal and
[u(x) will contain two arbitrary constants]. This means that each equal to m1, this reduces to
a0(D – m1)2 y = 0 or (D – m1)2 y = 0 (7)
d 2u du (since a0 ≠ 0)
a0 2
+ a1 + a2 u = 0 (3)
dx dx Let (D – m1)y = Y1(8)
Let y = v (x) represent a particular solution of the given Then, Eq. (7) becomes (D – m1)Y1 = 0. (9)
equation of Eq. (1). We have, then,
dY
d 2v dv Now, Eq. (9) is reduced to 1 − m1Y1 = 0, giving Y1 =
a0 2 + a1 + a2 v = F ( x )(4) C1em1x as the solution. dx
dx dx dy
Substituting in Eq. (8), − m1 y = c1e m1 x is a linear
Substituting y = u(x) + v(x) in Eq. (1), dx
equation. The general solution is given by ye–m1x = c2 + ∫c1em1x
d2 d × e − m1 x dx = c2 + c1 x
a0 2
(u + v ) + a1 (u + v ) + a2 (u + v )
dx dx
or y = c2 e m1 x + c1 xe m1 x = e m1 x (c2 + c1 x )
d 2u du d 2v dv where c1 and c2 are arbitrary constants.
= a0 2 + a1 + a2 u + a0 2 + a1 + a2 v
dx dx dx dx Case 3: Let the roots of (V) be complex. Let us assume the
= 0 + F(x) (by Eqs. (3) and (4)) roots as the conjugate pairs a ± ib. (The co-efficients a0, a1,
a2 being real, roots occur in conjugate pairs).
= F(x).
The general solution is y = c1e(α + iβ ) x + c2 e(α −iβ ) x
We infer that y = u(x) + v(x) is the general solution of
the Eq. (1). Thus, the general solution of Eq. (2) is the sum = c1eα x (cos β x + i sin β x ) + c2 eα x (cos β x − i sin β x )
of the general solution of the corresponding homogeneous
equation (2) and a particular solution of the given equation = ea x{(c1 + c2) cos b x + i(c1 – c2) sin b x}
(1). y = u(x) is called the complementary function of Eq. (2) = ea x{A1 cos b x + A2 sin b x).
and y = v(x) is called a particular integral of Eq. (1). The where A and B are arbitrary constants. We may now
general solution of Eq. (1) is given by y = u(x) + v(x). summarize the nature of the complementary function of
= [Complementary function] + [Particular integral] Eq. (1) as follows:
= CF + PI (in short). Roots of the Auxiliary Complementary Function of
Equation a0m 2 + a1m + a2 = 0 (1), or General Solution of (2)
To find the complementary function of Eq. (1) or to
obtain the general solution of the homogeneous equation Roots, real and distinct, say y = c1em1x + c2em2x
m1, m2
dy
(2): As y = emx is a solution of − my = 0, we assume y = Roots, real and equal, say y = (c1 + c2 x)em1x
dx each equals m1
emx (for some value of m) to be a solution of Eq. (2).
Roots, complex, say a ± ib y = ea x{c1cos bx + c2 sin bx}
d2 d
Then, a0 2 (e mx ) + a1 (e mx ) + a2 e mx must be equal Roots, complex and y = e a x[(c1 + c2x) cos bx + (c3 +
dx dx repeated, say m1 = m2 = a + c4 x) sin b x
to zero (or) emx{a0 m2 + a1 m + a2} = 0. ib and m3 = m4 = a – ib
Since emx cannot be equal to zero, a0 m2 + a1 m + a2 = 0 (5)
Eq. (5) is called the auxiliary equation corresponding to Example 23
(1) [or (2)]. Eq. (5) is quadratic in m and gives two values
Obtain the complementary function of the equation
for m, which may be real or complex.
d 2 y 7dy
Case 1: Let the roots of Eq. (5) be real and distinct, say m1 − + 6 y = x4 .
= dx 2 dx
and m2 (m1 ≠ m2). Then, y e=m1 x and y e m2 x are two distinct
solutions of (2) or y = C1e 1 + C2 e 2
mx mx
(6) Solution
(C1 and C2 are arbitrary constants) is the general solution d2 y dy
of (2) or the complementary function of (1). − 7 + 6 y = x4
dx 2 dx
Case 2: Let the roots of (5) be real and equal and each ⇒ (D – 7D + 6) y = x
2 4
Example 25 1
Let e kx = X1
Obtain the complementary function of the equation (D − k)
d2 y dy dX1
2
− 6 ⋅ + 10 y = e3 x . Then ( D − k ) X1 = e kx or − kX 1 = e kx
dx dx dx
This is a linear equation and the particular solu-
Solution tion of the above equation is xekx. Therefore, particular
d2 y dy 1
Given: − 6 ⋅ + 10 y = e3 x integral = xe kx .
dx 2 dx a0 ( k − m0 )
⇒ (D2 – 6D + 10)y = e3x 2. Suppose both the roots of the auxiliary equation are k.
Auxiliary equation is m2 – 6m + 10 = 0 Then, particular integral
1
6 ± 36 − 40 6 ± 2i = [e kx ]
m= = = 3±i a0 ( D − k ) 2
2 2
1 1
∴ The complementary function is given by yc = e3x(c1cos x = e kx
+ c2sin x). a0 ( D − k ) ( D − k )
1
To find a particular integral of Eq. (1) or to find a particular = [ xe kx ],
solution of the Eq. (1): a0 ( D − k )
1
d2 y dy Use the result in (1) . Now, let ( xe kx ) = X 2
a0 2
+ a1 + a2 y = F ( x ) D−k
dx dx
dX 2
We have, therefore, (D – k) X2 = xekx or − kX 2 = xe kx
We may write the above as (a0D2 + a1D + a2) y = F(x) or dx
f (D) y = F(x) where f (D) stands for (a0D2 + a1D + a2). which is a linear equation.
Particular integral y is that function of x independent of x 2 kx
Particular solution is X 2 = e or, particular inte-
arbitrary constants such that f (D) on y or f (D) y yields F(x). 2
1 x 2 kx
This is symbolically represented as y = {F ( x )}. gral in this case is given by y = e .
f ( D) 2
Example 26
Case 1: F(x) = ekx where k is a constant. Solve the differential equation:
We have D(ekx) = kekx, D2(ekx) = k2ekx … or, in general,
(D2 + 5D + 6)y = e–4x
Complementary function is c1e–3x + c2e–2x. Case 2: F(x) = sin kx or cos kx where k is a constant.
1 We have D{sin kx} = k cos kx
Particular integral = ⋅ e− 4 x
D 2 + 5D + 6 D2{sin kx} = – k2 sin kx
−4x
1 e
= e− 4 x = Similarly, D2{cos kx} = – k2 cos kx
(− 4) 2 + 5( − 4) + 6 2
If g(D2) is a polynomial in D2,
∴ General solution is
e− 4 x g(D2) {sin kx or cos kx} = g(–k2) sin kx or g(–k2) cos kx.
y= c1e −3 x + c2 e −2 x + .
2 1 1 1
Hence, sin kx = sin kx and cos kx
2
g( D ) g ( −k )
2 g( D 2 )
Example 27
1
Solve (3D2 – D – 10)y = 6e2x = cos kx, provided g(–k2) ≠ 0.
g ( −k 2 )
Solution We shall illustrate the above technique by considering
Given (3D2 – D – 10)y = 6e2x two examples.
Auxiliary equation 3m2 - m - 10 = 0
−5 Example 29
m = 2, .
3 Find the particular integral of the equation (D2 + 16)y
∴ Complementary function is = cos 3x.
−5
x
CF = c1e 2 x + c2 e 3 Solution
1 1 1 cos 3 x
PI = 6e 2 x PI = cos 3 x = cos 3 x =
3D 2 − D − 10 D2 + 16 − (3) + 16
2 7
1
= 6e 2 x Example 30
( D − 2) (3D + 5)
Find the particular integral of the equation (D2 – 5D + 6) y
1 1 1 1 2x
=6 e2 x = 6 e 1 = sin 3x.
D − 2 3D + 5 ( D − 2) 11
Solution
6 1 6
= e 2 x = xe 2 x 1
11 ( D − 2) 11 PI = sin 3 x
D 2 − 5D + 6
∴ General solution is
1
5
− x 6 = sin 3 x
y = c1e 2x
+ c2 e 3 + xe 2 x . −3 − 5 D + 6
2
11
1
Example 28 sin 3 x
−5 D − 3
Solve (D2 – 12D + 36)y = e6 x
5D − 3
Solution =− sin 3 x
(5 D + 3) (5 D − 3)
Given: (D 2 – 12D + 36) y = e6 x
Auxiliary equation is m2 – 12m + 36 = 0. (5 D − 3) 3 − 5D
= sin 3 x = sin 3 x
m2 – 12m + 36 = 0. − ( 25 D − 9)
2 25 × ( −9) − 9
m = 6, 6 1
Complementary function (CF) = (c1 + c2x)e6x = [(3 − 5 D ) sin 3 x ]
−234
1 1
PI = e6 x = e6 x −1
D2 − 12 D + 36 ( D − 6) 2 = [3 sin 3 x − 5 D(sin 3 x )]
234
x2 6x
= e −1
2! = [3 sin 3 x − 15 cos 3 x ]
234
∴ General solution is y = CF + PI
x2 6 x0 15 cos 3 x 3 sin 3 x
= (c1 + c2 x ) e6 x + e PI = −
2! 234 234
1 dz 1
Similarly, if we have to find [cos kx ]. ⇒ =
D2 + k2 dx x
1 dy dy dz
We write it as the real part of (eikx ) = ⋅
D + k2
2 dx dz dx
1 dy dy 1
= Real part of (eikx ) = ⋅
( D − ik )( D + ik ) dx dz x
x dy dy
= Real part of ( − i cos kx + sin kx ) =x
2k dz dx
x sin kx
= . d 2 y d dy d 1 dy
2k = = ⋅
dx 2 dx dx dx x dz
1 −x
sin kx = cos kx −1 dy 1 d dy
D2 +k 2 2k = +
x 2 dz x dx dx
1 x −1 dy 1 d dy dz
cos kx = sin kx = 2 +
D2 +k 2 2k
x dz x dz dz dx
d 2 y −1 dy 1 d 2 y
Example 31 = +
dx 2 x 2 dz x 2 dz 2
Solve the equation (D2 + 16) y = sin 4x.
d 2 y d 2 y dy d dy
Solution x2 2
= 2 − = − y
dx dz dz dz dz
Given: (D2 + 16) y = sin 4x
Auxiliary equation is m2 + 16 = 0 dy dy d2 y
Let =θ y ⇒ x = θ y, x 2 2 = θ (θ −1) y
m = ± 4i dz dx dx
−1
d3 y e2 z θ 3 + 9θ 2 + 28θ
Similarly x 3 = θ (θ −1)(θ − 2) y, and so on.
dx 3 1 + z
30 30
Then Eq. (1) is changed into a linear differential equation.
e 2 z θ 3 + 9θ 2 + 28θ
We solve this by methods discussed earlier. = 1 − z
30 30
Example 32
e2 z 28 2 z
d2 y dy = z− e
Solve x 2 2 + 3 x − 3 y = 0 30 (30) 2
dx dx
y = CF + PI
Solution
e2 z 28 2 z
Let x = ez or z = log x = C1e–z + e–z(C2 cos z + C3 sin z) + z− e
30 (30) 2
dy d2 y
Then x = θ y; x 2 2 = θ (θ −1) y C1 1 x 2 log x 28 2
dx dx = + (C2 cos(log x ) + C3 sin(log x )) + − x
x x 30 900
The above equation becomes
[θ (θ −1) + 3θ − 3] y = 0 Example 34
θ 2 + 2θ − 3 y = 0 d2 y dy
Solve ( 2 x − 1) + 2( 2 x − 1) − 100 y = 32( 2 x − 1) 2
dx 2 dx
Auxiliary equation is m2 + 2m – 3 = 0
⇒ (m + 3)(m – 1) = 0 Solution
⇒ m = –3, 1 Let 2x – 1 = u
∴ y = c1e–3z + c2ez du
2=
= c1x–3 + c2x. dx
dy dy du dy
Example 33 = ⋅ =2
dx du dx du
d3 y 2
2 d y + 8 dy + 2 = 2 log .
Solve x 3 + 6 x y x x d 2 y d dy d dy
dx 3 dx 2 dx = = 2
dx 2 dx dx dx du
Solution
d dy du d2 y
Put x = ez or z = log x. Then =2 ⋅ = 22 2
du du dx du
dy d2 y ∴ The given equation becomes
x = θ y, x 2 2 = θ (θ −1) y,
dx dx
d2 y dy
22 u 2 + 2 ⋅ 2u − 100 y = 32u 2
d3 y du 2 du
x3 = θ (θ −1)(θ − 2) y
dx 3 d2 y dy
u2 2 + u − 25 y = 8u 2
The given equation becomes du du
dy d2 y
[q(q – 1) (q – 2) + 6q(q – 1) + 8q + 2]y = e2z ⋅ z Let u = e , u
z = 0; x 2 2 = θ (θ −1)
dx dx
(q 3 + 3q 2 + 4q + 2)y = e2z ⋅ z
AE = m3 + 3m2 + 4m + 2 = 0 [θ (θ −1) + θ − 25] y = 8e 2 z
(m + 1)(m2 + 2m + 2) = 0 θ 2 − 25 y = 8e 2 z
m = – 1 or m = –1 ± i
AE = m2 – 25 = 0 ⇒ m = ±5
CF = C1e–z + e–z (C2cos z + C3sin z )
1 CF = C1e–5z + C2e5z
PI = ⋅ e2 z z
3 2
θ + 3θ + 4θ + 2 1 1 −8 2 z
PI = ⋅ 8e 2 z = 8.e 2 z 2 = e
θ 2 − 25 2 − 25 21
1
= e2 z z
(θ + 2 + 3(θ + 2) 2 + 4(θ + 2) + 2
)3 8 2z
y = CF + PI = C1e5z + C2e5z – e
21
1 8 2
= e2 z ⋅z = C1u −5 + C2 u 5 − u where u = (2x – 1).
θ3 + 9θ 2 + 28θ + 30 21
2 ∞
k = ∫ e − st dt + ∫ e − st ⋅ tdt
3. (a) Let k be a constant L {k } =
s 1 2
2 ∞ ∞
1 e − st e − st e − st
(b) L {1} = , s > 0 = ]∫ + t ⋅ ]∫ − ∫ ⋅ dt
s −s 1
−s 2 2 −s
n!
{ }
∞
4. L t n = ,s>0 e −2 s e − s 2e −2 s 1 e − st
s n +1 =− + + +
s s s s −s 2
s
5. L {cos at} = ,s>0 −e −2 s e − s e −2 s 1 −2 s
s2 + a2 = + +2 + 2e
s s s s
a
6. L {sin at} = ,s>0 e −2 s 1 e − s
s2 + a2 = 1 + + .
s s s
s
7. L {cosh at} = ,s> a
s2 − a2 Example 39
a Find the Laplace transform of the function
8. L {sinh at} = ,s> a
e − st e − st
2 ∞ ∞ − st
s2 − a2 e
f (t) = sin 2t, 0 < t < p = ∫ +t ⋅ ∫ −∫ ⋅ dt
− s 1 − s 2 2 − s
n!
9. L{t n ⋅ e at } = , n∈ Z+
( s − a) n +1 = 0, t > p
∞
1 Solution
10. L f (t ) = ∫ F ( s) ds ∞
t s
L{f (t)} = ∫ e − st f (t )dt
0
Example 37 π ∞
Find the Laplace transform of the function = ∫ e − st ⋅ sin 2tdt + ∫ e − st ⋅ 0 dt
0 π
f ( x ) = 5e 2 x + 7e −3 x π
= ∫ e − st sin 2tdt
Solution 0
L{ f ( x )} = L(5e 2 x + 7e −3 x ) e − st
= [− s sin 2t − 2 cos 2t ]]π0
= 5L(e2x) + 7L(e–3x) s2 + 4
1 1 2(1 − e −π s )
L{ f (t )} = 5 ⋅ + 7⋅ = .
s−2 s+3 s2 + 4
5 7
= + ⋅
s−2 s+3 Example 40
Example 37 Find the Laplace transform of the function f (t) = (sin t + cos t)2
Find L{f (t)} where Solution
f (t) = 0, 0 < t < 1 L{(sint + cost)2} = L{1 + sin2t} = L{1} + L{sin2t} =
= 1, 1 < t <2 1 2
+ 2
= t, t > 2. s s +4
u
du.
s 0
s NOTE
6. L−1 2 = cos at This is also called as Heavisides unit function
s + a2
s 4. Periodic function: If f(t) is a periodic function with
7. L−1 2 = cos hat period a i.e., f (t + a) = f(t), then
s − a2
a
1 1
8. L−1 2
s − a2
= sin hat
a
∫ e − st f (t )dt
0
L{ f(t)} =
1 1 at 1 − e − sa
9. L−1 2
= e sin bt G ( s)
( s − a ) 2
+ b b 5. Using partial fractions: If F(s) is of the from
H ( s)
s−a
10. L−1 2 + 2
= e at cos bt where G and H are polynomials in S then break F(s)
( s − a ) b into partial fractions and manipulate term by term.
1 1 6. Heavisides expansion formula: Let F(s) and G(s)
11. L−1 2 2 2 = 3 (sin at – at cos at) be two polynomials in ‘s’ where F(s) has degree less
( s + a ) 2a
than that of G(s). If G(s) has n distinct zeros ar, r = 1,
s 1 2, 3, …., n
12. L−1 2 2 2 = t sin at
( s + a ) 2a i.e., G(s) = (s – a1)(s – a2)…(s – an), then
To find the inverse Laplace transform we use the following F ( S ) n F (α r ) α t
methods.
L−1 =∑ e r
G ( S ) r =1 G ′(α r )
1. Using the following properties
(a) If L−1{F ( s)} = f (t ), than L−1{F ( s − a)} = e at f (t ) Transform of Special Functions
(b) If L−1{F ( s)} = f (t )) and f (0) = 0; then 7. Bessel function:
d x2 x4 x6
(i) L−1{sF ( s)} = ( f (t )) J0 (x) = 1 -
+ 2 2 − 2 2 2 +
dt 2 2 ⋅4 2 ⋅ 4 ⋅6
dn
(ii) L−1{s n F ( s)} = n ( f (t )) if f (0) = f 1(0) = 1
dt then L {J0 (x)} =
f (n-1)(0) = 0 s2 +1
1 −2t −1 1
L−1 = e L 5/ 2 t t
x
( s + 2) = ∫ f 2 ( x ) f1 (t − x )dx = ∫ sin 2 x (t − x )dx
5 / 2
s
0 0 4
5 3
−1 t t
t2 4t 2 e −2t t 1
= e −2t
5
=
3 π
=
40∫ x sin 2 xdx − ∫ x 2 sin 2 xdx
40
Γ
2 t
t x 1
= − cos 2 x + sin 2 x
e 2 −3s e −3 s
4 2 4 0
\ L−1 5/ 2
= e 2 L−1
( s + 2) ( s + 2)
5 / 2
t
1 x2 x 1
− − cos 2 x + sin 2 x + cos 2 x
4 4 2 2 4 0
= (t − 3)3/ 2 e −2( t − 4 ) ⋅ H (t − 3)
3 π 1
= (1 − t sin 2t − cos 2t )
(when expressed in terms of Heaviside’s unit step function) 16
Exercises
d2 y 4. Find the solution of tan y sec2 x dx + tan x sec2ydy = 0
1. The order and degree of the DE 2 = n2y respectively π
dx when x = y = .
are
⋅
4
(A) 1, 2 (B) 1, 1 (A) tan x tan y = 1
(C) 2, 2 (D) 2, 1 (B) cot x tan y = 1
2. The differential equation whose solution is y = mx + (C) tan x cot y = 1
4 (D) cot x cot y = 1
, where ‘m’ is parameter is
m 5. The general solution of the DE, (ex + 1)ydy = (y + 1)
dy dy dy
2 exdx is
x − y + 4
(A) = 0. (A) log (ex + 1) - log (y + 1) + c = 0
dx dx dx
2 (B) log (ex + 1) = y - log (y + 1) + c
dy dy
dx − dx + 4 = 0.
(B) (C) log (ex - 1) + log (y + 1) + c = 0
ex
dy
x − y + 4 = 0.
(D) log =c
(C)
dx y +1
2
dy dy dy
(D) x + + 4 = 0. 6. Solve = | x|
dx dx dx
3. If y = c1 log x + c2 log c3 + c4 ex + c5 is the general solu- x2 x2
tion of a homogeneous linear differential equation, then y=
(A) + c (B) y= + x+c
2 2
the order of the equation is
−x | x | x | x|
(A) 2 (B) 3 (C) y= + c (D)
y= +c
(C) 4 (D) 5 2 2
2
3 3
(C) tan-1 (3x + y + 1) = x + c π +2
y cosec x - x =
(A)
2
2 2x − y +1
(D) tan-1 = x+c cosec x π +2
3 3 (B) +x=
y 2
dy x − y
9. The general solution of = is
dx x + y π −2
y cosec x + x =
(C)
x2 + xy + y2 = k (B)
(A) x2 - y2 = k 2
x - 2xy - y = k (D)
(C) 2 2
x2y2 = k cosec x π +2
(D) −x=
dy x − 2 y +1 y 2
10. The general solution of = is
dx 2 x − 4 y + 3 dy
(A) x2 - 4xy - 6y = c 18. The general solution of x + y = y2logx is
dx
(B) x2 - 4xy + 4y2 + 2x - 6y = c
y = log x + cx (B)
(A) y = x + c log x
(C) x2 + 4xy + 4y2 + 2x - 6y = c
(D) x2 + 4xy - x + 6y = c 1 1
(C) = 1 + cx (D) = 1 + cx + log x
11. The solution of the differential equation 2xy dy + (x2 + y y
y2 +1)dx = 0 is dy
(A) x3 + xy2 + 3x = c 19. Consider the differential equation cos y + 3 x 2 sin y
dx
(B) x3 + 3xy2 + x = c = x2.
x3 To convert the above equation into linear form the sub-
(C) + xy2 + x = c
3 stituted variable is
(D) 3x2 + y2 + 2x = c (A) z = cos y (B) z = cosec y
12. The general solution of yexydx + (xexy + 2y)dy = 0 is (C) z = sin y (D) z = sec y
(A) ex + y2 = c (B) exy + y2 = c 20. The solution of (aD + bD + c) y = 0 whose auxiliary
2
22. The general solution of the differential equation 30. Solve the equation
d4x d2x d2 dy
+ 13 + 36 x = 0 is ______. 3x 2 2 + x − y = x 2 .
dt 4 dt 2 dx dx
x = (c1 + c2t) cos2t + (c3 + c4t) sin 3t
(A) y = C1x–3 + C2 x–1 + x3/7
(A)
x = c1e2t + c2e-2t + c3e3t + c4e-3t
(B) y = C1x3 + C2x + x2/7
(B)
x = (c1 + c2t) e2t + (c3 + c4t) e3t
(C) y = C1x1/3 + C2x-1 + x/7
(C)
x = c1 cos 2t + c2 sin 2t + c3cos 3t + c4 sin 3t
(D)
(D) y = C1x–1/3 + C2x + x2/7
23. The particular integral of (D2 - 4D + 3)y = e3x is 31. Laplace transform of 2sin2 2t = ______.
xe3 x 1 1 s
(A) (B)
e3x (A) + 2 (B)
2 s s + 16 s + 16
2
1
(C) e3x (D) xe2x 1 1 1 1
2 (C) − 2 (D) + 2
s s + 16 s s + 16
24. The particular integral of (D3 - 4D2)y = 6 is 32. The Laplace transform of (t + 1) is ______.
3
3 2
(A)
x2 (B) x 6 − 6 s + 3s 2 − s3 6 + 6 s + 3s 2 + s3
4 (A) 3
(B)
s s
3 −x 2
(C) - x2 (D) 6(1 + s + s 2 + s3 ) 6 + 6 s + 3s 2 + s3
4 4 (C) (D) 4
s 4
25. The particular of integral of (D2 + 3D + 2)y = cos 2x is s
33. The value of L {sinh 3tcos 3t} ______.
3 sin 2 x − cos 2 x
(A) 3 sin 2x - cos 2x (B) s 2 + 18 s 2 + 18
20 (A) 4 (B)
s + 81 s 4 + 324
cos 2 x − 3 sin 2 x cos x − sin 2 x
(C) (D)
10 40 3( s 2 − 18) 3( s 2 + 18)
(C) 4 (D)
2 6. The particular integral of (D2 - D) y = x2 - 2x + 4 is s + 324 s 4 − 324
(A) x3 - 8x + 4 (B) -x3 + 4x - 4 34. The value of L{t2cos 3t} is ______.
x3 −x 3 s 2 − 27 2 s( s 2 − 27)
(C) + 8x - 4 (D) - 4x - 4 (A) 2 (B)
3 3 ( s + 9) 4 ( s 2 + 9)3
0
16 t 8 t7 (t - x)dx is _______.
(C) (D)
35 π
⋅
105 π
⋅
1 t
(A) cos t (B) sin t
8 4 + 2s 2 2
4 1. The value of L-1 − is ______.
3s − 2 16 s 2 − 25 t t
(C)
t sin (D)
t cos
8 5t 5t 2 2
(A) sin h − cos h
3 4 4
47. Solve (D4 - 16)y = 1, y = y′ = y″ = y″′ = 0.
8 5t 5t
(B) e 2 / 3t − sin h − cos h −1
3 4 4 y=
(A) − [cos h 2t + sin h 2t ]
16
8 1 5t 1 5t
(C) e 2 / 3t − sin h − cos h
3 5 4 8 4 1
(B) y= (1 − cos h 2t + cos 2t )
(D) None of these 32
1
42. The inverse Laplace transform of is _____. −1 1
s2 − 8s + 20 (C) y= + (cos h 2t − sin t )
16 32
e 2t e 4t
(A) sin 2t (B) sin 2t
−1 1
2 2 y=
(D) + (cos h 2t + cos 2t )
(C)
e4t sin 2t (D)
e4t sin 4t 16 32
Answer Keys
Exercises
1. D 2. A 3. B 4. A 5. B 6. D 7. C 8. B 9. C 10. B
11. C 12. B 13. B 14. B 15. B 16. A 17. B 18. D 19. C 20. C
21. B 22. D 23. A 24. C 25. B 26. D 27. A 28. B 29. B 30. D
31. C 32. D 33. C 34. B 35. D 36. A 37. C 38. A 39. D 40. A
41. C 42. B 43. A 44. C 45. C 46. B 47. D 48. B
CHAPTER HIGHLIGHTS
Fourier Series where a0 an, bn are called Fourier coefficients and these are
obtained by
Periodic Function A function f(x) is said to be periodic π
if f(x + a) = f(x) for all x. The least value of a is called the 1
π −∫π
a0 = f ( x )dx,
period of f(x).
Example: sin x, cos x are periodic functions with period π
2p. 1
π −∫π
an = f ( x ) cos nxdx for n = 1, 2, 3,…
NOTES π
1
π −∫π
1. f (x) and g (x) are periodic functions with period k bn = f ( x ) sin nxdx for n = 1, 2, 3,…
then af (x) + bg (x) is also a periodic function with
period k.
2. If f (x) is a periodic function with period k, then the
k SOLVED EXAMPLES
period of f (ax) is .
a
Example 1
3. If the periods of functions f (x), g(x) and h(x) are a, b,
c, respectively, then the period of f (x) + g (x) + h (x) is Obtain the Fourier series expansion of f (x) = ex in (0, 2p).
the lcm of a, b and c.
Solution
2π
Euler’s Formula for the 1
Fourier Coefficients
a0 =
π ∫ f ( x )dx
0
Let f (x) is a periodic function whose period is 2p and is 2π
1
integrable over a period. Then f (x) can be represented by
trigonometric series.
=
π ∫ e x dx
0
∞ 2π
a0 1 x 1
2 n∫=1
f ( x) = + ( an cos nx + bn sin nx ) = e = (e 2π − 1) (1)
π 0 π
π π
2 2
π ∫0
e ax a0 =
f ( x )dx, an = ∫ f ( x ) cos nxdx
= [a cos bx + b sin bx ] π0
a2 + b2
2π Case 2: If f (x) is an odd function, then the Fourier series of
1 ex an odd function contains only sine terms, and is known as
∴ an = (cos nx + n sin nx )
π 1 + n
2
0 Fourier sine series.
∞
1 e 2π 1 f ( x ) = ∑ bn sin nx,
= (cos 2π n) −
π 1 + n 2
1 + n2 n =1
π
2π 2π 2
π ∫0
1 1 where bn = f ( x ) sin nxdx
bn =
π ∫ f ( x ) sin nxdx =
π ∫ ex sin nxdx
0 0
L
2 2 x sin nx
π 2 nπ
L ∫0
=− ∫ dx where bn = f ( x ) sin dx.
π0 8 n L
2 2 − x cos nx
π π
cos nx Half Range Expansion
=
π 8n
− ∫ n
dx
In the pervious examples we define the function f (x) with
n 0 0
the period 2L.
−4 Suppose f(x) is not periodic function and defined in half
= (π cos nπ ) the interval say (0, L) of lengths L. such expansions are
8π n2
known as half range expansions or half range Fourier series.
−1 In particular a half range expansion containing only cosine
= (con nπ ), n = 1, 2, 3, …
2n 2 series of f(x) in the interval (0, L) in a similar way half range
( −1) n +1 Fourier sine series contains only sine terms. To find the
= Fourier series of f(x) which is neither periodic nor even nor
2n2
odd we obtain Fourier cosine series and Fourier sine series
∞
( −1) n +1 of f(x) as follows. We define a function g(x) such that g (x) =
∴ f ( x) = ∑ 2
cos nx f (x) in the interval from (0, L) and g (x) is an even function
n =1 2n
in (–L, L) and is periodic with period 2L and g(x) is obtained
1 cos 2 x cos 3 x by previous methods which are discussed earlier. Similarly
= cos x − + 2 − .
2 22 3 we can obtain a fourier sine series as follows. Assume f (x) =
h(x) in (0, L) and h (x) is an odd function in the interval (–L,
L) with period 2L and evaluate h (x) by pervious methods
Function of Any Period (P = 2L) which are discussed earlier.
If the function f (x) is of period P = 2L has a Fourier series,
then f (x) can be expressed as, Example 3
a0 ∞
nπ nπ If f (x) = 1 – x in 0 < x < 1 find Fourier cosine series and
f ( x) = + ∑ ( an cos x + bn sin x) Fourier sine series.
2 n =1 L L
where the Fourier coefficients are as follows: Solution
1
L Given f (x) = 1 – x in 0 < x < 1 since f (x) is neither periodic
L −∫L
a0 = f ( x )dx nor even nor odd function.
Let us assume g(x) = f (x) = 1 – x in 0 < x < 1
L
1 nπ = 1 + x in – 1< x < 0
L −∫L
an = f ( x ) cos xdx
L
\ g (x) is even function in (–1, 1)
L
1 nπ a0 ∞ nπ x
bn = ∫
L −L
f ( x ) sin
L
xdx ∴ g( x) = + ∑ an cos
2 n =1 L
L 1
Fourier Series of Even and Odd Functions Let f (x) be an 2
L ∫0
a0 = f ( x )dx = 2 ∫ f ( x )dx (here L = 1)
even function in (–L, L), then the Fourier series is 0
a0 ∞ nπ x 1
x2
1
f ( x) = + ∑ an cos 1
= 2 ∫ (1 − x )dx = 2 x − = × 2 = 1
2 n =1 L 2
0 0 2
L L
1 2 L
Where a0 = ∫
L −L
f ( x )dx = ∫ f ( x )dx
L0 an =
2
∫ f ( x ) cos
nπ x
dx
L0 L
L
2 nπ x 1
an = ∫
L −L
f ( x ) cos
L
dx = 2 ∫ (1 − x ) cos nπ xdx
0
Let f (x) be an odd function in (–L, L) then Fourier series is sin nπ xx
1
∞
nπ x sin nπ x
1
f ( x ) = ∑ bn sin = 2 (1 − x ) − ∫ ( −1) nπ
n =1 L nπ 0 0
∂u ∂ 2u x X ′ Y + y XY ′ = 0
= c2 2
∂t ∂x ⇒ xX ′Y = – y XY ′
∂ 2 u ∂u X′ Y′
⇒ c2 − =0 ⇒ x = − y (3)
∂x 2 ∂t X Y
In Eq. (3), as LHS is a function of x alone and RHS is a
Comparing it with Eq. (1), we have function of y alone, they are equal only if each of them is
A = c2, B = 0 and C = 0 equal to some constant
\ B2 – 4AC = 02 − 4 × c2 × 0 = 0 X′ Y′
\ x = − y = k (say ) (4)
\ One dimensional heat equation is parabolic. X Y
Similarly, it can be easily observed that Where k is a constant
2. One-dimensional wave equation:
X′
∂2 y ∂2 y From Eq. (4), x =k ⇒ xX ′ = kX
= c 2 is hyperbolic (B2 − 4AC > 0) and X
∂t 2 ∂x 2 dX
⇒ x = kx
3. The Laplace equation: dx
∂2u ∂2u dX dx
2 + 2 = 0 is elliptic (B2 − 4AC < 0) ⇒ =k⋅
∂x ∂y X x
Integrating on both sides we have, Which is a linear equation with its auxiliary equation being
dx dx m – k = 0 ⇒ m = k
∫ X = k∫ x Hence its solution is X = C1 ekx(5)
⇒ log X = k log x + log C1 −T ′
⇒ log X = log xkC
Again from Eq. (3), +3= k
T
⇒ X = C1 xk(5) T′
⇒ = 3− k
Y′ T
Again from Eq. (4), − y =k
Y ⇒ T ′ = (3 – k)T
⇒ –yY ′ = kY ⇒ T ′ – (3 – k) T = 0 (6)
dY Which is a linear equation with its auxiliary equation being
⇒ y = −kY
dy m – (3 – k) = 0
dY dy ⇒ m = 3 – k
⇒ = −k
Y y \ The solution of Eq. (6) is T = C2 e(3 – k)t(7)
Substituting Eqs. (5) and (7) in Eq. (1), we get the general
Integrating on both sides, solution of given PDE (2) as
dY dy u = X . T = (C1ekx) (C2e(3 – k)t)
∫ Y = −k ∫ y
= C1C2ekx+(3-k)t
⇒ log Y = –k log y + log C2 \ u = cekx+(3-k)t; where c = c1 c2
⇒ log Y = log y -k C2 \ u(x, t) = cekx+(3-k)t(8)
⇒ Y = C2 y -k(6)
Given u(0, t) = 4et
Substituting Eqs. (5) and (6) in Eq. (2), we get the solution \ From Eq. (8), u(0, t) = ce(3-k)t = 4et
of Eq. (1) as, Comparing on both sides, we get
z = (C1xk) (C2y -k) C = 4, 3 – k = 1
= C1C2 xk y -k ⇒ C = 4; k = 2
k Substituting these in Eq. (8), we get the required solution of
x
∴ z = C where C = C1C2 . Eq. (2) as
y u(x, t) = 4e2x+t
Example 10 One Dimensional Diffusion Equation The diffusion
Solve the PDE ux + ut = 3u; u(0, t) = 4et by the method of equation is a partial differential equation that describes
separation of variables. density fluctuations in a material undergoing diffusion.
The partial differential equation representing the one
Solution dimensional diffusion equation is
Let u = X(x). T(t)(1) ∂u ∂ 2u
be the solution of the PDE =D 2
∂t ∂x
ux + ut = 3u (2) where u(x, t) is the density of the diffusing material at time t
∂u and D is diffusion coefficient
u = XT ⇒ ux = = X ′T and ut
∂x Example 11
∂u
= = XT ′ Find the solution of the one dimensional diffusion equation
∂t
∂u ∂ 2u
Substituting these in Eq. (2), we get = D 2 on the interval x ∈ [ 0, L ] with initial
∂t ∂x
X ′T + XT ′ = 3XT condition
Dividing throughout by XT, we have u (x, 0) = f ( x ), ∀x ∈ [ 0, L ]
X′ T′ and Dirichlet’s boundary conditions
+ =3
X T u(0, t) = u(L, t) = 0 ∀t > 0
X ′ −T ′ Solution
⇒ = +3= k (3), (say)
X T We will solve the one-dimensional diffusion equation
X′ ∂u ∂2u
From Eq. (3), =k ⇒ X ′ = kX = D 2 (1)
X ∂t ∂x
⇒ X ′ – kX = 0 (4) by the method of separation of variables.
∂u Example 50
1. L = su( x, s) − u( x, 0)
∂t Solve the wave equation of a stretched string given by
∂2u ∂2u ∂2u
2. L 2 = s 2 u( x, s) − su( x, 0) − ut ( x, 0) = 9 satisfying the boundary conditions u (x, 0) =
∂t ∂t 2 ∂x 2
Solution d 2u d 2u s2
or s2u (x, s) - su (x, 0) - ut (x, 0) = 9 ⋅
or 2 − u = 0
Taking Laplace transform on both sides of the equation dx 2 dx 9
∂u ∂ 2u (1)
=2 2
∂t ∂x ∞
∂u ∂ 2u
Also u (0, s) = ∫ F (t )e − st dt = F and u (x, s) = 0 as x → ∞
L = 2L 2 0
∂t ∂x
\ The general solution Eq. of (1) is u (x, s)
d 2u
s u - u (x, 0) = 2 ⋅ 2 s
x
−s
x
dx = c1e 3 + c2 e 3
∂2u s
or − u = -4sin 2px as u (x, 0) = 8sin 2px and u (x, s) = 0 as x → ∞ ⇒ c1 = 0
∂x 2 2
The general solution of the above equation is u and u (0, s) = F (s) = c2
sx
−
4 sin 2π x Hence, u (x, s) = F (s) e 3
= Ae( S / 2 ) X + Be −( S / 2 ) X −
s
−( 2π ) 2 −
2 − sx
∴ u( x, t ) = L−1 e 3 F ( s)
or
8 sin 2π x
u = Ae S /2 X
+ Be − S /2 X
+ (1) x x
8π 2 + s F t − 3 , t > 3
But u (0, t) = 0 = u (4, t) = as L− {F ( s) = F (t )},
\ u (0, s) = 0, u (4, s) = 0 x
0, t <
3
\ From Eq. (1), we have A + B = 0
when expressed in terms of Heaviside’s unit step function.
−4 8 sin 8π u(x, t)
and 0 = Ae + Be
4
s/2 s/2 +
8π 2 + s
x x
= F t − ⋅ H t − .
⇒ Ae =0
−4
+ Be
4
s/2 s/2
3 3
Exercises
4 8
sin 2 x sin 3 x
(C) sin x + + +
π3 π3 2 3
(C) (D) ⋅ ⋅
4
( n − 1)π
(A) cos cos nx
n −1 1 1 1
1+
(A) + + +
2( n − 1)π 2 4 6
(B) cos cos nx
n −1 1 1 1
1 − + − +
(B)
1 ∞ 2 4 6
(C) + ∑ cos nx
2 n =1
1 1 1 1
1+ + + +
(C)
( n − 1)π ( n + 1) 3 5 7 9
(D) cos n − 1 − cos n + 1
1 1 1 1
3. If f (x) = x in (-p, p) then the value of bn is ______.
3 1− + − +
(D)
3 5 7 9
2
(A) 3 [6 + n2π 2 ]
n Direction for questions 7 and 8:
2 cos n π Let f (x) = x2 in the interval (-p, p).
(B) 3 [6 − n2π 2 ]
n 7. The Fourier series of f (x) is _______.
2 π2
(C) 3 [6 − n2π 2 ] −
(A)
n 3
cos nπ π2 ∞
(D) 3 [6 + n2π 2 ] (B) + ∑ cos nx
n 3 n =1
4. Find the Fourier series of f (x) which is defined as ∞
follows: ∑ cos nx
(C)
f (x) = 2 0 < x < 1 n=1
= 3 1 < x < 2 π 2 ∞ ( −1) n
∞ (D) + ∑ 4 2 cos nx
πx 3 n =1 n
∑ sin
(A)
2
n =1 π2
8. The value of
is ⋅
5 2 ∞ 1 ( 2n − 1)π x 12
(B) − ∑ sin
2 π n =1 2n − 1 2 1 1 1
1 + + − +
(A)
2 3 4
2 ∞ 1 πx
(C) ∑ sin 1 1 1
π n =1 2n − 1 2 1 + 2 + 2 + 2 +
(B)
2 3 4
(D) None of these
1 1 1
1+
(C) + + +
Direction for questions 5 and 6: 2 3 4
π −x 1 1 1
Let f (x) = in the interval (0, 2p) 1−
(D) + − +
2 22 32 4 2
(A) density of and diffusion coefficient. 24. A string is stretched between two fixed points follows
(B) diffusion and density coefficient. ∂2 y ∂2 y
(C) viscosity and diffusion coefficient. the equation 2 = a 2 2 (t > 0, x > 0) satisfying the
∂t ∂x
(D) diffusion and viscosity coefficient. boundary conditions y(x, 0) = 0, x > 0 and y(0, t) = t
22. Which of the following pair can be represented by
the same partial differential equation? (Except pos- Lt y(x, t) = 0, t ≥ 0, Find y(x, t) in terms of Heaviside’s
x→∞
sibly a change in the constant multiplying the partial unit step function.
derivatives)
(A) The one-dimensional wave equation and the one (A) (t - x) H(t - x)
-dimensional heat equation. x x
(B) The one-dimensional wave equation and the two (B) t − a H t − a
-dimensional Laplace equation.
(C) (t - xa) H(t - xa)
(C) The one-dimensional heat equation and the two-
dimensional Laplace equation. (D) None of these
(D) The one-dimensional heat equation and the one- 25. The one dimensional wave equation is ______.
dimensional diffusion equation.
∂u ∂u
23. Solution of the one dimensional heat equation (A) = c
∂t ∂x
∂u ∂ 2 u
= x > 0, t > 0 satisfying the boundary condi-
∂t ∂x 2 ∂2u ∂ 2u
(B) 2 = c 2 2
tion u(0, t) = 1, u( x ,0) = 0 is ______. ∂t ∂x
x x ∂2u ∂2u
erf
(A) erf
(B) 2 t (C) 2 = c 2 2
2 t ∂x ∂t
1 x ∂2u ∂u
erf
(C) erf
(D) (D) 2 = c 2
2 t t ∂t ∂x
(A) elliptic k α )x k α )x
C cos( kt )[C1e(
(A) + C2 e −( ]
(B) parabolic
k α )x k α )x
(C) hyperbolic Ce kt [C1e(
(B) + C2 e −( ]
(D) None of these
k k
∂u Ce kt C1 cos
(C) x + C2 cos − x
5. The solution of the partial differential equation α α
∂t
∂ 2u k k
= α 2 is of the form [GATE, 2016] C sin( kt ) C1 cos
(D) x + C2 cos − x
∂x α α
Answer Keys
Exercises
1. A 2. D 3. B 4. B 5. C 6. D 7. D 8. D 9. B 10. B
11. C 12. B 13. B 14. B 15. C 16. C 17. A 18. B 19. B 20. C
21. A 22. D 23. A 24. B 25. B
CHAPTER HIGHLIGHTS
1 0 2
0
Row and Column Equivalence Matrices 0 1 2 0
0 1 0
3
Example: P = ; Q = 0 0 0 1
Row Equivalence Matrix If B is a matrix obtained by 0 0 4
1
0 0 0 0
applying a finite number of elementary row operations
successively on matrix A, then matrix B is said to be row 0 0 0
0
equivalent to A (or a row equivalent matrix of A). are echelon matrices. The number of non-zero rows (i.e.,
value of P and Q) are 3 and 2 respectively. The value of i and
Column Equivalence Matrix If B is obtained by applying a j are tabulated below
finite number of elementary column operations successively
i 1 2 3 i 1 2
on matrix A, then matrix B is said to be column equivalent P: Q:
to A (or a column equivalent matrix of A ). j 1 2 3 j 2 4
1 3 4 Normal form of a Matrix
Example: A = 2 5 −2 By means of elementary transformations, every matrix ‘A’
1 4 −3 of order m × n and rank r (> 0) can be reduced to one of the
following forms.
1 3 4 Ir 0 Ir
R2 − 2 R1 , R3 − R1 ∼ 0 −1 −10 = B (say )
1. 2. [Ir/0] 3. [Ir] 4. 0
0 0
0 1 −7 and these are called the normal forms. Ir is the unit matrix
B is a row equivalent matrix of A. of order ‘r’.
−7 −4
⇒ l = -1 and l = 6 are the eigen values. Eigen vector theorem.
corresponding to l = -1 is obtained as follows:
Solution
2 4 1 0 x1 0 The characteristic equation of
+ 1 =
3 3 0 1 x2 0 4 2
A= is | A − λ I | = 0
3 4 x1 0 −7 − 4
⇒ =
3 4 x2 0 4 − λ 2
⇒ =0
⇒ 3x1 + 4x2 = 0 −7 −4 − λ
4 ⇒ (4 - l)(- 4 - l) + 14 = 0
3x1 + 4x2 = 0 ⇒ x1 = − x2
3 ⇒ -16 - 4l + 4l + l2 + 14 = 0
\ Eigen vector corresponding to l = -1 is, ⇒ l2 - 2 = 0 (1)
−4 −4 By Cayley-Hamilton theorem, the matrix A satisfies its
x x
X = 1 = 3 2 = x2 3 characteristic equation (1).
x2 x \ A2 - 2I = O
2 1
1 0 0 0
Similarly eigen vector corresponding to l = 6 is obtained where I = and O =
as follows: 0 1 0 0
⇒ A2 = 2I(2)
2 4 1 0 x1 0
− 6 = Now A16 = (A2)8 = (2I)8 (From Eq. (2))
3 3 0 1 x2 0
1 0
= 28 I 8 = 256 I = 256
−4 4 x1 0 0 1
⇒ =
3 −3 x2 0
256 0
⇒ - 4x1 + 4x2 = 0 and 3x1 - 3x2 = 0 \ A16 =
0 256
⇒ x1 = x2
Eigen vector corresponding to l = 6 is, Example 12
x x 1 2 0 3
X = 1 = 2 = x2 .
x
2 2x 1 If A = 0 4 −5 ; then find the value of the
0 1 0
Example 10 matrix polynomial 3A9 - 18A8 + 39A7 - 32A6 + 12A5
Find the eigen values of the matrix - 26A4 + 16A3 + 24A2 - 50A + 40I.
6 2 2
Solution
A = 2 3 1 .
2 1 3 The characteristic equation of
2 0 3
Solution A = 0 4 −5 is | A - l I | = 0
Characteristic equation of the given matrix is A − λ = 0 0 1 0
6−λ 2 2 2−λ 0 3
⇒ 2 3−λ 1 =0 ⇒ 0 4−λ −5 = 0
2 1 3−λ 0 1 −λ
Exercises
1. Which of the following is false? a11 a12 a13
(A) Every diagonal matrix is a square matrix.
(B) Every unit matrix is a scalar matrix. 6. If D = a21 a22 a23 , then which of the following is
(C) Every square matrix is a diagonal matrix. a31 a32 a33
(D) Every scalar matrix is a diagonal matrix. true? (Here, Aij is the cofactor of the element aij)
1 a12 a1n
a11 A11 + a21 A12 + a23 A32 = D
(A)
a21 2 a2 n
2. If the trace of the matrix is 55 a11 A11 + a12 A12 + a13 A13 = D
(B)
a21 A12 + a23 A32 + a12 A21 = D
(C)
an1 an 2 n a12 A21 + a21 A12 + a31 A13 = D
(D)
then the value of n is
(A) 10 (B) 11 2 3 −3
(C) 9 (D) Cannot be determined 7. The determinant value of 1 −2 2 is
3. Which of the following statement is/are false? 7 4 −4
(A) AT . BT always defined for square matrices of same (A) 0 (B) 10
order. (C) -10 (D) 15
(B) AT ⋅ B is defined for matrices of the same order.
(C) tr(AT) + tr(BT) is always defined for matrices A, B of n! ( n + 1)! ( n + 2)!
same order. 8. The value of ( n + 1)! ( n + 2 )! ( n + 3)! is
(D) AT + BT is always defined for matrices A, B of same ( n + 2)! ( n + 3)! ( n + 4)!
order.
(A) 2n! (n + 1)!
4. Consider the following statements about two square
(B) 2n! (n + 1)! (n + 2)!
matrices A and B of the same order:
(C) (2n)! (n + 1)! (n + 2)!
P: (A + B)2 = A2 + 2AB + B2
(D) 2n! (n + 3)!
Q: (A + B) (A - B) = A2 - B2
xC xC x +1 C
Then, 0 1 1
(A) both P and Q are true. 9. If f(x) = 2 x C1 x
2 C2 2( x +1) C , then f (200) is
2
(B) both P and Q are false
6 x C2 6 x C3 6( x +1) C3
(C) both P and Q are true if A and B commute
(D) P is true but Q is false. (A) 200 (B) -200
(C) 0 (D) -2001
2 1 2 −2 x 3 x
2 3+i −1
5. If 1 0 1 x −2 0 = I3 ×3, then x =
2 2 1 2 −2 x − x 10. The determinant 3 − i 0 −1 + i is
−1 −1 − i 1
(A) -1 (B) 1
1 (A) purely imaginary (B) zero
(C) (D) 2
2 (C) real (D) 10
x y z 1 2 3 4 5
2
2x y 3z 3 4 5 6
11. If A = , then |A| = ______.
x y z 1 8. I. If A = 3 4 5 6 7 , then A-1 is symmetric.
2 2 2 4 5 6 7 8
5 6 7 8 9
(A) 10xyz (B) 1
1 II. If a non-singular matrix A is symmetric, then A-1 is
(C) 0 (D) (x3 + y3 + z3 - 3xyz) also symmetric.
2
Which of the following is correct?
12. If the elements of a row or column of a given square
matrix is multiplied by 2, then the value of determinant (A) Both I and II true. (B) Both I and II false.
is ______ times the original determinant. (C) I is true, II is false. (D) I is false, II is true.
1 19. A is a third order matrix. If the value of the square of the
(A) (B) 1
2 determinant of the matrix of co-factors of A is 28561,
(C) 2 (D) 4 then |A| equals
(A) 25 (B) ±13
13. If A is a square matrix of order k and det(kA) = 27
(C) 120 (D) ±169
det(A), then k = ______.
(A) 9 (B) 1 20. If A is a square matrix of order 3, then the product of A
(C) 2 (D) 3 and its transpose is
(A) unit matrix. (B) zero matrix.
14. If A and B are two square matrices of order 4 such that
(C) identity matrix. (D) symmetric matrix.
|A| = -2 and |B| = 5, then |4AB| is
(A) -80 (B) -160 21. If A and B are two skew symmetric matrices of the
(C) -2560 (D) -256 same order then AB is skew symmetric if and only if
15. I. (a - b), (b - c), (c - a) are factors of the determinant (A) AB + BA = O (B) AB - BA = O
1 1 1 (C) AB + BA = I (D) AB - BA = I
a b c . 1 2 3
22. Rank of the matrix A = is
a2 b2 c2 4 5 6
II. If the elements of a determinant are functions of x (A) 1 (B) 2
and its two rows or columns become identical (i.e., (C) 3 (D) 4
determinant equals zero,) when we substitute x = k, 2 −1 −3
then (x - k) is a factor of the determinant.
2 3. The rank of the matrix −4 2 6 is
Which of the following is correct? −10 5 15
(A) Both I and II true. (B) Both I and II false.
(C) I is true, II is false (D) I is false, II is true (A) 0 (B) 1
16. A lower triangular matrix A = (aij)n × n is singular if and (C) 2 (D) 3
only if 1
(A) aii = 0 for all i = 1, 2, … n
24. If A = (1 2 3) and B = 2 then r(AB) is
(B) aii = 0 for atleast one i = 1, 2, … n 3
(C) aii ≠ 0 for all i = 1, 2, … n
(D) aii ≠ 0 for atleast one i, i = 1, 2, … n (A) 0 (B) 1
(C) 2 (D) 4
2 −1 0
25. Which of the following matrix is row echelon form?
17. Inverse of the matrix 1 2 3 is
−4 1 −1 1 0 −1 2
0 1 0 3 0 1 2
1 0 −1
−5 −11 9 −5 −1 −3 (A) (B)
0 0 1 −2
0 1 0
−1 −2 2 (B)
(A) −11 −2 −6 0 0 0 0
−3 −6 5 9 2 5
1 0 0 0 0 0 0 1
5 11 9 5 1 3 0
2 −1 3 0 0 1
0
(C)
1 2 −2 11 2 6
(D) (C) (D)
0 0 −1 4 0 1 0 0
3 6 −5 9 −2 −5
0 0 0 0 1 0 0 2
26. Which of the following set of vectors are linearly 35. The sum and product of the eigen values of the matrix
dependent? 2 0 −1
(A) (2, 3, 3), (3, -1, 3), (4, -2, 5) 0 4 −2
(B) (3, 4, -1), (-1, 3, 1), (-2, -7, -2) is respectively
(C) (2, 1, 4), (1, -2, 2), (-3, 1, -6) 1 3 −5
(D) (1, 3, -5), (-5, -1, 3), (4, -2, -2) (A) 0, 24 (B) 1, -24
27. The system of equations 2x - y + 3z = 9, x + y + z = 0 (C) 2, 20 (D) 4, -24
and x - y + z = 0 has/is 2 0 1
36. The eigen values of a matrix A = 0 2 p are 1, 2,
(A) unique solution.
(B) infinite solutions.
(C) only zero solution. 1 0 q
(D) inconsistent. and 3. Then the values of p and q are ______.
28. The system of equations 6x + 7y + 8z = 1, 13x + 14y + (A) p = 0, q = 0
15z = 2 and x + 2y + 3z = 2 is (B) p = any real number, q = 2
(A) consistent with unique solution. (C) p = 2, q = 0
(B) consistent with infinite solutions. (D) p = 2, q = 2
(C) inconsistent. 0 −1 2 −3
1 0 4 6
(D) None of these
37. The eigen values of the matrix is
29. The value of l for which the following system of equa- −2 −4 0 5
tion does not have a solution is 3 −6 −5 0
x+y+z=6
(A) real only (B) imaginary
4x + ly - lz = 0
3x + 2y - 4z = -8 (C) zero only (D) imaginary or zero
(A) 3 (B) -3 38. The number of linearly independent eigen vectors of
(C) 0 (D) 1 5 2
−2 1 is ______.
30. If the number of variables in the linear homogeneous
system AX = O is n, then the system will have exactly (A) 0 (B) 1 (C) 2 (D) infinite
one solution X = O, if the rank of the matrix A is 39. Which of the following is an eigen vector for the matrix
(A) 1 (B) < n
1 4
(C) ≤ n (D) n ?
31. If the equations 2x - y - z = 0, kx - 3y + 2z = 0 and -3x 2 −1
+ 2y + kz = 0 have a non-zero solution, then the value 1 −1
of k is (A)
(B)
3 1
(A) 2 (B) 1
(C) 7 (D) Both 1 and 7 3 −2
(C)
(D)
32. The system of equations a + 3y + 5z = 0, 2x - 4ay + a 1
−2
z = 0, -4x + 18y + 7z = 0 has only trivial solution if a is
6 −6 2 −2
(A) -1 or -3 (B) 1 or -3 −6 5 −4
(C) not equal to 1, -3 (D) not equal to -1 and 3 40. For a matrix A = , X = 2 is an eigen
2 −4 1 −1
2 −1 0
33. The eigen values of 0 1 −1 is vector. The corresponding eigen value is ______.
(A) -2 (B) 1
0 0 3
(C) 2 (D) 13
(A) 0, 0, 0 (B) 0, 1, 0 41. Let A be a 2 × 2 square matrix with l1 = -2 and l2 =
(C) 2, 1, 3 (D) -2, -1, -3
−4 6
34. The
characteristic roots of the inverse of the matrix -3 as its eigen values and x1 = , x 2 = as its
2
2 1 −4 7
eigen vectors then A is given by
3 1 are
1
1 0 2 4 −6
2 2 4 −5 (B)
(A) 7 −9
(A) -1, -1, 5 (B) 1, 1, 5
1 1 −2 6 2 6
(C) 1, 1, (D) -1, -1, 7 −3 (D)
(C) −4 3
5 5
2 5 4 3 −1 −1
42. Consider the matrix A = 0 1 0 let B = A-1, then 45. For the matrix A = −1 3 −1 , consider the fol-
0 −3 −2 1 1 −1
B = ______ lowing statements
−1 1 2 (P) The characteristic equation of A is l3 - 5l2 + 4l = 0
(A) [A2 - A - 4I] (B) [A - A - 4I] (Q) A-1 exists
4 4
(R) The matrix A is diagonalizable
1 −1 2
(C) [A2 + A - 4I] (D) [A - A + 4I] Which of the above statements are TRUE?
4 4 (A) P, Q and R
2 3 (B) P and R but not Q
43. If A = , then A =
15
(C) P and Q but not R
4 6
(D) Q and R but not P
(A) 814A (B) 815A
(C) 8 A 16
(D) 15A 46. If P is a modal matrix and D is a spectral matrix of a
diagonalizable matrix A, then which of the following
2 0 0 relations is NOT TRUE among A, P and D?
4 4. If A = 3 6 7 , then the value of the matrix poly- (A) PD = AP (B) DP-1 = P-1A
9 0 1 (C) A2P = PD2 (D) DP = PA
nomial 2A10 - 18A9 + 40A8 - 25A7 + 9A6 - 20A5 + 13A4 47. If A is a 3 × 3 square matrix with eigen values 0, 2, 3
- 9A3 + 20A2 - 10A is ______. with P as its modal matrix, then the eigen values of the
2 0 0 4 0 0 matrix P-1 AP are _______.
3 6 7 6 12 14 (A) 0, 2, 3
(A)
(B)
(B) 0, 4, 6
9 7 1 18 0 2
1 1
(C) 0, ,
1 0 0 0 0 0 2 3
0 1 0 0 0 0
(C)
(D)
1 1
(D) 1, ,
0 0 1 0 0 0 2 3
(A) Both [S] and [D] are symmetric. 13. Let A = [aij], 1 ≤ i, j ≤ n with n ≥ 3 and aij = i ⋅ j. The
(B) Both [S] and [D] are skew-symmetric. rank of A is [GATE, 2015]
(C) [S] is skew-symmetric and [D] is symmetric. (A) 0 (B) 1
(D) [S] is symmetric and [D] is skew-symmetric. (C) n - 1 (D) n
9 5 14. For what value of p the following set of equations will
7. The eigen vales of matrix are [GATE, 2012] have no solution?
5 8 2x + 3y = 5
(A) -2.42 and 6.86 (B) 3.48 and 13.53 3x + py = 10
(C) 4.70 and 6.86 (D) 6.86 and 9.50 [GATE, 2015]
8. What is the minimum number of multiplications 15. The smallest and largest eigen values of the following
involved in computing the matrix product PQR? 3 −2 2
Matrix P has 4 rows and 2 columns, matrix Q has
matrix are: 4 −4 6 [GATE, 2015]
2 rows and 4 columns, a matrix R has 4 rows and 1
2 −3 5
column?
[GATE, 2013] (A) 1.5 and 2.5 (B) 0.5 and 2.5
3 2 1 1 (C) 1.0 and 3.0 (D) 1.0 and 2.0
2 1
9. Given the matrices J = 2 4 2 and K = 2
16. The two eigen values of the matrix have a
1 2 6 −1 1 p
the product KT JK is. [GATE, 2014] ratio of 3 : 1 for p = 2. What is another value of p for
0 1 2 3 which the eigen values have the same ratio of 3 : 1?
1 0 3 0 [GATE, 2015]
1 0. The determinant of matrix is. (A) -2 (B) 1
2 3 0 1
7 14
3 0 1 2 (C) (D)
3 3
[GATE, 2014]
17. Consider the following linear systems:
6 0 4 4 x + 2y - 3z = a
−2 14 8 18
11. The rank of the matrix is. 2x + 3y + 3z = b
14 −14 0 −10 5x + 9y - 6z = c
[GATE, 2014] This system is consistent if a, b and c satisfy the
equation [GATE, 2016]
1 2. The sum of Eigen values of the matrix [M] is.
(A) 7a - b - c = 0 (B) 3a + b - c = 0
215 650 795
(C) 3a - b + c = 0 (D) 7a - b + c = 0
Where [M] = 655 150 835 [GATE, 2014]
18. If the entries in each column of a square matrix M add
485 355 550
up to 1, then an eigen value of M is [GATE, 2016]
(A) 915 (B) 1355 (A) 4 (B) 3
(C) 1640 (D) 2180 (C) 2 (D) 1
Answer Keys
Exercises
1. C 2. A 3. C 4. C 5. B 6. B 7. A 8. B 9. C 10. C
11. C 12. C 13. D 14. C 15. A 16. B 17. B 18. D 19. B 20. D
21. A 22. B 23. B 24. B 25. A 26. C 27. A 28. C 29. A 30. D
31. D 32. C 33. C 34. C 35. B 36. B 37. B 38. B 39. B 40. D
41. B 42. A 43. A 44. B 45. B 46. D 47. A
CHAPTER HIGHLIGHTS
☞ Probability ☞ Statistics
☞ Some special continuous distributions ☞ Hypothesis testing
m Example 2
probability of happening of E is P(E) = . Since the
n If a card is drawn from a pack of cards, find the probability
event does not occur (n - m) times, the probability of non-
of getting a queen.
occurrence of E is P ( E ).
n−m m Solution
P( E ) = = − n = 1 − P( E )
n n When a card is drawn, the number of possible outcomes is
52. The number of favourable outcomes of getting a queen
Therefore, P ( E ) + P ( E ) = 1. card is 4.
4 1
= =
The required probability .
NOTES 52 13
1. Probability [P(E)] of the happening of an event E is Example 3
known as the probability of success and the probabil-
A bag contains 5 green balls and 4 red balls. If 3 balls are
ity [ P ( E )] of the non-happening of the event is the
picked from it at random, then find the odds against the
probability of failure. three balls being red.
2. If P(E) = 1, the event is called a certain event and if
P(E) = 0 the event is called an impossible event. Solution
3. Instead of saying that the chance of happening of an
m The total number of balls in the bag = 9. Three balls can be
event is , we can also say that the odds in favour
n selected from 9 balls in 9C3 ways.
of the event are m to (n - m) or the odds against the
Three red balls can be selected from 4 red balls in 4C3 ways.
event are (n - m) to m.
Probability of picking three red balls
4C 4 1 20
Addition Theorem of Probability = =
=
9C
3
= ; P( E )
3 84 21 21
If A and B are two events, then
P(A ∪ B) = P(A) + P(B) - P(A ∩ B) Odds against the three balls being red are
This result follows from the corresponding result in set 20 1
= P= ( E ) : P( E ) : = 20 : 1.
theory. If n (X) represents the number of elements in set X, 21 21
n (X ∪ Y) = n (X) + n (Y) - n (X ∩ Y).
Example 4
Example: If a die is rolled, what is the probability that the When two dice are rolled together, find the probability of
number that comes up is either even or prime? getting at least one 4.
A = The event of getting an even number = {2, 4, 6}
Solution
B = The event of getting a prime = {2, 3, 5}
Let E be the event that at least one dice shows 4. E be the
A ∪ B = {2, 3, 4, 5, 6}
event that no dice shows 4. The number of favourable out-
A ∩ B = {2} 25
comes of E is 5 × 5 = 25 ⋅ P ( E ) =
3 3 5 1 36
P ( A) = , P ( B) = , P ( A ∪ B) = and P ( A ∩ B ) = . We
6 6 6 6 25 11
\ P( E ) = 1 − P( E ) = 1 − = .
can verify that P(A ∪ B) = P(A) + P(B) - P(A ∩ B). 36 36
Example 5
SOLVED EXAMPLES
When two dice are rolled together find the probability that
Example 1 total score on the two dice will be 8 or 9.
When a cubical dice is rolled, find the probability of getting Solution
an even integer. When two dice are rolled, the total number of outcomes
Solution = 6 × 6 = 36.
When a dice is rolled, the number of possible out comes is Favourable outcomes for getting the sum 8 or 9 are {(2, 6),
6. The number of favourable outcomes of getting an even (6, 2), (3, 5), (5, 3), (4, 4), (3, 6), (6, 3), (4, 5), (5, 4)}, i.e.,
integer is 3. the total number of favourable outcomes = 9.
3 1 9 1
The required probability= = . The required probability = = .
6 2 36 4
Example 6 n( A ∩ B)
If two cards are drawn simultaneously from a pack of cards, A n( A ∩ B) n( S ) P ( A ∩ B)
what is the probability that both will be jacks or both are \ P = = = .
B n ( B ) n ( B ) P ( B)
queens? n( S )
Solution NOTES
Here two events are mutually exclusive, P(J ∪ Q) = P(J) 1. This definition is also valid for infinite sample spaces.
4
C2 2. The conditional probability of B given A is denoted by
+ P(Q). Probability of drawing two jacks is P ( J ) =
52
C2 B B P ( A ∩ B)
P and P = .
4
A A P ( A)
C2
Probability of drawing two queens is P (Q ) = 52
C2
Multiplication Theorem
P(J ∪ Q) = P(J) + P(Q) Let A and B be two events of certain random experiment
such that A occurs only when B has already occurred. Then,
4
C2 4
C2 4
C2 2 A
= + = 2 ⋅ = . for the conditional event , the total possible outcomes are
52
C2 52 C2 52
C2 221 B
the outcomes favourable to the event B and its favourable
outcomes are the outcomes favourable to both A and B.
Example 7
A n( A ∩ B)
When two cards are drawn from a pack of cards, find the So, P =
B n( B )
probability that the two cards will be kings or blacks.
n( A ∩ B) n( S ) 1
= × = P ( A ∩ B) ×
Solution n( S ) n( B) P ( B)
4C
That is, P ⋅ P ( B) = P ( A ∩ B )
The probability of drawing two kings = 2 A
52 C
2 B
26 C
2 This is called the multiplication theorem on probability.
The probability of drawing two black cards is = 52 C
2 Example 8
2C
2 A letter is selected at random from the set of English alpha-
The probability of drawing two black kings is 52 C
2 bet and it is found to be a vowel. What is the probability
\ The required probability that it is ‘e’?
4C 26 C 2C 55 Solution
2 2 2
= 52 C
+ 52 C
− 52 C
= . Let A be the event that the letter selected is ‘e’ and B be the
2 2 2 221
event that the letter is a vowel. Then, A ∩ B = {e} and B =
{a, e, i o, u}
Conditional Probability 1
Let S be a finite sample space of a random experiment and A P ( A ∩ B ) 26 1
So, P = = = .
A, B are events, such that P(A) > 0, P(B) > 0. If it is known B P ( B) 5 5
that the event B has occurred, in light of this we wish to 26
compute the probability of A, we mean conditional proba- Independent Events In a random experiment, if A, B are
bility of A given B. The occurrence of event B would reduce
events such that P(A) > 0, P(B) > 0 and if P = P(A) or
A
the sample space to B, and the favourable cases would now
be A ∩ B. B B
P = P(B) (conditional probability equals to uncondi-
A⋅∩⋅B⋅(new⋅favourable⋅set) A
tional probability) then we say A, B are independent events.
A Sample
space If A, B are independent, P(A ∩ B) = P(A) P(B).
B
B⋅(new⋅sample⋅space) Example 9
Two coins are tossed one after the other and let A be the
Notation The conditional probability of A given B is de- event of getting tail on second coin and B be the event of
noted by P A . A
getting head on first coin, then find P .
B B
K 1. p(x) ≥ 0;
Where A is an arbitrary event of S. 2. ∑p(x) = 1.
Solution
Example 14
x
10 10 1 −10
A family consists of five children. If the random variable P ( X ≤ 10) = ∫ f ( x ) dx = ∫ e dx
(X) represents the number of boys in that family then, 0 0 10
1. Find the expected value E(X) of X. −x
10
Y n
X y1 y2 y3 … yn ∑ j =1Pxy ( x i, y )
j
∑ i =1 Pxy ( x i , y j )
m
Py(y1) Py(y2) Py(y3) … Py(yn)
From the above table, it can be easily observed that the mar- 1.
The conditional probability mass (density) function
ginal probability mass functions of X and Y namely Px(x) x
and PY(y) respectively can be obtained from the joint prob- f X of X , given Y = y is defined as
y
Y
ability mass function Pxy (x, y) as
x f ( x, y )
F X = XY , where fY(y) ≠ 0 and
Px ( xi ) = ∑
n
P (x , y j ), for i = 1, 2, … , m y fY ( y )
j =1 xy i Y
2.
The conditional probability mass (density) function
And
y y
f Y of Y , given X = x is defined as f Y
X
x x
Py ( y j ) = ∑ j =1 Pxy ( xi , y j ) for j = 1, 2, 3, … , n
m
X
f XY ( x, y )
= where f X ( x ) ≠ 0.
•• Pxy (xi, yj) ≥ 0 ∀ i, j f X ( x)
∑ i =1 ∑ j =1 Pxy ( xi , y j ) = 1
m n
•• Independent Random Variables
•• The cumulative joint distribution function of the two Two discrete (continuous) random variables X and Y defined
dimensional random variable (X, Y) is given by Fxy (x, y) on the same sample space with joint probability mass (den-
= P(X ≤ x, Y ≤ y). sity) function PXY(x, y) are said to be independent, if and
only if,
PXY(x, y) = PX(x) PY(y)
Joint Probability Density Function
Where PX(x) and PY(y) are the marginal probability mass
Let X and Y are two continuous random variables on the
(density) functions of the random variables X and Y
same sample space S with fx(x) and fy(y) as the probabil-
respectively.
ity density functions respectively. Then a function fxy(x, y)
is called the joint probability density function of the two NOTE
dimensional random variable (X, Y) if the probability that If the random variables X and Y are independent then
the point (x, y) will lie in the infinitesimal rectangular region Pxy (a ≤ X ≤ b, c ≤ Y ≤ d) = Px (a ≤ X ≤ b) Py(c ≤ Y ≤ d)
of area dx dy is fxy(x, y) dx dy,
That is,
Statistics
1 1 1 1 Statistics is basically the study of numeric data. It includes
P x − dx ≤ X ≤ x + dx, y − dy ≤ Y ≤ y + dy
2 2 2 2 methods of collection, classification, presentation, analysis
and inference of data. Data as such is qualitative or quan-
= fXY (x, y) dx dy
titative in nature. If one speaks of honesty, beauty, colour,
∞ ∞ etc., the data is qualitative while height, weight, distance,
•• ∫−∞ ∫−∞ f XY ( x, y)dxdy = 1 marks, etc are quantitative.
•• The marginal probability density functions fX(x) and fY(y) of The present course aims to systematically study statistics
the two continuous random variables X and Y are given by, of quantitative data. The quantitative data can be divided
∞ ∞ into three categories
f x ( x) = ∫ f XY ( x, y )dy and f y ( y ) = ∫ f XY ( x, y )dx
−∞ −∞ 1. Individual series
•• The cumulative joint distribution function FXY(x, y) of the 2. Discrete series and
two-dimensional random variable (X, Y) (where X and Y 3. Continuous series
are any two continuous random variables defined on the
same sample space) is given by, Individual Series
x y
Examples:
FXY ( x, y ) = ∫
−∞ ∫−∞
f XY ( x, y )dxdy. 1. Heights of 8 students
5.0, 4.9, 4.5, 5.1, 5.3, 4.8, 5.1, 5.3 (in feet)
2. The weight of 10 students
Conditional Probability Functions 46, 48, 52, 53.4, 47, 56.8, 52, 59, 55, 52 (in kgs)
of Random Variables
Discrete Series
Let X and Y be two discrete (continuous) random variables Example:
defined on the same sample space with joint probability x : Number of children in a family
mass (density) function fXY(x, y), then f : Number of families
Measures of Dispersion 2
Σxi2 Σxi
Range Alternatively σ = − is a useful formula for
n n
The range of a distribution is the difference between the
computational purpose.
greatest and the least values observed.
Some Important Results Based on Range Some Results Based on SD
1. Range is a crude measure of dispersion as it is based 1. The square of standard deviation is termed as
only on the value of extreme observations. variance.
2. It is also very easy to calculate. 2. SD is the least mean square deviation.
3. It does not depend on the frequency of items. 3. If each item is increased by a fixed constant the SD
does not alter or SD is independent of change of
Quartile Deviation (QD) origin.
Q3 − Q1 4. Standard deviation depends on each and every data
QD = item.
2
5. For a discrete series in the form a, a + d, a + 2d, . . . (AP),
Individual Series The numbers are first arranged in
n2 − 1
ascending or descending order, then we find the quartiles the standard deviation is given by SD = d ,
Q1 and Q3 as 12
where n is number of terms in the series.
Q1 → size of [(n + 1)/4]th item
Q3 → size of [3(n + 1)/4]th item
Co-efficient of Variation (CV)
The first quartile (or the lower quartile) Q1 is that value of
SD
the variable, which is such that one-quarter of the observa- Co-efficient of variation (CV) is defined as, CV = ×100.
tions lies below it. The third quartile Q3 is that value of the AM
variable, which is such that three-quarters of the observa- This is a relative measure, which helps in measuring the
tions lie below it. consistency. Smaller the co-efficient of variation, greater is
the consistency.
Mean Deviation (MD)
It is defined as the arithmetic mean of the deviation from ori- Example 17
gin, which may be either mean or median or mode. For the individual series, compute the mean, median and
Individual Series mode 8, 11, 14, 17, 20, 23, 26, 29.
| x1 − A | + | x2 − A | + + | xn − A | Solution
MD =
n Σxi 8 + 11 + + 29
Mean = x = = = 18.5
where x1, x2, …, xn are the n observations and A is the mean n 8
or median or mode. Median: As the numbers are in ascending order and the
Some Results Based on MD number 17 and 20 being middle terms.
1. Mean deviation depends on all items. 17 + 20 37
Median = = = 18.5
2. By default, mean deviation is to be computed about 2 2
mean. Mode: As no term can be regarded as ‘most often found’,
3. Mean deviation about the median is the least. mode is not-defined. However using empirical formula,
| a−b| Mode = 3 median - 2 mean
4. Mean deviation of two numbers a and b is . = 3(18.5) - 2(18.5) = 18.5.
2
Standard Deviation (SD) Example 18
Standard deviation is referred to as root mean squared devi- The arithmetic mean of 8, 14, x, 20 and 24 is 16; then find x.
ation about the mean.
Solution
Individual Series
8 + 14 + x + 20 + 24
x= = 16
( x1 − x ) 2 + ( x2 − x ) 2 + + ( xn − x ) 2 5
SD (s) =
n
⇒ x = 80 − 66 = 14.
where x1, x2, …, xn are n observations with mean as x.
Critical Region and Critical Value 3. Two tailed test: If the alternative hypothesis H1 is
Consider the area under the probability curve of the sam- of not equal to type (For example, H1 : m ≠ m0 or H1
pling distribution of the test statistic which follows some : σ 12 ≠ σ 22), then the critical region lies on both sides
known distribution. The area under the probability curve (right and left tails) of the probability curve of the test
α
is divided into two regions, namely the region of rejection statistic S* such that the critical region of area lies
where null hypothesis is rejected and the region of accept- 2
ance where null hypothesis is accepted. α
on the right tail and the critical region of area lies
2
Critical Region (or) the Region of Rejection on the left tail as shown in the figure.
(or) the Significant Region Acceptance
Region of area
The region under the probability curve of the sampling dis- 1− α
tribution of the test statistic, where the null hypothesis (H0) Critical Region
of area α /2 Critical Region
is rejected is called the critical region. of area α /2
•• The area of the critical region is equal to the level of sig-
nificance a.
Test of Hypothesis (significance) Standard Error = SE = s* Test Statistic = Z Expansions for Notations
Sample mean (x) and σ x−µ m = Population mean
population mean (m) n σ* x = Sample mean
Means of two samples x1 − x 2 s = Population standard
σ 12 σ 22
( x1 and x2) + deviation
n1 n2 σ*
Sample standard deviation σ s −σ n = Sample size
(s) and population standard 2n σ * s = Sample standard
Difference deviation (s) deviation
between Sample standard deviations (s1 σ 12 σ 22 s1 − s2 P = Population proportion
and s2) +
2n1 2n2 σ* Q=1-P
P1Q1 P2Q2 P1 − P2
Two sample proportions (P1, P2) +
n1 n2 σ*
Example 21 1.9
= × 7 = 2.333
The dean of an engineering college claims that the aver- 5.7
age attendance of students in the final semester of B.Tech is
Critical region: As the alternative hypothesis is of ≠ type,
72.5% with a standard deviation of 5.7%. To test this claim,
the test should be a two tailed test, where the critical region
the attendance of a random sample of 49 students of final
lies on both sides of the curve as shown in the figure.
semester of B.Tech were examined, which showed the aver-
age as 74.4%. Can the claim be accepted or not at 1% level a = 0.01
of significance? α
⇒ = 0.005
2
Solution \ P(Z ≤ −Zα /2 ) = 0.05
Here population mean = m0 = 72.5 ⇒ − Zα / 2 = −2.575
Simple mean = x = 74.4 ⇒ Zα / 2 = 2.575
Population standard deviation = s
= 5.7 Critical Region
Critical Region
1
Level of significance = α = = 0.01
100
Sample size = n = 49
Null hypothesis H0: m = m0 = 72.5
Alternative hypothesis H1: m ≠ m0 (= 72.5) −Zα / 2 = − 2.575 Z = 0 Zα / 2 = + 2.575
Level of significance, a = 0.01
Test statistic, Decision: It can be easily observed that the value of the test
statistic lies between -Za/2 and Za/2 i.e., The test statistic is
x − µ0 74.4 − 72.5 not in the critical region.
Z= =
σ 5.7 Hence we accept the null hypothesis H0
n 49 \ The claim of the dean can be accepted.
x1 − x2
Z=
Zα = 1.645
σ 12 σ 22
+
n1 n2
Here a = 0.05 x1 − x2
\ P(Z ≥ Za) = 0.05 i.e., Z =
s12 s22
⇒ P(Z ≤ −Za) = 0.05 +
n1 n2
⇒ −Za = −1.645
⇒ Za = 1.645 (∵ The standard deviations s1 and s2 are given and the stand-
ard deviations of populations are unknown)
\ The critical region is to the right of Za
(i.e., to the right of Z = 1.645) under the standard normal (71 − 76) −5
= =
curve. 92 14 2 9
Decision: As the value of the test statistic Z = 2.5 is greater + 4 + 4
36 49
than that of Za = 1.645, the test statistic lies in the critical
region. Z = −2
on them (or k values are already available), then the number F-Distribution
of degrees of freedom can be obtained by n - k. A random variable F is said to follow the F-distribution
The number of degrees of freedom is denoted by n. with (n1, n2) degrees of freedom, if its probability density
Example: If x1, x2, . . . .xn are the observations given, then function f(F) is given by
the number of degrees of freedom for the mean x is n
ν1
ν1
(∵ we use all values to find x ) ν 1 2 2 −1
F
The number of degrees of freedom for the variance is n - 1
f (F ) = ν 2 ,F > 0
(∵ The variance depends on the mean) (ν 1 +ν 2 ) / 2
ν ν ν
β 1 , 2 1 + 1 F
2 2 ν 2
Student’s t-Distribution (or) t-Distribution
A random variable t is said to follow the t-distribution with The graph of F-distribution is given below.
n = n - 1 degrees of freedom, n being the sample size, if its
probability density function is given by f(F) F-distribution
curve
− (υ +1) / 2
1 t2
f (t ) = 1 + ; −∞ < t < ∞
1 υ υ α
υβ ,
2 2
α α
Chi-square Distribution
If a random variable X follows chi-square distribution
−tα = t1−α O tα t (denoted as c2-distribution or c2 (n)), then the probability
density function of X is given by
ν 1
−1 (ii) t0.05 with degrees of freedom J = 17 = The value
1 in the t-table at the intersection point of the column
f (χ 2 ) = e − χ / 2 ( χ 2 ) 2 , 0 ≤ χ 2 < ∞
2
1 Solution
=
F1−0.95 ( 24,15) (i) χ 02.05 with degrees of freedom J = 16
= The value in the c2 - table at the intersection point of
1 the column headed by a = 0.05 and the row headed by
=
F0.05 ( 24,15) J = 16
1 = 26.296
= = 0.4367.
2.29 (ii) χ 02.01 with degrees of freedom J = 21
= The value in the c2 - table at the intersection point
Example 28 of the column headed by a = 0.01 and the row headed
by J = 21
Find the values of
= 38.932
(i) χ 02.05 with degrees of freedom J = 16
Similarly, we have
(ii) χ 02.01 with degrees of freedom J = 21 (iii) χ 02.10 with degrees of freedom J = 4
(iii) χ 02.10 with degrees of freedom J = 4 = 7.779.
Example 29 Solution
The highest temperature in the month of June at a certain Population mean = m0 = 40
place is normally distributed with mean 40°C. The high- 43 + 37 + 35 + 39 + 38
est temperatures in June during the last five years are 43°C, Sample mean = x =
5
37°C, 35°C, 39°C and 38°C. From this data, can we con-
clude that the average highest temperature in June during the 192
∴x = = 38.4
last five years is less than the normal highest temperature? 5
(Test at 0.05 level of significance) Sample size v n = 5
O t
tα = −1.895
tα 20.7729
with υ = 16 =
2 17.7502
= t0.01 at υ = 16 F = 1.1703
= 2.583 Critical region:
t Here u1 = n1 - 1 = 13 - 1 = 12
\ The critical region is to the left of − α = -2.583 and to
2 And u2 = n2 - 1 = 9 - 1 = 8
tα
the right of = 2.583. \ The critical value is
2
Fa(u1, u2) = F0.05(12, 8) = 3.28
t
Decision: As the test statistic t = 2.0175 lies between − α And the critical region is to the right of 3.28
2
tα As F = 1.1703 < Fa(u1, u2) = 3.28,
and , we accept the null hypothesis.
2 we conclude that the two random samples are drawn from
\ There is no significant difference between the means of two normal populations with the same variance.
marks scored by the students of sections A and B.
Non-Parametric Tests
Example 31
The variances of two samples of sizes 9 and 13 are 15.778
Goodness of Fit Test
and 19.175 respectively. Test whether the two samples be To determine, if a population follows a specified theoretical
regarded as drawn from normal populations with the same distribution such as binomial, Poisson or normal distribution,
variance at 5% level of significance. c2 test can be used. c2 test, which is based on how good a fit
is there between the observed frequencies (Oi from the sam-
Solution ple) and the expected frequencies (Ei from the theoretical
distribution) is known as ‘goodness of fit test’.
Always the sample having higher variance will be taken as
Let a distribution be given. Let Oi and Ei (I = 1, 2, 3... n)
the first sample in this test of hypothesis.
be the observed and expected frequencies of the ith class (or
\ Variance of first sample = s12 n n
= 19.175 cell) such that ∑ Oi = ∑ Ei = N = Total Frequency.
i =1 i =1
Variance of second sample = s22
Test statistic (OR) Statistic for test of ‘goodness of fit’
= 15.778
n
Sample size of first sample = n1 = 13 (Oi − Ei ) 2
= χ2 = ∑
Sample size of second sample = n2 = 9 i =1 Ei
where n is the number of class intervals or cells, in
n 13 the given frequency distribution and c2 is a random
σˆ12 = 1 s12 = × 19.175
n1 − 1 (13 − 1) variable which is very closely approximated with c2-
distribution with degrees of freedom n.
σˆ12 = 20.7729
Degrees of Freedom for Goodness of Fit Test Based on the
n 9 theoretical distribution given, the degrees of freedom are as
σˆ 22 = 2 s22 = × 15.778
n
2 − 1 ( 9 − 1) given below.
1. For uniform distribution, n = n - 1
σˆ 22 = 17.7520 2. For binomial and Poison distribution, n = n - 2
3. For normal distribution, n = n - 3
Null hypothesis: H0: σˆ12 = σˆ 22
NOTES
Alternative hypothesis: H1: σˆ12 ≠ σˆ 22 Given data should satisfy the following conditions.
Level of significance: ∝ = 0.05 1. Sample size (OR) the number of sample observa-
tions, N should be more than 50(N ≥ 50)
Test statistic: 2. If individual frequencies (Oi and/or Ei) is/are small
The test statistic is say less than 10, then combine neighbouring frequen-
cies in such a way that, they will be ≥ 10.
σˆ12 3. The number of class or cells n should be neither
F=
σˆ 21
2
too small nor too large. Generally, 4 ≤ n ≤ 16.
In general, these tables arise in two kinds of problems. Various requirements for testing of hypothesis in these
two types of problems were described in the following
1. Test for Independence.
table.
2. Test for Homogeneity.
2. If χ > χα with n degrees of freedom, then accept the alternative hypothesis H1.
2 2
Example 33 Null hypothesis: H0: Heart problem and the drinking habits
Test the hypothesis with 1% level of significance that the are independent.
heart problem is independent of drinking (alcoholic drinks) Alternative hypothesis: H1: Heart problem and the drink-
habits from the following experimental data on 200 persons. ing habits are not independent.
Non Moderate Heavy Level of significance: α = 0.01
Drinkers Drinkers Drinkers
Test statistic:
Heart problem 25 40 35
r c (O − E ) 2
No Heart problem 50 30 20
X 2 = ∑∑
ij ij
i =1 j =1
Eij
Solution
Given contingency table is 2 3 (Oij − Eij ) 2
=∑ ∑ E
Non Moderate Heavy Row j =1
i =1 ij
Drinkers Drinkers Drinkers Total
Heart Problem 25 40 35 100 ( 25 − 37.5) 2 ( 40 − 35) 2 (35 − 27.5) 2
= + +
No Heart problem 50 30 20 100 37.5 35 27.5
Column Total 75 70 55
(50 − 37.5) 2 (30 − 35) 2 ( 20 − 27.5) 2
= + +
Total number of persons = N = 200 37.5 35 27.5
The expected frequency Eij is given by = 4.1667 + 0.7143 + 2.0454 + 4.1667 + 0.7143 + 2.0454
\ c2 = 13.8528
(ith row total) × ( jth column total)
Eij = Critical region:
Total number of personss (N )
Here α = 0.01
100 × 75 r = Number of rows in the table = 2
E11 = = 37.5,
200 c = Number of columns in the table = 3
100 × 70 \ Degrees of freedom = n = (r - 1) (c - 1)
E12 = = 35,
200 = (2 - 1) (3 - 1) = 2
100 × 55 n=2
E13 = = 27.5,
200 ∴ χα2 with υ degrees of freedom.
100 × 75 = χ 02.01 with 2 degrees of freedom
E21 = = 37.5,
200 = 9.210.
100 × 70 Decision: As the test statistic = c2 = 13.8528 > χα2
E22 = = 35,
200 (= 9.210), reject the null hypothesis.
100 × 55 i.e., accept the alternative hypothesis.
E23 = = 27.5
200 \ Heart problem and the drinking habits are not independent.
Exercises
1. If eight unbiased coins are tossed together, then the 2. If A and B are two mutually exclusive and exhaustive
probability that the number of heads exceeds the num- events and the probability that the non-occurrence of A
ber of tails is 3
is , then the probability of occurrence of B is
4
31 1
(A) (B) 1 1
128 2 (A) (B)
4 2
93 57 3 1
(C) (D) (C) (D)
256 256 4 16
on them and the remaining cards have the letter ‘N’ 20. If two events A and B are such that P ( A) = 0.4, P(B) =
printed on them. If three cards are picked up one after
the other at random, and placed on a table in that order, B
0.7 and P(A ∩ B) = 0.2, then P is
then what is the probability that the word formed will A∪ B
be ‘MAN’? 3 2
(A) (B)
5 1 5 5
(A) (B)
44 22 1 4
3 (C) (D)
3 4 5
(C) (D)
22 44 21. A cinema historian noted that for a brief period, all mov-
15. A and B pick a card at random from a well shuffled ies released were either directed by Nolan or starred
pack of cards, one after the other replacing it every time Bale. Also no movie directed by Nolan starred Bale.
till one of them gets a spade. The person who picks a The probability that a movie was directed by Nolan
spade is declared the winner. If A begins the game, then is 0.5, and the probability that a movie starred Bale is
the probability that B wins the game is 0.5. The probability that a movie is a hit if directed by
5 4 Nolan is 0.6, while the probability that a movie is a hit
(A) (B) given that Bale acted in it is 0.4. Given that a movie is
9 9
a hit, find the probability that it is directed by Nolan.
3 4 (A) 0.4 (B) 0.5
(C) (D)
7 7 (C) 0.6 (D) 0.7
16. A number is randomly chosen from the numbers 10 to 22. Probability mass function of a variate x is as follows:
99. It is observed that the sum of the digits of the num-
x 0 1 2 3 4
ber is ten. Find the probability that it is divisible by five.
1 1 P(X = x) k 2k 3k 4k 5k
(A) (B)
9 3
then P(x ≥ 3) =
1 2 1 4
(C) (D) (A) (B)
2 9 3 15
17. An unbiased coin is tossed a person gets `30 if the coin
shows head, and he loses `15 if the coin shows tail. If 3 5
(C) (D)
three coins are tossed, the probability that the person 5 7
gets `45 is 23. The expected number of trials required to open a door
3 1 using a bunch of n keys of which only one is the correct
(A) (B)
8 2 key is
n n −1
1 1 (A) (B)
(C) (D) 2 2
10 25
n +1
18. What is the probability of getting at least 6 heads when (C) (D) n
2
a coin is tossed 7 times if it is known that there are at
least 5 heads?
Direction for questions 24 and 25:
5 8
(A) (B) A variate x has the probability distribution as
29 29
x 4 8 12
9
(C) (D) None of these P(X = x) 1 3 1
29 3 15
5
3 6 1 24. Values of E(x) and E(x2) respectively are
19. If P(A) = , P(BC) = and P(A ∩ B) = , then find
5 7 4
104 160 102 150
Ac (A) , (B) ,
P c . 15 3 15 3
B
21 160 104 151
17 71 (C) , (D) ,
(A) (B) 3 5 15 3
60 120
2 5. The value of E[(3x + 2)2] is ________.
19 29 (A) 675.2 (B) 560.2
(C) (D)
60 60 (C) 134.56 (D) 567.2
(C) 1 - 7e 6
(D) 1 - 7e–6 4
4 8 4
15 15 5
32. The expected value of a random variable with uniform
distribution over the interval (2, 5) is 1 2
P(Y1 = yi)
3 3
1
(A) 2 (B) 2
2 Let X2 and Y2 be two discrete random variables with
1 1 joint probability mass function given as follows:
3 (D)
(C) 4
2 2 Y2
0 4 7 P(X2 = xi)
X2
1
33. If X is a continuous random variable with PDF f(x) = 1 3 1 3
4 -1
if - 2 ≤ x ≤ 2 and f(x) = 0 elsewhere, the mean of X is 7 14 14 7
______. 4 2 2 4
3
(A) 1 (B) 1.5 21 7 21 7
(C) 2 (D) 0 1 1 1
P(Y2 = yi)
34. If X is a uniformly distributed random variable in [1, 4] 3 2 6
3
Which of the following statements is TRUE about the
then P x > is
2 random variables X1, X2, Y1 and Y2?
(A) Only X1 and Y1 are independent. 49. If the standard deviation of 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11
(B) Only X2 and Y2 are independent. is M, then the standard deviation of 101, 102, 103, 104,
(C) X1 and Y1 are independent as well as X2 and Y2 are … and 111 is
independent. (A) M (B) 100 + M
(D) Neither X1 and Y1 are independent nor X2 and Y2 (C) 100 - M (D) M - 100
are independent. 50. If the standard deviation of 10, 20, 30, 40 and 50 is S,
39. If X and Y are two independent random variables with then the standard deviation of 20, 30, 40, 50 and 60 is
expectations 3 and 4 respectively. Then the expectation (A) S (B) S + 10
of X Y is (C) S - 10 (D) 10S
(A) 1 (B) 7 51. The arithmetic mean of five observations is 6.4 and the
(C) 12 (D) 16 variance is 8.24. If three of the observations are 3, 4, 8,
40. If X and Y are two independent random variables that then find the other two observations.
are uniformly distributed over the same interval [2, 5] (A) 6, 11 (B) 10, 7
(C) 8, 9 (D) 5, 12
11 11
then P X ≤ , Y ≥ is 52. The director of a sports academy claims that the aver-
4 3
age height of sports persons in their academy is more
1 2 than 170 cms. A random sample of 40 sports persons
(A) (B)
9 9 of that academy has an average height of 174 cm with
1 4 a standard deviation of 15 cm. Then the claim of the
(C) (D) director can be accepted with ______.
3 7
41. The mean of cubes of first 10 natural numbers is (A) both 1% as well as 5% levels of significance
(A) 305 (B) 300 (B) 5% level of significance but not with 1% level of
(C) 302.5 (D) 310 significance
(C) neither 5% nor 1% levels of significance
42. The mean of 25 observations was found to be 38. It was (D) no level of significance
later discovered that 23 and 38 were misread as 25 and
36, then the mean is 53. In a survey conducted on the increase in pay packages
(A) 32 (B) 36 of male and female managers across IT industry by tak-
(C) 38 (D) None of these ing a random sample of 32 male managers and another
random sample of 36 female managers, the following
43. If 3, 2 and 9 occur with frequencies 2, 5 and 3 respec- information was derived.
tively, then their arithmetic mean is
Average Standard
(A) 4.3 (B) 5 Sample Increases in Deviation of
(C) 6 (D) 4.8 Size Pay Package/ Increase in Pay
44. The median of first ten prime numbers is Annum Package/Annum
(A) 11 (B) 13 Male 32 20% 4%
(C) 12 (D) 10 Managers
Female 36 17% 3%
45. If the mean of a set of 12 observations is 10 and another Managers
set of 8 observations is 12, then the mean of combined
The standard error of the difference between the aver-
set is
age increase in pay packages is ______
(A) 12.6 (B) 10.8
(C) 12.8 (D) 10.6 3 3
(A) (B)
46. The mode of a distribution of 13 and its mean is 4 then 2 2
its median is 3 3
(A) 7 (B) 9 (C) (D)
2 4
(C) 8 (D) 11
47. Consider the non-decreasing series of the numbers, 1, 54. Match the following:
8, 8, 13, 14, 14, x, y, 18, 20, 31, 34, 38 and 40. If the List I List II
median of the series is 15, then the mode of the series is P. Standard error 1. Level of significance
(A) 14 (B) 16 Q. Type I error 2. Standard deviation of the sam-
pling distribution of a statistic
(C) 18 (D) Cannot be determined
R. Type II error 3. Accepting the null hypothesis
48. The standard deviation of 5, 5, 5, 5, 5, 5, 5, 13 is when it should be rejected
2 2 (B)
(A) 6 S. Size of a test 4. Rejecting the null hypothesis
when it should be accepted
(C) 5 (D) 7
each container. The probability that one of the ball is that there will be less than 4 penalties in a day is
Red and the other is Blue will be [GATE, 2011] ______. [GATE, 2014]
1 9 12. If {x} is a continuous, real valued random variable
(A) (B) defined over the interval (-∞, +∞) and its occurrence
7 49
is defined by the density function given as f(x) =
12 3
(C) (D)
49 7 1 1 x −a
2
⋅
where ‘a’ and ‘b’ are the statistical
2π * b e 2 b
6. In an experiment, positive and negative values are
equally likely to occur. The probability of obtaining at attributes of the random variable {x}. The value of the
most one negative value in five trials is a 1 x −a
2
1
[GATE, 2012] integral ∫ 2π * b
e 2 b dx is [GATE, 2014]
1 2 −∞
(A) (B)
32 32 (A) 1 (B) 0.5
3 6 π
(C) (D) p (D)
(C) ⋅
32 32 2
7. The annual precipitation data of a city is normally 13. Consider the following probability mass function
distributed with mean and standard deviation as (pmf) of a random variable X:
1000 mm and 200 mm, respectively. The probability
q if X = 0
that the annual precipitation will be more than 1200
mm is [GATE, 2012] p( x, q) = 1 − q if X =1
(A) <50% (B) 50% 0 othewise
(C) 75% (D) 100%
If q = 0.4, the variance of X is ________.
8. Find the value of l such that the function f (x) is a [GATE, 2015]
valid probability density function
14. The probability density function of a random vari-
F(x) = l(x - 1) (2 - x) for 1 ≤ x ≤ 2
able, x is
= 0 otherwise [GATE, 2013]
x
9. A fair (unbiased) coin was tossed four times in suc- f(x) = (4 - x2) for 0 ≤ x ≤ 2
4
cession and resulted in the following outcomes;
= 0 otherwise
(i) Head, (ii) Head (iii) Head (iv) Head. The probabil-
The mean, µx of the random variable is ______.
ity of obtaining a ‘Tail’ when the coin is tossed again
[GATE, 2015]
is [GATE, 2014]
15. X and Y are two random independent events. It is
1
(A) 0 (B) known that P(X ) = 0.40 and P ( X ∪ Y C ) = 0.7. Which
2
one of the following is the value of P ( X ∪ Y ) ?
4 1 [GATE, 2016]
(C) (D)
5 5 (A) 0.7 (B) 0.5
10. The probability density function of evaporation E on (C) 0.4 (D) 0.3
any day during a year in a watershed is given by 16. Probability density function of a random variable X is
1 given below
0 ≤ E ≤ 5 mm/day
F(E) = 5 0.25 if 1 ≤ x ≤ 5
0, otherwis f ( x) =
0 otherwise
The probability the E lies in between 2 and 4 mm/day P ( x ≤ 4) is ⋅ [GATE, 2016]
in a day in the watershed is (in decimal).
[GATE, 2014] 3 1
(A) (B)
11. A traffic office imposes on an average 5 number of 4 2
penalties daily on traffic violators. Assume that the
number of penalties on different days is independent 1 1
(C) (D)
and follows a Poisson distribution. The probability 4 8
17. If f(x) and g(x) are two probability density functions, (C) Mean of f (x) and g(x) are different; Variance of
x f (x) and g(x) are same.
a + 1 : −a ≤ x < 0 (D) Mean of f (x) and g(x) are different; Variance of
f (x) and g (x) are different.
f ( x) = x
− a + 1 : 0 ≤ x ≤ a 18. The spot speeds (expressed in km/h) observed at a
road section are 66, 62, 45, 79, 32, 51, 56, 60, 53 and
0 : otherwise
49. The median speed (expressed in km/h) is ______.
x (Note: answer with one decimal accuracy)
− a : −a ≤ x < 0 [GATE, 2016]
g( x) = x 1 9. Type II error in hypothesis testing is [GATE, 2016]
a :0 ≤ x ≤ a (A) acceptance of the null hypothesis when it is false
and should be rejected.
0 : otherwise
(B) rejection of the null hypothesis when it is true
Which one of the following statement is true? and should be accepted.
[GATE, 2016] (C) rejection of the null hypothesis when it is false
(A) Mean of f(x) and g(x) are same; Variance of f(x) and should be rejected.
and g(x) are same. (D) acceptance of the null hypothesis when it is true
(B) Mean of f(x) and g(x) are same; Variance of f(x) and should be accepted.
and g(x) are different.
Answer Keys
Exercises
1. C 2. C 3. B 4. A 5. D 6. D 7. C 8. D 9. A
10. (i) D (ii) A (iii) B 11. C 12. A 13. B 14. B 15. C 16. A 17. A
18. B 19. B 20. B 21. C 22. C 23. C 24. A 25. D 26. A 27. B
28. D 29. D 30. B 31. D 32. C 33. D 34. C 35. C 36. A 37. C
38. C 39. C 40. A 41. C 42. C 43. A 44. C 45. B 46. A 47. D
48. D 49. A 50. A 51. A 52. B 53. A 54. B 55. D 56. B 57. A
58. A 59. D 60. C 61. A
CHAPTER HIGHLIGHTS
value
SOLVED EXAMPLES 3 3
f ′( x ) = ⇒ f ′(0) =
2 1 + 3x 2
Example 1
π 11 −9 −9
If the number = 0.785398163 is approximated by , f ′′( x ) = 3
⇒ f ′′(0) =
4 14 4
then 4(1 − 3 x ) 2
(i) Find the number of digits upto which, this approxima- 27 1
tion is accurate. and f ′′′( x ) = × 5
8
(ii) Find the absolute and the percentage errors. (1 + 3 x ) 2
\ Substituting these in (2), we get
Solution
3 9 x2
π f ( x) = 1 + 3x ≈ 1 + × −
Given = 0.785398163 2 4 2!
4
π 11 3 9
Let x = = 0.785398163 (Exact value) and x = = 1 + × − x2
4 14 2 8
π x3
= 0.785714285 Approximate value of Truncation error = f ′′′(0)
4 3!
π 11 1 Max x3
(i) | x − x | = − ≤
4 14 3! 0 ≤ x ≤ 1
= | 0.785398163 - 0.785714285 |
27
= 3.16122 × 10-4 Max
0 ≤ x ≤1 8(1 + 3 x )
5 / 2
≤ 5 × 10-4
11 π 1 27
\ The approximation to is accurate upto three = (1)
3! 8
⋅
14 4
significant digits
= 0.5625.
(ii) Absolute error = | x − x |
11 π Methods for Finding the Real Roots
− = 3.16122 × 10 −4
=
14 4 (Zeros) of f (x) = 0
|x−x | The equation of the form f (x) = 0 is called an Algebraic
Percentage error = × 100 = 0.04%. or Transcendental according as f (x) is purely a polynomial
x
in x or contains some other functions such as exponential,
logarithmic and trigonometric functions etc.
Example 2
Examples:
Using the Taylor’s series expansion about x = 0, find a s econd
degree polynomial approximation to f ( x ) = 1 + 3x . Also x9 + 8x5 - 4x3 - 11x + 3 = 0 → Algebraic equation
1.
find the maximum error for this approximation when 10x4 - log(x2 - 3) + e-xsin x + tan2 x = 0 → Transcen
2.
x ∈ [0, 1]. dental equation
In this chapter, we obtain the solution of an equation f (x)
Solution
= 0, i.e., we mean to find the zeros of f (x).
We know that the Taylor’s series expansion of f(x) about x We shall now discuss few methods to find the real roots
= 0 is of both algebraic (with numerical coefficients) and tran-
x2 scendental equations.
f ( x ) = f (0) + xf ′(0) + f ′′(0) We first find an approximate value of the root of the given
2!
equation and then successively improve it to some desired
x3
+ f ′′′(0) + + ∞ (1) degree of accuracy.
3! We start with an initial approximate value, say x0, and
\ Considering the terms upto second degree, we have then find the better approximations successively x1, x2, x3 …,
x2 xn by repeating the same method.
f ( x ) ≈ f (0) + x f ′(0) + f ′′(0) (2) If the successive approximations at each step of a method
2!
approach the root more and more closely, we say that the
Here f ( x ) = 1 + 3 x ⇒ f (0) = 1 method converges.
NOTE
In bisection method, the convergence is very slow but Convergence of Regula Falsi Method
definite. The order of convergence is linear or of first order. The order of convergence of the method of false position is
greater than 1.
3
= = 0.6557.
4.5747 f(x 0)
Newton–Raphson Method 1 a
= (b − 1) xn + b −1
Let x0 be the approximate root of f (x) = 0 and let x1 = x0 + h b xn
be the correct root. Then f (x1) = 0
⇒ f (x0 + h) = 0 (1) Newton’s Iterative Formula to Find a
Expanding Eq. (1) using Taylor’s theorem, Reciprocal of a Number N
We get f (x0) + hf ′ (x0) + … = 0 The iterative formula is given by
− f ( x0 ) xi+1 = xi (2 - xiN)
⇒ h= ,
f ′( x0 )
Example 8
f ( x0 )
\ x1 = x0 − Find a real root of the equation -4x + cos x + 2 = 0, by
f ′( x0 )
Newton–Raphson method upto four decimal places assum-
Now x1 is the better approximation than x0. Proceeding ing x0 = 0.5
this way, the successive approximations x2, x3,…, xn+1 are
f ( xn ) Solution
given by xn +1 = xn −
f ′( xn ) Let f (x) = -4x + cos x + 2 and
This is called Newton–Raphson formula. f ′(x) = -4 - sin x
Example 9
Least Squares Approximation
Obtain the cube root of 14 using Newton–Raphson method,
with the initial approximation as 2.5. Least squares approximation method is one of the best
methods available for curve fitting.
Solution Let (x1, y1), (x2, y2), …, (xn, yn) be the pairs of observed
We know that, the iterative formula to find b
a is set of values of x and y. Let y = f (x) be the functional rela-
tionship sought between the variables x and y where f (x)
1 a consists of some unknown parameters. We need to find the
xn +1 = (b − 1) xn + b −1
b xn relationship y = f (x) by using the observed values.
Here b = 3 and a = 14 and let x0 = 2.5 Procedure
1 14 1. Find the residual di = yi - f (xi) (i = 1, 2, …, n) for
∴ x1 = 2 x0 + 2
3 x0 every pair of observed value yi and f(xi), the value of
1 14 the functional relation f (x) at x = xi
x1 = 2( 2.5) +
3 ( 2.5) 2 2. Find the sum of the squares of residuals corresponding
to all pairs of values of yi and f(xi) and let it be S
1 14 1
= 5 + = {5 + 2.24} = 2.413 n
3 6.25 3 i.e., S = ∑ ( yi − f ( x )) 2 .
i =1
1 14
x2 = 2( 2.413) + 3.
Find the values of the parameters in f(x) such that S is
3 ( 2.413) 2 minimum.
1 14
= 4.826 + . Fitting a Straight Line Let y = a + bx be a straight line to
3 5.822569 . be fitted to the data (x1 y1), (x2 y2), …, (xn yn).
1 \ Residual = di = yi - (a + bxi), i = 1, 2, …, n
= {4.826 + 2.4044} = 2.410
3 Sum of the squares of the residuals = S = ∑(yi - (a + bxi ))2
\ The approximate cube root of 14 is 2.41. Now we have to find the parameters a and b such that S is
minimum
Example 10
∂S
Find the reciprocal of 24 using Newton–Raphson method = ∑ 2[ yi − ( a + bxi )( −1)] and
∂a
with the initial approximation as 0.041.
∂S
= ∑ [2( yi − ( a + bxi ))( − xi )]
Solution ∂b
1 For S to be minimum,
The iterative formula to find is,
N
∂S ∂S
= 0 and =0
xn+1 = xn(2 - xnN) ∂a ∂a
⇒ ∑ [−2( yi − ( a + bxi ))] = 0 and The required values in the normal equations can be found
using the following table
∑ [2( yi − (a + bxi ))(− xi )] = 0 xi yi xi yi xi2
⇒ ∑ yi = n a + b∑ xi (1) 1 4 4 1
2 11 22 4
and ∑ xi yi = a∑ xi + b∑ xi 2 (2)
3 35 105 9
Eqs. (1) and (2) are known as normal equations.
4 100 400 16
By solving these equations, we get the values of ‘a’ and b
Solution ∑ xi = 0; ∑ yi = 76; ∑ xi yi = 91; ∑ xi2 = 10; ∑ xi3 = 0; ∑ xi4 = 34;
We have to fit the line y = a + bx to the given data. ∑ xi2 yi = 183
The normal equations are
Substituting these values in the normal equation we have
∑ yi = n a + b ∑ xi (1)
76 = 5a + b × 0 + c × 10
∑ xi yi = a ∑ xi + b ∑ xi2 (2) 91 = a × 0 + b × 10 + c × 0
NOTES and D4 y0 = D3 y1 - D3 y0
= (y4 - 3y3 + 3y2 - y1) - (y3 - 3y2 + 3y1 - y0)
1. Newton’s forward interpolation formula is
used to interpolate (estimate) the values y = \ D4 y0 = y4 - 4y3 + 6y2 - 4y1 + y0
f(x) near the beginning of the set of tabulated From the above discussion, one can observe that the coef-
values given or for estimating the value of y = ficients of yi s are nothing but the binomial coefficients with
f(x) to the left of the beginning. positive and negative signs alternatively.
2. Newton’s backward interpolation formula is \ Dn y0 = yn - nC1 yn-1 + nC2 yn-2 + … + (-1)ny0 .
used to interpolate (estimate) the values y =
f(x) near the end of the set of tabulated values Example 14
or for estimating the values of y = f(x) to the A function y = f (x) is given by the following table
right of the last tabulated value yn.
x 5 10 15 20 25
When the values x0, x1, x2, …, xn of x are not equally y = f(x) 31 42 51 62 76
spaced, then we can’t make use of Newton’s forward
as well as backward interpolation formulae. In such Using Newton’s forward interpolation formula, find the
situations, the following two interpolation formulae value of y at x = 7.
will be helpful.
Solution
3. Newton’s divided difference formula: If the function
First let us form the forward difference table:
y = f(x) takes the values y0, y1, y2, …, yn corresponding
to the values x0, x1, x2, …, xn (need not be equally x y = f(x) Dy D2y D3y D4y
spaced) of x, then the Newton’s divided difference 5 31
interpolation polynomial is given by 11
10 42 -2
y (x) = f (x) = y0 + (x - x0) [x0, x1] + (x - x0) (x - x1) 9 4
[x0, x1, x2] +…+ (x - x0)(x - x1)…(x - xn-1) [x0, x1, … 15 51 2 -3
11 -1
xn]
20 62 3
4. Lagrange’s interpolation formula: If the function y 14
= f (x) takes the values y0, y1, …, yn respectively at the 25 76
points x0, x1, x2, …, xn (need not be equally spaced)
of x, then the Lagrange’s interpolation polynomial is By Newton’s forward difference formula,
given by y (x) = f(x) = p( p −1) 2 p( p − 1)( p − 2) 3
yp = y0 + pD y0 + D y0 + D y0 +
( x − x1 )( x − x2 ) ( x − xn ) 2! 3!
y0 p( p − 1)( p − 2)( p − 3) 4
( x0 − x1 )( x0 − x2 ) ( x0 − xn ) ∆ y0 + (1)
4!
( x − x0 )( x − x2 ) ( x − xn )
+
y1 Here x0 = 5; h = 5 and x = 7
( x1 − x0 )( x1 − x2 ) ( x1 − xn )
x − x0 7 − 5
( x − x0 )( x − x1 ) ( x − xn −1 ) ∴p= = = 0.4
+ +
yn . h 5
( xn − x0 )( xn − x1 ) ( xn − xn −1 )
Substituting these in Eq. (1), we have
y(7) = 31 + (0.4) × 11 +
Example 13
(0.4)(0.4 − 1) (0.4)(0.4 − 1)(0.4 − 2)
If D denotes the forward difference operator, then show that × ( −2) +
2! 3!
D4 y0 = y4 - 4y3 + 6y2 - 4y1 + y0
(0.4)(0.4 − 1)(0.4 − 2)(0.4 − 3)
Hence find the general expression for Dny0 in terms of y0, y1, ×4+ × ( −3)
4!
…, yn that does not involve the difference operators.
= 31 + 4.4 + 0.24 + 0.256 + 0.1248 = 36.0208
Solution \ The value of y at x = 7 is 36.0208.
We know that Dy0 = y1 - y0 Example 15
D2 y0 = Dy1 - D y0 Following table shows the values of a function y = f(x) at 0,
= (y2 - y1) - (y1 - y0) 2, 5 and 9
= y2 - 2y1 + y0 x 0 2 5 9
D y0 = D y1 - D y0
3 2 2 y = f(x) 6 15 27 40
= (y3 - 2y2 + y1) - (y2 - 2y1 + y0) = y3 - 3y2 + 3y1 - y0 Using the Lagrange’s interpolation formula, find y(6).
Solution dy 1 ( 2 p − 1) 2
∴ = ∆y0 + ∆ y0
Given values of y = f (x) are dx h 2!
x 0 2 5 9 (3 p 2 − 6 p + 2) 3
+ ∆ y0 +
Y = f(x) 6 15 27 40 3!
By Lagrange’s interpolation formula, we have dy
∴ at x = x0
( x − x1 )( x − x2 )( x − x3 ) dx
f ( x) = y0
( x0 − x1 )( x0 − x2 )( x0 − x3 ) 1 1 1 1 1
= ∆y0 − ∆ 2 y0 + ∆3 y0 − ∆ 4 y0 + ∆50 …
( x − x0 ( x − x2 )( x − c3 ) h 2 3 4 5
+ y1
( x1 − x0 )( x1 − x2 )( x1 − x3 ) ( At i = x0; p = 0)
∴
( x − x0 )( x − x1 )( x − x3 ) d2 y d dy d dy dp 1
+ y2
And = = =
( x2 − x0 )( x2 − x1 )( x2 − x3 ) dx 2 dx dx dp dx dx h
( x − x0 )( x − x1 )( x − x2 )
+ y3 2 6p−6 3 12 p 2 − 36 p + 22 4 1
( x3 − x0 )( x3 − x1 )( x3 − x2 ) ∆ y0 + ∆ y0 + ∆ y0
3! 4! h
( x − 2)( x − 5)( x − 9) ( x − 0)( x − 5)( x − 9)
∴ f ( x) = ×6 + × d2 y 1 2 6p−6 3 12 p 2 − 36 p + 22
(0 − 2)(0 − 5)(0 − 9) ( 2 − 0))( 2 − 5)( 2 − 9) 2 = 2 ∆ y0 + ∆ y0 +
dx h 3! 1!
( x − 0)( x − 2)( x − 9) ( x − 0)( x − 2)( x − 5)
15 + × 27 + × 40
(5 − 0)(5 − 2)(5 − 9) (9 − 0)(9 − 2)(9 − 5) ∆ 4 y0
(6 − 2)(6 − 5)(6 − 9)
Now, y(6) = f(6) = d2 y
( −2) × ( −5) × ( −9)
2
dx at x = x0
(6 − 0)(6 − 5)(6 − 9) (6 − 0)(6 − 2)(6 − 9)
×6+ × 15 + × 1 2 11 4 5 5
2 x( −3) × ( −6) −5 x − 3 x( −4) = ∆ y0 − ∆ y0 + 12 ∆ y0 − 6 ∆ y0 +
3
h2
(6 − 0)(6 − 2)(6 − 5) 4 45 162 80
2+ × 40 = − + + = 30.5807. \ By using Newton’s forward interpolation formula,
9× 7× 4 5 7 5 21 the first and second derivatives of y = f (x) at x = x0 are
dy
given by
Numerical Differentiation dx x = x0
In numerical differentiation, we find the derivatives by using 1 1 1 1 1
the interpolation formulae. = ∆y0 − ∆ 2 y0 + ∆3 y0 − ∆ 4 y0 + ∆5 y0 …
h 2 3 4 5
1. Derivatives using newton’s forward difference
interpolation formula: We know that the Newton’s and
forward difference interpolation formula is
p( p −1) 2 p( p − 1)( p − 2) d2 y 1 2 11 4 5 5
y = y0 + pD y0 + D y0 + 2 = 2 ∆ y0 − ∆ y0 + 12 ∆ y0 − 6 ∆ y0 + …
3
2! 3! dx x = x0 h
D 3y0 +…
Differentiating both sides wrt p, 2. Derivatives using newton’s backward difference
interpolation formula: We know that the Newton’s
dy 2 p −1 2 3 p2 − 6 p + 2 3 backward difference interpolation formula is y = yn +
= ∆y0 + ∆ y0 + ∆ y0 +
dp 2! 3! p( p +1) 2 p( p + 1)( p + 2)
p∇ yn + ∇ yn + + ∇3yn …
x − x0 dp 1 2! 3!
As p = = =
h dx h Differentiating on both sides wrt p,
dy dy dp ( 2 p − 1) 2 dy 2 p +1 2 3 p2 + 6 p + 2 3
Now = = [∆y0 + ∆ y0 = ∇yn + ∇ yn + ∇ yn +
dx dp dx 2! dx 2 3!
3 p2 − 6 p + 2 1 x − xn dp 1
+ ∆y0 + ] As p = =
3! h h dx h
dy dy dp Solution
∴ =
dx dp dx As we have to find the first two derivatives of y = f (x) at x
= 7, (end point of the given data), we use the derivatives’
2 p +1 2 3 p2 + 6 p + 2 3 1 formulae obtained from Newton’s backward interpolation
= ∇yn + ∇ yn + ∇ yn + formula.
2 ! 3 ! h
The backward difference table for the given data is
dy 1 2 p +1 2 3 p2 + 6 p + 2
∴ = ∇yn + ∇ yn + X y = f(x) Dy D2y D3y D4y D5y
dx h 2! 3!
2 4
4
3 ∇ yn + 3 8 3
7 -2
4 15 5 6
dy
dx 12 -8 2
x = xn 5 27 -3 8
9 0
1 1 1 1 1 6 36 -3
= ∇yn + ∇ 2 yn + ∇3 yn + ∇ 4 yn + ∇5 yn + ... 6
h 2 3 4 5
7 42
( At x = xn; p = 0)
∴
x y = f(x) First divided differences Second divided differences Third divided differences Fourth divided differences
5 100
294 − 100
= 97 202 − 97
7−5 = 21
7 294 10 − 5 27 − 21
900 − 294 =1
= 202 11 − 5
10 − 7 310 − 202 1−1
10 900 = 27 =0
11 − 7 33 − 27 13 − 5
1210 − 900 =1
= 310
11 1,210 11 − 10 13 − 7
409 − 310
= 33
2028 − 1210 13 − 10
= 409
13 2,028 13 − 11
y0 y1 y2 3h
yn−1 yn = ( y0 + 3 y1 + 3 y2 + y3 )
8
Similarly,
x
O x0 x1 x2 xn−1 xn
x6
3h
Geometrically, the curve y = f (x) is replaced by n straight ∫ f ( x )dx =
8
( y3 + 3 y4 + 3 y5 + y6 ) and so on.
line segments joining the points (x0, y0) and (x1, y1); (x1, y1) x3
and (x2, y2); …, (xn-1, yn-1) and (xn, yn). The area bounded by
the curve y = f (x) is then approximately equal to the sum of Adding all these integrals from x0 to xn where ‘n’ is a
the areas of n trapeziums as shown in the figure. multiple of 3, we get
xn
3h
∫ f ( x )dx = [(y0 + yn) + 3(y1 + y2 + y4 + y5 + …+ yn-2)
Simpson’s One-third Rule x0
8
[Three-point Quadrature] + 2(y + y + y + … + y )]
3 6 9 n-3
Substituting n = 2 in the Newton–Cotes quadrature formula
taking the curve through (x0, y0), (x1, y1) and (x2, y2) as a 3
The above rule is called Simpson’s rule which is
parabola, so that the differences of order higher than 2 8
becomes zero, we get applicable only when ‘n’ is a multiple of 3.
x0 + 2 h
1 Example 18
∫ f ( x )dx = 2h y0 + ∆y0 + ∆ 2 y0
6 2
x0
Evaluate: ∫ 1 + x 2 dx taking h = 0.2 using
h 0
= ( y0 + 4 y1 + y2 )
3 (i) Trapezoidal rule and
Similarly, 1
(ii) Simpson’s rd rule
x0 + 4 h
h 3
∫ f ( x )dx =
3
( y2 + 4 y3 + y4 )
Solution
x0 + 2 h
Here, a = 0, b = 2, h = 0.2
x0 + nh
h
∫ f ( x )dx =
3
( yn − 2 + 4 yn −1 + yn )
So, n =
b−a 2−0
= = 10
x0 + ( n − 2 ) h h 0.2
1 1.0198 1.077 1.1661 1.2806 1.414 1.562 1.7204 1.8867 2.059 2.236
y = 1 + x2 y0 y1 y2 y3 y4 y5 y6 y7 y8 y9 y10
(i)
By Trapezoidal rule, By Simpon’s three-eighth rule,
π /2
2 3h
∫ 1+
x 2 dx
h
= [( y0 + y10 ) + 2( y1 + y2 + + y9 ) ] ∫ e cos x dx =
8
[(y0 + y3) + 3(y1 + y2)]
0
2 0
3 π
= × [(2.718 + 1) + 3(2.3774 + 1.6487)]
0.2 8 6
= [(1 + 2.236) + 2(1.0198 + 1.077 + 1.1661 +
2 π
= [(3.718) + (12.0783)] = 3.10159.
⋅
2 Single-step Methods
∫ + x 2 dx
1 1.
Taylor’s series method
0
2.
Picard’s method of successive approximation
h ( y0 + y10 ) + 4( y1 + y3 + y5 + y7 + y9 ) Multi-step Methods
=
3 + 2( y2 + y4 + y6 + y8 ) 1.
Euler’s method
0.2 2.
Modified euler’s method
=
[(1 + 2.236) + 4(1.0198 + 1.1661 + 1.414
3 3.
Runge–Kutta method
+ 1.7204 + 2.059) + 2(1.077 + 1.2806 + 1.562 + 4.
Predictor–Corrector methods
[Milne’s and Adam’s]
1.8867)]
∴ In general, Solution
h2 f (x, y) = 1 + xy
y n+1 = yn + hyn′ + y ″ +… will be the iterative formula. x0 = 0, y0 = 1
2! n
x
Example 22 The first approximation y1 = y0 + ∫ f ( x, y0 )dx
dy x0
Given = x - y2 with the initial condition y(0) = 1 x x
dx = 1 + ∫ 1 + xdx
1+ ∫ 1 + xy0 dx
Find y (0.1) using Taylor series method with step size 0.1. x0 0
x2
Solution y1 = 1 + x +
2
f (x, y) = x - y2 (0.1) 2
x = 0.1, x0 = 0, y0 = 1, h = 0.1 At x = 0.1, y1 = 1 + (0.1) + = 1.105
2
y′ = x - y2 ⇒ y′(0) = x0 - y02 = - 1; The second approximation y2,
y″ = 1 - 2yy′ ⇒ y″(0) = 1 - 2y0y0′ x
dy x0
Given the differential equation = f (x, y)(1) x
dx
⇒ y3 = 1 + ∫ (1 + xy2 )dx.
Integrate Eq. (1) from x0 to x, we get
0
x x
x
⎛ ⎡ x 2 x 3 x 4 ⎤⎞
∫ dy = ∫ f ( x, y)dx 1+ ∫ ⎜⎝1 + x ⎢⎣1 + x + + + ⎥⎟ dx
x0 x0
0
2 3 8 ⎦⎠
x
x
⎛ x3 x 4 x5 ⎞
⇒ y( x ) − y( x0 ) = ∫ f ( x, y)dx 1+ ∫ ⎜⎝1 + x + x 2 + + + ⎟ dx
x0
0
2 3 8⎠
x
x2 x3 x 4 x5 x6
⇒ y( x ) = y( x0 ) + ∫ f ( x, y)dx (2) =1+x+ + + + +
x0 2 3 8 15 48
Put y = y0, we get the first approximation, At x = 0.1,
x
(0.1) 2 (0.1)3 (0.1) 4 (0.1)5 (0.1)6
y3 = 1 + (0.1) +
∫
yn = y0 + f ( x, yn −1 )dx.
2
+
3
+
8
+
15
+
48
x0
(0.1) 2 (0.1)3 (0.1) 4 (0.1)5 (0.1)6
Example 23 + + + +
2 3 8 15 48
dy = 1 + 0.1 + 0.005 + 0.0003 + 0.0000125 +
Given = 1 + xy and y (0) = 1. Evaluate y (0.1) by Picard’s
dx 0.0000006 + 0.00000002
method upto three approximations. y3 = 1.105313
Exercises
1. Three of the roots of the equation x4 + lx3 + mx2 + nx + 3. If the equation x6 + 5x5 + 11x4 + 34x2 + 20x + 24 = 0 has
24 = 0 are 3, 1 and -2. Which of the following could be exactly four non-real roots, then the number of negative
the value of l + m - n? roots is
(A) 0 (B) 1 (A) 1 (B) 0
(C) 2 (D) 3 (C) 3 (D) 2
4. A student finds, by trial, two negative and one posi-
2. If one of the roots of the equation x3 + 5x2 - 12x - 36 = tive root(s) of the equation x5 + 5x4 + 2802x + 3024
0 is thrice another root, then the third root is = 103x3 + 329x2. How many non-real roots does the
(A) -6 (B) 3 equation have?
89 (A) 0 (B) 1
−
-2 (D)
(C)
13 (C) 2 (D) 4
5. If the equation 3x4 - 13x3 + 7x2 + 17x + a - 10 = 0 has 15. The Newton’s iterative formula to find the value of 3
N
exactly three positive roots, then a can be is
(A) 11 (B) 4 N
(A) xi +1 = 2 xi − 2
(C) 13 (D) 12 xi
6. If two of the roots of the equation x3 + 3x2 - 10x - 24 = 1 N
0 are such that one is twice the other, then the third root xi +1 = xi − 2
(B)
3 xi
is
(A) -4 (B) -3 1 N
xi +1 = 2 xi + 2
(C)
(C) -2 (D) 3 3 xi
7. If 2.236146 is an approximation to 5, then the rela- 1 N
xi +1 = 2 xi − 2
(D)
tive error is 3 xi
(A) 3.4883 × 10-5 (B) 4.8383 × 10-5
16. Find the second approximation to the cube root of 24
(C) 8.3483 × 10 -4
(D) 5.8438 × 10-4
correct to three decimal places using Newton’s iterative
8. The least number of terms required to be considered formula.
1 (A) 2.695 (B) 2.885
in the Taylor’s series approximation of f ( x ) =
(2 + x) (C) 3.001 (D) None of these
about x = 0 such that the truncation error is at most 5 × 17. The Newton’s iterative formula to find the value of
10-4 for x ∈ [0, 1] is 1
(A) 3 (B) 5 is
N
(C) 6 (D) 8 (A) xi+1 = xi(2 + xiN)
9. Let f (x) = x3 - x - 5 = 0. By bisection method first two (B) xi+1 = xi(2 - xiN)
approximations x0 and x1 are 1.5 and 2.25 respectively, (C) xi+1 =xi2(2 + xiN)
then x2 is (D) None of these
(A) 1.625 (B) 1.875
18. Find the reciprocal of 22 using Newton–Raphson
(C) 1.999 (D) None of these
method.
10. Find the fourth approximation of the root of the equa- (A) 0.0454545 (B) 0.4504504
tion x3 + x - 11 = 0, between 2 and 3, using Bisection (C) 0.54054 (D) None of these
method. 19. If the first approximation of the root of x3 - 3x - 5 = 0
(A) 1.925 (B) 2.832 is (x0 =) 2, then find x1 by Newton–Raphson method.
(C) 2.5215 (D) 2.0625
(A) 2.2806 (B) 2.2790
11. The absolute error bisection method is (C) 2.3333 (D) None of these
1 20. Find the first approximation of the real root by Newton–
(A) 2n (B) |b - a|
2n Raphson method for x4 + x3 - 7x2 - x + 5 = 0 by taking
1 x0 = 2.
(C) (D) |b - a|2n
|b−a| (A) 2.066 (B) 2.981
(C) 2.819 (D) None of these
12. If the first two approximations x0 and x1 to a root of x3
- x - 4 = 0 are 1.666 and 1.780 respectively, then find 21. If y = 2.6 + 0.7x is a line that fits the data:
x2 by Regula–Falsi method. x -2 -1 0 1 2
(A) 1.974 (B) 1.794 y 1 2 3 K 4
(C) 1.896 (D) None of these
Then the value of K is
13. Find the second approximation to the root of the equa- (A) 3 (B) 5
tion 2x - 5 = 3sin x between (2, 3) using the method of (C) 6 (D) 7
false position.
(A) 2.2523 (B) 2.012 22. If a curve y = abx is fitted to the following data, then the
value of ‘b’ is
(C) 2.8804 (D) None of these
x -2 -1 0 1 2
14. For N = 28 and x0 = 5.5, the first approximation to N
by Newton’s iteration formula is y 11 13 20 25 34
23. For a set of 5 pairs of values (x0, y0), (x1, y1), (x2, y2) (x3, 3
1
y3) and (x4, y4) of (x, y), if D3y0 = 4 and D3y1 = 10, then 29. Find the value of ∫ 1 + x 2 dx taking four intervals by
2
the value of ∇4y4 is trapezoidal rule and also find the error when compared
(A) 4 (B) 6 to its exact value.
(C) 10 (D) 14 (A) 0.1759, 0.000004 (B) 0.1826, 0.04
24. The value of y(1.5) computed from the following data (C) 0.1953, 0.004 (D) 0.1423, -0.0004
using Newton’s forward interpolation formula is 1
x2
x 1 2 3 4 5
30. Find ∫ 1 + 8 x3 dx using Trapezoidal rule by taking 4
0
y 6 7 12 21 34 strips.
(A) 6 (B) 6.5 (A) 0.0911 (B) 0.9011
(C) 0.1901 (D) None of these
(C) 7 (D) 7.5 1.5
dx
25. The Lagrange’s interpolation polynomial correspond- 31. The estimate of ⋅
∫ x
obtained using Simpson’s rule
ing to the pairs of values of x and y given in the follow- 0.5
ing table is with three point function evaluation exceeds the exact
value by
x 1 3 4 7 (A) 0.235 (B) 0.068
y 36 16 9 72 (C) 0.024 (D) 0.012
6
(A) x3 - 6x2 + 9x + 36
(B) x3 - 6x2 + 18x - 45
32. Find the value of ⋅
∫ x log xdx taking 4 strips by
2
(C) 3x3 + 4x2 - 5x + 27 1
Simpson’s rd rule upto four decimals.
(D) x3 - 7x2 + 5x + 37 3
(A) 21.8901 (B) 22.8661
dy d2 y (C) 23.6581 (D) None of these
26. The values of and at x = 5 from the following
dx dx 2 π /2
table respectively are 1
33. Evaluate ∫ sin x dx by Simpson’s
3
rule using six
0
x 0 1 2 3 4 5 intervals.
y = f(x) 1 4 9 16 21 28 (A) 0.97768 (B) 0.98869
(C) 0.99968 (D) None of these
(A) 15.00 and 26.45
(B) 13.73 and 23.33 34. Find the maximum error in evaluating the above when
compared to its exact value.
(C) 17.13 and 31.42
(A) 0.000032 (B) 0.00032
(D) 21.64 and 43.00 (C) 0.00000032 (D) 0.0032
dy 3
1 3
27. The value of at x = 12 from the table given below is 35. Evaluate ∫ dx by using Simpson’s rule by
dx 2+ x 2 8
0
x 10 15 20 25 taking 3 strips.
y = f(x) 354 332 291 260 (A) 0.507 (B) 0.5007
(C) 0.7839 (D) None of these
-3.9427
(A) (B) 4.6125 dy
-0.4652
(C) (D) 1.3549 36. If = 1 - 3xy2, y (0) = 0, then by Taylor’s method y
dx
dy (0.1) =
28. The magnitude of error when at x = 2 is found by
dx (A) 0.02 (B) 0.001
1 (C) 0.05 (D) 0.1
Newton’s forward interpolation formula for y = y =
x dy
using the following data is 37. If = 2x + y, y (0) = 1, the Picard’s approximate of y
dx
upto second degree terms is
x 2 4 6 8
x2
y = 1/x 0.5000 0.2500 0.1667 0.1250 (A) 1 + x + x2 (B) 1 + x +
2
(A) 0.0005 (B) 0.0025 x2
(C) 1 - x - (D) None of these
(C) 0.0125 (D) 0.0625 2
38. If y0 = 1, f (x0, y0) = 1.2, f (x1, y1(0)) = 1.9312, h = 0.3, by 42. For the above problem find y(0.8) using Milne’s correc-
modified Euler’s formula y1 = (1) tor formula.
(A) 1.4696 (B) 1.2015 (A) 0.3046 (B) 0.4036
(C) 1.325 (D) 1.525 (C) 0.436 (D) None of these
39. Using Euler’s modified method, find a solution of the
dy Direction for questions 43 and 44:
equation = x + y with y(0) = 1 at y(0.2).
dx 43. Find using the Adams–Bashforth corrector formula y
(A) 1.3902 (B) 1.2309 dy 1
(0.4), for the differential equation = xy , given
(C) 1.3092 (D) None of these dx 2
dy y(0.1) = 1.01, y(0.2) = 1.022, y (0.3) = 1.023, y01 = 0, y11
40. Find k1, by Runge–Kutta method of fourth order if
dx = 0.0505, y21 = 0.1022, y31 = 0.1535.
= 2x + 3y2 and y(0.1) = 1.1165, h = 0.1.
(A) 0.3993 (B) 0.9393 (A) 1.5418 (B) 1.0410
(C) 0.3939 (D) None of these (C) 1.4100 (D) None of these
44. For the above differential equation find y(0.5) using
Direction for questions 41 and 42: Adams–Bashforth predictor formula.
dy (A) 1.00463 (B) 1.06463
41. Find y(0.8) by Milne’s predictor formula, given =
dx (C) 1.00599 (D) None of these
x - y2, y2 = 0.0795, y (0.6) = 0.1762, y0 = 0.0000, y11 =
0.1996, y21 = 0.3937, y31 = 0.5689, h = 0.2 45. The Runge–Kutta methods has the error of order _____.
(A) 0.9304 (B) 0.4930 (A) 1 (B) 3
(C) 0.3049 (D) None of these (C) 5 (D) 2
7. In Newton–Raphson iterative method, the initial 9. For step-size, Dx = 0.4, the value of following integral
guess value (xinj) is considered as zero while finding 1
using Simpson’s rule is ________.
the roots of the equation: f (x) = -2 + 6x - 4x2 + 0.5x3. 3
The correction, Dx, to be added to xinj in the first itera- 0.8
tion is ______. [GATE, 2015] ∫ (0.2 + 25 x − 200 x 2 + 675 x3 − 900 x 4 + 400 x5 )dx
8. The quadratic equation x - 4x + 4 = 0 is to be solved
2 0
Answer Keys
Exercises
1. D 2. B 3. D 4. A 5. B 6. D 7. A 8. C 9. B 10. D
11. B 12. B 13. C 14. A 15. C 16. B 17. B 18. A 19. C 20. A
21. A 22. C 23. B 24. A 25. D 26. B 27. A 28. D 29. D 30. A
31. D 32. B 33. C 34. B 35. C 36. D 37. A 38. A 39. B 40. C
41. C 42. A 43. B 44. B 45. C