from sympy.core.symbol import symbols
from sympy.solvers.solveset import nonlinsolve
Since the example in the documentation is broken, how do I solve a non-linear system of equations numerically in SymPy?