solve_poly_system([y**2 - x**3 + 1, y*x], x, y)