Numerical Methods
Numerical Methods
Numerical Methods
NUMERICAL METHODS
(For B.A. and B.Sc. VIth Semester students of All Colleges affiliated to
Universities in Kumaun)
As per Kumaun University Syllabus
By
Dedicated
to
Lord
Krishna
Authors & Publishers
B rief C ontents
Dedication...............................................................(iII)
Brief Contents.........................................................(IV)
(iv)
N-3
Chapter-1
The Calculus of Finite Differences
Comprehensive Problems 1
Problem 1: Evaluate the following :
2x
(i) ∆ (ii) ∆ cot 2 x
( x + 1) !
(iii) ∆ sinh ( a + bx ) (iv) ∆ cosh ( a + bx )
∆2 x3
(v) ∆ tan ax (vi)
E x3
(vii) ∆4 ae x (viii) ∆ tan−1 ax.
2x 2 x +1 2x
Solution: (i) ∆ = − , taking the interval of differencing as unity
( x + 1) ! ( x + 2) ! ( x + 1) !
2x 2 2x 2 − x − 2 − x . 2 x
= − 1 = = ⋅
( x + 1) ! x + 2 ( x + 1) ! x + 2 ( x + 2) !
= − cosec 2 x +1
sin [ ax + a − ax ] sin a
= = ⋅
cos ax cos a ( x + 1) cos ax cos a ( x + 1)
∆2 x3 ( E − I )2 x3 ( E2 + I − 2 E ) x3 E2 x3 + Ix3 − 2 E x3
(vi) = = =
3 3 3
E x E x E x E x3
( x + 2)3 + x3 − 2 ( x + 1)3
=
( x + 1)3
x3 + 8 + 6 x2 + 12 x + x3 − 2 ( x + 1)3
=
( x + 1)3
6 ( x + 1) 6
= = ⋅
3
( x + 1) ( x + 1)2
(vii) ∆ ae x = a ∆ e x = a [ e x + h − e x ]
= ae x [ e h − 1] = ae x ( e − 1), for h = 1.
Now ∆2 ae x = ∆ [ ∆ ae x ]
= ∆ [ ae x ( e − 1)] = a ( e − 1) [ e x +1 − e x ]
= a ( e − 1) e x ( e − 1) = a ( e − 1)2 e x .
= 2 cos a + bx + sin
b b
2 2
π
sin a + bx + +
b b
= 2 sin
2 2 2
b + π
sin a + bx +
b
= 2 sin ⋅
2 2
N-5
b + π
∆2 sin ( a + bx ) = ∆ 2 sin sin a + bx +
b
And
2 2
= 2 sin
b sin a + b ( x + 1) + b + π − sin a + bx + b + π
2 2 2
2b + π
⋅ 2 cos a + bx +
b b
= 2 sin sin
2 2 2
2
2b + π
= 2 sin cos a + bx +
b
2 2
2
b + π
= 2 sin sin a + bx + 2
b
⋅
2 2
= − 2 sin a + bx + sin
b b
2 2
b + π
cos a + bx +
b
= 2 sin ⋅
2 2
b + π
∆2 cos ( a + bx ) = ∆ 2 sin cos a + bx +
b
And
2 2
= 2 sin
b cos a + b ( x + 1) + b + π − cos a + bx + b + π
2 2 2
2b + π
⋅ 2 sin a + bx +
b b
= 2 sin − sin
2 2 2
2
b + π
= 2 sin cos a + bx + 2
b
⋅
2 2
(iii) We have
cos ( a + bx + b ) cos ( a + bx )
= −
sin ( a + bx + b ) sin ( a + bx )
(i) ∆ [ f ( x ) + g( x )] = ∆ f ( x ) ± ∆ g( x ),
(ii) ∆ [ a f ( x )] = a ∆ f ( x ), a being a constant,
(iii) ∆ [ f ( x ) g( x )] = f ( x + 1)∆ g( x ) + g( x ) ∆ f ( x )
f ( x ) g ( x ) ∆ f ( x ) − f ( x ) ∆ g( x )
(iv) ∆ = ⋅
g ( x ) g( x ) g ( x + 1) (Kanpur 2012)
∴ ∆ ( x2 ) = ( x + h)2 − x2 = 2hx + h2 ;
∆ ( e x ) = e x + h − e x = e x ( e h − 1) ;
= a ∆3 x3 = a . 3 ! h3 = 6 ah3 .
Problem 5(ii): If f ( x ) = e ax , show that f(0 )and its leading differences form a geometrical progression .
Solution: We have f ( x ) = e ax .
∴ f (0 ) = e0 = 1
∆f (0 ) = f (0 + h) − f (0 ) = f ( h) − f (0 ) = e ah − 1
∆2 f (0 ) = f (0 + 2h) − 2 f (0 + h) + f (0 ) = f (2h) − 2 f ( h) + f (0 )
= e2 ah − 2e ah + 1 = ( e ah − 1)2
∆3 f (0 ) = f (0 + 3h) − 3 f (0 + 2h) + 3 f (0 + h) − f (0 )
= f (3h) − 3 f (2h) + 3 f ( h) − f (0 )
N-7
= e3 ah − 3e2 ah + 3e ah − 1
= ( e ah − 1)3
… … … … … … …
… … … … … … …
Thus f (0 ) and its leading differences form a G.P. whose first tem is 1 and common ratio is
( e ah − 1).
cos a x + + b ⋅
ah h
(i) ∆ sin ( ax + b ) = 2 sin
2 2
Solution: We have
∇ f ( x ) = f ( x ) − f ( x − h) ;
∆ E − 1 f ( x ) = ∆ f ( x − h) = f ( x ) − f ( x − h) ;
E −1 ∆ f ( x ) = E −1 [ f ( x + h) − f ( x )] = E −1 f ( x + h) − E −1 f ( x )
= f ( x ) − f ( x − h) ;
( I − E −1 ) f ( x ) = f ( x ) − E −1 f ( x ) = f ( x ) − f ( x − h).
= ( E3 − 3E2 + 3E − I ) f ( x ).
N-8
Problem 9: Find the successive differences of f ( x ) = ab cx and sum the first n differences.
Solution: Here f ( x ) = ab cx .
∴ ∆ f ( x ) = ab c ( x + h) − ab cx = ab cx . ( b ch − 1).
Again ∆2 f ( x ) = ∆ ∆ f ( x ) = ∆ [ ab c x . ( b c h − 1)]
= a ( b ch − 1) ∆ b cx = a ( b ch − 1) [ b c ( x + h) − b cx ]
= a ( b ch − 1) . b cx ( b ch − 1) = ab cx ( b ch − 1)2
and so on.
In general, ∆n f ( x ) = ab cx ( b ch − 1)n.
= ab cx ( b ch − 1) + ab cx ( b ch − 1)2 + … + ab cx ( b ch − 1)n
= ab cx ( b ch − 1) [1 + ( b ch − 1) + ( b ch − 1)2 + … + ( b ch − 1)n−1 ]
1 − ( b ch − 1)n ( b ch − 1)n − 1
= ab cx ( b ch − 1) cx ch
= ab ( b − 1) ⋅ ⋅
ch
1 − ( b − 1) b ch − 2
Solution: Let the required function be f ( x ) and the interval of differencing h = 1. Then
(i) As ∆ f ( x ) = a (i. e., a constant), therefore f ( x ) will be a polynomial of first degree and hence
f ( x ) = ax + b, where b is a constant.
(ii) As ∆ f ( x ) = ax + b …(1)
i. e., a polynomial of first degree, so f ( x ) will be a polynomial of second degree.
Let f ( x ) = Ax2 + Bx + C.
or ∆ f ( x ) = 2 Ax + ( A + B ) …(2)
Comparing (1) and (2), we get
2 Ax + ( A + B ) ≡ ax + b
1 1
giving A= a and B=b− a.
2 2
N-9
π . cos x + π
1 1
= k . 2 sin
2 2
= 2k ( − sin x ). …(2)
Comparing (1) and (2), we get
1
− 2k = 1 or k=− ⋅
2
1
Thus f ( x) = − sin x.
2
(iv) Here ∆ f ( x ) = e a + bx . …(1)
so that ∆ f ( x ) = k e a + b ( x +1) − k . e a + bx
= k . e a + bx ( e b − 1). …(2)
Comparing (1) and (2), we get
1
k ( e b − 1) = 1 or k= ⋅
b
e −1
1
∴ f ( x) = ⋅ e a + bx
eb − 1
e a + bx
or f ( x) = ⋅
eb − 1
Let f ( x) = k . 5 x and h =1
x +1
so that ∆ f ( x) = k 5 − k 5 = k 5 x [5 − 1] = 4k 5 x .
x
…(2)
Comparing (1) and (2), we get
1
4k 5 x = 5 x giving k = ⋅
4
1 x
Hence f ( x) = 5 .
4
∆2
(iii) [3 x(4) − 5 x(−2) ]
∆ x2
(iv) Find the relation between α, β and γ in order that α + βx + γ x2 may be expressible in one term in the
factorial notation. (Kanpur 2009)
2 2
(v) Find the value of E x when the values of x vary by a constant increment of 2.
(Gorakhpur 2013)
Solution: We have
∆
(i) [3 x(4) − 2 x(3) − 3 x(2) + x − 5] = 12 x(3) − 6 x(2) − 6 x (1) + 1
∆ x
∆2 ∆ ∆ ∆
(iii) [3 x (4) − 5 x(−2) ] = [3 x(4) − 5 x(−2) ] = [12 x (3) + 10 x(−3) ]
∆ x2 ∆ x ∆ x ∆ x
= 36 x (2) − 30 x(−4).
= ( a + bx )2 − ab − b2 x
= ( a2 − ab ) + (2ab − b2 ) x + b2 x2 .
∴ α + βx + γ x2 = ( a2 − ab ) + (2ab − b2 ) x + b2 x2 .
= ( x + 2 + 2)2 = ( x + 4)2
= x2 + 8 x + 16.
Problem 12: A second degree polynomial passes through (0, 1), (1, 3), (2, 7) and (3, 13). Find the
polynomial.
Solution: Let us consider the 2nd degree polynomial in factorial notation i. e., let
f ( x ) = Ax(2) + Bx(1) + C.
(1, 3) ⇒ x = 1, f ( x ) = 3;
(2, 7) ⇒ x = 2, f ( x ) = 7
and (3, 13) ⇒ x = 3, f ( x ) = 13.
Thus the difference table for the given values of x and f ( x ) is as follows :
x f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x )
0 1
2
1 3 2
4 0
2 7 2
6
3 13
∴ ∆ f ( x ) = 2 Ax(1) + B = 2 Ax + B.
Hence ∆ f (0 ) = B.
But ∆ f (0 ) = 2 (from table)
so that B =2.
Also ∆2 f ( x ) = 2 A.
∴ ∆2 f (0 ) = 2 A.
∴ 2 A = 2 giving A = 1.
Again for x = 0, f ( 0 ) = C giving C = 1, (from table).
Thus f ( x ) = 1 . x(2) + 2 . x (1) + 1 = x ( x − 1) + 2 x + 1
= x2 + x + 1.
= (1 + ∆ ) x f (0 )
x ( x − 1) 2
= f (0 ) + x ∆ f (0 ) + ∆ f (0 )
1. 2
x ( x − 1)
=1+ x .2 + ⋅2
2
= 1 + 2 x + x ( x − 1)
= x2 + x + 1.
N-12
Problem 13: A third degree polynomial passes through the points (0, − 1), (1, 1), (2, 1) and (3, 2). Find the
polynomial.
x f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x )
0 −1
2
1 1 −2
0 −1
2 1 −3
−3
3 −2
= 3 Ax ( x − 1) + 2Bx + C.
∴ ∆ f (0 ) = C and from table ∆ f (0 ) = 2.
Comparing the two values, we get C = 2.
Now ∆2 f ( x ) = 6 Ax(1) + 2B = 6 Ax + 2B.
∴ 2B = − 2 giving B = − 1.
Again ∆3 f ( x ) = 6 A.
∴ ∆3 f (0 ) = 6 A
= (1 + ∆ ) x f (0 )
x ( x − 1) 2 x ( x − 1) ( x − 2) 3
= f (0 ) + x ∆ f (0 ) + ∆ f (0 ) + ∆ f (0 )
1. 2 1. 2 . 3
x ( x − 1) x ( x − 1) ( x − 2)
= −1+ x .2 + ⋅ ( − 2) + ⋅ ( − 1)
2 6
1 3
=− ( x + 3 x2 − 16 x + 6).
6
2
∆ Ux ∆2U x
Problem 14: What is the difference between and ⋅ If U x = x3 and interval of
E Ux 2
E Ux
differencing is one, find expressions for both.
2
∆ Ux
Solution: In the first forward difference operator ∆ and the shifting operator E are to
E Ux
be operated on U x and then we are to find the square of the quotient of the two expressions thus
∆2 U x
found while in second forward difference operator ∆2 and shifting operator E2 are to be
2
E Ux
operated on U x and then we are to find the quotient of the two expressions thus obtained. For
U x = x3 ,
2 2 2
∆ Ux ∆ x3 ( x + 1)3 − x3
=
=
E Ux E x3 ( x + 1)
3
2
x3 x6 2 x3
= 1 − =1+ − ⋅
( x + 1)3 ( x + 1)6 ( x + 1)3
∆2 U x ∆2 x3 ∆ [( x + 1)3 − x3 ]
While, = =
E2 U x E2 x3 [ x + 2]3
( x + 2)3 − 2 . ( x + 1)3 + x3
=
( x + 2)3
x3 + 8 + 3 . 2 x2 + 3 . 4 x − 2 ( x3 + 1 + 3 x2 + 3 x ) + x3
=
x3 + 8 + 6 x2 + 12 x
6x + 6
=
x3 + 6 x2 + 12 x + 8
6 ( x + 1)
= ⋅
x3 + 6 x2 + 12 x + 8
N-14
x : 1 2 3 4 5 6 7
y : 2 4 8 ? 32 64 128
Explain why the result differs from 16. (Kanpur 2008; Kumaun 09)
(ii) Estimate the production of cotton in the year 1935 from the data given below :
Production f ( x ) : 17 ⋅1 13 ⋅ 0 14 ⋅ 0 9⋅6 — 12 ⋅ 4 18 ⋅ 2
(In millions)
Solution: (i) Here six values of x and y are given, therefore we can assume the given function
expressible as a polynomial of degree five.
Hence ∆6 y = 0 i. e., ∆6 f ( x ) = 0
or ( E − I )6 f ( x ) = 0
or [E6 − 6 C1 E5 + 6 C2 E4 − 6 C3 E 3 + 6 C4 E2 − 6 C5 E + 6 C6 ] f ( x ) = 0
or E6 f ( x ) − 6 E5 f ( x ) + 15 E4 f ( x ) − 20 E3 f ( x ) + 15 E2 f ( x )
− 6 E f ( x) + f ( x) = 0
or f ( x + 6) − 6 f ( x + 5) + 15 f ( x + 4) − 20 f ( x + 3) + 15 f ( x + 2)
− 6 f ( x + 1) + f ( x ) = 0.
For x = 1, we have
f (7) − 6 f (6) + 15 f (5) − 20 f (4) + 15 f (3) − 6 f (2) + f (1) = 0
or 128 − 6 × 64 + 15 × 32 − 20 f (4) + 15 × 8 − 6 × 4 + 2 = 0
or 128 − 384 + 480 − 20 f (4) + 120 − 24 + 2 = 0
or 20 f (4) = 730 − 408 = 322.
∴ f (4) = 16 ⋅1.
Explanation: From the given values we observe that the function is of the form y = 2 x . From
this we get 24 = 16 i. e., y = 16 for x = 4. This value differs from the calculated value as we have
assumed the function as a polynomial while actually the behaviour of the given values shows that
the function is y = 2 x and hence the difference.
(ii) As six observations are given, therefore we can fit a polynomial of fifth degree and hence the
sixth difference will be equal to zero.
∴ ∆6 f ( x ) = 0
or ( E − I )6 f ( x ) = 0
or [E6 − 6 C1 E5 + 6 C2 E4 − 6 C3 E 3 + 6 C4 E2 − 6 C5 E + 6 C6 ] f ( x ) = 0
or E6 f ( x ) − 6 E5 f ( x ) + 15 E4 f ( x ) − 20 E3 f ( x ) + 15 E2
f ( x) − 6 E f ( x) + f ( x) = 0
N-15
or f ( x + 6) − 6 f ( x + 5) + 15 f ( x + 4) − 20 f ( x + 3) + 15 f ( x + 2)
− 6 f ( x + 1) + f ( x ) = 0.
Taking x = 1931, we get
f (1937) − 6 f (1936) + 15 f (1935) − 20 f (1934) + 15 f (1933)
− 6 f (1932) + f (1931) = 0
or 18 ⋅ 2 − 6 × 12 ⋅ 4 + 15 × f (1935) − 20 × 9 ⋅ 6 + 15 × 14 − 6 × 13 + 17 ⋅1 = 0
or 15 f (1935) = 99 ⋅1 giving f (1935) = 6 ⋅ 60.
Solution: We have ∆u x = 2u x
∴ u x = 3 x for h = 1
∆u x = ∆3 x = 3 x + 1 − 3 x = 3 x (3 − 1) = 3 x ⋅ 2 = 2 ⋅ u x
Solution: We have u x = x ( x − 1) ( x − 2)
∆u x = ( x + 1) x ( x − 1) − x ( x − 1) ( x − 2)
= x ( x − 1) ( x + 1 − x + 2) = 3 x( x − 1).
x(3)
∴ ux = + C, where C is an arbitrary constant
3
1
= x ( x − 1) ( x − 2) + C.
3
Solution: (i) Converting the given function into factorial notation, we have
x3 + 4 x2 + 9 x + 12 = Ax (3) + Bx(2) + Cx(1) + D.
0 1 5
2 1 5 14 = C
0 2
3 1
7=B
0
1= A
N-16
Solution: Here five values of [ x, f ( x )] are known therefore we can assume that f ( x ) can be
represented by a polynomial of degree four.
Hence ∆5 f ( x ) = 0 or ( E − I )5 f ( x ) = 0
or [ E5 − 5 C1 E4 + 5 C2 E3 − 5 C3 E2 + 5 C4 E − 5 C5 ] f ( x ) = 0
or f ( x + ⋅ 5) − 5 f ( x + ⋅ 4) + 10 f ( x + ⋅ 3) − 10 f ( x + ⋅ 2)
+ 5 f ( x + ⋅1) − f ( x ) = 0,
because here the interval of differencing h = ⋅1.
At x = 0 (first tabulated value of x), we have
f (⋅5) − 5 f (⋅4) + 10 f (⋅3) − 10 f (⋅2) + 5 f (⋅1) − f (0 ) = 0 …(1)
Substituting the known values from the given data, (1) gives
⋅082 − 5 f (⋅4) + 10 × ⋅100 − 10 × ⋅111 + 5 f (⋅1) − ⋅135 = 0
or 5 f (⋅4) − 5 f (⋅1) = − ⋅163 …(3)
and (2) gives
⋅074 − 5 × ⋅ 082 + 10 f (⋅4) − 10 × ⋅100 + 5 × ⋅111 − f (⋅1) = 0
or 10 f (⋅4) − f (⋅1) = ⋅ 881. …(4)
Solving (3) and (4), we get
f (⋅1) = ⋅123 and f (⋅4) = ⋅ 090,
which are the required missing values.
N-17
or [E8 − 8 C1 E7 + 8 C2 E6 − 8 C3 E5 + 8 C4 E4 − 8 C5 E3
+ 8 C6 E2 − 8 C7 E + 8 C8 ] U0 = 0
or U8 − 8 C1 U7 + 8 C2 U6 − 8 C3 U5 + 8 C4 U4 − 8 C5 U3
+ 8 C6 U2 − 8 C7 U1 + U0 = 0
+ 8 C4 U4 = 0.
= 2 ∆2 ( x2 + 2 x + 1) + ∆ ( x2 + 2 x + 1) − 1 ( x2 + 2 x + 1)
= 2 ∆2 x2 + ∆ ( x2 + 2 x ) − x2 − 2 x − 1 [∵ ∆2 (2 x + 1) = 0, ∆ 1 = 0 ]
= 2∆ [( x + h)2 − x2 ] + [( x + h)2 + 2 ( x + h) − x2 − 2 x ] − x2 − 2 x − 1
= 2 ∆ (2hx + h2 ) + 2hx + h2 + 2h − x2 − 2 x − 1
= 4h ∆x + 2hx + h2 + 2h − x2 − 2 x − 1 [∵ ∆ h2 = 0 ]
= 4h ( x + h − x ) + 2hx + h2 + 2h − x2 − 2 x − 1
= 5h2 + 2hx + 2h − x2 − 2 x − 1 .
= ( ∆ + 1) [2∆ ( x2 + 2 x + 1) − 1 ( x2 + 2 x + 1)]
= ( ∆ + 1) [2 {( x + h)2 + 2( x + h) + 1 − x2 − 2 x − 1} − x2 − 2 x − 1]
= ( ∆ + 1) [2(2hx + h2 + 2h) − x2 − 2 x − 1]
N-18
= (∆ + 1)(4hx + 2h2 + 4h − x2 − 2 x − 1)
= 4h ∆ x − ∆ x2 − 2∆ x + 4hx + 2h2 + 4h − x2 − 2 x − 1
= 5h2 + 2hx + 2h − x2 − 2 x − 1.
(iii) We have ( E − 2) ( E − 1) (2 x / h + x )
= ( E − 2)[ E (2 x / h + x ) − (2 x / h + x )]
= ( E − 2)[2( x + h) / h + x + h − 2 x / h − x ]
= E (2 x / h + h) − 2(2 x / h + h)
= 2( x + h) / h + h − 2( x + h) / h − 2h = − h
(iv) We have ( E2 − 3E + 2) (2 x / h + x )
= E2 (2 x / h + x ) − 3E (2 x / h + x ) + 2(2 x / h + x )
Problem 23: Use the method of separation of symbols to prove the following identities :
1 2 1.3 4 1.3.5 6
(ii) Ux − ∆ U x −1 + ∆ Ux −2 − ∆ Ux −3 + …
8 8.16 8.16.24
1 1
= U x + 1 /2 − ∆U x + 1 /2 + ∆2U x + 1 /2 − …
2 4
(iii) U x = U x − 1 + ∆U x − 2 + ∆2U x − 3 + … + ∆n − 1 U x − n + ∆nU x − n.
= R.H.S.
1 1. 3 4 1. 3 . 5 6
(ii) We have U x − ∆2 U x −1 + ∆ U x −2 − ∆ U x −3 + …
8 8 . 16 8 . 16 . 24
1 1 . 3 4 −2 1 . 3 . 5 6 −3
= 1 − ∆2 E −1 + ∆ E − ∆ E + … U x
8 8 . 16 8 . 16 . 24
N-19
− 1 − 3
2 −1 2 −1
2
= 1+ − 1 ∆ E + 2 2 ∆ E
2 4 2! 4
− 1 − 3 − 5
3
2 2 2 ∆2 E −1
+
3! 4 + … U x
−1 / 2 −1 / 2
∆2 E −1 ∆2
= 1 + U x = E1 /2 E + Ux
4 4
−1 / 2
∆2
= E1 /2 1 + ∆ + Ux [∵ E ≡ 1 + ∆ ]
4
−1
∆ ∆ ∆2
= E1 /2 1 + U x = E1 /2 1 − + − … U x
2 2 22
1 1
= U x +1 /2 − ∆ U x +1 /2 + ∆2 U x +1 /2 − …
2 4
(iii) We have U x − ∆nU x − n = U x − ∆nE − nU x = (1 − ∆nE − n )U x
∆n E n − ∆n 1 E n − ∆n
= 1 − U = Ux = Ux
n x n
E E E n I
1 E n − ∆n
= Ux [∵ I = E − ∆ ]
E n E − ∆
= E − n[ E n − 1 + ∆ E n − 2 + ∆2 E n − 3 + … + ∆n − 1 ] U x , [Dividing Nr. by Dr.]
−1 −2 2 −3 n −1 − n
= [E + ∆E +∆ E +… + ∆ E ]Ux
2 n −1
= U x − 1 + ∆U x − 2 + ∆ U x − 3 + … + ∆ U x − n.
U0 + U1 + U2 + … + U n = U0 + EU0 + E2 U0 + … + E n U0
E n +1 − 1
= (1 + E + E2 + … + E n ) U0 = U0 , summing the G.P.
E −1
(1 + ∆ )n +1 − 1
= U0
∆
1
= [(1 + n +1C1 ∆ + n +1C2 ∆2 + n +1C3 ∆3 + … + ∆n +1 ) − 1] U0
∆
1 n +1
= [ C1 ∆ U0 + n +1C2 ∆2 U0 + n +1C3 ∆3 U0 + … + ∆n +1 U0 ]
∆
= n+1C1 U0 + n+1C2 ∆ U0 + n+1C3 ∆2 U0 + … + ∆n U0 .
N-20
= U n + x C1 ∆ E −1 U n + x +1C2 ∆2 E −2 U n + x +2 C3 ∆3 E −3 U n + …
= [1 + x C1 ∆ E −1 + x +1C2 ∆2 E −2 + x +2 C3 ∆3 E −3 + … ] U n
= [1 − ∆ E −1 ]− x U n
−x −x
∆ E − ∆
= 1 − Un = Un
E E
−x
=
1
U n = E x U n = U n+ x = U x + n .
E
(i) U0 + x
C1∆U1 + x C2 + … = U x + x C1 ∆2 U x − 1 + x C2 ∆4 U x − 2 + … (Kumaun 2013)
(ii) U x − U x +1 + U x +2 − U x +3 + …
1
2
=
1 2 1. 3 1 ∆4 U x −5 /2
U x −1 /2 − ∆ U x −3 /2 +
2 8 2! 8
1. 3 . 5 1
3
− ∆6 U x −7 /2 + … ⋅
3 ! 8
= [1 + x C1 ∆2 E −1 + x C2 ∆4 E −2 + … ] U x
x x
2 −1 x E + ∆2 E2 − E + 1
= (1 + ∆ E ) U x = Ux =
Ux
E E
1
= [1 + E ( E − 1)] x U x = E − x (1 + ∆E ) x U x
Ex
= (1 + ∆E ) x U0 = (1 + x C1 ∆E + x C2 ∆2 E2 + … ) U0
= U0 + x C1 ∆ U1 + x C2 ∆2U2 + …
N-21
(ii) We have
1
2
1 2 1. 3 1
U x −1 /2 − ∆ U x −3 /2 + ∆4 U x −5 /2
2 8 2 ! 8
1. 3 . 5 1
3
− ∆6 U x −7 /2 + …
3! 8
1 3
2
1 −1 / 2 1 1 2 −3 / 2 1
= E Ux − ∆ E U x + 2 2 ∆4 E −5 /2 U x
2 2 4 2! 4
1 3 5
3
1
− 2 2 2 ∆6 E −7 /2 U x + …
3! 4
− 1 − 3
2
2 2 1 2 −1
= E −1 /2 1 + − ∆2 E −1 +
1 1 1
∆ E
2 2 4 2! 4
− 1 − 3 − 5
3
2 2 2 1 2 −1
+ ∆ E + … U x
3! 4
−1 / 2
=
1 −1 / 2
E 1 + 1 ∆2 E −1 Ux
2 4
−1 / 2 −1 /2
1 −1 /2 4E + ∆2 1 −1 /2 4 + 4 ∆ + ∆2
= E Ux = E Ux
2 4E 2 4E
1 −1 / 2
= E . 2E1 /2 {(2 + ∆ )2 }−1 /2 U x = (2 + ∆ )−1 U x
2
= (1 + E )−1 U x = [1 − E + E2 − E3 + E4 − E5 + … ] U x
= U x − U x +1 + U x + 2 − U x + 3 + U x + 4 − U x + 5 + …
= {(1 + x ) + x ∆ }n U0 = {1 + x (1 + ∆ )}n U0
Solution: We have ∆ g ( x − 1) = g ( x ) − g ( x − 1)
∴ f ( x − 1) ∆ g ( x − 1) = f ( x − 1) g ( x ) − f ( x − 1) g ( x − 1)
⇒ ∆ [ f ( x − 1) ∆ g ( x − 1)] = ∆ [ f ( x − 1) g ( x )] − ∆ [ f ( x − 1) g ( x − 1)] …(1)
= f ( x ) g ( x + 1) − f ( x − 1) g ( x ) − f ( x ) g ( x ) + f )( x − 1) g ( x − 1)
N-22
= [ f ( x ) g ( x + 1) − f ( x − 1) g ( x )] − [ f ( x ) g ( x ) − f ( x − 1) g ( x − 1)]
= ∇ [ f ( x ) g ( x + 1)] − ∇ [ f ( x ) g ( x )]
= ∇ [ f ( x ) g ( x + 1) − f ( x ) g ( x )]
= ∇ [ f ( x ) { g ( x + 1) − g ( x )}]
= ∇ [ f ( x ) ∆ g ( x )].
Now from (1), we have
∆[ f ( x − 1) ∆ g ( x − 1)] = ∆ [ f ( x − 1) g ( x ) − f ( x − 1) g ( x − 1)]
= ∆ [ f ( x − 1){ g ( x ) − g ( x − 1)}] = ∆ [ f ( x − 1) ∇ g ( x )]
1 ∆2 y ∆3 y ∆4 y
Problem 26(ii): Show that Dy = [ ∆y − + − + …].
h 2 3 4
(Garhwal 2010)
Problem 27: Express the function 2 x3 − 3 x2 + 3 x − 10 in factorial notation, the interval of differencing
being unity. Also find its successive differences. (Kanpur 2010)
3 2 (3) (2) (1)
Solution: Let 2 x − 3 x + 3 x − 10 = Ax + Bx + Cx + D.
The coefficients A, B, C, D are easily obtained by the method of synthetic division. The procedure
is given below :
1 2 −3 3 − 10 = D
0 2 −1
2 2 −1 2=C
0 4
3 2 3=B
2= A
= 6 x ( x − 1) + 6 x + 2 = 6 x2 + 2,
∆2 f ( x ) = 6. 2 x(1) + 6 = 12 x + 6,
and ∆3 f ( x ) = 12.
N-23
Problem 28: If
x : 1 2 3 4 5
y : 2 5 10 20 30
x y = f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x )
1 2
2 5 2
5 3
3 10 5 −8
10 −5
4 20 0
10
5 30
Problem 29: If f (0 ) = − 3, f (1) = 6, f (2) = 8 and f (3) = 12, prepare forward difference table to find
∆3 f (0 ).
∆n 0 n +1 = n [ ∆n −1 0 n + ∆n 0 n ]
∆n −1 0 n = ( n − 1) [ ∆n − 2 0 n −1 + ∆n −1 0 n −1 ]
∆n − 2 0 n −1 = ( n − 2) [ ∆n − 3 0 n − 2 + ∆n − 2 0 n − 2 ]
… … … … …
… … … … …
∆2 03 = 2 [ ∆ 02 + ∆2 02 ]
∆ 02 = 1 [ ∆0 01 + ∆1 01 ].
N-24
n ( n − 1) 2
= ( −1)n 1 − n∆ + ∆ un [neglecting higher order differences
1. 2
as un is a polynomial of second degree]
n ( n − 1) 2
= ( −1)n un − n ∆ un + ∆ un
2
n2 − n 2
= ( −1)n ( an2 + bn + c ) − n∆ ( an2 + bn + c ) + ∆ ( an2 + bn + c )
2
n2 − n
= ( −1)n ( an2 + bn + c ) − n { a ∆ n2 + b ∆ n} + ⋅ a ∆2 n2
2
= ( −1)n ( an2 + bn + c ) − n { a ( n + 1)2 − an2 }
n2 − n
− bn ( n + 1 − n) + ⋅ a ∆ {( n + 1)2 − n2 }
2
n2 − n
= ( −1)n ( an2 + bn + c ) − n (2an + a + b ) + ⋅ a ∆ (2n + 1)
2
n2 − n
= ( −1)n ( an2 + bn + c ) − n (2an + a + b ) + ⋅ a {2 ( n + 1) − 2n}
2
= ( −1)n [( an2 + bn + c ) − n (2an + a + b ) + a ( n2 − n)]
= ( −1)n [ c − 2an]
= R. H. S.
N-25
Problem 32: Taking fifth order differences of U x to be constant and given U0 , U1 , … , U6 , prove that
1 25 ( c − b ) + 3 ( a − c )
U = c+
21 2 256
2
Solution: We have U = U5 /2 = E 5 /2 U0
21
2
5 5 5 5 5
− 1 − 1 − 2
5 2 2 2 2 2 2 3
= U0 + ∆ U0 + ∆ U0 + ∆ U0
2 2! 3!
5 5 5 5
− 1 − 2 − 3
2 2 2 2 4
+ ∆ U0
4!
5 5 − 1 5 − 2 5 − 3 5 − 4
2 2 2 2 5
+ 2 ∆ U0
5!
5 15 2 5 3 5 4 3 5
= U0 + ∆ U0 + ∆ U0 + ∆ U0 − ∆ U0 + ∆ U0
2 8 16 128 256
5 15 5
= U0 + ( E − 1) U0 + ( E − 1)2 U0 + ( E − 1)3 U0
2 8 16
5 3
− ( E − 1)4 U0 + ( E − 1)5 U0
128 256
5 15 5
= U0 + (U1 − U0 ) + ( E 2 − 2 E + 1) U0 + ( E 3 − 3E 2 + 3E − 1) U0
2 8 16
5
− ( E 4 − 4 E 3 + 6 E 2 − 4 E + 1) U0
128
3
+ [ E 5 − 5E 4 + 10 E 3 − 10 E 2 + 5E − 1] U0
256
5 15 5
= U0 + (U1 − U0 ) + (U2 − 2U1 + U0 ) + (U3 − 3U2 + 3U1 − U0 )
2 8 16
5
− (U4 − 4U3 + 6U2 − 4U1 + U0 )
128
3
+ (U5 − 5U4 + 10U3 − 10U2 + 5U1 − U0 )
256
= U0 1 − +
5 15 5 5 3 5 15 + 15 + 5 + 15
− − − + U1 −
2 8 16 128 256 2 4 16 32 256
+ U2
15 15 15 15 5 + 5 + 15 + U − 5 − 15 + 3 U
− − − + U3 4 5
8 16 64 128 16 32 128 128 256 256
N-26
3 25 75 75 25 3
= U0 − U1 + U2 + U3 − U4 + U5
256 256 128 128 256 256
3 25 75
= (U0 + U5 ) − (U1 + U4 ) + (U2 + U3 )
256 256 128
3 25 75 3a 25b 1 11
= a− b+ c= − + + c
256 256 128 256 256 2 128
1 3 ( a − c ) + 25 ( c − b )
= c+ ⋅
2 256
Solution: We have
∴ ∆3 u x = 0 ⇒ ∆2 u x +1 − ∆2 u x = 0 ⇒ ∆2 u x = ∆2 u x +1 .
u x +(1 /2) = u x + ( u x +1 − u x ) − ∆2 u x + ∆2 u x
1 1 1 1
2 8 2 2
1 1 ∆2 u x ∆2 u x +1
= ux + ( u x +1 − u x ) − +
2 8 2 2
1 1 2
= ( u x + u x +1 ) − ( ∆ u x + ∆2 u x +1 ).
2 16
2. We have ∇ f ( x ) = f ( x ) − f ( x − h).
3. If f ( x ) = 2 x3 − 3 x2 + 3 x − 10,
4. We have ∆3 (1 − x ) (1 − 2 x ) (1 − 3 x )
= ∆3 (1 − 6 x + 11x2 − 6 x3 )
= ( −6) . 3 ! . 1 = − 36.
5. If ∆ U x = x ( x − 1) = x (2), then
x(3) x ( x − 1) ( x − 2)
Ux = + C= + C.
3 3
6. See article 10 part ( f ), remark.
= a . 1 . n ! . 1 + b . 0 = a n !.
9. We have ∆3 f (0 ) = ( E − I )3 f (0 )
= ( E3 − 3E2 + 3E − I ) f (0 ) = E3 f (0 ) − 3E2 f (0 ) + 3E f (0 ) − If (0 )
= 12 − 3 × 8 + 3 × 6 − ( −3) = 12 − 24 + 18 + 3 = 9.
Then ∆ f ( x ) = f ( x + h) − f ( x ) = c − c = 0.
2. We have (1 + ∆ ) (1 − ∇) = EE −1 = 1.
4. We have ∆2 f ( a + h) = ( E − I )2 f ( a + h)
= ( E2 − 2E + I ) f ( a + h)
= E2 f ( a + h) − 2Ef ( a + h) + I f ( a + h)
= ( I − 2 E − 1 + E − 2 ) f ( a + 2 h)
= f ( a + 2h) − 2 f ( a + h) + f ( a).
N-28
6. We have ∆ [ f ( x ) g ( x )] = f ( x + h) g ( x + h) − f ( x ) g ( x )
= f ( x + h) g ( x + h) − f ( x + h) g ( x ) + f ( x + h) g ( x ) − f ( x ) g ( x )
= f ( x + h) [ g ( x + h) − g ( x )] + g ( x ) [ f ( x + h) − f ( x )]
= E f ( x) . ∆ g ( x) + g ( x) . ∆ f ( x)
f ( x) g ( x) . ∆ f ( x) − f ( x) . ∆ g ( x)
7. We have ∆ = ⋅
g ( x) g ( x) . E g ( x)
8. We have ∆3 = ( E − I )3 = E3 − 3E2 + 3E − I .
9. If f ( x ) = x4 − 12 x3 + 24 x2 − 30 x + 9,
then ∆4 f ( x ) = 1 . 4 ! . 1 = 24, if the interval of differencing is 1.
True or False
1. We know that both E and ∆ are linear operators.
2. We know that E = 1 + ∆ or ∆ = E − 1.
∆ [ f ( x ) g ( x )] = E f ( x ) ∆ g ( x ) + g ( x ) ∆ f ( x )
= f ( x + 1) ∆ g ( x ) + g ( x ) ∆ f ( x ).
Hence, the given statement is true.
4. We have ∆ ( e x ) = e x + h − e x = e x ( e h − 1).
ex x ex
∴ ∆ =e and ∆ + C = e x .
h h
e − 1 e − 1
ex
So, if ∆ f ( x ) = e x , then f ( x ) = + C.
h
e −1
❍❍❍
N-29
Chapter-2
Interpolation with Equal Intervals
Comprehensive Problems 1
Problem 1: What do you mean by Interpolation ? What are the underlying assumptions for the validity
of the various methods used for interpolation ?
Problem 3: By constructing a difference table, find the 7th term as well as the general term of the sequence
0, 0, 2, 6, 12, 20, …
1 0
0
2 0 2
2 0
3 2 2
4 0
4 6 2
6 0
5 12 2
8
6 20
Problem 4: In an examination the number of candidates who secured marks between certain limits were as
follows :
Marks : 0 – 19 20 – 39 40 – 59 60 – 79 80 – 99
No. of candidates : 41 62 65 50 17
Estimate the number of candidates getting marks less than 70. (Kumaun 2008, 10, 12)
f ( x)
Marks less than No. of candidates securing ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x )
x marks less than x
20 41
62
40 103 3
65 − 18
60 168 − 15 0
50 − 18
80 218 − 33
17
100 235
Obviously, the number of candidates getting marks less than 70 = f (70 )i. e., f ( x )when x = 70.
By Newton-Gregory forward difference formula, we have
u ( u − 1) 2
f ( a + uh) = f ( a) + u ∆ f ( a) + ∆ f ( a)
2!
u ( u − 1) ( u − 2) 3 u ( u − 1) ( u − 2) ( u − 3) 4
+ ∆ f ( a) + ∆ f ( a) + … .
3! 4!
Problem 5: Find the polynomial which fits the data in the following table :
x : 3 5 7 9 11
y = f ( x) : 6 24 58 108 174
N-31
Solution: The forward difference table for the given data is as given below.
x y = f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x )
3 6
18
5 24 16
34 0
7 58 16
50 0
9 108 16
66
11 174
Here, the first tabulated value of the argument x is 3 and the interval of differencing h is 2. We have
x − 3
f ( x ) = f [3 + ( x − 3)] = f 3 + ⋅2
2
x − 3 1 x − 3 x − 3
= 1 + ∆+ ⋅ − 1 ∆2 f (3)
2 2 ! 2 2
[ ∵ Third and higher differences of f ( x ) are all zero]
x − 3 1 x − 3 x − 5 2
= f (3) + ∆f (3) + ⋅ ∆ f (3)
2 2 2 2
x − 3 ( x − 3) ( x − 5)
= 6 + ⋅ 18 + ⋅ 16
2 8
= 6 + 9 ( x − 3) + 2 ( x − 3) ( x − 5) = 2 x2 − 7 x + 9
Solution: We have
f ( x ) = E x f (0 ) [ ∵ Here, the interval of differencing h is 1]
= (1 − ∇)− x f (0 ) [ ∵ E −1 = 1 − ∇ or E = (1 − ∇)−1 ]
x ( x + 1) 2 x ( x + 1) ( x + 2) 3
= 1 + x ∇ + ∇ + ∇ + … f (0 )
1 . 2 1 . 2 . 3
x ( x + 1) 2 x ( x + 1) ( x + 2) 3
= f (0 ) + x ∇f (0 ) + ∇ f (0 ) + ∇ f (0 ) +
2! 3!
x ( x + 1) ( x + 2) … ( x + n − 1) n
…+ ∇ f (0 ) + … …(1)
n!
N-32
Now ∇r f (0 ) = (1 − E −1 )r f (0 )
r
E − 1
= r −r
f (0 ) = ( E − 1) E f (0 )
E
= ( E − 1)r f [0 + ( − r ) . 1] = ∆r f ( − r ).
∴ ∇f (0 ) = ∆f ( −1), ∇2 f (0 ) = ∆2 f ( −2),
∇2 f (0 ) = ∆3 f ( −3), … .
Problem7: The gamma function, Γ ( x) is tabulated for various values of x in the following table :
1⋅15 1⋅16 1⋅17 1⋅18 1⋅19
x :
Find Γ (1⋅1673).
x f ( x) = Γ ( x) ∆f ( x ) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x )
1⋅15 ⋅93304
− ⋅ 00324
1⋅19 ⋅92088
Problem 8: The following are the marks obtained by 492 candidates in a certain examination.
Marks : No. of candidates
0 – 40 210
40 – 45 43
45 – 50 54
50 – 55 74
55 – 60 32
60 – 65 79
Marks f ( x)
less than No. of candidates ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x ) ∆5 f ( x )
x securing marks
less than x
40 210
43
45 253 11
54 9
50 307 20 − 71
74 − 62 222
55 381 − 42 151
32 89
60 413 47
79
65 492
Here a = 40, h = 5.
By Newton-Gregory forward difference formula, we have
u ( u − 1) 2
f ( a + uh) = f ( a) + u ∆ f ( a) + ∆ f ( a)
2!
u ( u − 1) ( u − 2) 3 u ( u − 1) ( u − 2) ( u − 3) 4
+ ∆ f ( a) + ∆ f ( a)
3! 4!
u ( u − 1) ( u − 2) ( u − 3) ( u − 4) 5
+ ∆ f ( a) + … . …(1)
5!
N-34
(a) Obviously, the number of candidates who secured more than 48 but not more than 50 marks
= the number of candidates who secured 49 or 50 marks
= (the number of candidates who secured less than 51 marks)
− (the number of candidates who secured less than 49 marks)
= f (51) − f (49).
Hence, the number of candidates who secured less than 48 but not less than 45 marks
= f (48) − f (45)
= 280 − 253 = 27.
Problem 9: From the following table of yearly premiums for policies maturing at quinquennial ages,
interpolate the premium for policies maturing at the age of 12 years.
Age (years) x : 10 15 20 25 30 35
Premium f ( x ) : 3 ⋅ 54 3 ⋅ 22 2 ⋅ 91 2 ⋅ 60 2 ⋅ 31 2 ⋅ 04
x f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x ) ∆5 f ( x )
10 3 ⋅ 54
− ⋅ 32
15 3 ⋅ 22 ⋅01
− ⋅ 31 − ⋅ 01
20 2 ⋅ 91 ⋅00 ⋅03
− ⋅ 31 ⋅02 − ⋅ 05
25 2 ⋅ 60 ⋅02 − ⋅ 02
− ⋅ 29 ⋅00
30 2 ⋅ 31 ⋅02
− ⋅ 27
35 2 ⋅ 04
Problem 10: Estimate the sale for 1966 using the following table :
Year : 1931 1941 1951 1961 1971 1981
Sale in thousands : 12 15 20 27 39 52
(Lucknow 2010)
x − 1931
Solution: Let us introduce a new variable u given by u = ;
10
∴ u takes the values 0, 1, 2 , 3, 4.
The forward difference table is as under :
x u f ( u) ∆ f ( u) ∆2 f ( u) ∆3 f ( u) ∆4 f ( u) ∆5 f ( u)
1931 0 12
3
1941 1 15 2
5 0
1951 2 20 2 3
7 3 −10
1961 3 27 5 −7
12 −4
1971 4 39 1
13
1981 5 52
x : 0 ⋅1 0 ⋅2 0 ⋅3 0 ⋅4 0 ⋅5
x 1⋅10517 1⋅ 2140 1⋅ 34986 1⋅ 49182 1⋅ 64872
e :
x f ( x) = e x ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x )
⋅1 1⋅10517
⋅10883
⋅2 1⋅ 2140 ⋅02703
⋅13586 − ⋅ 02093
⋅4 1⋅ 49182 ⋅01494
⋅15690
⋅5 1⋅ 64872
u ( u − 1) 2
Thus e 0⋅243 = f (⋅1) + u ∆ f (⋅1) + ∆ f (⋅1)
2
u ( u − 1) ( u − 2) 3 u ( u − 1) ( u − 2) ( u − 3) 4
+ ∆ f (⋅1) + ∆ f (⋅1) + 0
6 24
1
or e 0⋅243 = 1⋅10517 + ⋅10883 × 1⋅ 43 + × ⋅ 02703 × 1⋅ 43 × ⋅ 43
2
1
+ × ( − ⋅ 02093) × 1⋅ 43 × ⋅ 43 × ( − ⋅ 57)
6
1
+ × ⋅ 02977 × 1⋅ 43 × ⋅ 43 × ( − ⋅ 57) × ( − 1⋅ 57)
24
∴ From (1),
1
e 0⋅411 = 1⋅10517 + ⋅10883 × 3 ⋅11 + × ⋅ 02703 × 3 ⋅11 × 2 ⋅11
2
1
+ × ( − ⋅ 02093) × 3 ⋅11 × 2 ⋅11 × 1⋅11
6
1
+
× ⋅ 02977 × 3 ⋅11 × 2 ⋅11 × 1⋅11 × ⋅11
24
= 1⋅10517 + ⋅ 3384613 + ⋅ 0886867 − ⋅ 0254087 + ⋅ 0008042 = 1⋅ 5083.
Problem 12: Given log x for x = 310, 320, 330, 340, 350 and 360 according to the logarithmic table.
Find the value of log 3375.
Solution: Obtaining values of log x for given values of x from logarithmic table, the difference
table is as under :
x f ( x ) = log x ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x ) ∆5 f ( x )
310 2 ⋅ 49136
⋅01378
360 2 ⋅ 5563
Here a = 310, h = 10 ;
x − a 337 ⋅ 5 − 310
∴ u= = = 2 ⋅ 75.
h 10
2 ⋅ 75 × 1⋅ 75 2
∴ log 337 ⋅ 5 = f ( a) + 2 ⋅ 75 × ∆ f ( a) + ∆ f ( a)
2!
2 ⋅ 75 × 1⋅ 75 × ⋅ 75 3 2 ⋅ 75 × 1⋅ 75 × ⋅ 75 × ( − ⋅ 25) 4
+ ∆ f ( a) + ∆ f ( a)
3! 4!
2 ⋅ 75 × 1⋅ 75
= 2 ⋅ 49136 + 2 ⋅ 75 × ⋅ 01378 + × ( − ⋅ 00042),
2
leaving further terms as their value is very small
= 2 ⋅ 49136 + ⋅ 03792 − ⋅ 00101 = 2 ⋅ 52829.
N-39
Problem 13: The following table gives the amount of a chemical dissolved in water :
Temp : 10° 15° 20° 25° 30° 35°
Solubility : 19 ⋅ 97 21⋅ 51 22 ⋅ 47 23 ⋅ 52 24 ⋅ 65 25 ⋅ 89
x° f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x ) ∆5 f ( x )
(in degree)
10 19 ⋅ 97
1⋅ 54
15 21⋅ 51 − ⋅ 58
⋅96 ⋅67
20 22 ⋅ 47 ⋅09 − ⋅ 68
1⋅ 05 − ⋅ 01 ⋅72
25 23 ⋅ 52 ⋅08 ⋅04
1⋅13 ⋅03
30 24 ⋅ 65 ⋅11
1⋅ 24
35 25 ⋅ 89
Problem 14: Find the value of sin 30 ° 15′ 30 ′ ′ from the following table :
Angle x ° : 30 31 32 33 34
sin x ° : ⋅5000 ⋅5150 ⋅5299 ⋅5446 ⋅5592
x f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x )
30 ⋅5000
⋅0150
31 ⋅5150 − ⋅ 0001
⋅0149 − ⋅ 0001
33 ⋅5446 − ⋅ 0001
⋅0146
34 ⋅5592
− ⋅ 0001 31 31
− 1 − 2
31
+ ×
3! 120 120 120
⋅0002 31 31
− 1 − 2 − 3 + 0
31 31
+
4 ! 120 120 120 120
= ⋅ 50387 = ⋅ 5039.
x: 45° 50 ° 55° 60 °
x f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x )
45° 0 ⋅ 7071
0 ⋅ 0589
50 ° 0⋅7660 −0 ⋅ 0057
0⋅0532 −0 ⋅ 0007
60 ° 0⋅8660
Problem 16: From the following table, for what value of x, y is minimum? Also find this value of y.
x : 3 4 5 6 7 8
x f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x )
3 0.205
0.035
4 0.240 − 0.016
0.019 0.001
5 0.259 − 0.016 − 0.001
0.003 0.001
6 0.260 − 0.015 0.000
– 0.012 0.001
7 0.250 − 0.014
– 0.026
8 0.224
Here a = 3 and h = 1.
By Newton’s forward interpolation formula, we get
u ( u − 1) ( u − 2) 2 u ( u − 1) 3
f ( x ) = f ( a + uh) = f ( a) + u∆ f ( a) + ∆ f ( a) + ∆ y,
2! 3!
neglecting the higher order differences
u( u − 1)
f ( x ) = f ( a + uh) = 0 ⋅ 205 + u(0 ⋅ 035) − ( − 0 .016) …(1)
2
Differentiating (1) w.r.t. u, we get
2u − 1
f ′( x ) = hf ′ ( a + uh) = 0 ⋅ 035 + ( − 0 ⋅ 016)
2
For maximum or minimum, we must have f ′( x ) = 0
2u − 1
0 ⋅ 035 + ( − 0 ⋅016) = 0
2
or u = 2 ⋅6875
∴ x = a0 + uh = 3 + 2 ⋅ 6875 × 1 = 5 ⋅ 6875
Hence f ( x ) is minimum at x = 5.6875.
Putting u = 26865
. in (1), we get the minimum value of f ( x )
1
f ( x ) = 0 ⋅ 205 + (2 ⋅ 6875)(0 ⋅ 035) + (2 ⋅ 6875)(2 ⋅ 6875 − 1)( − 0 ⋅ 016)
2
= 0 ⋅ 2628.
True or False
1. See article 1.
2. See article 1.
3. See article 1.
❍❍❍
N-44
Chapter-3
Interpolation with Unequal
Intervals of the Argument
Comprehensive Problems 1
Problem 1 (a): Given that f (0 ) = 8, f (1) = 68, f (5) = 123, construct a divided difference table,
using the table determine the value of f (2). (Kumaun 2012)
x f ( x) ∆| f ( x ) ∆| 2 f ( x )
0 8
68 − 8
= 60
1−0
1 68 13 ⋅ 75 − 60
= − 9 ⋅ 25
123 − 68 5−0
= 13 ⋅ 75
5 −1
5 123
x log x = f ( x ) ∆| f ( x ) ∆| 2 f ( x ) ∆| 3 f ( x )
2 ⋅30103
⋅17609
3 ⋅47712 − ⋅ 0333883
⋅075925 ⋅0082846
5 ⋅62897 ⋅008035
⋅108065
7 ⋅84510
N-45
× ( − ⋅ 0333883) + (4 ⋅ 7 − 2) (4 ⋅ 7 − 3) (4 ⋅ 7 − 5) × ⋅ 0082846
= 0 ⋅ 611812 .
1
Problem 3: Show that nth order divided difference f ( x0 , x1, x2 , … , x n ) for the function ux = is equal
x
( − 1)n
to .
x0 x1 x2 … x n (Kumaun 2013)
Solution: We have
f ( x1 ) − f ( x0 ) 1 1 1 1 x0 − x1
f ( x0 , x1 ) = = − =
x1 − x0 x1 − x0 x1 x0 x1 − x0 x0 x1
( − 1) ( x1 − x0 ) ( − 1)
= = …(1)
( x1 − x0 ) x0 x1 x0 x1
f ( x1, x2 ) − f ( x0 , x1 )
Again f ( x0 , x1, x2 ) =
x2 − x0
1 ( −1) ( −1)
= − , [Using (1)]
x2 − x0 x1 x2 x0 x1
( −1) x0 − x2 ( − 1)2 ( x2 − x0 )
= =
x2 − x0 x0 x1 x2 ( x2 − x0 ) x0 x1 x2
( − 1)2
= …(2)
x0 x1 x2
1 ( −1)2 ( − 1)2
= − , [(using (2)]
x3 − x0 x1 x2 x3 x0 x1 x2
( −1)2 ( x0 − x3 ) ( −1)3
= = .
( x3 − x0 ) x0 x1 x2 x3 x0 x1 x2 x3
x : 5 11 27 34 42
x f ( x) ∆| f ( x ) ∆| 2 f ( x ) ∆| 3 f ( x ) ∆| 4 f ( x )
5 23
146
11 899 40
1
1026
27 17315 69 0
1
2613
34 35606 100 0
1
4113
42 68510 a 0
1
b
3 c a′ 0
1
b′
3 c a′ ′
b′ ′
3 c
= c + ( x − 3) b ′ + ( x − 3)2 a′ ′ + ( x − 3)3 × 1
= − 13 + ( x − 3) . 2 + ( x − 3)2 . 6 + ( x − 3)3
Problenm 4(ii): The observed values of a function are respectively 168, 120, 72 and 63 at the four
positions 3, 7, 9 and10 of the independent variable. What is the best estimate you can give for the value of the
function at the position 6 of the independent variable ? (Purvanchal 2011; Bundelkhand 12)
x f ( x) ∆| f ( x ) ∆| 2 f ( x ) ∆| 3 f ( x )
3 168
120 − 168
= − 12
7−3
7 120 − 24 + 12
= −2
72 − 120 9−3
= − 24
9−7 5+2
=1
10 − 3
9 72 − 9 + 24
=5
63 − 72 10 − 7
= −9
10 − 9
10 63
Hence f ( x ) = f ( x0 ) + ( x − x0 ) ∆| f ( x0 ) + ( x − x0 ) ( x − x1 ) ∆| 2 f ( x0 )
+ ( x − x0 ) ( x − x1 ) ( x − x2 ) ∆| 3 f ( x0 )
= x3 − 21 x2 + 119 x − 27.
= 147.
Problem 5: The following table gives the normal weights of babies during the first twelve months of life :
Age in months ( x ) : 0 2 5 8 10 12
1 1
Weight in lbs f ( x ) : 7 10 15 16 18 21
2 4
(7 − 2) (7 − 5) (7 − 8) (7 − 10 ) (7 − 12) 15
f (7) = ×
(0 − 2) (0 − 5) (0 − 8) (0 − 10 ) (0 − 12) 2
(7 − 0 ) (7 − 5) (7 − 8) (7 − 10 ) (7 − 12) 41
+ ×
(2 − 0 ) (2 − 5) (2 − 8) (2 − 10 ) (2 − 12) 4
(7 − 0 ) (7 − 2) (7 − 8) (7 − 10 ) (7 − 12)
+ × 15
(5 − 0 ) (5 − 2) (5 − 8) (5 − 10 ) (5 − 12)
(7 − 0 ) (7 − 2) (7 − 5) (7 − 10 ) (7 − 12)
+ × 16
(8 − 0 ) (8 − 2) (8 − 5) (8 − 10 ) (8 − 12)
(7 − 0 ) (7 − 2) (7 − 5) (7 − 8) (7 − 12)
+ × 18
(10 − 0 ) (10 − 2) (10 − 5) (10 − 8) (10 − 12)
(7 − 0 ) (7 − 2) (7 − 5) (7 − 8) (7 − 10 )
+ × 21
(12 − 0 ) (12 − 2) (12 − 5) (12 − 8) (12 − 10 )
15 287 175 63 21
= − +5+ − +
128 384 12 16 32
= ⋅11718 − ⋅ 7473 + 5 + 14 ⋅ 583 − 3 ⋅ 937 + ⋅ 652 = 15 ⋅ 67 lbs.
Problem 6: Determine by Lagrange’s formula the percentage number of criminals under 35 years :
= −10 ⋅ 4 + 50 ⋅ 475 + 42 ⋅ 05 − 4 ⋅ 72
= 77 ⋅ 405 = 77%
N-49
x Ux ∆| U x ∆| 2 U x ∆| 3 U x
0 8
11 − 8
=3
1−0
1 11 19 − 3
=4
68 − 11 4−0
= 19
4 −1 9−4
=1
5−0
4 68 55 − 19
=9
123 − 68 5 −1
= 55
5−4
5 123
= 8 + 3 x + 4 x2 − 4 x + x3 − 5 x2 + 4 x = x3 − x2 + 3 x + 8.
Problem 8: Using Lagrange’s formula find the form of the function given by f ( x) :
x : 3 2 1 −1
f ( x) : 3 12 15 −21
(Kanpur 2008)
3 2
Solution: Proceed as in Ex. 10. Ans. x − 9 x + 17 x + 6.
Solution: Here we are required to find y0 i. e., the value of y for x = 0 and given y−3 , y−1 ,
y1 , y3 i. e., the values of y x for x = − 3, − 1, 1 and 3 i. e., x0 = − 3, x1 = − 1, x2 = 1 and x3 = 3.
Using Lagrange’s interpolation formula, we have
N-50
(0 + 1) (0 − 1) (0 − 3) (0 + 3) (0 − 1) (0 − 3)
y0 = ⋅ y− 3 + ⋅ y− 1
( −3 + 1) ( −3 − 1) ( −3 − 3) ( −1 + 3) ( −1 − 1) ( −1 − 3)
(0 + 3) (0 + 1) (0 − 3) (0 + 3) (0 + 1) (0 − 1)
+ ⋅ y1 + ⋅ y3
(1 + 3) (1 + 1) (1 − 3) (3 + 3) (3 + 1) (3 − 1)
1 9 9 1
=− y− 3 + y− 1 + y1 − y3
16 16 16 16
y−3 + +
1 1 1 1 1 1
=− y + + y − y
16 2 16 −1 2 16 1 16 3
1 1
= ( y1 + y−1 ) − [( y3 − y1 ) − ( y−1 − y−3 )]
2 16
x y = f ( x) ∆| f ( x ) ∆| 2 f ( x ) ∆| 3 f ( x )
– 1⋅941 1⋅9
0⋅1062699
0⋅0942507 0⋅0004869
0⋅0842459
1⋅248 2⋅2
+ ( x − x0 ) ( x − x1 ) ( x − x2 ) ∆| 3 f ( x0 )
Problem 11: If y35⋅0 = 1175, y35⋅5 = 1280, y39⋅5 = 2180, y40⋅5 = 2420, find y40 by
divided differences.
Solution: The divided difference table is as given below :
x yx ∆| y x ∆| 2 y x ∆| 3 y x
35 1175
1280 − 1175
= 210
35 ⋅ 5 − 35
40 ⋅ 5 2420
+ ( x − x0 ) ( x − x1 ) ( x − x2 ) y ( x0 , x1, x2 , x3 ).
Substituting the values from the table and taking x = 40, we get
y40 = 1175 + (40 − 35) × 210 + (40 − 35) (40 − 35 ⋅ 5) × 3 ⋅ 33
Problem 12: Apply Lagrange’s formula to find f (1⋅ 50 ), using the following values of the function
f ( x) :
x : 1⋅ 00 1⋅ 20 1⋅ 40 1⋅ 60 1⋅ 80 2 ⋅ 00
Solution: Using the Lagrange’s interpolation formula for the given values, we have
N-52
= 0 ⋅12693.
Problem 13: Apply Lagrange’s formula to find the cubic polynomial which includes the following values of
x and y x :
x : 0 1 4 6
yx : 1 −1 1 −1
( x − 1) ( x − 4) ( x − 6) ( x − 0 ) ( x − 4) ( x − 6)
yx = ⋅1 + ⋅ ( −1)
(0 − 1) (0 − 4) (0 − 6) (1 − 0 ) (1 − 4) (1 − 6)
( x − 0 ) ( x − 1) ( x − 6) ( x − 0 ) ( x − 1) ( x − 4)
+ ⋅1 + ⋅ ( −1)
(4 − 0 ) (4 − 1) (4 − 6) (6 − 0 ) (6 − 1) (6 − 4)
x3 − 11 x2 + 34 x − 24 x3 − 10 x2 + 24 x x3 − 7 x2 + 6 x x3 − 5 x2 + 4 x
= − + −
−24 15 −24 60
17 x3 29 x2 207 x 5
=− + + + ⋅
120 24 120 4
Problem 14: Given the following data, find f ( x) as a polynomial in powers of ( x − 5), by extending the
table to include argument x = 5 repeated as many times as may be necessary :
f (0 ) = 4, f (2) = 26, f (3) = 58, f (4) = 112 , f (7) = 466, f (9) = 922 .
Also find f ′ (5), f ′ ′(5) and f ′ ′ ′(5).
N-53
x f ( x) ∆| f ( x ) ∆| 2 f ( x ) ∆| 3 f ( x )
0 4
26 − 4
= 11
2−0
2 26 32 − 11
=7 11 − 7
3−0 =1
58 − 26 4−0
= 32
3−2
3 58 54 − 32
= 11 16 − 11
4−2 =1
112 − 58 7−2
= 54
4−3
4 112 118 − 54
= 16 22 − 16
7−3 =1
466 − 112 9−3
= 118
7−4
7 466 228 − 118
= 22 a − 22
9−4 =1
922 − 466 5−4
= 228
9−7
9 922 b − 228
=a
c − 922 5−7
=b
5−9 a′ − a
=1
5−7
5 c b′−b
= a′
5−9
a′ ′ − a′
b′ =1
5−9
5 c a′ ′
b ′′
5 c
From divided difference table we observe that the third divided differences of f ( x ) are constant,
therefore f ( x ) is a polynomial of degree 3. As f ( x ) is to be expressed as a polynomial in powers of
( x − 5) therefore argument x = 5 is to be repeated three times in the divided difference table.
In the extended divided difference table the values of c , b, b ′, a, a′, a′ ′ are obtained as follows :
I. Complete the column V by the third difference of f ( x ) which is constant and 1 in this
case.
a − 22
II. From column V, =1 giving a = 23
5−4
a′ − a a′ − 23
then =1 or =1 giving a′ = 21
5−7 −2
a′ ′ − a′ a′ ′ − 21
and =1 or =1 giving a′ ′ = 17.
5−9 −4
N-54
= c + ( x − 5) b ′ + ( x − 5)2 a′ ′ + ( x − 5)3 × 1
Problem 15: In the following table h is the height above the sea level and p the barometric pressure ;
calculate p when h = 5280.
p : 27 25 23 20
(5 − 2) (5 − 3) (5 − 4) (5 − 7) (5 − 1) (5 − 3) (5 − 4) (5 − 7)
f (5) = ×2+ ×4
(1 − 2) (1 − 3) (1 − 4) (1 − 7) (2 − 1) (2 − 3) (2 − 4) (2 − 7)
(5 − 1) (5 − 2) (5 − 4) (5 − 7) (5 − 1) (5 − 2) (5 − 3) (5 − 7)
+ ×8+ × 16
(3 − 1) (3 − 2) (3 − 4) (3 − 7) (4 − 1) (4 − 2) (4 − 3) (4 − 7)
(5 − 1) (5 − 2) (5 − 3) (5 − 4)
+ × 128
(7 − 1) (7 − 2) (7 − 3) (7 − 4)
2 32 128 128
=− + − 24 + +
3 5 3 15
= − ⋅ 66 + 6 ⋅ 4 − 24 + 42 ⋅ 66 + 8 ⋅ 53
= 57 ⋅ 59 − 24 ⋅ 66 = 32 ⋅ 93.
From the given values of function we observe that the function is of the form 2 x i. e., f ( x ) = 2 x
giving f (5) = 25 = 32 which differs from the interpolated value 32 ⋅ 93 due to the fact that for
applying Lagrange’s formula we have assumed the given function to be expressible as a
polynomial while actually the given function is of the form 2 x .
Problem 17: Four equidistant values u−1 , u0 , u1 and u2 being given, a value is interpolated by
Lagrange’s formula. Show that it may be written in the form
( y2 − 1) 2 x ( x2 − 1) 2
u x = yu0 + xu1 + y ∆ u−1 + ∆ u0 ,
3! 3!
Solution: We have
∆2 u−1 = ( E − 1)2 u−1= ( E 2 − 2 E + 1) u−1
= u1 − 2u0 + u−1 .
Similarly,
∆2 u0 = ( E 2 − 2E + 1) u0 = u2 − 2u1 + u0 .
x ( x − 1) ( x − 2) x ( x − 1) ( x − 2) ( x − 1) x ( x + 1)
=− u−1 + − ( x − 1) + + u
6 3 6 0
x ( x − 1) ( x − 2) ( x − 1) x ( x + 1) ( x − 1) x ( x + 1)
+ x − − u1 + u2
6 3 6
x ( x − 1) ( x − 2) ( x − 2) ( x − 1) ( x + 1)
=− u−1 + u0
6 2
( x − 2) x ( x + 1) ( x − 1) x ( x + 1)
− u1 + u2 .
2 6
…(1)
For u−1 , u0 , u1 , u2 as known values, by Lagrange’s formula, we have
( x − 0 ) ( x − 1) ( x − 2) ( x + 1) ( x − 1) ( x − 2)
ux = u−1 + u0
( −1 − 0 ) ( −1 − 1) ( −1 − 2) (0 + 1) (0 − 1) (0 − 2)
( x + 1) ( x − 0 ) ( x − 2) ( x + 1) ( x − 0 ) ( x − 1)
+ u1 + u2
(1 + 1) (1 − 0 ) (1 − 2) (2 + 1) (2 − 0 ) (2 − 1)
x ( x − 1) ( x − 2) ( x + 1)( x − 1) ( x − 2)
=− u−1 + u0
6 2
( x + 1) x ( x − 2) ( x + 1) x ( x − 1)
− u1 + u2 .
2 6
…(2)
From (1) and (2), we have
y ( y2 − 1) 2 x ( x2 − 1) 2
u x = y u0 + xu1 + ∆ u−1 + ∆ u0 .
3! 3!
Problem 18: Prove that the Lagrange’s formula can be put in the form
n φ ( x ) f ( xr )
Pn ( x ) = Σ
r = 0 ( x − xr ) φ′ ( xr )
n
where φ ( x ) = Π ( x − xr ).
r =0 (Rohilkhand 2011)
φ ( x ) f ( xr )
n 1
= Σ ( x − x ) ( x − x ) … ( x − x ) ⋅
r = 0 ( x − xr ) r 0 r 1 r n
n
Now φ ( x ) = Π ( x − xr ).
r =0
∴ φ′( xr ) = [ φ′( x )] x = x r
= ( xr − x0 ) ( xr − x1 ) … ( xr − xr −1 )( xr − xr +1 ) … ( xr − xn ).
n φ ( x ) f ( xr )
∴ Pn ( x ) = Σ ⋅
r = 0 ( x − xr ) φ′ ( xr )
N-57
xr , xr +1 of degree n − 1.
4. We have f (2) = 23 − 2 × 2 = 4, f (4) = 43 − 2 × 4 = 56.
f (4) − f (2) 56 − 4
Now f (2 , 4) = = = 26.
4−2 2
5. See article 2.
2. The nth divided difference can be expressed as the quotient of two determinants each of
order ( n + 1).
3. See article 4.
4. See article 7.
f ( b ) − f ( a) (1 / b2 ) − (1 / a2 )
5. We have f ( a, b ) = =
b−a b−a
a2 − b2 (a − b) (a + b) a+ b
= = =−
a b2 ( b − a)
2
a2 b2 ( b − a) a2 b2
True or False
1. We know that Lagrange’s formula can be used for both equal and unequal intervals.
2. We know that divided differences are symmetric functions of their arguments. So, the
value of any divided difference is independent of the order of the arguments.
f ( b ) − f ( a)
, then ∆| =
1 1
3. If f ( x) =
x b a b−a
(1 / b ) − (1 / a) a−b 1
= = =− ⋅
b−a ab ( b − a) ab
N-58
f ( x1 ) − f ( x 0 )
4. We have f ( x0 , x1 ) = ⋅
x1 − x 0
5. We know that the nth divided differences of a polynomial of degree n are constant.
6. We know that the divided differences can be expressed as the product of multiple integrals.
7. The given formula is Newton’s divided differences interpolation formula and not
Lagrange’s interpolation formula.
❍❍❍
N-59
Chapter-4
Central Difference Interpolation Formulae
Comprehensive Problems 1
Problem 1: Prove that
1 2
(i) δ = ∆ E −1 /2 = ∇ E1 /2 (ii) µ 2 = 1 + δ
4
(Rohilkhand 2011;
Gorakhpur 11; Bundelkhand 11; Purvanchal 12)
1 1
(iii) E 1 /2 = µ + δ (iv) E −1 /2 = µ − δ
2 2
(Rohilkhand 2011)
1
(v) √ (1 + δ2µ 2 ) = 1 + δ2 (vi) δ2 y0 = y1 − 2 y0 + y−1
2
(Kumaun 2012)
n n!
n k
(vii) δ y x = Σ ( −1) y x +(n /2)− k .
k =0 k ! (n − k ) !
= ( E 1 /2 − E −1 /2 ) f ( x ).
∴ δ = E 1 /2 − E −1 /2 = ( E − 1) E −1 /2
= ∆ E −1 /2 . [ ∵ E = 1 + ∆ or E − 1 = ∆ ]
Again δ = E1 /2 − E −1 /2 ⇒ δ = (1 − E −1 ) E 1 /2 .
∴ δ = ∇ E 1 /2 [∵ 1 − E −1 = ∇]
Hence δ = ∆ E −1 /2 = ∇ E 1 /2 .
(ii) We have δ = E 1 /2 − E −1 /2
1
and µ= ( E 1 /2 + E −1 /2 ).
2
1 1
∴ µ2 = ( E 1 /2 + E −1 /2 )2 = [( E 1 /2 − E −1 /2 )2 + 4 E 1 /2 E −1 /2 ]
4 4
1
= [ 4 + ( E 1 /2 − E −1 /2 )2 ]
4
1 1
= (4 + δ2 ) = 1 + δ2 .
4 4
N-60
= ( E − 2 E 1 /2 E −1 /2 + E −1 ) y0
= ( E − 2 + E −1 ) y0
= Ey0 − 2 y0 + E −1 y0
= y1 − 2 y0 + y−1.
(vii) We have δ n y x = ( ∆ E −1 /2 )n y x [∵ δ = ∆E −1 /2 ]
= ∆ n E − n /2 y x
= ∆ n y x −(n /2)
(ii) ∆ + ∇ = ∆ / ∇ − ∇ / ∆.
Solution: (i) We have δ = E 1 /2 − E −1 /2
1
and µ= ( E 1 /2 + E −1 /2 ).
2
1
∴ µδ= ( E 1 /2 + E −1 /2 ) ( E 1 /2 − E −1 /2 )
2
1
= [ ( E 1 /2 )2 − ( E −1 /2 )2 ]
2
1
= ( E − E −1 )
2
1
= [(1 − E −1 ) + ( E − 1)]
2
1
= ( ∇ + ∆ ). [∵ E = 1 + ∆ and E −1 = 1 − ∇ ]
2
1
Again µ δ = [(1 − E −1 ) + ( E − 1)], as proved above
2
1 1
= (1 − E −1 ) + ( E − 1)
2 2
1 1
= ( E − 1) E −1 + ∆ [∵ E − 1 = ∆ ]
2 2
1 1
=∆E −1 + ∆.
2 2
1 1 1
Hence µ δ = ( ∇ + ∆ ) = ∆E −1 + ∆.
2 2 2
(ii) We have ∆ + ∇ = ( E − 1) + (1 − E −1 )
= E − E −1 ...(1)
−1
∆ ∇ E −1 1− E
Also − = −
∇ ∆ 1 − E −1 E −1
E (1 − E −1 ) E −1 ( E − 1)
= −
−1 E −1
1− E
= E − E −1 ...(2)
Problem 3: If D, E, δ and µ be the operators with usual meanings and if hD = U , where h is the interval of
differencing , prove that
(i) E = eU (ii) e − U = 1 − ∇
h h2 2 h3 3
= 1 f ( x) + Df ( x ) + D f ( x) + D f ( x) + …
1! 2! 3!
hD ( hD)2 ( hD)3
= 1 + + + + … f ( x )
1! 2! 3!
= e hD f ( x ).
∴ E = e hD = e U , where hD = U .
∴ ( e U )−1 = E −1
or e −U = 1 − ∇ [∵ E −1 = 1 − ∇]
u yu ∆ yu ∆ 2 yu ∆ 3 yu
−1 10
−5
0 5 8
3 −9
1 8 −1
2
2 10
N-63
12516 − 12520
∴ u= = − 0 ⋅ 4.
10
We require the value of y for u = −0 ⋅ 4.
The difference table is given below :
u yu ∆ yu ∆ 2 yu ∆ 3 yu
−2 111⋅ 803399
0 ⋅ 044712
1 111⋅ 937483
(0 ⋅ 6) ( −0 ⋅ 4)
+ ( −0 ⋅ 000018)
2
(0 ⋅ 6) ( −0 ⋅ 4) ( −1⋅ 4)
+ ( −0 ⋅ 000001)
6
= 111⋅ 892806 − 0 ⋅ 017878 + ⋅ 000002 + negligible term
= 111⋅ 874930.
N-64
Problem 6: Apply (i) Stirling’s formula , (ii) Bessel’s formula (iii) Gauss forward formula,
(iv) Laplace-Everett’s formula to find the value of f (1⋅ 22) from the following table which gives the values of
1 2
1 t − t
f ( t) =
∫
√ (2 π ) 0
e 2 dt at intervals of t = 0 ⋅ 5 from t = 0 to t = 3 :
t f ( t)
0 ⋅ 00 0 ⋅ 000
0 ⋅ 50 0 ⋅191
1⋅ 00 0 ⋅ 341
1⋅ 50 0 ⋅ 433
2 ⋅ 00 0 ⋅ 477
2 ⋅ 50 0 ⋅ 494
3 ⋅ 00 0 ⋅ 499
Solution: Take 1⋅ 00 as the origin and 0 ⋅ 5 as the unit, thus the new values of t are
1⋅ 22 − 1⋅ 00
−2 , − 1, 0, 1, 2 , 3, 4. When t = 1⋅ 22 , the new value of t = = 0 ⋅ 44 .
0 ⋅5
Construct the following difference table :
t f ( t) ∆ f ( t) ∆2 f ( t) ∆3 f ( t) ∆4 f ( t) ∆5 f ( t) ∆6 f ( t)
−2 0 ⋅ 000
0 ⋅191
−1 0 ⋅191 −0 ⋅ 041
0 ⋅150 −0 ⋅ 017
0 0 ⋅ 341 −0 ⋅ 058 0 ⋅ 027
0 ⋅ 092 0 ⋅ 010 −0 ⋅ 016
1 0 ⋅ 433 −0 ⋅ 048 0 ⋅ 011 −0 ⋅ 001
0 ⋅ 044 0 ⋅ 021 −0 ⋅ 017
2 0 ⋅ 477 −0 ⋅ 027 −0 ⋅ 006
0 ⋅ 017 0 ⋅ 015
3 0 ⋅ 494 −0 ⋅ 012
0 ⋅ 005
4 0 ⋅ 499
t ( t 2 − 12 ) ∆3 f ( −1) + ∆3 f ( −2) t 2 ( t 2 − 12 ) 4
+ + ∆ f ( −2)
3! 2 4!
= 0 ⋅ 341 + (⋅44) (0 ⋅ 092) + (⋅22) (⋅56) (0 ⋅ 058) − (⋅24) (⋅44) (⋅56) (0 ⋅ 010 )
+ (⋅06) (⋅44) (⋅56) (0 ⋅ 027) − (2 ⋅ 44) (⋅06) (⋅44) (⋅56) (1⋅ 56) (⋅0032)
= ⋅ 389 approx.
N-66
u ( u2 − 1) 2 u ( u2 − 1) ( u2 − 4) 4
f ( t) = u f (0 ) + ∆ f ( −1) + ∆ f ( −2)
3! 5!
t ( t2 − 1) 2 t ( t2 − 1) ( t2 − 4) 4
+ t f (1) + ∆ f (0 ) + ∆ f ( − 1) ,
3 ! 5 !
where u = (1 − t).
(⋅56) (⋅562 − 1)
∴ f (0 ⋅ 44) = ⋅56 × 0 ⋅ 341 + ( − 0 ⋅ 058)
6
(⋅44) (⋅442 − 1)
+ ⋅44 × 0 ⋅ 433 + ( − 0 ⋅ 048)
3!
(⋅44) (⋅442 − 1) (⋅442 − 4)
+ (0 ⋅ 011)
5!
= 0 ⋅ 389 approximately.
Problem 7: The values of e − x for certain equidistant values of x are given in the following table. Find the
value of e − x when x = 1.7489 by Bessel’s as well as Stirling formula.
x e− x
1.72 0.1790661479
1.73 0.1772844100
1.74 0.1755204006
1.75 0.1737739435
1.76 0.1720448638
1.77 0.1703329888
1.78 0.1686381473
Problem 8: Use Gauss interpolation formula to find y41 with the help of the following data :
u yu ∆ yu ∆ 2 yu ∆ 3 yu ∆ 4 yu
−2 3678 ⋅ 2
−683 ⋅1
−1 2995 ⋅1 88 ⋅1
−595 ⋅ 0 −17
0 2400 ⋅1 71⋅1 9⋅9
−523 ⋅ 9 −7 ⋅1
1 1876 ⋅ 2 64 ⋅ 0
−459 ⋅ 9
2 1416 ⋅ 3
u yu ∆ yu ∆ 2 yu ∆ 3 yu
−1 0 ⋅ 2707
0 ⋅ 0320
0 0 ⋅ 3027 0 ⋅ 0039
0 ⋅ 0359 0 ⋅ 0010
1 0 ⋅ 3386 0 ⋅ 0049
0 ⋅ 0408
2 0 ⋅ 3794
Problem 10: From the following table find the value of f(⋅5437) by Gauss formula :
x f ( x)
⋅51 ⋅529244
⋅52 ⋅537895
⋅53 ⋅546464
⋅54 ⋅554939
⋅55 ⋅563323
⋅56 ⋅571616
⋅57 ⋅579816
u yu ∆ yu ∆2 yu ∆3 yu ∆4 yu ∆5 yu ∆6 yu
–3 ⋅529244
⋅008651
–2 ⋅537895 – ⋅000082
⋅008569 – ⋅000012
–1 ⋅546464 – ⋅000094 ⋅000015
⋅008475 ⋅000003 – ⋅000018
0 ⋅554939 – ⋅000091 – ⋅000003 ⋅000019
⋅008384 ⋅000000 ⋅000001
1 ⋅563323 – ⋅000091 – ⋅000002
⋅008293
2 ⋅571616 – ⋅000093 – ⋅000002
⋅008200
3 ⋅579816
δ f ( x ) = f x + h − f x − h
1 1
2 2
1 1
f x + h + f x − h , where h is the interval of differencing.
1
and µ f ( x) =
2 2 2
Second part: We have
L.H.S. = µ f ( x ) δ g ( x ) + µ g ( x ) δ f ( x )
1 1 1 1 1
= [ f ( x + h) + f ( x − h)] [ g( x + h) − g( x − h)]
2 2 2 2 2
1 1 1 1 1
+ [ g ( x + h) + g ( x − h)] [ f ( x + h) − f ( x − h)]
2 2 2 2 2
1 1 1 1 1 1
= f ( x + h) g ( x + h) − f ( x − h) g ( x − h)
2 2 2 2 2 2
1 1 1 1 1 1
+ g ( x + h) f ( x + h) − g ( x − h) f ( x − h),
2 2 2 2 2 2
the other terms in the two products cancel each other
N-70
1 1 1 1
= f (x + h) g ( x + h) − f ( x − h) g ( x − h)
2 2 2 2
= δ [ f ( x ) g ( x )].
Problem 12: Apply Bessel’s formula to find a polynomial of degree three or less which takes the following
values of the function u x :
x : 4 6 8 10
f ( x) : 1 3 8 20
x −6
Solution: Let u = .
2
The difference table is given below :
x u f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x )
4 −1 1
6 0 3 3
5 4
8 1 8 7
12
10 2 20
Bessel’s formula is
1 1 u ( u − 1) ( ∆2 y−1 + ∆2 y0 )
f ( x) = ( y0 + y1 ) + ( u − ) ∆y0 +
2 2 2 2
u − 1 u u − 1)
2
+ ∆3 y−1 + ...
6
u − 1 u ( u − 1)
1 1 u ( u − 1) (3 + 7) 2
= (3 + 8) + ( u − ) 5 + + ×4
2 2 2 2 6
11 5 5 (2u − 1) u ( u − 1)
= + (2u − 1) + ( u2 − u) +
2 2 2 3
2 3 3 2 17
= u + u + u+ 3
3 2 6
3 2
2 x − 6 3 x − 6 17 x − 6
= + + +3
3 3 2 2 6 2
1 3 17
= ( x − 6)3 + ( x − 6)2 + ( x − 6) + 3.
12 8 12
1 3 17
∴ f ( x) = ( x − 6)2 + ( x − 6)2 + ( x − 6) + 3.
12 8 12
N-71
θ u yu ∆ yu ∆2 yu ∆3 yu ∆4 yu ∆5 yu ∆6 yu
0° −3 0 ⋅ 0000
⋅0875
5° −2 0 ⋅ 0875 ⋅0013
⋅0888 ⋅0015
10° −1 0 ⋅1763 ⋅0028 ⋅0002
⋅0916 ⋅0017 − ⋅ 0002
15° 0 0 ⋅ 2679 ⋅0045 ⋅0000 ⋅0011
⋅0961 ⋅0017 ⋅0009
20° 1 0 ⋅ 3640 ⋅0062 ⋅0009
⋅1023 ⋅0026
25° 2 0 ⋅ 4663 ⋅0088
⋅1111
30° 3 0 ⋅ 5774
Stirling’s formula is
1 u2 2
yu = y0 + u . [ ∆y0 + ∆y−1 ] + ∆ y− 1
2 2
u ( u2 − 1) 1 3 u2 ( u2 − 1) 4
+ [ ∆ y−1 + ∆3 y−2 ] + ∆ y− 2
3! 2 4!
u ( u2 − 1) ( u2 − 4) 1 5 u2 ( u2 − 1) ( u2 − 4) 6
+ [ ∆ y−2 + ∆5 y−3 ] + ∆ y− 3 .
5! 2 6!
⋅2 (⋅2)2
∴ y⋅2 = 0 ⋅ 2679 + [⋅0961 + ⋅ 0916] + (⋅0045)
2 2
(⋅2) (⋅04 − 1) [⋅0017 + ⋅ 0017] (⋅04) (⋅04 − 1)
+ + (⋅0000 )
6 2 24
(⋅2) (⋅04 − 1) (⋅04 − 4) ( − ⋅ 0002 + ⋅ 0009)
+ ⋅
120 2
(⋅04) (⋅04 − 1) (⋅04 − 4)
+ (⋅0011)
720
= 0 ⋅ 2679 + ⋅ 01877 + ⋅ 00009 + negligible quantities
= 0 ⋅ 28676.
Thus the estimated value of tan 16° = 0 ⋅ 28676.
N-72
δ2
µ2 = 1 + ⋅
4
1
4. We know that σ = ⋅
δ
E 1 /2
5. We know that σ = ⋅
E −1
True or False
1. The given formula is Gauss’s forward formula.
2. We know that Bessel’s interpolation formula is the mean of Gauss’s forward formula and
the third formula due to gauss.
3. We know that ∆ − ∇ = δ2 . See Ex. 1 part (iii).
4. See article 1.
❍❍❍
N-73
Chapter-5
Solution of Algebraic and
Transcendental Equations
Comprehensive Problems 1
Problem 1: Find the root of x2 − 5 x + 2 = 0 correct to five decimal places which lies between 4 and 5 by
using Newton-Raphson method. (Kumaun 2007, 08, 11, 13; Kanpur 08; 11; Purvanchal 08)
= 22 ⋅ 25 − 22 ⋅ 5, which is −ve
and f (4 ⋅ 6) = (4 ⋅ 6)2 − 5 × 4 ⋅ 6 + 2 = 21⋅16 − 23 + 2
2 xn2 − 5 xn − xn2 + 5 xn − 2 x 2 −2
= = n ⋅
2 xn − 5 2 xn − 5
xn2 − 2
Thus, xn +1 = ⋅ …(2)
2 xn − 5
Problem 2: Find the root of log x = cos x , by Newton-Raphson method up to five decimal places.
(Lucknow 2009)
for x = 1⋅ 3 and f ( x ) = + ive for x = 1⋅ 4, therefore the root lies between 1⋅ 3 and 1⋅ 4.
Taking x0 = 1⋅ 3 and applying Newton-Raphson formula, we have
f ( x0 ) log e 1⋅ 3 − cos 1⋅ 3
x1 = x0 − = 1⋅ 3 −
f ′ ( x0 ) {1 / (1⋅ 3)} + sin 1⋅ 3
⋅11394 − ⋅ 26750
= 1⋅ 3 − = 1⋅ 3029.
⋅76923 + ⋅ 96356
Now taking x1 = 1⋅ 3029 as approximate root and repeating Newton-Raphson formula, we have
f ( x1 ) log e 1⋅ 3029 − cos (1⋅ 3029)
x2 = x1 − = 1⋅ 3029 − = 1⋅ 30295.
f ′ ( x1 ) 1 / (1⋅ 3029) + sin (1⋅ 3029)
Now taking x2 = 1⋅ 30295 as approximate root and repeating Newton-Raphson formula, we have
f ( x2 ) log e 1⋅ 30295 − cos (1⋅ 30295)
x3 = x2 − = 1⋅ 30295 −
f ′ ( x2 ) 1 / (1⋅ 30295) + sin (1⋅ 30295)
= 1⋅ 30295.
As x2 = x3 , therefore the root is 1⋅ 30295.
1
Problem 3: Find a root of x = + sin x near x = 1⋅ 5.
2
Now taking x1 = 1⋅ 4972 as approximate root and repeating Newton-Raphson formula, we have
f ( x1 ) 1⋅ 4972 − ⋅ 5 − sin (1⋅ 4972)
x2 = x1 − = 1⋅ 4972 −
f ′ ( x1 ) 1 − cos (1⋅ 4972)
1⋅ 4972 − ⋅ 5 − sin 85 ⋅ 82° ∵ 22 = π = 180 °
= 1⋅ 4972 −
1 − cos 85 ⋅ 82° 7
Problem 4: The equation x4 − 5 x3 − 12 x2 + 76 x − 79 = 0 has two roots close to x = 2 . Find the roots
correct upto four decimals.
2 1 −5 −12 76 −79
2 −6 −36 80
2 1 −3 −18 40 1 = f (2)
2 −2 −40
2 1 −1 −20 0 = f ′(2)
2 2
1
1 1 −18 = f ′ ′(2) ⇒ f ′ ′(2) = − 36
2
−2 f ( c ) −2 f (2)
Now x=c± =2± ⋅
f ′ ′ (c ) f ′ ′ (2)
− 2 × 1
∴ x=2± = 2 ± √ (⋅05555) = 2 ± ⋅ 2357.
− 36
Problem 5: Describe the method of false position to solve an equation f ( x) = 0. Workout two steps of this
method to solve x3 + 1⋅ 2 x2 − 4 x − 4 ⋅ 8 = 0.
N-76
Putting x2 = y, we have
y ( y − 5)2 = (2 y − 6)2 or y3 − 14 y2 + 49 y − 36 = 0
or y ( y2 + 49) = 14 y2 + 36
or y2 ( y2 + 49)2 = (14 y2 + 36)2
or z 3 − 98 z 2 + 1393 z − 1296 = 0
or z ( z 2 + 1393) = 98 z 2 + 1296
∴ | p1|8 = q1 ≈ 6818
1686433
| p2|8 = q2 ≈
6818
1679616
| p3|8 = q3 ≈ ⋅
1686433
Taking square roots three times, we get
| p1| = (6818)1 /8 = 3 ⋅ 0144433
Problem 8: Use Graeffe’s method to solve the equation given below, squaring three times :
x3 − 5 x2 − 17 x + 20 = 0. (Kumaun 2010; Gorakhpur 13)
Collecting all terms containing even powers of x to one side and all terms containing odd powers of
x to the other side and squaring both sides, we get from (1)
x2 ( x2 − 17)2 = 25 ( x2 − 4)2 …(2)
y ( y − 17)2 = 25 ( y − 4)2
− 25600000000 = 0
3 2
or t + (383842 − 6265009) t + (36833670000 − 800960000 ) t
− 25600000000 = 0
3 2
or t − 5881167 t + 36032710000 t − 25600000000 = 0 …(3)
Let the roots of (3) be q1 , q2 , q3 .
Then q1 = p18 = | p1|8 , q 2 = | p2|8 , q 3 = | p3|8 .
1 b
is xn +1 = xn + ⋅ …(1)
2 xn
Take f ( x ) = x2 − 30.
N-79
1 b 1 30
and x4 = x3 + = 5 ⋅ 4772 + = 5 ⋅ 477225.
2 x3 2 5 ⋅ 4772
Problem 10: Obtain the cube root of 12 to five decimal places by Newton-Raphson method.
Solution: Let x = (12)1 /3 .
Then x3 = 12 or x3 − 12 = 0.
Therefore the cube root of 12 is nothing but the root of the equation
f ( x ) ≡ x3 − 12 = 0. …(1)
We have f ( x ) = x3 − 12 and f ′ ( x ) = 3 x2 .
Also f (2) = 23 − 12 = 8 − 12 = − 4 which is −ve
and f (3) = 33 − 12 = 27 − 12 = 15 which is +ve.
x 3 − 12 3 xn3 − xn3 + 12
= xn − n =
3 xn2 3 xn2
2 xn3 + 12 1
= = 2 xn + 12 ⋅
3 xn2 3 2
xn
N-80
1
Thus, xn +1 = 2 xn + 12 ⋅ …(2)
3 2
xn
1
x1 = 2 x0 + 12 = 1 2 × 2 ⋅ 2 + 12
3 2 3 2
(2 ⋅ 2)
x0
=
1 4 ⋅ 4 + 12 = 1 (4 ⋅ 4 + 2 ⋅ 4793388) = 2 ⋅ 2931129.
3 4 ⋅ 84 3
Taking n = 1, we have from (2), the second approximation
1
x2 = 2 x1 + 12 = 1 2 × 2 ⋅ 2931129 + 12
2 3 2
3 x1 (2 ⋅ 2931129)
=
1 4 ⋅ 5862258 + 12
3 5 ⋅ 2583668
1
= (4 ⋅ 5862258 + 2 ⋅ 2820774) = 2 ⋅ 2894344.
3
Taking n = 2 , we have from (2), the third approximation
1 12 1 12
x3 = 2 x2 + = 2 × 2 ⋅ 2894344 +
3 2
x2 3 2
(2 ⋅ 2894344)
= 3 − ⋅116 = 2 ⋅ 8838.
Now taking x1 = 2 ⋅ 8838 as approximate root and again applying Newton-Raphson formula, we
obtain x2 = 2 ⋅ 8832 .
N-81
x
Problem 13: Find the real root of the equation e − 3 x = 0, using iteration method.
x
Solution: Let f ( x) ≡ e − 3 x = 0 so that f ′ ( x ) = e x − 3.
We have f (0 ) = e 0 − 3 × 0 = 1 = + ive
As f ( x ) is +ive for x = 0 and −ive for x = 1, therefore root lies between x = 0 and x = 1.
Again f (⋅5) = e ⋅5 − 3 × ⋅ 5 = 1⋅ 6487 − 1⋅ 5 = ⋅14872 = + ive
= ⋅ 6 + ⋅ 0187 = ⋅ 6187.
Repeating the formula taking approximate root as ⋅6187, we have
f ( x1 ) 1⋅ 8565 − 1⋅ 8561
x2 = x1 − = ⋅ 6187 −
f ′( x1 ) 1⋅ 8563 − 3
⋅0004
= ⋅ 6187 + = ⋅ 61904.
1⋅1435
Thus the approximate real root of the given equation is 0 ⋅ 61904.
N-82
Squaring, we get
x2 ( x2 + 17)2 = (8 x2 + 10 )2 .
Putting x2 = y, we get
y ( y + 17)2 = (8 y + 10 )2
or y3 − 30 y2 + 129 y − 100 = 0
or y3 + 129 y = 30 y2 + 100
or y ( y2 + 129) = 30 y2 + 100.
Squaring, we get
y2 ( y2 + 129)2 = (30 y2 + 100 )2 .
Putting y2 = z , we have
10641
| p2|4 = q2 ≈
642
10000
| p3|4 = q3 ≈ ⋅
10641
Taking square roots two times, we get
| p1| = (642)1 /4 = 5 ⋅ 03366
By actually substituting in (1), we see that p1 = 5 ⋅ 03366 satisfies (1) while p1 = − 5 ⋅ 03366 does
not satisfy (1).
So, p1 = 5 ⋅ 03366.
Similarly, we see that the admissible values of p2 , p3 are p2 = 2 ⋅ 01772 , p3 = 0 ⋅ 98459.
Hence, the approximate values of the roots of the given equation are
5 ⋅ 03366, 2 ⋅ 01772 , 0 ⋅ 98459.
1
Problem 15: The equation sin x = 5 x − 2 can be put as x = sin−1 (5 x − 2)and also as x = ( sin x + 2)
5
suggesting two iterative procedures for its solution. Which of these, if any, would succeed and which would fail
to give the root in the neighbourhood of 0 ⋅ 5.
Solution: As given sin x = 5 x − 2 ;
∴ x = sin−1 (5 x − 2).
We can apply the procedure again to get another approximation and so on.
1
Now let us consider x = (sin x + 2).
5
1
In this case taking x = ⋅ 5 = on (R.H.S.), the first approximation is given by
2
x1 =
1 sin 1 + 2 ⋅
5 2
1
Now sin is not known, hence this method fails even for first approximation. Hence we would
2
succeed from first one.
1 N
xn +1 = 2 xn + 2 [ Prove this formula before using it]
3 xn
N-84
Here N = 10.
1 10
Therefore, x1 = first approximation = 2 × 2 ⋅ 5 + = 2⋅2 ;
3 (2 ⋅ 5)2
1 10
x2 = second approximation = 2 × 2 ⋅ 2 + = 2 ⋅15537 ;
3 (2 ⋅ 2)2
1 10
= 2 × 2 ⋅15537 + = 2 ⋅15444 .
2
3 (2 ⋅ 15537 )
62 − 29 7
For x0 = 6, we get x1 = 6 − =6− = 5 ⋅ 417.
2×6 12
Now f (0 ) = 0 − ⋅ 33 = − ive.
and f (1) = 1 − ⋅ 33 = + ive.
∴ Root lies between 0 and 1. Taking x0 = 1 and applying Newton-Raphson method, we have
f ( x0 ) (1)4 − ⋅ 33 ⋅67
x1 = x0 − =1− =1−
f ′ ( x0 ) 4 × (1)3 4
= ⋅ 8325.
N-85
Now taking x1 = ⋅ 8325 as approximate root and again applying Newton-Raphson formula, we
have
f ( x1 ) (⋅8325)4 − ⋅ 33
x2 = x1 − = ⋅ 8325 −
f ′ ( x1 ) 4 × (⋅8325)3
⋅1503268
= ⋅ 8325 − = 0 ⋅ 7673.
2 ⋅ 307876
Problem 18: By using Newton-Raphson method, find the root of x4 − x − 10 = 0, which is nearer to
x = 2 , correct to three places of decimals. (Rohilkhand 2002; Agra 12; Bundelkhand 12;
Kumaun 12)
Solution: Here f ( x) = x4 − x − 10.
∴ f ′( x ) = 4 x3 − 1.
f ( xn ) x 4 − xn − 10 3 xn4 + 10
∴ xn +1 = xn − = xn − n = ⋅
f ′( xn ) 4 xn3 − 1 4 xn3 − 1
3 x 24 + 10 3 (1⋅ 856)4 + 10
x3 = = = 1⋅ 856.
3
4 x2 − 1 4 (1⋅ 856)3 − 1
3. The real roots of the equation f ( x ) = 0 are the abscissae of the points where the curve
y = f ( x ) crosses the x-axis.
7. Newton’s iterative formula for obtaining a−1 is xn+1 = xn (2 − axn ). See article 2, part (v).
8. Newton’s iterative formula for obtaining the square root of a is
1 a
xn +1 = xn + ⋅
2 xn
N-86
9. The method which does not require any prior information of the roots such as approximate
value of the root etc., is Graeffe’s root squaring method.
10. Let f ( x ) = x3 − x − 1. Then f (1) = 1 − 1 − 1 = − 1 and f (2) = 8 − 2 − 1 = 5.
Since f (1) and f (2) are of opposite signs, therefore at least one root or an odd number of
roots of f ( x ) = 0 lie between 1 and 2.
11. All the three numbers1, 2 , 5 satisfy the equation x3 − 8 x2 + 17 x − 10 = 0. Hence, the exact
roots of the equation x3 − 8 x2 + 17 x − 10 = 0 are 1, 2 , 5.
12. See article 2 (v).
13. Lagrange’s method is used to find any value of a function when the given values of the
independent variable are not equidistant, not for finding roots.
14. See article 2 (ii).
True or False
1. If f ( x ) is exactly divisible by x − a, then a is a root of the equation f ( x ) = 0.
Chapter-6
Numerical Differentiation
Comprehensive Problems 1
dy d2 y
Problem 1: Find and at x = 1 from the following table :
dx dx2
x : 1 2 3 4 5 6
Solution: Here the values of the argument are equally spaced and the derivatives are required at
x = 1, which is near the beginning of the table so we shall use Newton’s forward formula.
The forward difference table is given below :
x y = f ( x) ∆ f ( x) ∆2 f ( x ) ∆3 f ( x ) ∆4 f ( x ) ∆5 f ( x )
1 198669
96851
2 295520 − 2953
93898 − 938
3 389418 − 3891 39
90007 − 899 8
4 479425 − 4790 47
85217 − 852
5 564642 − 5642
79575
6 644217
x ( x − 1) ( x − 2) 3 x ( x − 1) ( x − 2) ( x − 3) 4
+ ∆ f ( a) + ∆ f ( a)
3! 4!
x ( x − 1) ( x − 2) ( x − 3) ( x − 4) 5
+ ∆ f ( a)
5!
N-88
( x2 − x ) 2 ( x3 − 3 x2 + 2 x ) 3
= f ( a) + x ∆ f ( a) + ∆ f ( a) + ∆ f ( a)
2 6
( x4 − 6 x3 + 11x2 − 6 x ) 4
+ ∆ f ( a)
24
( x5 − 10 x4 − 35 x3 − 50 x2 + 24 x ) 5
+ ∆ f ( a).
120
Differentiating with respect to x twice, we get
(2 x − 1) 2 (3 x2 − 6 x + 2) 3
h f ′ ( a + xh) = ∆ f ( a) + ∆ f ( a) + ∆ f ( a)
2 6
(4 x3 − 18 x2 + 22 x − 6) 4
+ ∆ f ( a)
24
(12 x2 − 36 x + 22) 4
and h2 f ′ ′ ( a + xh) = ∆2 f ( a) + ( x − 1) ∆3 f ( a) + ∆ f ( a)
24
11 5
and f ′ ′ (1) = − 2953 + ( −1) × ( − 938) + × 39 − × 8
12 6
dy d2 y
Hence, at x = 1, = 98007 and = − 1986.
dx dx2
x : 2 4 9 13 16 21 29
Solution: In this case the values of the argument are unequally spaced, so we may use Newton’s
divided difference formula.
x f ( x) ∆| f ( x ) ∆| 2 f ( x ) ∆| 3 f ( x ) ∆| 4 f ( x ) ∆| 5 f ( x )
2 57
644
4 1345 1765
12999 556
9 66340 7881 45
83928 1186 1
13 402052 22113 64
238719 2274 1
16 1118209 49401 89
633927 4054
21 4287844 114265
2119372
29 2124282
Since the fifth divided differences are constant, Newton’s divided difference formula is
f ( x ) = f ( a) + ( x − a) ∆| f ( a) + ( x − a) ( x − b ) ∆| 2 f ( a)
+ ( x − a)( x − b ) ( x − c ) ∆| 3 f ( a) + ( x − a) ( x − b ) ( x − c ) ( x − d ) ∆| 4 f ( a)
+ ( x − a) ( x − b ) ( x − c ) ( x − d ) ( x − e ) ∆| 5 f ( a)
[Upto fifth differences]
Differentiating with respect to x thrice, we have
f ′ ′ ′ ( x ) = 6 ∆| 3 f ( a) + [24 x − 6 ( a + b + c + d )] ∆| 4 f ( a)
+ [60 x2 − 24 x ( a + b + c + d ) + 6 ( ab + cd ) + 6 ( a + b ) ( c + d )
Problem 3: Find the first three derivatives of the function tabulated below at the point x = 2 ⋅ 5 :
Solution: In this case the values of the argument are unequally spaced, so we may use here
Newton’s divided difference formula.
N-90
x f ( x) ∆| f ( x ) ∆| 2 f ( x ) ∆| 3 f ( x )
1⋅ 5 3 ⋅ 375
6 ⋅ 71
1⋅ 9 6 ⋅ 059 5 ⋅ 90
12 ⋅ 61 1
2⋅5 13 ⋅ 625 7 ⋅ 60
22 ⋅ 49 1
3⋅2 29 ⋅ 368 10 ⋅ 00
40 ⋅ 49 1
4⋅3 73 ⋅ 907 13 ⋅ 40
76 ⋅ 67
+ ( x − a) ( x − b ) ( x − c ) ∆| 3 f ( a)
+ [3 x2 − 2 x ( a + b + c ) + ab + ac + bc ] ∆| 3 f ( a),
f ′ ′ ( x ) = 2 ∆|2 f ( a) + [6 x − 2 ( a + b + c )] ∆| 3 f ( a),
f ′ ′ ′ ( x ) = 6 ∆| 3 f ( a).
Here a = 1⋅ 5, b = 1⋅ 9, c = 2 ⋅ 5.
On putting x = 2 ⋅ 5, we get
f ′ (2 ⋅ 5) = 6 ⋅ 710 + [5 − 3 ⋅ 4] × 5 ⋅ 9
+ [3 × 6 ⋅ 25 − 5 (1⋅ 5 + 1⋅ 9 + 2 ⋅ 5) + 11⋅ 35] × 1
= 6 ⋅ 710 + 9 ⋅ 44 + ⋅ 600 = 16 ⋅ 750,
f ′ ′ (2 ⋅ 5) = 2 [5 ⋅ 9] + [15 ⋅ 00 − 2 × 5 ⋅ 9] × 1 = 15 ⋅ 00,
f ′ ′ ′ (2 ⋅ 5) = 6 × 1 = 6 ⋅ 00.
N-91
x : 3 5 11 27 34
(Purvanchal 2008)
Solution: Here the values of the argument x are not equally spaced, therefore we shall use
Newton’s divided difference formula. The divided difference table is given below :
x f ( x) ∆| f ( x ) ∆| 2 f ( x ) ∆| 3 f ( x ) ∆| 4 f ( x )
3 − 13
18
5 23 16
146 1
11 899 40 0
1026 1
27 17315 69
2613
34 35606
+ ( x − x0 ) ( x − x1 ) ( x − x2 ) ∆| 3 f ( x0 )
+ ( x − 3) ( x − 5) ( x − 11) ∆| 3 f (3)
= − 13 + ( x − 3) × 18 + ( x − 3) ( x − 5) × 16 + ( x − 3) ( x − 5) ( x − 11) × 1
∴ f ′ ( x ) = 18 + 16 (2 x − 8) + (3 x2 − 38 x + 103).
dy
Problem 5: Find at x = 1 from the following table :
dx
dy d2 y
Problem 7: Find and of y = x1 /3 at x = 50 from the following table :
dx dx2
x: 50 51 52 53 54 55 56
1 /3
y= x : 3⋅6840 3⋅7084 3⋅7325 3⋅7563 3⋅7798 3⋅8030 3⋅8259
Solution: Since the derivatives are required at x = 50, which is near the beginning of the table,
therefore we shall use Newton’s forward formula. The difference table is given below :
x y = x1 /3 ∆ y ∆2 y ∆3 y
50 3⋅6840
0⋅0244
51 3⋅7084 – 0⋅0003
0⋅0241 0
52 3⋅7325 – 0⋅0003
0⋅0238 0
53 3⋅7563 – 0⋅0003
0⋅0235 0
54 3⋅7798 – 0⋅0003
0⋅0232 0
55 3⋅8030 – 0⋅0003
0⋅0229
56 3⋅8259
N-93
Problem 8: Find the first two derivatives of f ( x) at x = 1 from the following table :
x : −2 −1 0 1 2 3 4
f ( x) : 104 17 0 −1 8 69 272
Problem 9: Using divided differences, find the value of f ′ (8), given that
f (6) = 1⋅ 556, f (7) = 1⋅ 690, f (9) = 1⋅ 908, f (12) = 2 ⋅158 .
Solution: Divided difference formula is given by
f ( x ) = f ( x 0 ) + ( x − x 0 ) ∆| f ( x 0 ) + ( x − x 0 ) ( x − x1 ) ∆| 2 f ( x 0 )
+ ( x − x 0 ) ( x − x1 ) ( x − x 2 ) ∆| 3 f ( x 0 )
+ ( x − x 0 ) ( x − x1 ) ( x − x 2 ) ( x − x 3 ) ∆| 4 f ( x 0 )
+ ( x − x 0 ) ( x − x1 ) ( x − x 2 ) ( x − x 3 ) ( x − x 4 ) ∆| 5 f ( x 0 ) + … …(1)
x f ( x) ∆| f ( x ) ∆| 2 f ( x ) ∆| 3 f ( x )
6 1⋅ 556
⋅134
7 1⋅ 690 − ⋅ 0083
⋅109 ⋅00051
9 1⋅ 908 − ⋅ 0052
⋅083
12 2 ⋅158
+ x 0 x1 + x1 x 2 + x 2 x 0 ] ∆| 3 f ( x 0 ).
N-94
Putting x 0 = 6, x1 = 7, x 2 = 9, x = 8 and the values of the divided differences from the table, we
get
f ′ (8) = ⋅134 + [2 × 8 − 6 − 7] ( − ⋅ 0083)
+ [3 × 64 − 2 × 8 (6 + 7 + 9) + 6 × 7 + 7 × 9 + 9 × 6] (⋅00051)
= ⋅134 − ⋅ 0249 + (192 − 352 + 159) (⋅00051)
= ⋅134 − ⋅ 0249 − ⋅ 00051
= ⋅134 − ⋅ 0254 = 0 ⋅10859.
x: 0 1 3 4 5 7 9
Problem 11: Show that the expressions given below are approximations to the third derivative of y x :
∆3 y0 + x − ∆4 y0
3
(i)
2
Problem 12: Use Stirling’s formula to find the first derivatives of the function y = 2e x − x − 1
tabulated below at the point x − 0 .6.
x : 0.4 0.5 0.6 0.7 0.8
y : 1.5836 1.7974 2.0442 2.3275 2.6510
N-95
Solution: Here we require the derivative at the point x = 0 .6 which lies near the middle of the
table, so we use Stirling’s formula.
Take x = 0 .6, x0 = 0 .6, h = 1
x − x0 0 .6 − 0 .6
∴ u= = =0
h 1
The Stirling’s formula is given by
u4 − u2 4
+ ∆ y− 2 + … …(1)
24
The difference table is given below :
x u yu ∆ yu ∆2 yu ∆3 yu ∆4 yu
0.4 –2 1.5836
0.2138
0.2468 0.0035
0.2833 0.0037
0.3235
0.8 2 2.6510
1 1 1 4
and y0 ′ = ( y1 − y−1 ) − (δ2 y1 − δ2 y−1 ) + (δ y1 − δ4 y−1 ) − … ⋅
2h 6 30
(Bundelkhand 2012)
2
δ
Solution: We know that µ 2 = 1 + ⋅
4
−1 / 2
δ2 δ2 3 4
∴ µ −1 = 1 + =1− + δ −…
4 8 128
hD hD −1 δ2 3 4 δ2 3 4
Now = µ = 1 − + δ − … 1 − 8 + 128 δ − …
δµ δ 24 640
δ2 1 4
=1− + δ −…
6 30
δµ δ2 1 4 µ δ3 δ5
or D= 1− + δ − … = δ − + − … ⋅
h 6 30 h 6 30
Operating these operators on y, we get
µ δ3 δ5
Dy = δ − 6 + 30 − … y
h
µ δ3 y δ5 y
⇒ y′= δy − + − …
h 6 30 …(1)
By definition, we have
1
δ = E 1 /2 − E −1 /2 and µ = ( E 1 /2 + E −1 /2 ).
2
1
∴ δµ = ( E − E −1 ) = µδ.
2
1 1
∴ (δµ ) y0 = ( E − E −1 ) y0 = ( y1 − y−1 )
2 2
…(2)
= µδy0
1 µ δ3 µ δ5
y0 ′ = µ δy0 − 6 y0 + 30 y0 − …
h
1 ( y − y ) − 1 (δ2 y − δ2 y ) + 1 (δ4 y − δ4 y ) − … ⋅
= −1 −1 −1
2h 1 6
1
30
1
N-97
True or False
1. See article 1.
2. See article 1.
3. In case the values of the argument are unequally spaced, and we desire to find the derivative
of the function at a point, we should use Newton’s divided difference formula.
4. See Ex. 9.
❍❍❍
N-98
Chapter-7
Numerical Integration
Comprehensive Problems 1
3
Problem 1: Calculate the approximate value of ∫ x4 dx by using
−3
Solution: Divide the range of integration ( − 3, 3) into six equal parts each of width
3 − ( − 3)
= = 1. Hence h = 1. The values of the function for each point of sub-division are given
6
below :
x y = x4
x0 = − 3 ( − 3)4 = 81
x0 + h = − 2 ( − 2)4 = 16
x 0 + 2h = − 1 ( − 1)4 = 1
x 0 + 3h = 0 04 = 0
x 0 + 4h = 1 14 = 1
x 0 + 5h = 2 24 = 16
x 0 + 6h = 3 34 = 81
Problem 2(i): Use Simpson’s rule dividing the range into ten equal parts, to show that
1 log (1 + x2 )
∫0 2
1+ x
dx = 0 ⋅173.
Solution: Dividing the range of integration (0, 1) into 10 equal parts by taking h = 0 ⋅1, we
log (1 + x2 )
compute the values of the function y = at each point of subdivision.
1 + x2
log (1 + x2 )
x x2 1 + x2 log (1 + x2 ) y=
1 + x2
0 0 1 0 0 = y0
1
By Simpson’s ‘ ’ rule, we have
3
1 log ( 1 + x2 ) h
∫0 1+ x 2
dx = [ y0 + y10 + 4 ( y1 + y3 + y5 + y7 + y9 )
3
+ 2 ( y2 + y4 + y6 + y8 )]
⋅1
= [0 + ⋅ 3465735 + 4 (⋅0098518 + ⋅ 0790621 + ⋅1785148
3
+ ⋅ 2676349 + ⋅ 3278048 + 2 (⋅0377122 + ⋅1279482
+ ⋅ 2260916 + ⋅ 301644)]
⋅1
= [⋅3465735 + 3 ⋅ 4514736 + 1⋅ 386792]
3
⋅1
= (5 ⋅184839) = ⋅17282793
3
= ⋅1728 correct to four decimal places
= ⋅173 correct to three decimal places.
N-100
6 dx
Problem 2(ii): Using Simpson’s one-third rule, find ∫0 (1 + x )2
⋅
(Bundelkhand 2011)
Solution: Dividing the range of integration (0, 6) into 6 equal parts by taking h = 1, we compute
the values of the function y = 1 / (1 + x )2 at each point of subdivision.
x 0 1 2 3 4 5 6
(1 + x )2 1 4 9 16 25 36 49
1
By Simpson’s ‘ ’ rule, we have
3
6 dx h
∫0 (1 + x )2
= [ y0 + y6 + 4 ( y1 + y3 + y5 ) + 2 ( y2 + y4 )]
3
1
= [1 + ⋅ 02041 + 4 (⋅25 + ⋅0625 + ⋅02778) + 2 (⋅11111 + ⋅ 04)]
3
1
= [1⋅ 02041 + 1⋅ 36112 + ⋅ 30222] = 0 ⋅ 8945833
3
= 0 ⋅ 8946 correct to four decimal places.
1 x2 1
Problem 2(iii): Find the value of log 2 from ∫ dx, using Simpson’s ‘ ’ rule, by dividing the
0 1 + x3 3
range of integration into four equal parts. Also find the error.
Solution: Dividing the range of integration (0, 1) into four equal parts by taking h = ⋅ 25, we
compute the values of the function y = x2 / (1 + x3 ) at each point of sub-division.
1
By Simpson’s ‘ ’ rule, we have
3
1 x2 h
∫0 1+ x 3
dx =
3
[ y0 + y4 + 4 ( y1 + y3 ) + 2 y2 ]
⋅25
= [0 + ⋅ 5 + 4 (⋅0615384 + ⋅ 3956044) + 2 (⋅2222222)]
3
⋅25
= [⋅5 + 1⋅ 8285712 + ⋅ 4444444]
3
⋅25
= (2 ⋅ 7730156) = 0 ⋅ 2310846
3
= 0 ⋅ 23108 correct to five decimal places.
N-101
1
= (0 ⋅ 6931471) = 0 ⋅ 23105 correct to five decimal places.
3
Therefore the error = 0 ⋅ 23105 − 0 ⋅ 23108 = − ⋅ 00003.
1
Remark: We have log 2 = 0 ⋅ 2310846 ⇒ log 2 = 0 ⋅ 6932538.
3
Also, from tables, log e 2 = 0 ⋅ 6931471.
Therefore, the error = 0 ⋅ 6932538 − 0 ⋅ 6931471 = 0 ⋅ 0001067.
12 dx
Problem 3(i): Evaluate ∫ by numerical methods.
5 x (Rohilkhand 2011)
1
Solution: Let y = f ( x ) = ⋅
x
Dividing the range of integration (5, 12) into seven equal parts by taking h = 1, we compute the
values of the function y = 1 / x at each point of sub-division.
x: 5 6 7 8 9 10 11 12
12 1 h
∫5 x
dx =
2
[( y0 + y7 ) + 2 ( y1 + y2 + y3 + y4 + y5 + y6 )]
1
= [⋅2 + ⋅ 083333 + 2 (⋅166667 + ⋅142857 + ⋅125
2
+ ⋅111111 + ⋅1 + ⋅ 090909) ]
1
= (⋅283333 + 1⋅ 473088) = ⋅ 87821.
2
Also by actual integration
12 1 12
∫5 x
dx = [log x ]12
5
= log 12 − log 5 = log
5
= log 2 ⋅ 4
= ⋅ 87547.
∴ The error = ⋅ 87547 − ⋅ 87821 = − ⋅ 00274.
21
Problem 3(ii): Evaluate ∫ dx by Simpson’s rule with 4 strips (i. e., 5 ordinates)
1 x
and 8 strips (i. e., 9 ordinates) respectively. Determine the error by direct integration.
Solution: Let y = f ( x ) = 1 / x.
N-102
Case I: The range of integration (1, 2)is divided into four equal parts. In this case, h = 1 / 4 = ⋅ 25.
x 1 5 /4 3 /2 7 /4 2
y 1 ⋅8 ⋅666667 ⋅571429 ⋅5
y0 y1 y2 y3 y4
1
In this case by Simpson’s ‘ ’ rule,
3
21 h
∫
1 x
dx = [( y0 + y4 ) + 4 ( y1 + y3 ) + 2 y2 ]
3
1
= [1 + ⋅ 5 + 4 (⋅8 + ⋅ 571429) + 2 (⋅666667)]
12
1
= (1⋅ 5 + 5 ⋅ 485716 + 1⋅ 333334)
12
1
= (8 ⋅ 31905) = 0 ⋅ 693254. …(1)
12
Case II: The range of integration (1, 2) is divided into eight equal parts. In this case,
h = 1 / 8 = ⋅125.
When x = 1, y =1= y ;
0
when x = 1⋅125, y = ⋅ 888889 = y1 ;
when x = 1⋅ 25, y = ⋅ 8 = y2 ;
when x = 1⋅ 375, y = ⋅ 727273 = y3 ;
when x = 1⋅ 5, y = ⋅ 666667 = y4 ;
when x = 1⋅ 625, y = ⋅ 615385 = y5 ;
when x = 1⋅ 75, y = ⋅ 571429 = y6 ;
when x = 1⋅ 875, y = ⋅ 533333 = y7 ;
and when x = 2, y = ⋅ 5 = y8 .
1
By Simpson’s ‘ ’ rule, we have
3
21 h
∫1 x
dx = [ y0 + y8 + 4 ( y1 + y3 + y5 + y7 ) + 2 ( y2 + y4 + y6 )]
3
1/8
= [1 + ⋅ 5 + 4 (⋅888889 + ⋅ 727273 + ⋅ 615385 + ⋅ 533333)
3
+ 2 (⋅8 + ⋅ 666667 + ⋅5571429)]
1
= (1⋅ 5 + 11⋅ 05952 + 4 ⋅ 076192)
24
= 0 ⋅ 693155 …(2)
2 1
By actual integration
∫1 x
dx = [log x ]12 = log 2 = 0 ⋅ 693147.
7 dx
Problem 3(iii): Use Simpson’s rule to prove that loge 7 is approximately 1⋅ 9587 using ∫ ⋅
1 x
Solution: Divide the range (1, 7) into six equal parts each of width 1. Hence h = 1. The values
of y for each point of sub-division are given below :
x: x0 = 1 x1 = 2 x2 = 3 x3 = 4 x4 = 5 x5 = 6 x6 = 7
y: y0 = 1 y1 = 1 / 2 y2 = 1 / 3 y3 = 1 / 4 y4 = 1 / 5 y5 = 1 / 6 y6 = 1 / 7
1 dx
Problem 4(i): Show that ∫ = log 2 = 0 ⋅ 69315.
0 1+ x (Lucknow 2008; Purvanchal 09)
Solution: To use Simpson’s ‘1/3’ rule, divide the range of integration (0, 1) into ten equal
parts each of width 1/10. Hence h = 1 / 10 = ⋅1, n = 10.
The values of y at each point of sub-division are given below :
1
x y=
1+ x
x0 = 0 1
x 0 + h = ⋅1 1 / 1⋅1 = ⋅ 9090909
x 0 + 2h = ⋅ 2 1 / 1⋅ 2 = ⋅ 8333333
x 0 + 3h = ⋅ 3 1 / 1⋅ 3 = ⋅ 7692307
x 0 + 4h = ⋅ 4 1 / 1⋅ 4 = ⋅ 7142857
x 0 + 5h = ⋅ 5 1 / 1⋅ 5 = ⋅ 6666666
x 0 + 6h = ⋅ 6 1 / 1⋅ 6 = ⋅ 6250000
x 0 + 7h = ⋅ 7 1 / 1⋅ 7 = ⋅ 5882352
x 0 + 8h = ⋅ 8 1 / 1⋅ 8 = ⋅ 5555555
x 0 + 9h = ⋅ 9 1 / 1⋅ 9 = ⋅ 5263157
x 0 + 10 h = 1 1 / 2 = ⋅ 5000000
1 dx h
∫0 1 + x = 3 [ y0 + y10 + 4 ( y1 + y3 + y5 + y7 + y9 )
+ 2 ( y2 + y4 + y6 + y8 )]
⋅1
= [1⋅ 5 + 4 × 3 ⋅ 4595391 + 2 × 2 ⋅ 7281745]
3
⋅1
= × 20 ⋅ 794505 = ⋅ 6931501
3
= ⋅ 69315.
1 dx
The actual value of
∫0 1+ x
= [log (1 + x )]1
0
10 dx
Problem 4(ii): Calculate (upto 4 places of decimal ) ∫ by dividing the range into eight equal
2 1+ x
parts.
Solution: Divide the range of integration (2, 10 ) into eight equal parts each of width
10 − 2
= 1. Hence h = 1. The values of y for each point of sub-division are given below :
8
1
x y=
1+ x
x0 = 2 1/3
x0 + h = 3 1/4
x 0 + 2h = 4 1/5
x 0 + 3h = 5 1/6
x 0 + 4h = 6 1/7
x 0 + 5h = 7 1/8
x 0 + 6h = 8 1/9
x 0 + 7h = 9 1 / 10
x 0 + 8h = 10 1 / 11
1 1 1 1 1
+ 4 + + +
1
+ 2 + +
1 1 1 1
= +
3 3 11 4 6 8 10 5 7 9
1
= [⋅4242424 + 4 (⋅6416666) + 2 (⋅4539682)]
3
1
= × 3 ⋅ 8988453 = 1⋅ 2996151 = 1⋅ 2996.
3
0⋅7 1 /2
Problem 5(i): Evaluate ∫ x e − x dx approximately by using a suitable formula.
0⋅5
Solution: Divide the range of integration (⋅5, ⋅ 7) into 4 equal parts each of width
⋅7 − ⋅ 5
= ⋅ 05. Hence h = ⋅ 05. The values of y for each point of sub-division are given below :
4
x y = x1 /2 e − x
x 0 = ⋅ 50 ⋅4288818
x 0 + h = ⋅ 55 ⋅4278774
x 0 + 2h = ⋅ 60 ⋅4251076
x 0 + 3h = ⋅ 65 ⋅4208867
x 0 + 4h = ⋅ 70 ⋅4154730
⋅05
= [⋅8443548 + 3 ⋅ 3950564 + ⋅ 8502152]
3
= ⋅ 0848271.
4
Problem 5(ii): Evaluate ∫ e x
dx, by Simpson’s rule, given that
0
= 54 ⋅ 60 − 1 = 53 ⋅ 60.
π /2 sin x 1
Problem 6(i): Calculate ∫ e dx correct to four decimal places by Simpson’s ‘ ’ rule, dividing
0 3
the range of integration (0, π / 2) into four equal parts.
π
Solution: Dividing the range of integration (0, π / 2) into four equal parts by taking h = , we
8
compute the values of the function y = esin x at each point of sub-division.
x sin x y = esin x
x0 = 0 0 y0 = 1
x0 + h = π / 8 0.3826834 y1 = 1⋅ 4662138
x 0 + 2h = π / 4 0.7071067 y2 = 2 ⋅ 028115
x 0 + 3h = 3 π / 8 0.9238795 y3 = 2 ⋅ 5190442
x 0 + 4h = π / 2 1 y4 = 2 ⋅ 7182818
1
By Simpson’s ‘ ’ rule, we get
3
π /2 sin x h
∫0 e dx =
3
[ y0 + y4 + 4 ( y1 + y3 ) + 2 y2 ]
π /8
= [1 + 2 ⋅ 7182818 + 4 (1⋅ 4662138 + 2 ⋅ 5190442)
3
+ 2 (2 ⋅ 028115)]
π
= [23 ⋅ 7155438] = 3 ⋅1043574
24
= 3 ⋅1044, correct to four decimal places.
π /2 3
Problem 6(ii): Calculate ∫ e sin x dx correct to four decimal places by Simpson’s ‘ ’ rule, dividing
0 8
the range of integration (0, π / 2) into 3 equal parts.
Solution: Dividing the range of integration (0, π / 2) into 3 equal parts by taking h = π / 6, we
compute the values of the function y = esin x at each point of sub-division.
x 0 π /6 π /3 π /2
sin x 0 0 ⋅5 0 ⋅ 8660254 1
3
By Simpson’s ‘ ’ rule we have
8
π / 2 sin x 3h
∫0 e dx =
8
[ y0 + y3 + 3 ( y1 + y2 )]
3 π
= ⋅ [1 + 2 ⋅ 7182818 + 3 (1⋅ 6487213 + 2 ⋅ 377427)]
8 6
π
= [3 ⋅ 7182818 + 12 ⋅ 078445]
16
π
= [15 ⋅ 796727] = 3 ⋅10168
16
= 3 ⋅1017 correct to four decimal places.
1⋅6 x
Problem 7: Find by Weddle’s rule the value of the integral I = ∫ dx by taking 12
0⋅4 sinh x
sub-intervals.
1⋅ 6 − 0 ⋅ 4
Solution: Here h = = 0 ⋅1.
12
x
Let y= ⋅ We make the following calculations :
sinh x
x sinh x y
0 ⋅4 ⋅4107523 ⋅9738229 = y0
0 ⋅5 ⋅5210953 ⋅9595173 = y1
0 ⋅6 ⋅6366535 ⋅9424278 = y2
0 ⋅7 ⋅7585837 ⋅9227722 = y3
0 ⋅8 ⋅8881059 ⋅9007934 = y4
0 ⋅9 1⋅ 0265167 ⋅8767514 = y5
1⋅ 0 1⋅1752012 ⋅8509181 = y6
1 ⋅1 1⋅ 3356475 ⋅8235705 = y7
1⋅ 2 1⋅ 5094614 ⋅7949855 = y8
1⋅ 3 1⋅ 6983824 ⋅7654342 = y9
By Weddle’s rule
1⋅6 x 3h
∫0⋅4 sinh x
dx =
10
[ y0 + y12 + 5 ( y1 + y5 + y7 + y11
+ y2 + y4 + y8 + y10 + 2 y6 + 6 ( y3 + y9 )]
3 (⋅1)
= [1⋅ 647346 + 16 ⋅ 821514 + 3 ⋅ 3733844
10
+ 1⋅ 7018362 + 10 ⋅129238]
3
= (33 ⋅ 673319) = 1⋅ 0101996.
100
π /2
Problem 8: Calculate an approximate value of ∫ sin x dx by the
0
(i) Trapezoidal rule,
1
(ii) Simpson’s ‘ ’ rule, using 11 ordinates. (Purvanchal 2011)
3
Solution: First we divide the range of integration into ten equal parts by taking the interval of
differencing h = π / 20 and then we compute the values of the function f ( x ) = sin x for each point
of sub-division. These computed values are as shown in the following table.
x y = sin x x y = sin x
x0 = 0 ⋅00000 x 0 + 6h = 6 π / 20 ⋅80902
x 0 + h = π / 20 ⋅15643 x 0 + 7h = 7 π / 20 ⋅89101
x 0 + 2h = 2 π / 20 ⋅30902 x 0 + 8h = 8 π / 20 ⋅95106
x 0 + 3h = 3 π / 20 ⋅45399 x 0 + 9h = 9 π / 20 ⋅98769
x 0 + 4h = 4 π / 20 ⋅58779 x 0 + 10 h = 10 π / 20 1⋅ 00000
x0 + 5h = 5 π / 20 ⋅70711
6 dx
Problem 9: Evaluate ∫ by using
0 1 + x2
Solution: Divide the range of integration (0, 6) into six equal parts each of width
6−0
= 1.
6
Hence h = 1. The values of y for each point of sub-division are given below :
1
x y=
1 + x2
x0 = 0 1 / 1 = 1⋅ 0000000
x0 + h = 1 1 / 2 = ⋅ 5000000
x 0 + 2h = 2 1 / 5 = ⋅ 2000000
x 0 + 3h = 3 1 / 10 = ⋅1000000
x 0 + 4h = 4 1 / 17 = ⋅ 0588235
x 0 + 5h = 5 1 / 26 = ⋅ 0384615
x 0 + 6h = 6 1 / 37 = ⋅ 0270270
1
= [1⋅ 0270270 + 2 ⋅ 5538462 + ⋅ 517647] = 1⋅ 3661734.
3
By Simpson’s ‘3/8’ rule, we get
6 dx 3h
∫0 1 + x2
=
8
[ y0 + y6 + 3 ( y1 + y2 + y4 + y5 ) + 2 y3 ]
3
= [1⋅ 0270270 + 2 ⋅ 391855 + ⋅ 200000 ]
8
3
= × 3 ⋅ 618882 = 1⋅ 3570808.
8
N-110
3
= [1⋅ 0270270 + 2 ⋅ 6923075 + ⋅ 8588235]
10
= 1⋅ 3734474.
⋅
11
Problem 10: Using the data of the following table compute the integral ∫ xy dx, by Simpson’s rule :
0⋅5
x : 0 ⋅5 0 ⋅6 0 ⋅7 0 ⋅8 0 ⋅9 1⋅ 0 1 ⋅1
1
By Simpson’s ‘ ’ rule,
3
⋅
11 h
∫0⋅5 xy dx = 3 [⋅2402 + 1⋅ 02091 + 4 (⋅34014 + ⋅ 58096 + ⋅ 8658)
+ 2 (⋅4543 + ⋅ 71865)]
⋅1
= [1⋅ 26111 + 7 ⋅1476 + 2 ⋅ 3459] = ⋅ 358487.
3
Problem 11(i): Find, from the following table, the area bounded by the curve and the x-axis from
x = 7 ⋅ 47 to x = 7 ⋅ 52.
x 7 ⋅ 47 7 ⋅ 48 7 ⋅ 49 7 ⋅ 50 7 ⋅ 51 7 ⋅ 52
f ( x) 1⋅ 93 1⋅ 95 1⋅ 98 2 ⋅ 01 2 ⋅ 03 2 ⋅ 06
7⋅52
Solution: The required area = ∫ f ( x ) dx. …(1)
7⋅74
Problem 11(ii): A rocket is launched from the ground, its acceleration is registered during the first 80
seconds and is given in the table below :
t ( sec ) : 0 10 20 30 40 50 60 70 80
acc (m/sec ) 2
30 ⋅ 00 31⋅ 63 33 ⋅ 44 35 ⋅ 47 37 ⋅ 75 40 ⋅ 33 43 ⋅ 25 46 ⋅ 69 50 ⋅ 67
Here h = 10.
1
By Simpson’s ‘ ’ rule, we have
3
80 h
0 ∫ f dt = [30 ⋅ 00 + 50 ⋅ 67 + 4 (31⋅ 63 + 35 ⋅ 47 + 40 ⋅ 33 + 46 ⋅ 69)
3
+ 2 (33 ⋅ 44 + 37 ⋅ 75 + 43 ⋅ 25)]
10
= [80 ⋅ 67 + 616 ⋅ 48 + 228 ⋅ 88] = 3086 ⋅ 77.
3
Hence velocity at time t = 80 is 3086 ⋅ 77 m / sec.
Problem 12: Using the following table of values, approximate by Simpson’s rule the arc length of the
graph y = 1 / x between the points (1, 1) and (5, 1 / 5).
x : 1 2 3 4 5
1 + x 4
: 1⋅414 1⋅031 1⋅007 1⋅002 1⋅001
x4
∴ The arc length of the curve between the points (1, 1) and (5, 1 / 5)
5 1 + x4
=
∫1
x4
dx
h
= [1⋅ 414 + 1⋅ 001 + 4 (1⋅ 031 + 1⋅ 002) + 2 (1⋅ 007)],
3
1
by Simpson’s ‘ ’ rule, where h = 1
3
1
= (2 ⋅ 415 + 8 ⋅132 + 2 ⋅ 014) = 4 ⋅187.
3
Problem 13: Prove Simpson’s formula
b b−a
∫a f ( x ) dx =
6n
[ f ( x 0 ) + 4 f ( x1 ) + 2 f ( x2 ) + ... + f ( x2 n )],
loge 2 = 0 ⋅ 69315.
N-112
+ 4 { f ( x 0 + h) + f ( x 0 + 3h) + …}
+ 2 { f ( x 0 + 2h) + f ( x 0 + 4h) + …}] …(1)
b−a
Putting x 0 = a, n = 2n, x 0 + 2nh = x2 n = b, h = in (1), we get the required form of
2n
Simpson’s rule.
When n = 1, we should divide the whole range (1, 2) into 2 equal parts by three points
x 0 , x1, x 2 . Thus the above formula becomes
2 dx x2 1
∫1 x
=
∫ x0 f ( x ) dx where f ( x ) =
x
, x2 = 2, x0 = 1
b−a 2 − 1 1 1 1
= [ f ( x0 ) + 4 f ( x1 ) + f ( x2 )] = + 4⋅ +
6n 6 × 1 1 3 / 2 2
= ⋅ 69444.
When n = 2, divide the whole range (1, 2) into 4 equal parts by the points x 0 , x1, x 2 , x 3 , x 4 .
2 dx x4 b−a
We have
∫1 x
=
∫ x0 f ( x ) dx =
6n
[ f ( x 0 ) + 4 f ( x1 ) + 2 f ( x 2 ) + 4 f ( x 3 ) + f ( x 4 )]
1 1 1 1 1 1
= + 4⋅ + 2⋅ + 4⋅ + = ⋅ 69325.
12 1 5 /4 3 /2 7 / 4 2
2 dx 2
The actual value of
∫1 x
= [log e x ] = log e 2 = ⋅ 69315.
1
2⋅5 1
Problem 14: Prove that ∫1⋅5 U x dx =
24
[U1 + 22U2 + U3 ] ⋅
1 1
(ii)
∫0 f ( x ) dx =
12
[5 f (1) + 8 f (0 ) − f ( − 1)] ⋅
N-113
1
Solution: (i) Here we are to express ∫ f ( x ) dx in terms of f (0 ), f (1) and f (2).
0
( x − 1) ( x − 2) ( x − 0 ) ( x − 2) ( x − 0 ) ( x − 1)
f ( x) = f (0 ) + f (1) + f (2)
(0 − 1) (0 − 2) (1 − 0 ) (1 − 2) (2 − 0 ) (2 − 1)
1 2 1
= ( x − 3 x + 2) f (0 ) − ( x2 − 2 x ) f (1) + ( x2 − x ) f (2).
2 2
1 1
1 1 x3 3 2 x3
∴
∫0 f ( x ) dx =
2
f (0 )
3
− x + 2 x − f (1)
2 0 3
− x2
0
1
1 x3 x2
+ f (2) −
2 3 2
0
5 2 1
= f (0 ) + f (1) − f (2)
12 3 12
1
= [5 f (0 ) + 8 f (1) − f (2)].
12
1
(ii) Here we are to express
∫0 f ( x ) dx in terms of f (1), f (0 ) and f ( − 1).
1 1 2
+
2
f ( − 1)
∫0
( x − x ) dx
1 1 1
1 x3 x2 x3 1 x3 x2
= f (1) + − f (0 ) − x + f ( − 1) −
2 3 2 3 0 2 3 2
0 0
5 2 1
= f (1) + f (0 ) − f ( − 1)
12 3 12
1
= [5 f (1) + 8 f (0 ) − f ( − 1)] ⋅
12
f ⋅ Since we are known four pairs of values of x and f ( x ), therefore f ( x ) can be represented as
5
2
N-114
a polynomial in x of degree 3. Consequently third differences of f ( x ) are constant and fourth and
higher differences of f ( x ) are all zero.
1
We have x0 = − and h = 1.
2
f ( x ) = f − + x +
1 1
Now [Form f ( x0 + uh)]
2 2
= E x + (1 /2) f − = (1 + ∆ ) x + (1 /2) f −
1 1
2 2
x + 1 x − 1
1 1 1 2 2 ∆2 f − 1
= f − + x + ∆f − +
2 2 2 1⋅ 2 2
x + 1 x − 1 x − 3
2 2 2
∆3 f − ⋅
1
+
1⋅ 2 ⋅ 3 2
1 3 3 2 1 1
x − x − x + ∆ f − dx
3 3
+
6 2 4 8 2
2 2
1 x2 1 1 x3 1
= f − [ x]2 + ∆ f − + x + ∆2 f −
1 1
− x
2 0 2 2 2 2 2 3 4
0 0
2
1 3 1 x4 x3 x2 3
+ ∆ f − − − + x
6 2 4 2 8 8
0
= 2 f − + 3 ∆ f − +
1 1 13 2 1 1 3 1
∆ f − + ∆ f −
2 2 12 2 24 2
= 2 f − + 3 f − f − +
1 1 1 13 3 1 + f − 1
f −2f
2 2 2 12 2 2 2
1 5 3 + 3 f 1 − f − 1
+ f −3f
24 2 2 2 2
f − +
23 1 23 3
f
1 1 1 5
= f + f +
24 2 24 2 24 2 24 2
1 1
f − + 23 f + 23 f + f ⋅
1 3 5
=
24 2 2 2 2
Solution: We have ∆ f − = f − − f −
3 1 3
2 2 2
N-115
∆ f = f − f ⋅
1 3 1
and
2 2 2
1 /2
So here we are to obtain a formula for
∫−1 /2 f ( x ) dx in terms of
f − , f − , f , f ⋅
3 1 1 3
2 2 2 2
Since we are known four pairs of values of x and f ( x ), therefore f ( x ) can be represented as a
polynomial in x of degree 3. Consequently third differences of f ( x ) are constant.
3
Here x0 = − and h = 1.
2
f ( x ) = f − + x +
3 3
Now [Form f ( x 0 + uh)]
2 2
= E x + (3 /2) f −
3
2
= (1 + ∆ ) x + (3 /2) f −
3
2
x + 3 x + 1
3 3 3 2 2 ⋅ ∆2 f − 3
= f − + x + ∆ f − +
2 2 2 1⋅ 2 2
x + 3 x + 1 x − 1
2 2 2 ∆3 f − 3
+
1⋅ 2 ⋅ 3 2
= f − + x + ∆ f − + x2 + 2 x + ∆2 f −
3 3 3 1 3 3
2 2 2 2 4 2
1 3 3 2 1
x + x − x − ∆ f − ⋅
3 3 3
+
6 2 4 8 2
1 /2 3 1 /2 3 1 /2
f ( x ) dx = f − dx + ∆ f − x + dx
3
∴
∫−1 /2 ∫
2 −1 / 2 2 −1 / 2 ∫
2
1 2 3 1 /2 2
x + 2 x + dx
3
+
2
∆ f −
2 −1 / 2 ∫ 4
1 3 3 1 /2 3 3 2 1
x + x − x − dx
3
+
6
∆ f −
∫
2 −1 / 2 2 4 8
3 1 /2 3 1 /2 3
= 2 f − dx + 2 ∆ f −
2 0 ∫ 2 0 2
dx
∫
1 2 3 1 /2 2 3 3 1 /2 3 2 3
x + dx +2 ⋅ ∆ f −
1 3
+ 2⋅
2
∆ f −
2 0 ∫ 4 6
x − dx
2 0 2 8 ∫
0, if f ( − x ) = − f ( x )
∵ a
∫ f ( x ) dx = a
−a
∫
2 0 f ( x ) dx, if f ( − x ) = f ( x )
N-116
= f − + ∆ f − +
3 3 3 10 2 3 1 3 3
∆ f − − ∆ f −
2 2 2 24 2 24 2
= f − + f − − f − + f − 2 f − + f −
3 3 1 3 5 1 1 3
2 2 2 2 12 2 2 2
1 3
f − 3 f + 3 f − − f −
1 1 3
−
24 2 2 2 2
f − +
13 1 13 1
f
1 3 1 3
=− f − + f −
24 2 24 2 24 2 24 2
12 1 12 1
f − + f − f − − f
1 1 1 1 1 3 1 3
= f − + f +
24 2 24 2 24 2 24 2 24 2 24 2
=
1 1 1 1 f − 1 − f − 3 − f 3 − f 1
f − + f +
2 2 2
2 2 2 24 2
1 1 1 1 3 1
= f − + f + ∆ f − − ∆ f ⋅
2 2 2 24 2 2
b
Solution: Here we are to express ∫ f ( x ) dx in terms of f ( − a), f (0 ) and f ( b ).
−a
x2 + ax
+ f ( b ).
b ( b + a)
b
b 1 x3 bx2
∴
∫−a f ( x ) dx =
a (a + b)
f ( − a)
3
−
2
−a
b b
1 x3 x2 1 x3 ax2
− f (0 ) + (a − b) − abx + f (b) +
ab 3 2 − a b ( b + a) 3 2
−a
1 b3 b3 a3 ba2
= − + + f ( − a)
a (a + b) 3 2 3 2
1 b3 b2 a3 a2 ( a − b )
− + (a − b) − ab2 + − − a2 b f (0 )
ab 3 2 3 2
1 b3 ab2 a3 a3
+ + + − f (b)
b ( b + a) 3 2 3 2
N-117
1 1
= (2a3 + 3a2 b − b3 ) f ( − a) + ( a3 + 3a2 b + 3ab2 + b3 ) f (0 )
6a ( a + b ) 6ab
1
+ (2b3 + 3b2 a − a3 ) f ( b )
6 b ( b + a)
1 1
= ( a + b )2 (2a − b ) f ( − a) + ( a + b )3 f (0 )
6a ( a + b ) 6ab
( b + a)2 (2b − a)
+ f (b)
6 b ( b + a)
a+ b
= [ b (2a − b ) f ( − a) + ( a + b )2 f (0 ) + a (2b − a) f ( b )].
6ab
Problem 18: Prove that if f ( x) is a continuous function whose fifth differences are constant, then
1 8 5
∫−1 f ( x ) dx =
9
f (0 ) + [ f { √ (0 ⋅ 6)} + f { − √ (0 ⋅ 6)}] ⋅
9
Solution: Proceed as in problem 17.
f (3n). Since we are given four pairs of values of x and f ( x ) i. e.,(0, f (0 )), ( n, f ( n)), (2n, f (2n)) and
(3n, f (3n)), therefore f ( x ) can be represented as a polynomial in x of degree 3. Consequently
fourth and higher differences of f ( x ) are all zero.
We have x 0 = 0 and h = n.
n
Let I =
∫0 f ( x ) dx.
x−0
Put = u i. e., x = nu.
n
Then dx = n du.
1
∴ I =n
∫0 f ( nu) du
1
=n
∫0 E u f (0 ) du [∵ E u f (0 ) = f (0 + hu) = f (0 + nu) = f ( nu)]
1
=n
∫0 (1 + ∆ )u f (0 ) du
1 u ( u − 1) 2 u ( u − 1) ( u − 2) 3
∆ f (0 ) du
=n
∫0
f (0 ) + u ∆ f (0 ) +
1⋅ 2
∆ f (0 ) +
1⋅ 2 ⋅ 3
1 1 2 1 3
= n f (0 ) + ∆ f (0 ) − ∆ f (0 ) + ∆ f (0 )
2 12 24
1 1
= n f (0 ) + { f ( n) − f (0 )} − { f (2n) − 2 f (n) + f (0 ) }
2 12
= n 1 − −
1 1 1 1 + 1 + 1 f ( n) − 1 + 1 f ( 2n) + 1 f (3n)
− f (0 ) +
2 12 24 2 6 8 12 8 24
3 1
= n f (0 ) + {19 f ( n) − 5 f (2n) + f (3n)} ⋅
8 24
7. See article 6.
8. Newton Raphson method is used for finding successively better approximations to the
roots of real-valued function.
9. See article 5.
10. See article 3.
x0 + nh n2 n3 n2 ∆2 y0
+ … upto ( n + 1) terms .
∫x0 y dx = h ny0 +
2
∆ y0 +
3
−
2
2 !
See article 2.
2. The trapezoidal rule is applicable to any number of sub-intervals whether even or odd.
3. The formula obtained by putting n = 2 in the general quadrature formula is Simpson’s
one-third rule.
N-119
n
7. If Ckn are Cote’s numbers, then Σ Ckn = 1.
k= 0
True or False
1. We know that the process of computing the value of a definite integral from a set of
2. Simpson’s one-third rule requires the division of the whole range of integration into an
3. Trapezoidal rule is
xn h
∫x0 y dx =
2
[ y0 + 2 ( y1 + y2 + … + yn −1 ) + yn ].
Comprehensive Problems 1
Problem 1: Compute an approximate value of (1⋅ 04) 3⋅01.
= (3) (1) (⋅04) + (1)3 log (1) (⋅01) = 0 ⋅12. [Using (2)]
Problem 2: The deflection at the centre of a rod of length l and diameter d supported at its ends and loaded at
the centre with a weight w varies as wl3 d − 4 . What is the percentage increase in the deflection corresponding to
the percentage increase in w, l and d of 3, 2 and 1 respectively ?
w l3
Solution: Let D be the deflection of the rod at the centre, then D = k
d4
or log D = log k + log w + 3 log l − 4 log d.
δD δw δl δd
Differentiating, we get = +3 −4
D w l d
100 δ D 100 δ w 100 δ l δd
or = +3× − 4 × 100 = 3 + 3 × 2 − 4 × 1 = 5.
D w l d
Hence, the required percentage increase in the deflection is 5%.
Problem 3: Show that the acceleration due to gravity is reduced nearly by 1% at an altitude equal to
0 ⋅ 5% of earth’s radius. Given that at an external point x kilometers from the earth’s centre, such an
2
acceleration is given by g , where r is the radius of the earth.
r
x
Solution: Let a be the acceleration due to gravity, then we have a = g ( r / x)2 . Taking log , we
get log a = log g + 2 log r − 2 log x .
δa δx
∴ = − 2 , since g and r are constants
a x
δa δx Given δx × 100 = 0 ⋅ 5
i. e., × 100 = − 2 × 100 = − 2 × 0 ⋅ 5 = − 1.
a x x
Problem 4: In estimating the number of bricks in a pile which is measured to be 5m × 10 m × 5m, count of
bricks is taken as 100 bricks per m3 . Find the error in the cost when the tape is stretched 2% beyond its
standard length. The cost of bricks is ` 2,000 per thousand bricks.
Solution: Let the length, breadth and height of the pile be x, y and z respectively so that its
volume V = x y z .
Taking log, we get log V = log x + log y + log z
δV δ x δ y δ z
∴ = + +
V x y z
δ V 100 δ x 100 δ y 100 δ z
or 100 = + + =2+ 2+ 2=6
V x y z
6V 6 (5 × 10 × 5)
or δV = = = 15 m3 .
100 100
Number of bricks in δ V = 15 × 100 = 1500.
1500 × 2000
Thus the error in the cost = = ` 3000
1000
which results in a loss to the brick seller.
Problem 5: Find expression for the remainder term in Bessel's formula when terminated with a difference
of even order.
to compute the value of log (1⋅2) correct to seven decimal places and find the number of terms retained.
Solution: We have
1+ x x3 x5 x2 n −1
log e = 2 x + + + …+ + Rn( x )
1− x 3 5 2n −1
1+ x 1
Let =1⋅ 2 i. e., x = ⋅
1− x 11
N-122
2 n +1
1 2 1
Now ( ER )max. at x = is = .
11 2n + 1 11
which gives n ≥ 3.
Hence retaining the first three terms of the given series, we have
x3 x5
at x = = 0 ⋅1823215.
1
log e(1⋅ 2) = 2 x + +
3 5 11
Problem 7: Derive the expression for the remainder in Newton's forward formula
Solution: See article 8.
EA 0 .21356 × 10 −3
Relative error ER = =
True Value 2
= 0 .151 × 10 −3 .
2. See article 4.
2
3. True value of = 0 .666666
3
2
Approximate value of = 0 .667
3
EA .000334
Relative error ER = = = 0 .0005.
True value .666666
N-123
1. See article 4.
2. See article 4.
True or False
1. See article 3.
2. Here E A = 625483
. − 625 = 0.483
EA .483
and ER = = =.000772.
625.483 625.483
EP = ER × 100 = .077%.
❍❍❍
N-124
Chapter-9
Simultaneous Linear Algebraic Equations
Comprehensive Problems 1
Problem 1: Solve the following equations by factorization method i. e., by triangularization method :
x1 + x 2 + x 3 = 1, 3 x1 + x 2 − 3 x 3 = 5, x1 − 2 x 2 − 5 x 3 = 10.
1 1 1 x1 1
i. e., AX = B, where A = 3 1 −3 , X = x2 , B = 5 ⋅
1 −2 −5 x3 10
1 1 1
=3 1 −3 ⋅
1 −2 −5
1 0 0 1 1 1
∴ L = 3 1 0 ,U = 0 −2 −6 ⋅
3 0 0 3
1 1
2
N-125
Now AX = B
⇒ LUX = B [ ∵ A = LU ]
y1
⇒ LY = B, where UX = Y = y2
y3
1 0 0 y1 1
⇒ 3 1 0 y = 5
2
1 3 1 y3 10
2
y1 1
⇒ 3 y1 + y2 = 5
3 10
y1 + y 2 + y3
2
3
⇒ y1 = 1, 3 y1 + y2 = 5, y1 + y2 + y3 = 10
2
⇒ y1 = 1, y2 = 2, y3 = 6.
Again UX = Y
1 1 1 x1 1
⇒ 0 −2 −6 x2 =2
0 0 3 x3 6
x1 + x2 + x3 1
⇒ − 2 x − 6x =2
2 3
3 x3 6
⇒ x1 + x2 + x3 = 1, − 2 x2 − 6 x3 = 2 , 3 x3 = 6
⇒ x3 = 2 , x2 = − 7, x1 = 6.
Hence, the required solution of the given system of equations is
x1 = 6, x2 = − 7, x3 = 2 .
into the form LU, where L is unit lower triangular matrix and U is upper triangular matrix and hence solve
the system of equations
5 x − 2 y + z = 4, 7 x + y − 5z = 8, 3 x + 7 y + 4z = 10.
Solution: The given system of equations in matrix form AX = B is
5 −2 1 x 4
7 1 −5 y = 8 ⋅
3 7 4 z 10
Now let A = LU
5 −2 1 1 0 0 u11 u12 u13
i. e., 7 1
−5 = l21 1 0 0 u22 u23
3 7 4 l31 l32 1 0 0 u33
1 5 −2
0 0
1
19 32
7 −
∴ L= 1 0 and U = 0 5 5 ⋅
5
327
3 41 0 0
5 19 1 19
∵ AX = B and A = LU , therefore
LUX = B.
Now let
UX = Y so that LY = B
N-128
1 0 0 y1 4
7
i. e., 1 0 y2 8
5
3 41 y 10
1
3
5 19
7 12
giving y1 = 4 ; y1 + y 2 = 8 ⇒ y 2 =
5 5
3 41 3 41 12 46
and y1 + y 2 + y 3 = 10 ⇒ ⋅ 4 + ⋅ + y3 = 0 ⇒ y3 = ⋅
5 19 5 19 5 19
Now UX = Y implies
5 −2 1 4
x
0 19 − 32 y = 12
5 5 5
z
327 46
0 0
19 19
327 46 46
giving z= or z = ;
19 19 327
19 32 12 19 32 46 12 284
y− z= or y− × = or y = ;
5 5 5 5 5 327 5 327
284 46 122
and 5 x − 2 y + z = 4, or 5 x − 2 × + = 4 or x = ⋅
327 327 109
∴ The solution of the given system of equations is
122 284 46
x= , y= and z = ⋅
109 327 327
x − 3y − z = − 30
2x − y − 3z = −5
5x − y − 2z = 142 .
To eliminate x from the second and third equations with the help of first equation we perform
R2 → R2 − 2 R1 and R3 → R3 − 5R1 . So the given equations are equivalent to
x −3y − z = − 30
5y − z = 55
14 y + 3z = 292
or
x −3y − z = − 30
1
y − (1 / 5) z = 11 , by R2 → R2
5
14 y + 3z = 292.
N-129
Now to eliminate y from the third equation with the help of second equation we perform
R3 → R3 − 14 R2 . So the given equations are equivalent to
x −3y − z = − 30 …(1)
y − (1 / 5) z = 11 …(2)
(29 / 5) z = 138. …(3)
or x1 − x2 + x3 = 6
6 x2 − x3 = −9 , by R2 → R2 − 2 R1
5 x2 − 5 x3 = −20 R3 → R3 − 3R1
or
x1 − x2 + x3 = 6
by R2 → 1 / 6 R2
x2 − (1 / 6) x3 = ( −3 / 3) ,
R3 → 1 / 5 R3
x2 − x3 = −4
or
x1 − x2 + x3 = 6 …(1)
x2 − (1 / 6) x3 = − 3 /2, …(2)
− (5 / 6) x3 = − 5 /2 …(3)
by R3 → R3 − R2
To eliminate x from the second and third equations with the help of first equation we perform
R2 → R2 − 4 R1 and R3 → R3 − 2 R1 . So the given equations are equivalent to
2x + y + 4z = 12
− 7y − 14 z = −25
9y − 9z = 9
or
2x + y + 4z = 12
by R2 → ( −1 / 7) R2
y + 2z = (25 / 7) ,
R3 → (1 / 9) R3
y − z = 1
or
2x + y + 4z = 12 …(1)
y + 2z = (25 / 7) …(2)
− 3z = ( −18 / 7) , by R3 → R3 − R2 . …(3)
Solution: Let us put the given system of equations in the matrix form AX = B
2 −3 10 x 3
i. e., −1 4 2 y = 20 ⋅
5 2 1 z −12
Now let
A = LU
N-131
2 −3 10 0 0 u11 u12
1 u13
i. e., −1 4 2= 1 0 0 u22
l21 u23
5 2 1 l31 l32 1 0
0 u33
2 −3 10 u11 u12 u13
or −1 4 2 = l21 u11 l21 u12 + u22 l21 u13 + u23
5 2 1 l31 u11 l31 u12 + l32 u22 l31 u13 + l32 u23 + u33
1 0 0 y1 3
i. e., −1/2 1 0 y2 = 20
5 / 2 19 / 5 1 y3 −12
43 1012
giving y1 = 3, y2 = and y3 = − ⋅
2 10
Now UX = Y implies
2 −3 10 x 3
0 5 /2 7 y= 43 / 2
0 0 − 253 / 5 z − 1012 / 10
giving 2 x − 3 y + 10 z = 3,
5 43 253 −1012
y + 7z = and − z= ⋅
2 2 5 10
Hence z = 2 , y = 3 and x = − 4.
∴ The solution of the given system of equations is
x = − 4, y = 3 and z = 2 .
N-132
Solution: The given system of equations can be written in the form AX = B i. e.,
1 2 3 x1 14
2 5 2 x 18 ⋅
2 =
3 1 5 x3 20
←A→ ← X→ ←B→
Now let A = LU , where L is a unit lower triangular matrix and U is an upper triangular matrix.
Thus we have
1 2 3 1 0 0 u11 u12 u13
2 5 2 = l 1 0 0 u22 u23
21
3 1 5 l31 l32 1 0 0 u33
1 2 3 u11 u12 u13
or 2 5 2 =l u l21 u12 + u22 l21 u13 + u23 ⋅
21 11
3 1 5 l31 u11 l31 u12 + l32 u22 l31 u13 + l32 u23 + u33
y1
Now let UX = Y = y2 ⇒ LY = B i. e.,
y3
1 0 0 y1 14
2 1 0 y2 = 18 , giving
3 −5 1 y3 20
y1 = 14, 2 y1 + y2 = 18, 3 y1 − 5 y2 − y3 = 20
and y1 = 14, y2 = − 10, y3 = − 72.
N-133
Now UX = Y implies
1 2 3 x1 14
0
1 −4 x 2 = −10 ⋅ , giving
0 0 −24 x 3 −72
− 24 x 3 = − 72 i.e., x 3 = 3,
x 2 − 4 x 3 = − 10 ⇒ x 2 = 4 × 3 − 10 = 2
and x1 + 2 x 2 + 3 x 3 = 14 ⇒ x1 + 2 ⋅ 2 + 3 ⋅ 3 = 14 ⇒ x 3 = 1.
Hence the solution of the given equations is
x1 = 1, x 2 = 2 and x 3 = 3.
Solution: Since the requirement for iteration is satisfied by these equations, we solve each
equation for the unknown having the largest coefficient and the new equations are
1
x= (95 − 11 y + 4 z ), …(1)
83
1
y= (104 − 7 x − 13z ), …(2)
52
1
and z= (71 − 3 x − 8 y ). …(3)
29
Jacobi Iterative method: Starting with x = 0, y = 0, z = 0, we get
95 104 71
x (1) = = 1⋅14 ; y (1) = = 2 and z (1) = = 2 ⋅ 45.
83 52 29
These are first approximations. We proceed to obtain the second approximations as follows :
1 1
x (2) = [95 − 11 y (1) + 4 z (1) ] = (95 − 11 × 2 + 4 × 2 ⋅ 45) = ⋅ 998,
83 83
1 1
y (2) = [104 − 7 x (1) − 13z (1) ] = (104 − 7 × 1⋅14 − 13 × 2 ⋅ 45) = 1⋅ 234,
52 52
1 1
and z (2) = [71 − 3 x (1) − 8 y (1) ] = (71 − 3 × 1⋅14 − 8 × 2) = 1⋅ 779.
29 29
1 1
x (3) = [95 − 11 y (2) + 4 z (2) ] = (95 − 11 × 1⋅ 234 + 4 × 1⋅ 779) = 1⋅ 07,
83 83
1 1
y (3) = [104 − 7 x (2) − 13 z (2) ] = (104 − 7 × ⋅ 998 − 13 × 1⋅ 779) = 1⋅ 421,
52 52
1 1
and z (3) = [71 − 3 x (2) − 8 y (2) ] = (71 − 3 × ⋅ 998 − 8 × 1⋅ 234) = 2 ⋅ 005.
29 29
Continuing in this way further iterations can be obtained.
N-135
2 x + 10 y + z = 13,
2 x + 2 y + 10 z = 14.
N-136
Solution: To apply Jacobi iteration method we start with x = 0, y = 0, z = 0 as the initial values
of x, y, z .
The given system of equations can be written as
12 y z 12
x= − − ⇒ x (1) = (first approximation for x)
10 10 10 10
13 x z 13
y= − − ⇒ y (1) = (first approximation for y)
10 5 10 10
14 x y 14
and z= − − ⇒ z (1) = (first approximation for z).
10 5 5 10
12 y (1) z (1) 12 13 14 93
Again x (2) = − − = − − =
10 10 10 10 10 × 10 10 × 10 100
13 x (1) z (1) 13 12 14 92
y (2) = − − = − − =
10 5 10 10 10 × 5 10 × 10 100
14 x (1) y (1) 14 12 13 90
z (2) = − − = − − = ⋅
10 5 5 10 10 × 5 10 × 5 100
Repeating the process with x (2), y (2), z (2), we get
12 y (2) z (2) 12 92 90
x (3) = − − = − − = 1⋅ 018
10 10 10 10 100 × 10 100 × 10
13 x (2) z (2) 13 93 90
y (3) = − − = − − = 1⋅ 024
10 5 10 10 100 × 5 100 × 10
14 x (2) y (2) 14 93 92
z (3) = − − = − − = 1⋅122 .
10 5 5 10 100 × 5 100 × 10
Similarly x (4) = 1⋅ 2 − ⋅1024 − ⋅1122 = 0 ⋅ 9854
y (4) = 1⋅ 3 − ⋅ 2036 − ⋅1122 = ⋅ 9842
z (4) = 1⋅ 4 − ⋅ 2036 − ⋅ 2048 = ⋅ 9916.
Again x (5) = 1⋅ 2 − 0 ⋅ 09842 − 0 ⋅ 09916 = 1⋅ 00242
y (5) = 1⋅ 3 − 0 ⋅19708 − 0 ⋅ 09916 = 1⋅ 00376
z (5) = 1⋅ 4 − 0 ⋅19708 − 0 ⋅19884 = 1⋅ 00408.
True or False
1. See article 2(II), part (c).
2. See article 2(I), part (b).
3. We know that the linear equations expressed in the matrix form AX = B are consistent if
and only if rank A = rank [ A B].
❍❍❍
N-138
Chapter-10
Curve Fitting
Comprehensive Problems 1
Problem 1: Find the line of fit to the following data :
x : 0 5 10 15 20 25
y : 12 15 17 22 24 30
Solution: Here the number of values of x is 6 i. e. even and the values of x are equidistant at an
interval of 5 i. e., h = 5, so we take
{x − (10 + 15) / 2} x − 12 ⋅ 5
u= =
2⋅5 2⋅5
v = y − 20
Let the least square line of v on u be v = a + bu. Then the normal equations are
Σv = ma + bΣu
and Σvu = aΣu + bΣu2
x y u v uv u2
0 12 −5 −8 40 25
5 15 −3 −5 15 9
10 17 −1 −3 3 1
15 22 1 2 2 1
20 24 3 4 12 9
25 30 5 10 50 25
0 0 122 70
x : 0 1 2 3 4
Solution: Here the number of values of x is odd and the values of x are equidistant. So we take the
origin at the middle value 2 for the x series.
Let us take X = x − 2 and Y = y. Let the curve of best fit be Y = a + bX + cX2 . Then the normal
equations are
ΣY = na + b Σ X + c Σ X2
Σ XY = a Σ X + b Σ X2 + c Σ X3 .
2 2 3 4
Σ X Y = aΣ X + bΣ X + c Σ X
x y X Y XY X2 X2Y X3 X4
0 1 –2 1 –2 4 4 –8 16
2 1⋅3 0 1⋅3 0 0 0 0 0
Substituting the values from the table in the normal equations, we get
12 ⋅ 9 = 5a + 10 c,11⋅ 3 = 10 b and 33 ⋅ 5 = 10 a + 34c.
Solving these equations simultaneously, we get
a = 1⋅ 48, b = 1⋅13, c = 0 ⋅ 55.
Hence the curve of best fit is
Y = 1⋅ 48 + 1⋅13 X + 0 ⋅ 55 X2 .
or y = 1⋅ 42 − 1⋅ 07 x + 0 ⋅ 55 x2 .
Problem 3: Fit a second degree parabola to the following data taking x as the independent variables :
x : 1 2 3 4 5 6 7 8 9
y : 2 6 7 8 10 11 11 10 9
Σ XY = a Σ X + b Σ X2 + c Σ X3
Σ X2Y = a Σ X2 + b Σ X3 + c Σ X4
x y X Y XY X2 X2Y X3 X4
2 6 −3 −1 3 9 −9 −27 81
3 7 −2 0 0 4 0 −8 16
4 8 −1 1 −1 1 1 −1 1
5 10 0 3 0 0 0 0 0
6 11 1 4 4 1 4 1 1
7 11 2 4 8 4 16 8 16
8 10 3 3 9 9 27 27 81
9 9 4 2 8 16 32 64 256
0 11 51 60 −9 0 708
Substituting the values from the table in the normal equations, we get
11 = 9a + 0 + 60 c,
51 = 0 + 60 b, − 9 = 60 a + 0 + 708c
Solving these equations simultaneously, we get
a = 3, b = 0 ⋅ 85, c = − 0 ⋅ 27
Hence the curve of best fit is
Y = 3 + 0 ⋅ 85 X − 0 ⋅ 27 X2
= 3 + 0 ⋅ 85 x − 4 ⋅ 25 − 0 ⋅ 27 x2 + 2 ⋅ 7 x − 6 ⋅ 75
or y = − 1 + 3 ⋅ 55 x − 0 ⋅ 27 x2
x : 1 2 3 4 5
Solution: Changing the origin by u = x − 3, v = ( y − 1450 ) / 5, let the curve of best fit
be v = a + bu + cu2 .
Σ v = na + b Σ u + c Σ u2
Σ uv = a Σu + b Σ u2 + c Σ u3
Σ u2 v = a Σu2 + b Σ u3 + c Σ u4
x y u v u2 u4 uv u2 v
3 1390 0 –12 0 0 0 0
4 1625 1 35 1 1 35 35
Total 0 –2 10 34 411 73
Changing the origin back i. e, putting u = x − 3 and v = ( y − 1450 ) / 5, the required curve of best fit
is
1
( y − 1450 ) = − 11⋅ 4 + 41⋅1 ( x − 3) + 5 ⋅ 5 ( x − 3)2
5
or y = 1024 + 40 ⋅ 5 x + 27 ⋅ 5 x2 .
Problem 5: If P is a pull required to lift a load W by means of a pulley block, find a linear law of the form
P = mW + c connecting P and W, using the following data :
P : 12 15 21 25
W : 50 70 100 120
Solution: We have
P = mW + c.
The normal equations are
Σ P = 4c + m Σ W
ΣWP = c ΣW + m ΣW 2
N-142
W P W2 WP
50 12 2500 600
70 15 4900 1050
Substituting the values from the table in the normal equations, we get
73 = 4c + 340 m; 6750 = 340 c + 31800 m
Solving these equations simultaneously, we get
m = 0 ⋅1879, c = 2 ⋅ 2785
Hence the curve of best fit is
P = 2ɺ ⋅ 2759 + 0 ⋅1879W
Now, when W = 150 kg, then P = 30 ⋅ 4635 kg.
x : 1 2 3 4 5 6 7 8
x y Y = log y xY x2
1 1⋅0 0⋅0000 0⋅0000 1