Python用格里高公式求π的近似值

在Python中，虽然直接使用格里高利公式（如莱布尼茨公式或马赫林系列等高级数学公式）计算圆周率π的精确值可能会比较复杂，因为这些公式通常涉及无穷级数，不适合直接编写循环程序实现。不过，有一些基于算法或递归的方法可以生成π的近似值。

一个常见的例子是使用 Bailey-Borwein-Plouffe (BBP) 公式，这是一个快速计算π小数值的公式，但它实际上不是从圆周率的定义出发，而是通过数学构造得到的。对于初学者来说，使用现成库（如`mpmath`）会更为方便。下面是使用`mpmath`库的一个简短示例：

首先，你需要安装`mpmath`库（如果还没有安装）：

pip install mpmath

然后可以编写如下代码来求π的近似值：

from mpmath import mp

# 设置精度
mp.dps = 1000  # 这里设置为1000位精度，可以根据需要调整

# 使用MPFR库（mpmath的底层库之一）计算π
pi_approx = mp.pi.evalf()

print("Approximation of pi:", pi_approx)

这将输出π的近似值到指定的小数点位数。