Function log1p 

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

ln(1 + x), without the loss of digits near 0 that ln(1 + x) has. log1p(-1) is -∞.