Python Calculate log(1+x)/x for x near 0
<p>Is there a way to correctly calculate the value of log(1+x)/x in python for values of x close to 0? When I do it normally using np.log1p(x)/x, I get 1. I somehow seem to get the correct values when I use np.log(x). Isn't log1p supposed to be more stable?</p>

<p>So I found one answer to this. I used a library called decimal.</p>

<pre><code>from decimal import Decimal
x = Decimal('1e-13')
xp1 = Decimal(1) + x
print(xp1.ln()/x)
</code></pre>

<p>This library seems to be much more stable than numpy.</p>

python numerical-stability