pub fn differentiate_with_stencils<F: Fn(f64) -> f64 + Sync>( x: f64, f: F, d: u128, stencils: Vec<i128>, ) -> f64