for such problems, broyden gives a generalization of the one-dimensional newton's method, replacing the derivative with the jacobian . the jacobian matrix is determined iteratively, based on the secant equation in the finite-difference approximation: