Interpolation Formulae

Download as pdf or txt
Download as pdf or txt
You are on page 1of 11

Interpolation Formulae

2.1 Introduction :

Suppose, 𝑦 = 𝑓(𝑥) is a given function. Then, 𝑓(𝑥0 ), 𝑓(𝑥1 ), 𝑓(𝑥2 ), ……….,


𝑓(𝑥𝑛 ) are the corresponding values for the values of independent
variable 𝑥, i.e., 𝑥0 , 𝑥1 , 𝑥2 , ………., 𝑥𝑛 respectively with respect to the given
function. When, we have to find out some intermediate terms or values
of a given function with given set of values (both arguments and
entries), we apply some special processes which are called
Interpolations.

Interpolations are applied for both with Equal and Unequal Intervals.

2.2 Newton-Gregory Formula for Forward Interpolation with Equal


Intervals :

Suppose, 𝑦 = 𝑓(𝑥) is a given function.

Suppose, 𝑓(𝑎), 𝑓(𝑎 + ℎ), 𝑓(𝑎 + 2ℎ), ………. , 𝑓(𝑎 + 𝑛ℎ) are (𝑛 + 1)
corresponding values for (𝑛 + 1) equidistant values 𝑎, 𝑎 + ℎ, 𝑎 +
2ℎ,………. , 𝑎 + 𝑛ℎ of the independent variable 𝑥 respectively, the
interval of differencing is ℎ.

Now, we consider, 𝑓(𝑥) is a nth degree polynomial in 𝑥, such that


𝑓(𝑥) = 𝐴0 + 𝐴1 (𝑥 − 𝑎) + 𝐴2 (𝑥 − 𝑎)(𝑥 − 𝑎 − ℎ)
+ 𝐴3 (𝑥 − 𝑎)(𝑥 − 𝑎 − ℎ)(𝑥 − 𝑎 − 2ℎ) +……………………
+ 𝐴𝑛 (𝑥 − 𝑎)(𝑥 − 𝑎 − ℎ)……(𝑥 − 𝑎 − ̅̅̅̅̅̅̅
𝑛 − 1 ℎ) …………. (1)
where 𝐴0 , 𝐴1 , 𝐴2 , … … … . , 𝐴𝑛 are constants to be determined.
Putting successively 𝑥 = 𝑎, 𝑎 + ℎ, 𝑎 + 2ℎ, … … … , 𝑎 + 𝑛ℎ in (1), we have,
𝑓(𝑎) = 𝐴0 ( 𝑤ℎ𝑒𝑛 𝑥 = 𝑎 ) ……………………… (2)

𝑓(𝑎 + ℎ) = 𝐴0 + 𝐴1 ℎ ( 𝑤ℎ𝑒𝑛 𝑥 = 𝑎 + ℎ )
 𝑓(𝑎 + ℎ) = 𝑓(𝑎) + 𝐴1 ℎ
 𝐴1 ℎ = 𝑓(𝑎 + ℎ) − 𝑓(𝑎)
𝑓(𝑎+ℎ)−𝑓(𝑎)
 𝐴1 =

∆𝑓(𝑎)
 𝐴1 = ………………………. (3)

𝑓(𝑎 + 2ℎ) = 𝐴0 + 𝐴1 . 2ℎ + 𝐴2 . 2ℎ. ℎ


∆𝑓(𝑎)
 𝑓(𝑎 + 2ℎ) = 𝑓(𝑎) + . 2ℎ + 𝐴2 . 2ℎ2

 𝑓(𝑎 + 2ℎ) = 𝑓(𝑎) + 2[𝑓(𝑎 + ℎ) − 𝑓(𝑎)] + 𝐴2 . 2ℎ2


 𝑓(𝑎 + 2ℎ) = 2𝑓(𝑎 + ℎ) − 𝑓(𝑎) + 𝐴2 . 2ℎ2
 𝐴2 . 2ℎ2 = 𝑓(𝑎 + 2ℎ) − 2𝑓(𝑎 + ℎ) + 𝑓(𝑎)
 𝐴2 . 2ℎ2 = ∆2 𝑓(𝑎)
∆2 𝑓(𝑎) ∆2 𝑓(𝑎)
 𝐴2 = = …………………. (4)
2ℎ2 2 ! ℎ2

Proceeding in this way, we get

∆3 𝑓(𝑎)
𝐴3 =
3 ! ℎ3
∆4 𝑓(𝑎)
𝐴4 =
4 ! ℎ4

……………………
……………………
∆𝑛 𝑓(𝑎)
𝐴𝑛 =
𝑛 ! ℎ𝑛

Now, substituting these values of 𝐴0 , 𝐴1 , 𝐴2 , ……….. , 𝐴𝑛 in (1), we get


∆𝑓(𝑎) ∆2 𝑓(𝑎)
𝑓(𝑥) = 𝑓(𝑎) + (𝑥 − 𝑎) + (𝑥 − 𝑎)(𝑥 − 𝑎 − ℎ)
ℎ 2 ! ℎ2
∆3 𝑓(𝑎)
+ (𝑥 − 𝑎) (𝑎 − ℎ)(𝑥 − 𝑎 − 2ℎ) +……………….
3 ! ℎ3
∆ 𝑛 𝑓(𝑎)
+ (𝑥 − 𝑎)(𝑥 − 𝑎 − ℎ)……..…(𝑥 − 𝑎 − ̅̅̅̅̅̅̅
𝑛 − 1 ℎ) .……… (5)
𝑛 ! ℎ𝑛

Then, we put

𝑥 = 𝑎 + ℎ𝑢
𝑥−𝑎
 𝑢=

Then, the equation (5) becomes

∆𝑓(𝑎) ∆2 𝑓(𝑎)
𝑓(𝑎 + ℎ𝑢)=𝑓(𝑎) + 𝑢ℎ. + 𝑢ℎ(𝑢ℎ − ℎ)
ℎ 2 ! ℎ2
∆3 𝑓(𝑎)
+ 𝑢ℎ(𝑢ℎ − ℎ)(𝑢ℎ − 2ℎ) + ……………….
3 ! ℎ3
∆ 𝑛 𝑓(𝑎)
+ 𝑢ℎ(𝑢ℎ − ℎ)(𝑢ℎ − 2ℎ) … . . (𝑢ℎ − ̅̅̅̅̅̅̅
𝑛 − 1 h)
𝑛 ! ℎ𝑛

𝑢(𝑢−1) 𝑢(𝑢−1)(𝑢−2)
𝑓(𝑎 + ℎ𝑢)=𝑓(𝑎)+𝑢∆𝑓(𝑎)+ ∆2 𝑓(𝑎)+ ∆3 𝑓(𝑎)+
2! 3!
𝑢(𝑢−1)(𝑢−2)……..(𝑢−𝑛+1)
…………+ ∆𝑛 𝑓(𝑎) ………..…. (6)
𝑛!

This formula is known as Newton-Gregory Formula for Forward


Interpolation with equal interval.

2.3 Worked out Examples :

Example (1) : Given that sin 45°=0.7071, sin 50°=0.7660,


sin 55°=0.8192, sin 60°=0.8660. Find the value of sin 52°.
Solution : Let us consider the function as
𝑦 = 𝑓(𝑥) = 𝑠𝑖𝑛𝑥°
Here, 𝑎 = 45, ℎ = 5 and 𝑥 = 𝑎 + ℎ𝑢 = 52
𝑎 + ℎ𝑢 = 52
 45 + 5𝑢 = 52
 5𝑢 = 52 − 45 = 7
 𝑢 = 1.4

The difference table is

Argument Entry 1st Difference 2nd Difference 3rd Difference


𝑥 𝑓(𝑥) ∆𝑓(𝑥) ∆2 𝑓(𝑥) ∆3 𝑓(𝑥)

45 0.7071
0.0589
50 0.7660 −0.0057
0.0532 −0.0007
55 0.8192 −0.0064
0.0468
60 0.8660

Now applying Newton-Gregory Forward Interpolation Formula (Upto


third differences), we get
𝑢(𝑢−1) 𝑢(𝑢−1)(𝑢−2)
𝑓(𝑎 + ℎ𝑢) = 𝑓(𝑎) + 𝑢∆𝑓(𝑎) + ∆2 𝑓(𝑎)+ ∆3 𝑓(𝑎)
2! 3!
1.4(1.4−1)
∴ 𝑓(52) = 0.7071+1.4 x 0.0589 + x (−0.0057)
2!
1.4(1.4−1)(1.4−2)
+ x (−0.0007)
3!

⇒ 𝑠𝑖𝑛52° = 0.7071 + 0.08246−0.0016+0.00004


⇒ 𝑠𝑖𝑛52° = 0.7880 (approximately)
Hence, 𝑠𝑖𝑛52° = 0.7880 (approximately)

Example (2) : If 𝑙𝑥 represents the number of persons living at age 𝑥


In a life table, find as accurately as the data will permit 𝑙𝑥 for values of
𝑥 = 35, 45 𝑎𝑛𝑑 47. Given,
𝑙20 = 512, 𝑙30 = 439, 𝑙40 = 346, 𝑙50 = 243.

Solution : Given that


𝑙20 = 512, 𝑙30 = 439, 𝑙40 = 346, 𝑙50 = 243.

(i) First, we have to find out the value, when 𝑥 = 35 i.e., 𝑙35
Here, 𝑎 = 20, ℎ = 10 and 𝑥 = 𝑎 + ℎ𝑢 = 35
𝑎 + ℎ𝑢 = 35
 20 + 10𝑢 = 35
 10𝑢 = 35 − 20 = 15
 𝑢 = 1.5

The difference table is

Argument Entry 1st Difference 2nd Difference 3rd Difference


𝑥 𝑙𝑥 ∆𝑙𝑥 ∆ 2 𝑙𝑥 ∆ 3 𝑙𝑥

20 512
−73
30 439 −20
−93 10
40 346 −10
−103
50 243

Now applying Newton-Gregory Forward Interpolation Formula (Upto


third differences), we get
𝑢(𝑢−1) 𝑢(𝑢−1)(𝑢−2)
𝑓(𝑎 + ℎ𝑢) = 𝑓(𝑎) + 𝑢∆𝑓(𝑎) + ∆2 𝑓(𝑎)+ ∆3 𝑓(𝑎)
2! 3!
1.5(1.5−1)
∴ 𝑓(35) = 𝑙35 = 512 + 1.5 x (−73) + x (−20)
2!
1.5(1.5−1)(1.5−2)
+ x 10
3!

⇒ 𝑙35 = 512 −109.5 − 7.5 − 0.625


⇒ 𝑙35 = 394.4
Hence, 𝑙35 = 394, approximately

(ii) Secondly, we have to find out the value, when 𝑥 = 42 i.e., 𝑙42
Here, 𝑎 = 20, ℎ = 10 and 𝑥 = 𝑎 + ℎ𝑢 = 42
𝑎 + ℎ𝑢 = 42
 20 + 10𝑢 = 42
 10𝑢 = 42 − 20 = 22
 𝑢 = 2.2

The difference table is

Argument Entry 1st Difference 2nd Difference 3rd Difference


𝑥 𝑙𝑥 ∆𝑙𝑥 ∆ 2 𝑙𝑥 ∆ 3 𝑙𝑥

20 512
−73
30 439 −20
−93 10
40 346 −10
−103
50 243

Now applying Newton-Gregory Forward Interpolation Formula (Upto


third differences), we get
𝑢(𝑢−1) 𝑢(𝑢−1)(𝑢−2)
𝑓(𝑎 + ℎ𝑢) = 𝑓(𝑎) + 𝑢∆𝑓(𝑎) + ∆2 𝑓(𝑎)+ ∆3 𝑓(𝑎)
2! 3!
2.2(2.2−1)
∴ 𝑓(42) = 𝑙42 = 512 + 2.2 x (−73) + x (−20)
2!
2.2(2.2−1)(2.2−2)
+ x 10
3!

⇒ 𝑙42 = 512 −160.6 − 26.4 + 0.88 = 325.88


Hence, 𝑙42 = 326, approximately

(iii) Finally, we have to find out the value, when 𝑥 = 47 i.e., 𝑙47
Here, 𝑎 = 20, ℎ = 10 and 𝑥 = 𝑎 + ℎ𝑢 = 47
𝑎 + ℎ𝑢 = 47
 20 + 10𝑢 = 47
 10𝑢 = 47 − 20 = 27
 𝑢 = 2.7

The difference table is


Argument Entry 1st Difference 2nd Difference 3rd Difference
𝑥 𝑙𝑥 ∆𝑙𝑥 ∆2 𝑙𝑥 ∆3 𝑙𝑥

20 512
−73
30 439 −20
−93 10
40 346 −10
−103
50 243

Now applying Newton-Gregory Forward Interpolation Formula (Upto


third differences), we get
𝑢(𝑢−1) 𝑢(𝑢−1)(𝑢−2)
𝑓(𝑎 + ℎ𝑢) = 𝑓(𝑎) + 𝑢∆𝑓(𝑎) + ∆2 𝑓(𝑎)+ ∆3 𝑓(𝑎)
2! 3!
2.7(2.7−1)
∴ 𝑓(47) = 𝑙47 = 512 + 2.7 x (−73) + x (−20)
2!
2.7(2.7−1)(2.7−2)
+ x 10
3!
⇒ 𝑙47 = 512 −197.1 − 45.9 + 5.355
⇒ 𝑙47 = 274.355
Hence, 𝑙47 = 274, approximately

Example (3) : The values of 𝑓(𝑥) for 𝑥 = 0, 1, 2, … . . , 6 are given by


𝑥 : 0 1 2 3 4 5 6
𝑓(𝑥) : 5 7 10 15 25 49 99
Estimate the value of 𝑓(2.2) using only five of the given values.

Solution : Here, last five values of 𝑓(𝑥) for 𝑥 = 2, 3, 4, 5, 6 are taken into
consideration so that 2.2 occurs in the beginning of the table.
Here, 𝑎 = 2, ℎ = 1 and 𝑥 = 𝑎 + ℎ𝑢 = 2.2
𝑎 + ℎ𝑢 = 2.2
 2 + 𝑢 = 2.2
 𝑢 = 2.2 − 2 = 0.2
The difference table is

1st 2nd 3rd 4th


Argument Entry
Difference Difference Difference Difference
𝑥 𝑓(𝑥) 2 3
∆𝑓(𝑥) ∆ 𝑓(𝑥) ∆ 𝑓(𝑥) ∆4 𝑓(𝑥)

2 10
5
3 15 5
10 9
4 25 14 3
24 12
5 49 26
50
6 99
Now applying Newton-Gregory Forward Interpolation Formula (Upto
fourth differences), we get
𝑢(𝑢−1) 𝑢(𝑢−1)(𝑢−2)
𝑓(𝑎 + ℎ𝑢) = 𝑓(𝑎) + 𝑢∆𝑓(𝑎) + ∆2 𝑓(𝑎)+ ∆3 𝑓(𝑎)
2! 3!
𝑢(𝑢−1)(𝑢−2)(𝑢−3)
+ ∆4 𝑓(𝑎)
4!
0.2(0.2−1) 0.2(0.2−1)(0.2−2)
∴ 𝑓(2.2) = 10 + 0.2 x 5 + x5+ x9
2! 3!
0.2(0.2−1)(0.2−2)(0.2−3)
+ x3
4!
0.2(−0.8)(−1.8)
 𝑓(2.2) = 10 + 1−0.4 + x9
6
0.2(−0.8)(−1.8)(−2.8)
+ x3
24

 𝑓(2.2) = 11− 0.4 + 0.432 − 0.1008 = 10.9312

Hence, 𝑓(2.2) = 10.93 approximately

Example (4) : The following statistical data are given for 250 students
who secured the marks in Mathematics in Higher Secondary
Examination of a Junior College :

Marks Secured Number of Students

30…..40 20
40…..50 25
50…..60 75
60…..70 100
70…..80 30

How many students secured more than 45 marks ?

Solution : Given data are considered as


Marks Secured Number of Students

Less than 40 20
Less than 50 45
Less than 60 120
Less than 70 220
Less than 80 250

Here, 𝑎 = 40, ℎ = 10 and 𝑥 = 𝑎 + ℎ𝑢 = 55


𝑎 + ℎ𝑢 = 55
 40 + 10𝑢 = 55
 10𝑢 = 55 − 40 = 15
 𝑢 = 1.5
The difference table is

Marks Number
1st 2nd 3rd 4th
secured of
Difference Difference Difference Difference
less than students 2 3
𝑥 𝑓(𝑥) ∆𝑓(𝑥) ∆ 𝑓(𝑥) ∆ 𝑓(𝑥) ∆4 𝑓(𝑥)

40 20
25
50 45 50
75 −25
60 120 25 −70
100 −95
70 220 −70
30
80 250
Now applying Newton-Gregory Forward Interpolation Formula (Upto
fourth differences), we get

𝑢(𝑢−1) 𝑢(𝑢−1)(𝑢−2)
𝑓(𝑎 + ℎ𝑢) = 𝑓(𝑎) + 𝑢∆𝑓(𝑎) + ∆2 𝑓(𝑎)+ ∆3 𝑓(𝑎)
2! 3!
𝑢(𝑢−1)(𝑢−2)(𝑢−3)
+ ∆4 𝑓(𝑎)
4!
1.5(1.5−1) 1.5(1.5−1)(1.5−2)
∴ 𝑓(55) = 20 + 1.5 x 25 + x 50 + x (−25 )
2! 3!

1.5(1.5−1)(1.5−2)(1.5−3)
+ x (−70 )
4!

⇒ 𝑓(55) = 20 + 37.5 + 18.75 + 1.5625 – 1.64062


⇒ 𝑓(55) = 76.17

Hence, Number of students who secured less than 55 is 76.

You might also like