pub fn bisect( v0: f64, v1: f64, expected: f64, tolerance: f64, f: impl Fn(f64) -> f64, ) -> Result<f64>