Broyden’s second method.
Updates inverse Jacobian by an optimal formula. There is NxN matrix multiplication in every iteration.
The best norm |F(x)|=0.003 achieved in ~20 iterations.
Recommended.
scipy.optimize.broyden1
scipy.optimize.broyden3