a, b, c = var('a b c')
eq = a * x^2 + b * x + c
s = solve(eq, x)
show(s)

x1 = s[0].rhs(); show(x1)

x2 = s[1].rhs(); show(x2)

eq.subs(x = x1)

_.expand()

eq.subs(x = x2).expand()

s2 = solve(x^3 + 4 * x^2 - 3 * x + 1, x); show(s2)

latex(s2)

n(s2[0].rhs()), n(s2[1].rhs()), n(s2[2].rhs())

show(n(s2[1].rhs()))

y = var('y')
solve([x + y == 3, x - y == 1], [x, y])

solve([x + y == 2, 2 * x + 2 * y == 4], [x, y])

solve([x + y == 2, 2 * x + 2 * y == 5], [x, y])

