Solving nonlinear simultaneous equations with a generalization of Brent's method

We describe an implementation of a generalization of Brent's method for solving systems of nonlinear equations. Some important features of the algorithm, like step control, discretization of derivatives and stopping criteria, are discussed. In particular we give numerical experiences which show that a stopping criterion proposed by D. Gay is efficient.