Если вы предпочитаете sympy, вы можете использовать nsolve .
>>> nsolve([x+y**2-4, exp(x)+x*y-3], [x, y], [1, 1])
[0.620344523485226]
[1.83838393066159]
Первый аргумент - это список уравнений, второй - список переменных, а третий - начальный Угадай.