We have h=xi−xi−1=0.3 and s=−1/3, therefore, make a table with first to fourth divided differences. For example, calculate the first value in the First divided differences column (where n=1 is the number of the divided difference:
d1=h⋅nf(1.3)−f(1.0)=0.3⋅10.62008−0.76519=−0.48371.
Consequently, for the second, third and fourth divided difference n=2,3,4.
The table will look like this:
The table is adapted from Burden & Faires (2004)
We will use the underlined numbers in the Stirling's formula. The Stirling’s formula to approximate f(1.5) with x0=1.6 is
f(1.5)==0.45540+(−31)(20.3)((−0.54895)+(−0.57861))++(−31)2(0.3)2(−0.04944)++21(−31)((−31)2−1)(0.3)3(0.06588+0.06807)++(−31)2((−31)2−1)(0.3)4(0.00183)==0.51182. References
R. Burden, J. Faires (2004). Numerical Analysis. 8 edition.
Comments