Example 2.  Consider the function  [Graphics:Images/Halley'sMethodMod_gr_175.gif], which has a root at  [Graphics:Images/Halley'sMethodMod_gr_176.gif].   
2 (a).  Use the Newton-Raphson formula to find the root.    Use the starting value  [Graphics:Images/Halley'sMethodMod_gr_177.gif]
2 (b).  Use Halley's formula to find the root.    Use the starting value  [Graphics:Images/Halley'sMethodMod_gr_178.gif]

Solution 2 (a).

[Graphics:../Images/Halley'sMethodMod_gr_179.gif]


[Graphics:../Images/Halley'sMethodMod_gr_180.gif]

Form the Newton-Raphson iteration function  g(x).

[Graphics:../Images/Halley'sMethodMod_gr_181.gif]



[Graphics:../Images/Halley'sMethodMod_gr_182.gif]

[Graphics:../Images/Halley'sMethodMod_gr_183.gif]

[Graphics:../Images/Halley'sMethodMod_gr_184.gif]

[Graphics:../Images/Halley'sMethodMod_gr_185.gif]


We start the iteration with  [Graphics:../Images/Halley'sMethodMod_gr_186.gif]  and carry 100 digits in the computations, by telling Mathematica the precision of [Graphics:../Images/Halley'sMethodMod_gr_187.gif] by issuing the command p[0] = N[[Graphics:../Images/Halley'sMethodMod_gr_188.gif],100].  Next, a short program is written to compute the first seven terms in the iteration:

[Graphics:../Images/Halley'sMethodMod_gr_189.gif]



[Graphics:../Images/Halley'sMethodMod_gr_190.gif]
[Graphics:../Images/Halley'sMethodMod_gr_191.gif]
[Graphics:../Images/Halley'sMethodMod_gr_192.gif]
[Graphics:../Images/Halley'sMethodMod_gr_193.gif]
[Graphics:../Images/Halley'sMethodMod_gr_194.gif]
[Graphics:../Images/Halley'sMethodMod_gr_195.gif]
[Graphics:../Images/Halley'sMethodMod_gr_196.gif]
[Graphics:../Images/Halley'sMethodMod_gr_197.gif]


Since the root is known to be exactly   [Graphics:../Images/Halley'sMethodMod_gr_198.gif]  we can have Mathematica list the error  [Graphics:../Images/Halley'sMethodMod_gr_199.gif]  at each step in the iteration:

[Graphics:../Images/Halley'sMethodMod_gr_200.gif]



[Graphics:../Images/Halley'sMethodMod_gr_201.gif]
[Graphics:../Images/Halley'sMethodMod_gr_202.gif]
[Graphics:../Images/Halley'sMethodMod_gr_203.gif]
[Graphics:../Images/Halley'sMethodMod_gr_204.gif]
[Graphics:../Images/Halley'sMethodMod_gr_205.gif]
[Graphics:../Images/Halley'sMethodMod_gr_206.gif]
[Graphics:../Images/Halley'sMethodMod_gr_207.gif]


Looking at the error, we see that the number of accurate digits is doubling at each step in the computations, hence convergence is proceeding quadratically.

Verify the convergence rate.  At the simple root  [Graphics:../Images/Halley'sMethodMod_gr_208.gif]  we can explore the ratio [Graphics:../Images/Halley'sMethodMod_gr_209.gif].

[Graphics:../Images/Halley'sMethodMod_gr_210.gif]


k

[Graphics:../Images/Halley'sMethodMod_gr_211.gif]

[Graphics:../Images/Halley'sMethodMod_gr_212.gif]

[Graphics:../Images/Halley'sMethodMod_gr_213.gif]

[Graphics:../Images/Halley'sMethodMod_gr_214.gif]

[Graphics:../Images/Halley'sMethodMod_gr_215.gif]

[Graphics:../Images/Halley'sMethodMod_gr_216.gif]

[Graphics:../Images/Halley'sMethodMod_gr_217.gif]

[Graphics:../Images/Halley'sMethodMod_gr_218.gif]

[Graphics:../Images/Halley'sMethodMod_gr_219.gif]

[Graphics:../Images/Halley'sMethodMod_gr_220.gif]

[Graphics:../Images/Halley'sMethodMod_gr_221.gif]

[Graphics:../Images/Halley'sMethodMod_gr_222.gif]

[Graphics:../Images/Halley'sMethodMod_gr_223.gif]

[Graphics:../Images/Halley'sMethodMod_gr_224.gif]

[Graphics:../Images/Halley'sMethodMod_gr_225.gif]

[Graphics:../Images/Halley'sMethodMod_gr_226.gif]

[Graphics:../Images/Halley'sMethodMod_gr_227.gif]

[Graphics:../Images/Halley'sMethodMod_gr_228.gif]

[Graphics:../Images/Halley'sMethodMod_gr_229.gif]

[Graphics:../Images/Halley'sMethodMod_gr_230.gif]

[Graphics:../Images/Halley'sMethodMod_gr_231.gif]

[Graphics:../Images/Halley'sMethodMod_gr_232.gif]

[Graphics:../Images/Halley'sMethodMod_gr_233.gif]

[Graphics:../Images/Halley'sMethodMod_gr_234.gif]

[Graphics:../Images/Halley'sMethodMod_gr_235.gif]

[Graphics:../Images/Halley'sMethodMod_gr_236.gif]

[Graphics:../Images/Halley'sMethodMod_gr_237.gif]

[Graphics:../Images/Halley'sMethodMod_gr_238.gif]

[Graphics:../Images/Halley'sMethodMod_gr_239.gif]

[Graphics:../Images/Halley'sMethodMod_gr_240.gif]

Therefore, the Newton-Raphson iteration is converging quadratically.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004