Write `quadratic(a, b, c)` returning a sorted tuple of two real roots of a*x^2 + b*x + c. Assume two real roots exist. Return as a list `[smaller, larger]`.