I want to solve (many times) a large (say 100x100) non-linear (cubic) system of algebraic equations.
I am using scipy.optimize.root non-linear solver to which I give
