pub fn integrate_2d_symmetric( f: &dyn Fn(&f64, &f64) -> f64, x_min: &f64, x_max: &f64, num_points: u128, ) -> f64