pub fn inverse_langevin(y: &f64) -> f64