Function lambertw 

pub fn lambertw<F>(x: F) -> F
where
    F: Float,

W0(x), the principal branch of the Lambert W function: the w of at least -1 with w e^w = x. A real x below -1/e gives NaN, except within half an ulp of -1/e, where W0 is -1.