In [2]:
# solve a classic quadratic equation, display the results
from sympy import *
init_printing(use_latex=True)
a,b,c,x  = symbols('a b c x')
solve(a*x**2+b*x+c,x)
Out[2]:
$$\begin{bmatrix}\frac{- b + \sqrt{- 4 a c + b^{2}}}{2 a}, & - \frac{b + \sqrt{- 4 a c + b^{2}}}{2 a}\end{bmatrix}$$