Python log（1 + x）/ x for x near 0を計算する

Question

xの値が0に近い場合、pythonのlog（1 + x）/ xの値を正しく計算する方法はありますか？通常、np.log1p（x）/ xを使用すると、1になります。np.log（x）を使用すると、何とか正しい値を得ているようです。 log1pはより安定しているとは思われませんか？

Answer 1

私はこれに1つの答えを見つけました。私はdecimalというライブラリを使いました。

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

このライブラリはnumpyよりはるかに安定しているようです。

Answer 2

np.log1p(1+x) 

これはあなたにlog(2+x)を与えます。それをnp.log1p(x)に変更してください。