Example 3.  Use the Pseudo-Newton's method to solve the nonlinear system  

        [Graphics:Images/BroydenMethodMod_gr_89.gif]     Solution 3.

First, enter the coordinate functions [Graphics:../Images/BroydenMethodMod_gr_90.gif] and construct the vector function  [Graphics:../Images/BroydenMethodMod_gr_91.gif]  using Mathematica, and then find the Jacobian matrix [Graphics:../Images/BroydenMethodMod_gr_92.gif],  and approximate Jacobian [Graphics:../Images/BroydenMethodMod_gr_93.gif].  

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

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

 

Second, graph the curves [Graphics:../Images/BroydenMethodMod_gr_96.gif] and [Graphics:../Images/BroydenMethodMod_gr_97.gif] using Mathematica.  The points of intersection are the solutions we seek.  

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


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

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

(i)  Use the Pseudo-Newton method to find a numerical approximation to the solution near  [Graphics:../Images/BroydenMethodMod_gr_101.gif].  

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


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

Compare with the Newton-Raphson method.  

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


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

Observation.  The results for the subroutines NewtonSystem and PseudoNewtonSystem are very close.  Obviously, the Jacobian and its approximation agree to several significant digits.

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


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

(ii)  Use the Pseudo-Newton method to find a numerical approximation to the solution near  [Graphics:../Images/BroydenMethodMod_gr_108.gif].  

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


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

Compare with the Newton-Raphson method.  

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


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

Observation.  The results for the subroutines NewtonSystem and PseudoNewtonSystem are very close.  Obviously, the Jacobian and its approximation agree to several significant digits.

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


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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005