Python Decimal类型精度错误

我有一个相当令人困惑的问题,我怀疑这与科学记数法和小数精度有关。以下是我的部分代码:

    def atan(x):
        # Calculate arctan(1/x)
        x = Decimal(x)
        current_value = Decimal(0)
        divisor = 1
        x_squared = x * x
        current_term = 1 / x


        while True:
            current_value += current_term

            divisor += 2
            current_term = (-current_term / x_squared) / divisor
            print(current_term)

            # The issue
            if current_term == Decimal(0):
                break

        return current_value

    print(atan(5))

这是基于公式 atan(1/x) = 1/x - 1/(3x^3) + 1/(5x^5) - ...

然而,我发现 current_term 每一次循环迭代都会变小,其值将变为4E-80000。自从我设置了小数精度 getcontext().prec 到20,当前术语甚至不应支持这些值。我想不知怎么的 当前_时间 不是十进制类型,而是科学记数法/浮点类型,但python告诉它仍然是十进制类型。

arctan(1/5)的正确值约为0.1973955。我得到的值是0.1973545,从第5位开始是错误的。即使我手动中断循环,由于某些原因,该值仍然是错误的。 如有任何解决此问题的帮助,我们将不胜感激。