scipy.optimize.broyden1

scipy.optimize.broyden1(F, xin, iter=10, alpha=0.10000000000000001, verbose=False)

Broyden’s first method.

Updates Jacobian and computes inv(J) by a matrix inversion at every iteration. It’s very slow.

The best norm |F(x)|=0.005 achieved in ~45 iterations.