New quasi-Newton method for solving systems of nonlinear equations

We propose a new Broyden method for solving systems of nonlinear equations, which uses the first derivatives, but is more efficient than the Newton method (measured by the computational time) for larger dense systems. The new method updates QR or LU decompositions of nonsymmetric approximations of the Jacobian matrix, so it requires O(n ²) arithmetic operations per iteration in contrast with the Newton method, which requires O(n ³) operations per iteration. Computational experiments confirm the high efficiency of the new method.