Example 4.  Show how Newton's method is in reality a fixed point iteration scheme.  Use the system    
        [Graphics:Images/NewtonSystemMod_gr_288.gif]    

Solution 4.

First, enter the coordinate functions [Graphics:../Images/NewtonSystemMod_gr_289.gif] and construct the vector function [Graphics:../Images/NewtonSystemMod_gr_290.gif] using Mathematica, and then find the Jacobian matrix [Graphics:../Images/NewtonSystemMod_gr_291.gif].  

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


[Graphics:../Images/NewtonSystemMod_gr_293.gif]
[Graphics:../Images/NewtonSystemMod_gr_294.gif]

[Graphics:../Images/NewtonSystemMod_gr_295.gif]
[Graphics:../Images/NewtonSystemMod_gr_296.gif]
[Graphics:../Images/NewtonSystemMod_gr_297.gif]

Construct the Newton-Raphson fixed point function [Graphics:../Images/NewtonSystemMod_gr_298.gif].

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



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


Use fixed point iteration and the above fixed point function [Graphics:../Images/NewtonSystemMod_gr_301.gif] to find a numerical approximation to the solution near [Graphics:../Images/NewtonSystemMod_gr_302.gif].  

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



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

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

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

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

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

A fixed point satisfies the equation  [Graphics:../Images/NewtonSystemMod_gr_309.gif].  Our last approximation is stored in  [Graphics:../Images/NewtonSystemMod_gr_310.gif],  check it out.  

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


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

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

Accuracy is determined by the tolerance and number of iterations.  How accurate was the solution?

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


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

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

Compare these values with those obtained with the Newton-Raphson method for finding a numerical approximation to the solution near [Graphics:../Images/NewtonSystemMod_gr_317.gif].  

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



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

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

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

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

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

A solution to the system satisfies [Graphics:../Images/NewtonSystemMod_gr_324.gif].  Our last approximation is stored in [Graphics:../Images/NewtonSystemMod_gr_325.gif], check it out.  

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


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

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

Accuracy is determined by the tolerance and number of iterations.  How accurate was the solution?

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


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

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

This Newton-Raphson solution is the same as the one obtained using fixed point iteration.

Use fixed point iteration and the above fixed point function [Graphics:../Images/NewtonSystemMod_gr_332.gif] to find a numerical approximation to the solution near [Graphics:../Images/NewtonSystemMod_gr_333.gif].  

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



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

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

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

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

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

A fixed point satisfies the equation  [Graphics:../Images/NewtonSystemMod_gr_340.gif].  Our last approximation is stored in  [Graphics:../Images/NewtonSystemMod_gr_341.gif],  check it out.  

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


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

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

Compare these values with those obtained with the Newton-Raphson method for finding a numerical approximation to the solution near [Graphics:../Images/NewtonSystemMod_gr_345.gif].  

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



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

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

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

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

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

A solution to the system satisfies [Graphics:../Images/NewtonSystemMod_gr_352.gif].  Our last approximation is stored in [Graphics:../Images/NewtonSystemMod_gr_353.gif], check it out.  

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


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

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

This Newton-Raphson solution is the same as the one obtained using fixed point iteration.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004