Example 1.  Use the value  [Graphics:Images/NewtonImprovedMod_gr_124.gif]  and compare Methods A,B and C for finding the double root  [Graphics:Images/NewtonImprovedMod_gr_125.gif]  of the equation  [Graphics:Images/NewtonImprovedMod_gr_126.gif].

Solution 1.

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

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

 

 

First, we will compute the iterations for each method, and afterward a table comparing the methods is given.  For the method C, all the iterations in the linear search are included.  

Using formula (2), the standard Newton-Raphson method.  

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

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

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

 

 

Using formula (9) Method A, the accelerated Newton-Raphson method with  m=2.  

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

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

 

 

Using formula (11) Method B, the modified Newton-Raphson method.

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

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

 

 

Using formula (13) Method C, the adaptive Newton-Raphson method.

The details for obtaining  [Graphics:../Images/NewtonImprovedMod_gr_136.gif] are:

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

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

 

 

Since  [Graphics:../Images/NewtonImprovedMod_gr_139.gif]  we reject  [Graphics:../Images/NewtonImprovedMod_gr_140.gif]  and set  [Graphics:../Images/NewtonImprovedMod_gr_141.gif]  and continue the iteration using formula (13).  The subroutine makes all these choices automatically.
        

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

  

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

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

 

 

Summary of the above results.

 

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

 

    Observe in this example that the standard Newton-Raphson method converges linearly and methods A and B converge quadratically.  The reader can use formulas (12) to verify that [Graphics:../Images/NewtonImprovedMod_gr_146.gif] is the order of the root [Graphics:../Images/NewtonImprovedMod_gr_147.gif].  The new method C is almost as fast as methods A and B.

 

Reminder.  The goal of this investigation is to show how the adaptive Newton-Raphson method is superior to the standard Newton-Raphson method, because of the limitations of Methods A and B.

 

Caveat.  Why is it difficult to locate a multiple root.  Because the function values themselves are essentially "noise" when you get close to a multiple root.  

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


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

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


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


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

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

One way to tell is to graph their difference.  

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


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

 

 

 

The formula  [Graphics:../Images/NewtonImprovedMod_gr_156.gif]  should be considered the "true value."
So one should be suspect of the computation  
[Graphics:../Images/NewtonImprovedMod_gr_157.gif]  for values of  x  is near  x=1.

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


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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004