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

        [Graphics:Images/BroydenMethodMod_gr_242.gif]   
        [Graphics:Images/BroydenMethodMod_gr_243.gif]  
        [Graphics:Images/BroydenMethodMod_gr_244.gif]  
    
Hint.  There are four solutions.  Good starting vectors are  [Graphics:Images/BroydenMethodMod_gr_245.gif].  

Solution 6.

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

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

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

 

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

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


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

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

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

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

 

 

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


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

 

 

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

 

 

 

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

 

 

[Graphics:../Images/BroydenMethodMod_gr_263.gif]
[Graphics:../Images/BroydenMethodMod_gr_264.gif]


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

 

 

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

 

 

Use the subroutine to get the answer.

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


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

Compare Broyden's result with Newton-Raphson's result.

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


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

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

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


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

Compare Broyden's result with Newton-Raphson's result.

 

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


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

(iii)  Use the Newton-Raphson method to find a numerical approximation to the solution near  [Graphics:../Images/BroydenMethodMod_gr_276.gif].  

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


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

Compare Broyden's result with Newton-Raphson's result.

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


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

(iv)  Use the Newton-Raphson method to find a numerical approximation to the solution near  [Graphics:../Images/BroydenMethodMod_gr_281.gif].  

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


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

Compare Broyden's result with Newton-Raphson's result.

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


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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005