Exercise 7.  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 4D space:  

        [Graphics:Images/BroydenMethodMod_gr_286.gif]   
        [Graphics:Images/BroydenMethodMod_gr_287.gif]  
        [Graphics:Images/BroydenMethodMod_gr_288.gif]   
        [Graphics:Images/BroydenMethodMod_gr_289.gif]  
    
Hint.  There is one real solutions.  A good starting vectors is  [Graphics:Images/BroydenMethodMod_gr_290.gif] Solution 7.

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

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

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

 

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

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

 

 

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


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

 

 

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

 

 

 

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

 

 

[Graphics:../Images/BroydenMethodMod_gr_302.gif]
[Graphics:../Images/BroydenMethodMod_gr_303.gif]


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

 

 

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

 

 

Use the subroutine to get the answer.

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


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

Do you think that there might be other solutions ?  We can use Mathematica to explore the situation.

 

Exploration

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2005