CH 4 - Finite Differences
CH 4 - Finite Differences
Numerical Analysis I
Prepared By
Biruk Endeshaw Mekonnen
March 8, 2021
Chapter 4
Finite Differences
4.1 INTRODUCTION
The calculus of finite differences deals with the changes that take place in the value of the
dependent variable due to finite changes in the independent variable from this we study the
relations that exist between the values, which can be assumed by function, whenever the
independent variable changes by finite jumps whether equal or unequal.
4.2 Differences
Consider a function y = f (x) defined on (a, b). x and y are the independent and dependent
variables respectively. If the points x0 , x1 , . . . , xn are taken at equidistance i.e., xi = x0 + ih, i =
0, 1, 2, . . . , n, then the value of y, when x = xi , is denoted as yi , where yi = f (xi ). Here, the val-
ues of x are called arguments and the values of y are known as entries. The interval h is called
the difference interval. To determine the values of function y = f (x) for given intermediate
values of x, the following three types of differences are useful:
The following table shows how the forward differences of all orders can be formed.
x y ∆y ∆2 y ∆3 y ∆4 y
x0 y0
∆y0
x1 y1 ∆2 y0
∆y1 ∆3 y0
x2 y2 ∆2 y1 ∆4 y0
∆y2 ∆3 y1
2
x3 y3 ∆ y2
∆y3
x4 y4
Table 4.1: Forward differences
Example 4.1. Construct a forward difference table for the following values:
x 0 5 10 15 20
f (x) 7 11 14 18 24
Solution: Forward difference table for given data is:
x y ∆y ∆2 y ∆3 y ∆4 y
0 7
4
5 11 -1
3 2
10 14 1 -1
4 1
15 18 2
6
20 24
Solution: For x = 0, 1, 2, 3, 4, 5, we get the values of y are 1, 1, 9, 31, 73, 141. Therefore, differ-
ence table for these data is as:
x f (x) ∆y ∆2 y ∆3 y
0 1
0
1 1 8
8 6
2 9 14
22 6
3 31 20
42 6
4 73 26
68 6
5 141 32
100
6 241
x f (x) ∆y ∆2 y ∆3 y ∆4 y ∆5 y
0 3
9
1 12 60
69 -10
2 81 50 -259
119 269 755
3 200 -219 496
-100 227
4 100 8
-92
5 8
Solution: We know that ∆ex =ex+h − ex = ex (eh − 1), where h is the interval of differencing.
x
Therefore, ex = eh1−1 ∆ex = ∆ ehe−1
ex
Hence, required function is given by eh −1
.
x y ∇y ∇2 y ∇3 y ∇4 y
x0 y0
∇y1
x1 y1 ∇2 y2
∇y2 ∇3 y3
x2 y2 ∇2 y3 ∇4 y4
3
∇y3 ∇ y4
2
x3 y3 ∇ y4
∇y4
x4 y4
Table 4.2: Backward differences
Hence, ∇5 f (6) = 0.
Example 4.6. Construct the backward difference table for y = log x given that:
x 10 20 30 40 50
y 1 1.3010 1.4771 1.6021 1.6990
and find the values of ∇3 log 40 and ∇4 log 50.
x y ∇y ∇2 y ∇3 y ∇4 y
10 1
0.3010
20 1.3010 -0.1249
0.1761 0.0738
30 1.4771 -0.0511 -0.0508
0.1250 0.0230
40 1.6021 -0.0281
0.0969
50 1.6990
It is noted in Table 4.3 that all odd differences have fraction suffices and all the even differ-
ences are with integer suffices.
i.e., ∇[af (x + h) + bg(x + h)] = [af (x + h) + bg(x + h)] − [af (x) + bg(x)]
= a[f (x + h) − f (x)] + b[g(x + h) − g(x)]
= a∇f (x + h) + b∇g(x + h)
There are few relations defined between these operators. Some of them are:
1. ∆ = E − 1 or E = 1 + ∆
2. ∇ = 1 − E −1 or E = (1 − ∇)−1
3. δ = E 1/2 − E −1/2
4. ∆ = E∇ = ∇E = δE 1/2
6. µ = 12 (E 1/2 + E −1/2 )
Proof:
Solution:
Solution:
2 5x + 12 2 5x + 12
∆ =∆
x2 + 5x + 6 (x + 2)(x + 3)
2 2 3 2 3
=∆ + =∆ ∆ +∆
x+2 x+3 x+2 x+3
1 1 1 1
=∆ 2 − +3 −
x+3 x+2 x+4 x+3
1 1
= −2∆ − 3∆
(x + 2)(x + 3) (x + 3)(x + 4)
1 1 1 1
= −2 − −3 −
(x + 3)(x + 4) (x + 2)(x + 3) (x + 4)(x + 5) (x + 3)(x + 4)
4 6
= +
(x + 2)(x + 3)(x + 4) (x + 3)(x + 4)(x + 5)
2(5x + 16)
=
(x + 2)(x + 3)(x + 4)(x + 5)
Example 4.9. Evaluate ∆n eax+b ; where the interval of differencing taken to be unity.
n!hn
n1
∆ = (−1)n
x x(x + h) . . . (x + nh)
Solution:
n 1 n−1 1 n−1 1 1
∆ =∆ ∆ =∆ −
x x x+h x
n−1 −h
=∆
x(x + h)
n−2 1
= (−h)∆ ∆
x(x + h)
n−2 1 1
= (−1)∆ ∆ −
x x+h
n−2 1 1 1 1
= (−1)∆ − − −
x+h x x + 2h x + h
n−2 2 1 1
= (−1)∆ − −
x + h x x + 2h
−2h2
n−2
= (−1)∆
x(x + h)(x + 2h)
2!h2
2 n−2
= (−1) ∆
x(x + h)(x + 2h)
3!h3
3 n−3
= (−1) ∆
x(x + h)(x + 2h)(x + 3h)
..
.
n!hn
= (−1)n
x(x + h) . . . (x + nh)
h i
∆f (x)
Example 4.12. Show that ∆ log f (x) = log 1 + f (x)
Example 4.13. Sum the following series 13 + 23 + 33 + . . . + n3 using the calculus of finite
differences.
S = u0 + u1 + u2 + . . . + un−1
= (1 + E + E2 + E3 + . . . + En−1 )u0
n
(1 + ∆)n − 1
E −1
= u0 = u0
E−1 ∆
1 n(n − 1) 2 n(n − 1)(n − 2) 3 n
= 1 + n∆ + ∆ + ∆ + . . . + ∆ − 1 u0
∆ 2! 3!
n(n − 1) n(n − 1)(n − 2)
=n+ ∆u0 + ∆u0 + . . .
2! 3!
We know ∆u0 = u1 − u0 = 23 − 13 = 7.
∆2 u0 = u2 − 2u1 + u0 = 33 − 2(2)3 + 13 = 12.
Example 4.14. Given that u0 = 3, u1 = 12, u2 = 81, u3 = 200, u4 = 100, u5 = 8. Find the value of
∆5 u0 .
∆5 u0 = (E − 1)5 u0
= (E 5 − 5E 4 + 10E 3 − 10E 2 + 5E − 1)u0
7. Evaluate: