pub fn divergence( fx: impl Fn(f64, f64, f64) -> f64, fy: impl Fn(f64, f64, f64) -> f64, fz: impl Fn(f64, f64, f64) -> f64, x: f64, y: f64, z: f64, h: f64, ) -> f64