pub fn find_root_or_minimum<T: MyFloat>( a: &T, b: &T, f: impl Fn(&T) -> T, epsilon: Fpt, ) -> DualResult<T>