pub fn log1p<F: Float>(x: F) -> F
ln(1 + x), without the loss of digits near 0 that ln(1 + x) has. log1p(-1) is -∞.
ln(1 + x)
log1p(-1)