pub fn lambert_W0(x: f64) -> f64
This computes the principal branch of the Lambert W function, W_0(x).