Example 2.  Use Newton's method to find the roots of the cubic polynomial  [Graphics:Images/SecantMethodMod_gr_120.gif].  
2 (a) Fast Convergence.  Investigate quadratic convergence at the simple root  [Graphics:Images/SecantMethodMod_gr_121.gif],  using the starting value  [Graphics:Images/SecantMethodMod_gr_122.gif]
2 (b) Slow Convergence.  Investigate linear convergence at the double root  [Graphics:Images/SecantMethodMod_gr_123.gif],  using the starting value  [Graphics:Images/SecantMethodMod_gr_124.gif]

Solution 2.

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

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

Graph the function.

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

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

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

The secant iteration formula  [Graphics:../Images/SecantMethodMod_gr_130.gif]  is

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

[Graphics:../Images/SecantMethodMod_gr_132.gif]
[Graphics:../Images/SecantMethodMod_gr_133.gif]

2 (a) Fast Convergence.  Investigate quadratic convergence at the simple root  [Graphics:../Images/SecantMethodMod_gr_134.gif],  using the starting value  [Graphics:../Images/SecantMethodMod_gr_135.gif]

First, do the iteration one step at a time.  
Type each of the following commands in a separate cell and execute them one at a time.

[Graphics:../Images/SecantMethodMod_gr_136.gif]
[Graphics:../Images/SecantMethodMod_gr_137.gif]

[Graphics:../Images/SecantMethodMod_gr_138.gif]
[Graphics:../Images/SecantMethodMod_gr_139.gif]

[Graphics:../Images/SecantMethodMod_gr_140.gif]
[Graphics:../Images/SecantMethodMod_gr_141.gif]

[Graphics:../Images/SecantMethodMod_gr_142.gif]
[Graphics:../Images/SecantMethodMod_gr_143.gif]

[Graphics:../Images/SecantMethodMod_gr_144.gif]
[Graphics:../Images/SecantMethodMod_gr_145.gif]

[Graphics:../Images/SecantMethodMod_gr_146.gif]
[Graphics:../Images/SecantMethodMod_gr_147.gif]

[Graphics:../Images/SecantMethodMod_gr_148.gif]
[Graphics:../Images/SecantMethodMod_gr_149.gif]

Notice that convergence is fast and the sequence is converging to the simple root  [Graphics:../Images/SecantMethodMod_gr_150.gif]  

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

[Graphics:../Images/SecantMethodMod_gr_152.gif]
[Graphics:../Images/SecantMethodMod_gr_153.gif]
[Graphics:../Images/SecantMethodMod_gr_154.gif]
[Graphics:../Images/SecantMethodMod_gr_155.gif]
[Graphics:../Images/SecantMethodMod_gr_156.gif]
[Graphics:../Images/SecantMethodMod_gr_157.gif]
[Graphics:../Images/SecantMethodMod_gr_158.gif]
[Graphics:../Images/SecantMethodMod_gr_159.gif]
[Graphics:../Images/SecantMethodMod_gr_160.gif]

[Graphics:../Images/SecantMethodMod_gr_161.gif]
[Graphics:../Images/SecantMethodMod_gr_162.gif]
[Graphics:../Images/SecantMethodMod_gr_163.gif]

At the simple root  [Graphics:../Images/SecantMethodMod_gr_164.gif], the order of convergence is known to be [Graphics:../Images/SecantMethodMod_gr_165.gif].  We can explore the ratio  [Graphics:../Images/SecantMethodMod_gr_166.gif]  for  k  sufficiently large.

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

k

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

 

2 (b) Slow Convergence.  Investigate linear convergence at the double root  [Graphics:../Images/SecantMethodMod_gr_202.gif],  using the starting value  [Graphics:../Images/SecantMethodMod_gr_203.gif]

First, do the iteration one step at a time.  
Type each of the following commands in a separate cell and execute them one at a time.

[Graphics:../Images/SecantMethodMod_gr_204.gif]
[Graphics:../Images/SecantMethodMod_gr_205.gif]

[Graphics:../Images/SecantMethodMod_gr_206.gif]
[Graphics:../Images/SecantMethodMod_gr_207.gif]

[Graphics:../Images/SecantMethodMod_gr_208.gif]
[Graphics:../Images/SecantMethodMod_gr_209.gif]

[Graphics:../Images/SecantMethodMod_gr_210.gif]
[Graphics:../Images/SecantMethodMod_gr_211.gif]

[Graphics:../Images/SecantMethodMod_gr_212.gif]
[Graphics:../Images/SecantMethodMod_gr_213.gif]

[Graphics:../Images/SecantMethodMod_gr_214.gif]
[Graphics:../Images/SecantMethodMod_gr_215.gif]

[Graphics:../Images/SecantMethodMod_gr_216.gif]
[Graphics:../Images/SecantMethodMod_gr_217.gif]

Notice that convergence is slow, but the sequence is converging to  the double root  [Graphics:../Images/SecantMethodMod_gr_218.gif]  

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

[Graphics:../Images/SecantMethodMod_gr_220.gif]
[Graphics:../Images/SecantMethodMod_gr_221.gif]
[Graphics:../Images/SecantMethodMod_gr_222.gif]
[Graphics:../Images/SecantMethodMod_gr_223.gif]
[Graphics:../Images/SecantMethodMod_gr_224.gif]
[Graphics:../Images/SecantMethodMod_gr_225.gif]
[Graphics:../Images/SecantMethodMod_gr_226.gif]
[Graphics:../Images/SecantMethodMod_gr_227.gif]
[Graphics:../Images/SecantMethodMod_gr_228.gif]
[Graphics:../Images/SecantMethodMod_gr_229.gif]
[Graphics:../Images/SecantMethodMod_gr_230.gif]
[Graphics:../Images/SecantMethodMod_gr_231.gif]
[Graphics:../Images/SecantMethodMod_gr_232.gif]
[Graphics:../Images/SecantMethodMod_gr_233.gif]
[Graphics:../Images/SecantMethodMod_gr_234.gif]
[Graphics:../Images/SecantMethodMod_gr_235.gif]
[Graphics:../Images/SecantMethodMod_gr_236.gif]
[Graphics:../Images/SecantMethodMod_gr_237.gif]
[Graphics:../Images/SecantMethodMod_gr_238.gif]
[Graphics:../Images/SecantMethodMod_gr_239.gif]
[Graphics:../Images/SecantMethodMod_gr_240.gif]
[Graphics:../Images/SecantMethodMod_gr_241.gif]
[Graphics:../Images/SecantMethodMod_gr_242.gif]
[Graphics:../Images/SecantMethodMod_gr_243.gif]
[Graphics:../Images/SecantMethodMod_gr_244.gif]
[Graphics:../Images/SecantMethodMod_gr_245.gif]
[Graphics:../Images/SecantMethodMod_gr_246.gif]
[Graphics:../Images/SecantMethodMod_gr_247.gif]
[Graphics:../Images/SecantMethodMod_gr_248.gif]
[Graphics:../Images/SecantMethodMod_gr_249.gif]
[Graphics:../Images/SecantMethodMod_gr_250.gif]

[Graphics:../Images/SecantMethodMod_gr_251.gif]
[Graphics:../Images/SecantMethodMod_gr_252.gif]
[Graphics:../Images/SecantMethodMod_gr_253.gif]

Compare our result with Mathematica's built in numerical root finder.

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

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

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

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

This can also be done with Mathematica's built in symbolic solve procedure.

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

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

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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004