pub fn solve_quadratic(a: Real, b: Real, c: Real) -> (Real, Real)