pub fn calculate_single_integral_simpson<E: Fn(f64) -> f64>( equation: E, first_integral_begin: f64, first_integral_end: f64, first_integral_step: f64, ) -> Result<f64, Error>