def test_solve_triangualted():
    f_1 = x**2 + y + z - 1
    f_2 = x + y**2 + 1 - z
    f_3 = x + y + z**2 - 1

    a, b = sqrt(2) - 1, -sqrt(2) - 1

    assert solve_triangulated([f_1, f_2, f_3], x, y, z) == \
        [(0, 0, 1), (0, 1, 0), (1, 0, 0)]

    dom = QQ.algebraic_field(sqrt(2))

    assert solve_triangulated([f_1, f_2, f_3], x, y, z, domain=dom) == \
        [(0, 0, 1), (0, 1, 0), (1, 0, 0), (a, a, a), (b, b, b)]
dataset/ETHPy150Open sympy/sympy/sympy/solvers/tests/test_polysys.py/test_solve_triangualted/swapped operands of `-`