pub fn forward_diff_const<const N: usize, F>( x: &[F; N], f: &'_ dyn Fn(&[F; N]) -> Result<F, Error>, ) -> Result<[F; N], Error>where F: Float + FromPrimitive,