Function richardson_derivative 

pub fn richardson_derivative<T, F>(
    f: F,
    x: T,
    h: T,
    order: usize,
) -> CoreResult<T>
where
    T: Float + StableComputation,
    F: Fn(T) -> T,

Stable numerical differentiation using Richardson extrapolation