Example 3.  Given  [Graphics:Images/NevilleAlgorithmMod_gr_291.gif]  and the nodes [Graphics:Images/NevilleAlgorithmMod_gr_292.gif] compute the entries in Aitken's tableau and Neville's tableau for evaluating  [Graphics:Images/NevilleAlgorithmMod_gr_293.gif]  at  [Graphics:Images/NevilleAlgorithmMod_gr_294.gif]  
Use the recursive subroutines.

Solution 3.

Enter the function and nodes and point of evaluation.  

[Graphics:../Images/NevilleAlgorithmMod_gr_295.gif]


[Graphics:../Images/NevilleAlgorithmMod_gr_296.gif]

 

 

3 (a).  Compute the values for Neville's tableau.

[Graphics:../Images/NevilleAlgorithmMod_gr_297.gif]

Neville's tableau is

[Graphics:../Images/NevilleAlgorithmMod_gr_298.gif]

 

 

    The accuracy of the computations in Neville's tableau  improves as you go down the table and to the right. Without knowing where [Graphics:../Images/NevilleAlgorithmMod_gr_299.gif] lies in relation to the nodes  [Graphics:../Images/NevilleAlgorithmMod_gr_300.gif] it is difficult to determine whether the diagonal approximations are converging faster or the bottom row.  For illustration, consider the sequence [Graphics:../Images/NevilleAlgorithmMod_gr_301.gif] on the diagonal.

[Graphics:../Images/NevilleAlgorithmMod_gr_302.gif]


[Graphics:../Images/NevilleAlgorithmMod_gr_303.gif]

 

 

The Neville polynomials  [Graphics:../Images/NevilleAlgorithmMod_gr_304.gif]  are converging to  [Graphics:../Images/NevilleAlgorithmMod_gr_305.gif].

[Graphics:../Images/NevilleAlgorithmMod_gr_306.gif]


[Graphics:../Images/NevilleAlgorithmMod_gr_307.gif]

[Graphics:../Images/NevilleAlgorithmMod_gr_308.gif]

 

 

 

3 (b).  Compute the values for Aitken's tableau.

[Graphics:../Images/NevilleAlgorithmMod_gr_309.gif]


[Graphics:../Images/NevilleAlgorithmMod_gr_310.gif]

 

 

    The accuracy of the computations in Aitken's tableau  improves as you go down the table and to the right. Without knowing where [Graphics:../Images/NevilleAlgorithmMod_gr_311.gif] lies in relation to the nodes  [Graphics:../Images/NevilleAlgorithmMod_gr_312.gif] it is difficult to determine whether the diagonal approximations are converging faster or the bottom row.  For illustration, consider the sequence [Graphics:../Images/NevilleAlgorithmMod_gr_313.gif] on the diagonal.

[Graphics:../Images/NevilleAlgorithmMod_gr_314.gif]


[Graphics:../Images/NevilleAlgorithmMod_gr_315.gif]

The Aitken polynomials  [Graphics:../Images/NevilleAlgorithmMod_gr_316.gif]  are converging to  [Graphics:../Images/NevilleAlgorithmMod_gr_317.gif].

 

 

[Graphics:../Images/NevilleAlgorithmMod_gr_318.gif]


[Graphics:../Images/NevilleAlgorithmMod_gr_319.gif]

[Graphics:../Images/NevilleAlgorithmMod_gr_320.gif]

 

 

Observation.  

It is easy verify that the entries on the diagonal are the same  
[Graphics:../Images/NevilleAlgorithmMod_gr_321.gif], [Graphics:../Images/NevilleAlgorithmMod_gr_322.gif], [Graphics:../Images/NevilleAlgorithmMod_gr_323.gif],[Graphics:../Images/NevilleAlgorithmMod_gr_324.gif], and [Graphics:../Images/NevilleAlgorithmMod_gr_325.gif].

[Graphics:../Images/NevilleAlgorithmMod_gr_326.gif]


[Graphics:../Images/NevilleAlgorithmMod_gr_327.gif]

 

 

Observation.  

It is easy verify that the entries on the diagonal are the same  [Graphics:../Images/NevilleAlgorithmMod_gr_328.gif], [Graphics:../Images/NevilleAlgorithmMod_gr_329.gif], [Graphics:../Images/NevilleAlgorithmMod_gr_330.gif],[Graphics:../Images/NevilleAlgorithmMod_gr_331.gif], [Graphics:../Images/NevilleAlgorithmMod_gr_332.gif], and [Graphics:../Images/NevilleAlgorithmMod_gr_333.gif].

[Graphics:../Images/NevilleAlgorithmMod_gr_334.gif]


[Graphics:../Images/NevilleAlgorithmMod_gr_335.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005