pub fn calculate_bounds_emg( mu: f64, sigma: f64, lambda: f64, step_size: f64, target: f64, lower_start: f64, upper_start: f64, n_steps: Option<usize>, ) -> (f64, f64)