Example 5.  Given  [Graphics:Images/NevilleAlgorithmMod_gr_386.gif]  and the nodes [Graphics:Images/NevilleAlgorithmMod_gr_387.gif] compute the entries in Aitken's tableau and Neville's tableau for evaluating  [Graphics:Images/NevilleAlgorithmMod_gr_388.gif]  at  [Graphics:Images/NevilleAlgorithmMod_gr_389.gif]  
Use the rearrangement of the nodes [Graphics:Images/NevilleAlgorithmMod_gr_390.gif] and the recursive subroutines.

Solution 5.

Enter the function and nodes and point of evaluation.  

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


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

 

 

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

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

Neville's tableau is

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

 

 

Remark.  Notice the ordering of the nodes in the second column.

    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_395.gif] lies in relation to the nodes  [Graphics:../Images/NevilleAlgorithmMod_gr_396.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_397.gif] on the diagonal.

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


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

 

 

Note.  The convergence is slightly faster than in Example 3 (a) where the nodes were chosen in ascending order  [Graphics:../Images/NevilleAlgorithmMod_gr_400.gif]  and the  entries on the diagonal were:  

 

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

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

The Neville polynomials  [Graphics:../Images/NevilleAlgorithmMod_gr_403.gif]  are converging to  [Graphics:../Images/NevilleAlgorithmMod_gr_404.gif].

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


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

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

 

 

 

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

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


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

 

 

Remark.  Notice the ordering of the nodes in the second column.

    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_410.gif] lies in relation to the nodes  [Graphics:../Images/NevilleAlgorithmMod_gr_411.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_412.gif] on the diagonal.

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


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

 

 

Note.  The convergence is slightly faster than in Example 3 (b) where the nodes were chosen in ascending order  [Graphics:../Images/NevilleAlgorithmMod_gr_415.gif]  and the  entries on the diagonal were:  

 

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

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

The Aitken polynomials  [Graphics:../Images/NevilleAlgorithmMod_gr_418.gif]  are converging to  [Graphics:../Images/NevilleAlgorithmMod_gr_419.gif].

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


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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005