Added support for solving with Broyden's method
Optimizations for jacobian calculations