More Background

    A drawback of Newton's method is the requirement that the Jacobian matrix be computed, which requires the calculation of    partial derivatives per iteration.  One obvious improvement is to use a finite difference approximation for the Jacobian matrix.  For two dimensions this is

        
where h is small.  Notice that this will require    function evaluations.  

Exploration

(c) John H. Mathews 2005