fn real_roots(a: f64, b: f64, c: f64, d: f64) -> Vec<f64>
The real roots of a x³ + b x² + c x + d, degenerating to the quadratic and linear cases when the leading coefficients vanish.
a x³ + b x² + c x + d