fn depressed_cubic_dominant(g: f64, h: f64) -> f64
Dominant root of depressed cubic x^3 + gx + h = 0.
Section 2.2 of Orellana and De Michele.