Implementation of Broyden's root finding function to numerically compute the root of a system of nonlinear equations
Broyden(system, vars, x = 0, tol = .Machine$double.eps^0.4, maxiters = 10000)
A list of functions
A character string list of variables that appear in the functions
A starting vector
The tolerance for the function specifying how precise it will be
maximum number of iterations.