log1p(x) computes log(1+x), with precision more accurate than standard log() function for small value of x.