Function gradient 

pub fn gradient<F: Fn([f32; 3]) -> f32>(
    field: &F,
    p: [f32; 3],
    eps: f32,
) -> [f32; 3]