scipy.optimize.broyden3

scipy.optimize.broyden3(F, xin, iter=10, alpha=0.40000000000000002, verbose=False)

Broyden’s second method.

Updates inverse Jacobian by an optimal formula. The NxN matrix multiplication is avoided.

The best norm |F(x)|=0.003 achieved in ~20 iterations.

Recommended.