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

Solution 2.

[Graphics:../Images/Newton'sMethodMod_gr_117.gif]

[Graphics:../Images/Newton'sMethodMod_gr_118.gif]

Graph the function.

[Graphics:../Images/Newton'sMethodMod_gr_119.gif]

[Graphics:../Images/Newton'sMethodMod_gr_120.gif]

[Graphics:../Images/Newton'sMethodMod_gr_121.gif]

The Newton-Raphson iteration formula  g[x]  is

[Graphics:../Images/Newton'sMethodMod_gr_122.gif]

[Graphics:../Images/Newton'sMethodMod_gr_123.gif]
[Graphics:../Images/Newton'sMethodMod_gr_124.gif]

2 (a) Fast Convergence.  Investigate quadratic convergence at the simple root  [Graphics:../Images/Newton'sMethodMod_gr_125.gif],  using the starting value  [Graphics:../Images/Newton'sMethodMod_gr_126.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/Newton'sMethodMod_gr_127.gif]
[Graphics:../Images/Newton'sMethodMod_gr_128.gif]

[Graphics:../Images/Newton'sMethodMod_gr_129.gif]
[Graphics:../Images/Newton'sMethodMod_gr_130.gif]

[Graphics:../Images/Newton'sMethodMod_gr_131.gif]
[Graphics:../Images/Newton'sMethodMod_gr_132.gif]

[Graphics:../Images/Newton'sMethodMod_gr_133.gif]
[Graphics:../Images/Newton'sMethodMod_gr_134.gif]

[Graphics:../Images/Newton'sMethodMod_gr_135.gif]
[Graphics:../Images/Newton'sMethodMod_gr_136.gif]

[Graphics:../Images/Newton'sMethodMod_gr_137.gif]
[Graphics:../Images/Newton'sMethodMod_gr_138.gif]

[Graphics:../Images/Newton'sMethodMod_gr_139.gif]
[Graphics:../Images/Newton'sMethodMod_gr_140.gif]

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

[Graphics:../Images/Newton'sMethodMod_gr_142.gif]

[Graphics:../Images/Newton'sMethodMod_gr_143.gif]
[Graphics:../Images/Newton'sMethodMod_gr_144.gif]
[Graphics:../Images/Newton'sMethodMod_gr_145.gif]
[Graphics:../Images/Newton'sMethodMod_gr_146.gif]
[Graphics:../Images/Newton'sMethodMod_gr_147.gif]
[Graphics:../Images/Newton'sMethodMod_gr_148.gif]
[Graphics:../Images/Newton'sMethodMod_gr_149.gif]
[Graphics:../Images/Newton'sMethodMod_gr_150.gif]
[Graphics:../Images/Newton'sMethodMod_gr_151.gif]
[Graphics:../Images/Newton'sMethodMod_gr_152.gif]
[Graphics:../Images/Newton'sMethodMod_gr_153.gif]

At the simple root  [Graphics:../Images/Newton'sMethodMod_gr_154.gif]  we can explore the relationship  [Graphics:../Images/Newton'sMethodMod_gr_155.gif]  for  k  sufficiently large.

This will be done by investigating the ratio  [Graphics:../Images/Newton'sMethodMod_gr_156.gif]  for  k  sufficiently large.

 

[Graphics:../Images/Newton'sMethodMod_gr_157.gif]

k

[Graphics:../Images/Newton'sMethodMod_gr_158.gif]

[Graphics:../Images/Newton'sMethodMod_gr_159.gif]

[Graphics:../Images/Newton'sMethodMod_gr_160.gif]

[Graphics:../Images/Newton'sMethodMod_gr_161.gif]

[Graphics:../Images/Newton'sMethodMod_gr_162.gif]

[Graphics:../Images/Newton'sMethodMod_gr_163.gif]

[Graphics:../Images/Newton'sMethodMod_gr_164.gif]

[Graphics:../Images/Newton'sMethodMod_gr_165.gif]

[Graphics:../Images/Newton'sMethodMod_gr_166.gif]

[Graphics:../Images/Newton'sMethodMod_gr_167.gif]

[Graphics:../Images/Newton'sMethodMod_gr_168.gif]

[Graphics:../Images/Newton'sMethodMod_gr_169.gif]

[Graphics:../Images/Newton'sMethodMod_gr_170.gif]

[Graphics:../Images/Newton'sMethodMod_gr_171.gif]

[Graphics:../Images/Newton'sMethodMod_gr_172.gif]

[Graphics:../Images/Newton'sMethodMod_gr_173.gif]

[Graphics:../Images/Newton'sMethodMod_gr_174.gif]

[Graphics:../Images/Newton'sMethodMod_gr_175.gif]

[Graphics:../Images/Newton'sMethodMod_gr_176.gif]

[Graphics:../Images/Newton'sMethodMod_gr_177.gif]

[Graphics:../Images/Newton'sMethodMod_gr_178.gif]

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

 

Evaluate the quantity  [Graphics:../Images/Newton'sMethodMod_gr_180.gif] at the root  [Graphics:../Images/Newton'sMethodMod_gr_181.gif].

 

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

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

Which is close to the computed value  [Graphics:../Images/Newton'sMethodMod_gr_184.gif]  

 

 

2 (b) Slow Convergence.  Investigate linear convergence at the double root  [Graphics:../Images/Newton'sMethodMod_gr_185.gif],  using the starting value  [Graphics:../Images/Newton'sMethodMod_gr_186.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/Newton'sMethodMod_gr_187.gif]
[Graphics:../Images/Newton'sMethodMod_gr_188.gif]

[Graphics:../Images/Newton'sMethodMod_gr_189.gif]
[Graphics:../Images/Newton'sMethodMod_gr_190.gif]

[Graphics:../Images/Newton'sMethodMod_gr_191.gif]
[Graphics:../Images/Newton'sMethodMod_gr_192.gif]

[Graphics:../Images/Newton'sMethodMod_gr_193.gif]
[Graphics:../Images/Newton'sMethodMod_gr_194.gif]

[Graphics:../Images/Newton'sMethodMod_gr_195.gif]
[Graphics:../Images/Newton'sMethodMod_gr_196.gif]

[Graphics:../Images/Newton'sMethodMod_gr_197.gif]
[Graphics:../Images/Newton'sMethodMod_gr_198.gif]

[Graphics:../Images/Newton'sMethodMod_gr_199.gif]
[Graphics:../Images/Newton'sMethodMod_gr_200.gif]

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

[Graphics:../Images/Newton'sMethodMod_gr_202.gif]

[Graphics:../Images/Newton'sMethodMod_gr_203.gif]
[Graphics:../Images/Newton'sMethodMod_gr_204.gif]
[Graphics:../Images/Newton'sMethodMod_gr_205.gif]
[Graphics:../Images/Newton'sMethodMod_gr_206.gif]
[Graphics:../Images/Newton'sMethodMod_gr_207.gif]
[Graphics:../Images/Newton'sMethodMod_gr_208.gif]
[Graphics:../Images/Newton'sMethodMod_gr_209.gif]
[Graphics:../Images/Newton'sMethodMod_gr_210.gif]
[Graphics:../Images/Newton'sMethodMod_gr_211.gif]
[Graphics:../Images/Newton'sMethodMod_gr_212.gif]
[Graphics:../Images/Newton'sMethodMod_gr_213.gif]
[Graphics:../Images/Newton'sMethodMod_gr_214.gif]
[Graphics:../Images/Newton'sMethodMod_gr_215.gif]
[Graphics:../Images/Newton'sMethodMod_gr_216.gif]
[Graphics:../Images/Newton'sMethodMod_gr_217.gif]
[Graphics:../Images/Newton'sMethodMod_gr_218.gif]
[Graphics:../Images/Newton'sMethodMod_gr_219.gif]
[Graphics:../Images/Newton'sMethodMod_gr_220.gif]
[Graphics:../Images/Newton'sMethodMod_gr_221.gif]
[Graphics:../Images/Newton'sMethodMod_gr_222.gif]
[Graphics:../Images/Newton'sMethodMod_gr_223.gif]
[Graphics:../Images/Newton'sMethodMod_gr_224.gif]
[Graphics:../Images/Newton'sMethodMod_gr_225.gif]
[Graphics:../Images/Newton'sMethodMod_gr_226.gif]
[Graphics:../Images/Newton'sMethodMod_gr_227.gif]
[Graphics:../Images/Newton'sMethodMod_gr_228.gif]
[Graphics:../Images/Newton'sMethodMod_gr_229.gif]
[Graphics:../Images/Newton'sMethodMod_gr_230.gif]
[Graphics:../Images/Newton'sMethodMod_gr_231.gif]

At the double root  [Graphics:../Images/Newton'sMethodMod_gr_232.gif]  we can explore the relationship  [Graphics:../Images/Newton'sMethodMod_gr_233.gif]  for  k  sufficiently large.

This will be done by investigating the ratio  [Graphics:../Images/Newton'sMethodMod_gr_234.gif]  for  k  sufficiently large.

 

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

k

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

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

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

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

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

[Graphics:../Images/Newton'sMethodMod_gr_241.gif]

[Graphics:../Images/Newton'sMethodMod_gr_242.gif]

[Graphics:../Images/Newton'sMethodMod_gr_243.gif]

[Graphics:../Images/Newton'sMethodMod_gr_244.gif]

[Graphics:../Images/Newton'sMethodMod_gr_245.gif]

[Graphics:../Images/Newton'sMethodMod_gr_246.gif]

[Graphics:../Images/Newton'sMethodMod_gr_247.gif]

[Graphics:../Images/Newton'sMethodMod_gr_248.gif]

[Graphics:../Images/Newton'sMethodMod_gr_249.gif]

[Graphics:../Images/Newton'sMethodMod_gr_250.gif]

[Graphics:../Images/Newton'sMethodMod_gr_251.gif]

[Graphics:../Images/Newton'sMethodMod_gr_252.gif]

[Graphics:../Images/Newton'sMethodMod_gr_253.gif]

[Graphics:../Images/Newton'sMethodMod_gr_254.gif]

[Graphics:../Images/Newton'sMethodMod_gr_255.gif]

[Graphics:../Images/Newton'sMethodMod_gr_256.gif]

[Graphics:../Images/Newton'sMethodMod_gr_257.gif]

[Graphics:../Images/Newton'sMethodMod_gr_258.gif]

[Graphics:../Images/Newton'sMethodMod_gr_259.gif]

[Graphics:../Images/Newton'sMethodMod_gr_260.gif]

[Graphics:../Images/Newton'sMethodMod_gr_261.gif]

[Graphics:../Images/Newton'sMethodMod_gr_262.gif]

[Graphics:../Images/Newton'sMethodMod_gr_263.gif]

[Graphics:../Images/Newton'sMethodMod_gr_264.gif]

[Graphics:../Images/Newton'sMethodMod_gr_265.gif]

[Graphics:../Images/Newton'sMethodMod_gr_266.gif]

[Graphics:../Images/Newton'sMethodMod_gr_267.gif]

[Graphics:../Images/Newton'sMethodMod_gr_268.gif]

[Graphics:../Images/Newton'sMethodMod_gr_269.gif]

[Graphics:../Images/Newton'sMethodMod_gr_270.gif]

[Graphics:../Images/Newton'sMethodMod_gr_271.gif]

[Graphics:../Images/Newton'sMethodMod_gr_272.gif]

[Graphics:../Images/Newton'sMethodMod_gr_273.gif]

[Graphics:../Images/Newton'sMethodMod_gr_274.gif]

[Graphics:../Images/Newton'sMethodMod_gr_275.gif]

[Graphics:../Images/Newton'sMethodMod_gr_276.gif]

[Graphics:../Images/Newton'sMethodMod_gr_277.gif]

[Graphics:../Images/Newton'sMethodMod_gr_278.gif]

[Graphics:../Images/Newton'sMethodMod_gr_279.gif]

[Graphics:../Images/Newton'sMethodMod_gr_280.gif]

[Graphics:../Images/Newton'sMethodMod_gr_281.gif]

[Graphics:../Images/Newton'sMethodMod_gr_282.gif]

[Graphics:../Images/Newton'sMethodMod_gr_283.gif]

[Graphics:../Images/Newton'sMethodMod_gr_284.gif]

[Graphics:../Images/Newton'sMethodMod_gr_285.gif]

[Graphics:../Images/Newton'sMethodMod_gr_286.gif]

[Graphics:../Images/Newton'sMethodMod_gr_287.gif]

[Graphics:../Images/Newton'sMethodMod_gr_288.gif]

[Graphics:../Images/Newton'sMethodMod_gr_289.gif]

[Graphics:../Images/Newton'sMethodMod_gr_290.gif]

[Graphics:../Images/Newton'sMethodMod_gr_291.gif]

[Graphics:../Images/Newton'sMethodMod_gr_292.gif]

[Graphics:../Images/Newton'sMethodMod_gr_293.gif]

[Graphics:../Images/Newton'sMethodMod_gr_294.gif]

[Graphics:../Images/Newton'sMethodMod_gr_295.gif]

[Graphics:../Images/Newton'sMethodMod_gr_296.gif]

[Graphics:../Images/Newton'sMethodMod_gr_297.gif]

[Graphics:../Images/Newton'sMethodMod_gr_298.gif]

[Graphics:../Images/Newton'sMethodMod_gr_299.gif]

[Graphics:../Images/Newton'sMethodMod_gr_300.gif]

[Graphics:../Images/Newton'sMethodMod_gr_301.gif]

[Graphics:../Images/Newton'sMethodMod_gr_302.gif]

[Graphics:../Images/Newton'sMethodMod_gr_303.gif]

[Graphics:../Images/Newton'sMethodMod_gr_304.gif]

[Graphics:../Images/Newton'sMethodMod_gr_305.gif]

[Graphics:../Images/Newton'sMethodMod_gr_306.gif]

[Graphics:../Images/Newton'sMethodMod_gr_307.gif]

[Graphics:../Images/Newton'sMethodMod_gr_308.gif]

[Graphics:../Images/Newton'sMethodMod_gr_309.gif]

[Graphics:../Images/Newton'sMethodMod_gr_310.gif]

[Graphics:../Images/Newton'sMethodMod_gr_311.gif]

[Graphics:../Images/Newton'sMethodMod_gr_312.gif]

[Graphics:../Images/Newton'sMethodMod_gr_313.gif]

[Graphics:../Images/Newton'sMethodMod_gr_314.gif]

[Graphics:../Images/Newton'sMethodMod_gr_315.gif]

[Graphics:../Images/Newton'sMethodMod_gr_316.gif]

[Graphics:../Images/Newton'sMethodMod_gr_317.gif]

[Graphics:../Images/Newton'sMethodMod_gr_318.gif]

[Graphics:../Images/Newton'sMethodMod_gr_319.gif]

[Graphics:../Images/Newton'sMethodMod_gr_320.gif]

[Graphics:../Images/Newton'sMethodMod_gr_321.gif]

[Graphics:../Images/Newton'sMethodMod_gr_322.gif]

[Graphics:../Images/Newton'sMethodMod_gr_323.gif]

[Graphics:../Images/Newton'sMethodMod_gr_324.gif]

[Graphics:../Images/Newton'sMethodMod_gr_325.gif]

[Graphics:../Images/Newton'sMethodMod_gr_326.gif]

[Graphics:../Images/Newton'sMethodMod_gr_327.gif]

[Graphics:../Images/Newton'sMethodMod_gr_328.gif]

[Graphics:../Images/Newton'sMethodMod_gr_329.gif]

[Graphics:../Images/Newton'sMethodMod_gr_330.gif]

[Graphics:../Images/Newton'sMethodMod_gr_331.gif]

[Graphics:../Images/Newton'sMethodMod_gr_332.gif]

[Graphics:../Images/Newton'sMethodMod_gr_333.gif]

[Graphics:../Images/Newton'sMethodMod_gr_334.gif]

[Graphics:../Images/Newton'sMethodMod_gr_335.gif]

[Graphics:../Images/Newton'sMethodMod_gr_336.gif]

[Graphics:../Images/Newton'sMethodMod_gr_337.gif]

 

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

 

[Graphics:../Images/Newton'sMethodMod_gr_338.gif]

[Graphics:../Images/Newton'sMethodMod_gr_339.gif]

[Graphics:../Images/Newton'sMethodMod_gr_340.gif]

[Graphics:../Images/Newton'sMethodMod_gr_341.gif]

 

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

 

[Graphics:../Images/Newton'sMethodMod_gr_342.gif]

[Graphics:../Images/Newton'sMethodMod_gr_343.gif]

[Graphics:../Images/Newton'sMethodMod_gr_344.gif]

[Graphics:../Images/Newton'sMethodMod_gr_345.gif]

[Graphics:../Images/Newton'sMethodMod_gr_346.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004