from sympy.solvers.solveset import nonlinsolve
x, y, z = symbols('x, y, z', real=True)
Since the example in the documentation is broken, how do I solve a non-linear system of equations numerically in SymPy?