如何提高此python代码中的计算精度

我尝试了100次迭代,并且工作正常:

x = 0
y = 0
vecx = []
vecy = []

for i in range(100):
    new_x = 1.3 + 0.3*x + 0.6*x*y - 0.6*y - y**2
    new_y = 0.1 - 0.7*x + 0.5*(x**2) - 0.8*x*y + 0.1*y - 0.6*(y**2)
    
    x, y = new_x, new_y  # Update x and y for the next iteration        

    vecx.append(x)
    vecy.append(y)

    print('x', x)
    print('y', y)

    

您可以使用 decimal Python模块:

from decimal import Decimal, getcontext

# Set the precision for decimal calculations
getcontext().prec = 50  # You can adjust the precision as needed

x = Decimal('0')
y = Decimal('0')
vecx = []
vecy = []

for i in range(5000):
    x_new = Decimal('1.3') + Decimal('0.3') * x + Decimal('0.6') * x * y - Decimal('0.6') * y - y**2
    y_new = Decimal('0.1') - Decimal('0.7') * x + Decimal('0.5') * (x**2) - Decimal('0.8') * x * y + Decimal('0.1') * y - Decimal('0.6') * (y**2)
    x, y = x_new, y_new
    vecx.append(x)
    vecy.append(y)

您应该尝试使用十进制库进行精确的浮点数字计算。 from decimal import Decimal 并以字符串形式提供数字 y = Decimal('1.005')