Module

for

The Newton Polynomial

     

Background.

    We have seen how to expand a function  [Graphics:Images/NewtonPolyMod_gr_1.gif] in a Maclaurin polynomial about [Graphics:Images/NewtonPolyMod_gr_2.gif] involving the powers [Graphics:Images/NewtonPolyMod_gr_3.gif] and a Taylor polynomial about [Graphics:Images/NewtonPolyMod_gr_4.gif] involving the powers [Graphics:Images/NewtonPolyMod_gr_5.gif].  These polynomials have a single "center" [Graphics:Images/NewtonPolyMod_gr_6.gif].  Polynomial interpolation can be used to construct the polynomial of  degree  [Graphics:Images/NewtonPolyMod_gr_7.gif]  that passes through the n+1 points  [Graphics:Images/NewtonPolyMod_gr_8.gif],  for  [Graphics:Images/NewtonPolyMod_gr_9.gif].  If multiple "centers"    [Graphics:Images/NewtonPolyMod_gr_10.gif]  are used, then the result is the so called Newton polynomial.  We attribute much of the founding theory to Sir Isaac Newton (1643-1727).

 

Theorem (Newton Polynomial).  Assume that  [Graphics:Images/NewtonPolyMod_gr_11.gif] and  [Graphics:Images/NewtonPolyMod_gr_12.gif] for  [Graphics:Images/NewtonPolyMod_gr_13.gif]  are distinct  values.  Then

    [Graphics:Images/NewtonPolyMod_gr_14.gif],
    
where [Graphics:Images/NewtonPolyMod_gr_15.gif] is a polynomial that can be used to approximate  [Graphics:Images/NewtonPolyMod_gr_16.gif],

    [Graphics:Images/NewtonPolyMod_gr_17.gif]  

and we write  

    [Graphics:Images/NewtonPolyMod_gr_18.gif].

The Newton polynomial goes through the [Graphics:Images/NewtonPolyMod_gr_19.gif] points  [Graphics:Images/NewtonPolyMod_gr_20.gif],  i.e.

    [Graphics:Images/NewtonPolyMod_gr_21.gif]    for   [Graphics:Images/NewtonPolyMod_gr_22.gif].  

The remainder term  [Graphics:Images/NewtonPolyMod_gr_23.gif] has the form

    [Graphics:Images/NewtonPolyMod_gr_24.gif],

for some value [Graphics:Images/NewtonPolyMod_gr_25.gif] that lies in the interval [Graphics:Images/NewtonPolyMod_gr_26.gif].  The coefficients  [Graphics:Images/NewtonPolyMod_gr_27.gif]  are constructed using divided differences.

Proof  Newton Polynomials  Newton Polynomials  

 

Definition. Divided Differences.

The divided differences for a function [Graphics:Images/NewtonPolyMod_gr_28.gif] are defined as follows:

    [Graphics:Images/NewtonPolyMod_gr_29.gif]  

The divided difference formulae are used to construct the divided difference table:

[Graphics:Images/NewtonPolyMod_gr_30.gif]

[Graphics:Images/NewtonPolyMod_gr_31.gif]

[Graphics:Images/NewtonPolyMod_gr_42.gif]

[Graphics:Images/NewtonPolyMod_gr_43.gif]

[Graphics:Images/NewtonPolyMod_gr_51.gif]

[Graphics:Images/NewtonPolyMod_gr_52.gif]

[Graphics:Images/NewtonPolyMod_gr_32.gif]

[Graphics:Images/NewtonPolyMod_gr_33.gif]

[Graphics:Images/NewtonPolyMod_gr_34.gif]

[Graphics:Images/NewtonPolyMod_gr_35.gif]

[Graphics:Images/NewtonPolyMod_gr_36.gif]

[Graphics:Images/NewtonPolyMod_gr_37.gif]

[Graphics:Images/NewtonPolyMod_gr_38.gif]

[Graphics:Images/NewtonPolyMod_gr_39.gif]

[Graphics:Images/NewtonPolyMod_gr_40.gif]

[Graphics:Images/NewtonPolyMod_gr_41.gif]

  

[Graphics:Images/NewtonPolyMod_gr_44.gif]

[Graphics:Images/NewtonPolyMod_gr_45.gif]

[Graphics:Images/NewtonPolyMod_gr_46.gif]

[Graphics:Images/NewtonPolyMod_gr_47.gif]

  

   

[Graphics:Images/NewtonPolyMod_gr_48.gif]

[Graphics:Images/NewtonPolyMod_gr_49.gif]

[Graphics:Images/NewtonPolyMod_gr_50.gif]

  

   

   

[Graphics:Images/NewtonPolyMod_gr_53.gif]

[Graphics:Images/NewtonPolyMod_gr_54.gif]

  

   

   

   

[Graphics:Images/NewtonPolyMod_gr_55.gif]

The coefficient [Graphics:Images/NewtonPolyMod_gr_56.gif] of the Newton polynomial  [Graphics:Images/NewtonPolyMod_gr_57.gif] is  [Graphics:Images/NewtonPolyMod_gr_58.gif]  and it is the top element in the column of the i-th divided differences.

The Newton polynomial of  degree  [Graphics:Images/NewtonPolyMod_gr_59.gif]  that passes through the  [Graphics:Images/NewtonPolyMod_gr_60.gif] points  [Graphics:Images/NewtonPolyMod_gr_61.gif],  for  [Graphics:Images/NewtonPolyMod_gr_62.gif]  is  

    [Graphics:Images/NewtonPolyMod_gr_63.gif]

Proof  Newton Polynomials  Newton Polynomials  

 

 

    The cubic curve in the figure below illustrates a Newton polynomial of degree n = 3.  
    
    

[Graphics:Images/NewtonPolyMod_gr_64.gif]

[Graphics:Images/NewtonPolyMod_gr_65.gif]


[Graphics:Images/NewtonPolyMod_gr_66.gif]


[Graphics:Images/NewtonPolyMod_gr_67.gif]


[Graphics:Images/NewtonPolyMod_gr_68.gif]


[Graphics:Images/NewtonPolyMod_gr_69.gif]
[Graphics:Images/NewtonPolyMod_gr_70.gif]


[Graphics:Images/NewtonPolyMod_gr_71.gif]
[Graphics:Images/NewtonPolyMod_gr_72.gif]


Theorem.  (Error Bounds for Newton Interpolation, Equally Spaced Nodes)  Assume that  [Graphics:Images/NewtonPolyMod_gr_73.gif]  defined on [Graphics:Images/NewtonPolyMod_gr_74.gif],  which contains the equally spaced nodes  [Graphics:Images/NewtonPolyMod_gr_75.gif].  Additionally, assume that    [Graphics:Images/NewtonPolyMod_gr_76.gif]  and the derivatives of  [Graphics:Images/NewtonPolyMod_gr_77.gif]  up to the order  [Graphics:Images/NewtonPolyMod_gr_78.gif]  are continuous and bounded on the special subintervals  [Graphics:Images/NewtonPolyMod_gr_79.gif], [Graphics:Images/NewtonPolyMod_gr_80.gif], [Graphics:Images/NewtonPolyMod_gr_81.gif], [Graphics:Images/NewtonPolyMod_gr_82.gif], and [Graphics:Images/NewtonPolyMod_gr_83.gif], respectively;  that is,

    [Graphics:Images/NewtonPolyMod_gr_84.gif],  

for  [Graphics:Images/NewtonPolyMod_gr_85.gif].  The error terms corresponding to these three cases have the following useful bounds on their magnitude  

(i).    [Graphics:Images/NewtonPolyMod_gr_86.gif][Graphics:Images/NewtonPolyMod_gr_87.gif]   is valid for  [Graphics:Images/NewtonPolyMod_gr_88.gif],  

(ii).    [Graphics:Images/NewtonPolyMod_gr_89.gif][Graphics:Images/NewtonPolyMod_gr_90.gif]   is valid for  [Graphics:Images/NewtonPolyMod_gr_91.gif],  

(iii).    [Graphics:Images/NewtonPolyMod_gr_92.gif][Graphics:Images/NewtonPolyMod_gr_93.gif]   is valid for  [Graphics:Images/NewtonPolyMod_gr_94.gif],  

(iv).    [Graphics:Images/NewtonPolyMod_gr_95.gif][Graphics:Images/NewtonPolyMod_gr_96.gif]   is valid for  [Graphics:Images/NewtonPolyMod_gr_97.gif],  

(v).    [Graphics:Images/NewtonPolyMod_gr_98.gif][Graphics:Images/NewtonPolyMod_gr_99.gif]   is valid for  [Graphics:Images/NewtonPolyMod_gr_100.gif].  

Proof  Newton Polynomials  Newton Polynomials  

 

Animations (Newton Polynomials  Newton Polynomials).  Internet hyperlinks to animations.

 

Computer Programs  Newton Polynomials  Newton Polynomials  

 

Algorithm (Newton Interpolation Polynomial).   To construct and evaluate the Newton polynomial of  degree  [Graphics:Images/NewtonPolyMod_gr_101.gif]  that passes through the n+1 points  [Graphics:Images/NewtonPolyMod_gr_102.gif],  for  [Graphics:Images/NewtonPolyMod_gr_103.gif]     

    [Graphics:Images/NewtonPolyMod_gr_104.gif]  

where  
    [Graphics:Images/NewtonPolyMod_gr_105.gif]  
and  
    [Graphics:Images/NewtonPolyMod_gr_106.gif]  

Remark 1.  Newton polynomials are created "recursively."

    [Graphics:Images/NewtonPolyMod_gr_107.gif].

Remark 2.  Mathematica's arrays are lists and the subscripts must start with  1  instead of  0.

 

Mathematica Subroutine (Newton Polynomial).

[Graphics:Images/NewtonPolyMod_gr_108.gif]


Example 1.   Form the Newton polynomials of degree  n = 1,2, 3, 4, and 5  for the function  [Graphics:Images/NewtonPolyMod_gr_109.gif]  over the interval  [Graphics:Images/NewtonPolyMod_gr_110.gif]  using equally spaced nodes selected from the following list  
[Graphics:Images/NewtonPolyMod_gr_111.gif]  
Solution 1 (a).
Solution 1 (b).
Solution 1 (c).
Solution 1 (d).
Solution 1 (e).
Solution 1 (f).
Solution 1 (g).

 

Example 2.  Error Analysis.  Investigate the error for the Newton polynomial approximations in Example 1.
Solution 2 (a).
Solution 2 (b).
Solution 2 (c).
Solution 2 (d).
Solution 2 (e).

 

Example 3.  Summary of the maximum error and the error bounds over the intervals  [Graphics:Images/NewtonPolyMod_gr_342.gif]  for the Newton polynomials in Example 2.
3 (a).  Find  [Graphics:Images/NewtonPolyMod_gr_343.gif].  
3 (b).  Find  [Graphics:Images/NewtonPolyMod_gr_344.gif].  
3 (c).  Find  [Graphics:Images/NewtonPolyMod_gr_345.gif].  
3 (d).  Find  [Graphics:Images/NewtonPolyMod_gr_346.gif].  
3 (e).  Find  [Graphics:Images/NewtonPolyMod_gr_347.gif].  
Solution 3.

 

Example 4.  Application to number theory.
4 (a).  Find the formula for the sum of the first  n  integers.
4 (b).  Find the formula for the sum of the squares of the first  n integers.  
Solution 4.

 

Various Scenarios and Animations for Newton Polynomials.

[Graphics:Images/NewtonPolyMod_gr_403.gif]

Example 5.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_404.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_405.gif], [Graphics:Images/NewtonPolyMod_gr_406.gif], and [Graphics:Images/NewtonPolyMod_gr_407.gif].  
Solution 5 (a).
Solution 5 (b).
Solution 5 (c).

 

Example 6.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_435.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_436.gif].  
Solution 6.

 

Example 7.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_446.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_447.gif].  
Solution 7.

 

Example 8.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_457.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_458.gif], and [Graphics:Images/NewtonPolyMod_gr_459.gif].  
Solution 8 (a).
Solution 8 (b).

 

Example 9.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_478.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_479.gif].  
Solution 9.

 

Example 10.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_491.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_492.gif].  
Solution 10.

 

Example 11.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_502.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_503.gif].  
Solution 11.

 

Example 12.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_513.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_514.gif].  
Solution 12.

 

Example 13.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_524.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_525.gif].  
Solution 13.

 

Example 14.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_535.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_536.gif].  
Solution 14.

 

Example 15.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_550.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_551.gif].  
Solution 15.

 

Example 16.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_565.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_566.gif].  
Solution 16.

 

Example 17.  Find the Newton polynomial approximation for  [Graphics:Images/NewtonPolyMod_gr_580.gif], on the interval [Graphics:Images/NewtonPolyMod_gr_581.gif].  
Solution 17.

 

Animations (Newton Polynomial Approximation  Newton Polynomial Approximation).  Internet hyperlinks to animations.

 

Old Lab Project (Newton Interpolation Polynomial  Newton Interpolation Polynomial).  Internet hyperlinks to an old lab project.  

 

Research Experience for Undergraduates

Newton Polynomials  Newton Polynomials  Internet hyperlinks to web sites and a bibliography of articles.  

 

Download this Mathematica Notebook Newton Polynomials

 

Return to Numerical Methods - Numerical Analysis

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004