I would like to use sympy non linear solver (nsolve) for solving a set of equations. One of these equations contain a custom python function that call a subroutine in