Exercise 9.  Observe that the subroutine ImprovedBroyden involves vector functions and is not dependent on the dimension.
Use the subroutine ImprovedBroyden to solve the nonlinear system in 3D space:  

        [Graphics:Images/BroydenMethodMod_gr_683.gif]   
        [Graphics:Images/BroydenMethodMod_gr_684.gif]  
        [Graphics:Images/BroydenMethodMod_gr_685.gif]  
    
Hint.  There are four solutions.  Good starting vectors are  [Graphics:Images/BroydenMethodMod_gr_686.gif].  

Solution 9.

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

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

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

 

Second, graph the surfaces  [Graphics:../Images/BroydenMethodMod_gr_692.gif],  [Graphics:../Images/BroydenMethodMod_gr_693.gif]  and  [Graphics:../Images/BroydenMethodMod_gr_694.gif]  using Mathematica.  The points of intersection are the solutions we seek.  

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


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

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

(i)  Use the improved Broyden's method to find a numerical approximation to the solution near  [Graphics:../Images/BroydenMethodMod_gr_698.gif].  

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


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

Compare with the standard Broyden's method.

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


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

(ii)  Use the improved Broyden's method to find a numerical approximation to the solution near  [Graphics:../Images/BroydenMethodMod_gr_703.gif].  

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


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

Compare with the standard Broyden's method.

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


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

(iii)  Use the improved Broyden's method to find a numerical approximation to the solution near  [Graphics:../Images/BroydenMethodMod_gr_708.gif].  

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


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

Compare with the standard Broyden's method.

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


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

(iv)  Use the improved Broyden's method to find a numerical approximation to the solution near  [Graphics:../Images/BroydenMethodMod_gr_713.gif].  

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


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

Compare with the standard Broyden's method.

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


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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005