pub fn richardson_derivative<T, F>( f: F, x: T, h: T, order: usize, ) -> Result<T, CoreError>where T: Float + StableComputation, F: Fn(T) -> T,
Stable numerical differentiation using Richardson extrapolation