from sympy import * x = Symbol('x') y = Symbol('y') print(solve([3 * x * x* x- y*y*y - 3, x* x + y*y - 1],[x, y]))