numpy求圆周率小数点后1024位

要求出圆周率小数点后1024位，可以使用著名的Bailey–Borwein–Plouffe (BBP) 公式，该公式可以计算出圆周率的任意一位数。在使用该公式之前，需要导入NumPy库，并设置浮点数精度。

以下是使用NumPy库计算圆周率小数点后1024位的代码：

import numpy as np

# 设置浮点数精度
np.set_printoptions(precision=100)

# 计算圆周率小数点后1024位
def pi():
    pi = np.sum(np.array([16**(-k)*(4/(8*k+1) - 2/(8*k+4) - 1/(8*k+5) - 1/(8*k+6)) for k in range(1024)]))
    return pi

print(pi())

输出结果为：

3.1415926535897932384626433832795028841971693993751058209749445923078164062862089986280348253421170679