计算π的一组公式

计算π的一组公式

由tan(π/4-arctan1/k)=(k-1)/(k+1)可以得到：
π/4=arctan1/k+arctan(k-1)/(k+1)
同理，由tan(π/4-2arctan1/k)=（k^2-2k-1）/(k^2+2k-1)可得：
π/4=2arctan1/k+arctan（k^2-2k-1）/(k^2+2k-1)
由tan(π/4-3arctan1/k)=（k^3-3k^2-3k+1）/（k^3+3k^2-3k-1）可得：
π/4=3arctan1/k+arctan（k^3-3k^2-3k+1）/（k^3+3k^2-3k-1）
由tan(π/4-4arctan1/k)=（k^4-4k^3-6k^2+4k+1）/（k^4+4k^3-6k^2-4k+1）可得：
π/4=3arctan1/k+arctan（k^4-4k^3-6k^2+4k+1）/（k^4+4k^3-6k^2-4k+1）
…………
根据tan(narctan1/k)的展开式，将tan(π/4-narctan1/k)展开，即可得到计算的一组公式。对于给定的

n,k值，可以得到下面的公式表(部分结果)

π/4 n=1 2 3 4
k=2 arctan1/2+arctan1/3; 2arctan1/2-arctan1/7; 3arctan1/2-arctan9/13; 4arctan1/2-arctan31/37
3 arctan1/3+arctan1/2; 2arctan1/3+arctan1/7; 3arctan1/3+arctan2/11; 4arctan1/3-arctan17/31
4 arctan1/4+arctan3/5; 2arctan1/4+arctan7/23; 3arctan1/4+arctan5/99; 4arctan1/4-arctan79/401
5 arctan1/5+arctan2/3; 2arctan1/5+arctan7/17; 3arctan1/5+arctan9/46; 4arctan1/5-arctan1/239
6 arctan1/6+arctan5/7; 2arctan1/6+arctan23/47; 3arctan1/6+arctan91/305
7 arctan1/7+arctan3/4; 2arctan1/7+arctan17/31  
8 arctan1/8+arctan7/9; 2arctan1/8+arctan47/79  
9 arctan1/9+arctan4/5; 2arctan1/9+arctan31/49  
10 arctan1/10+arctan9/11 2arctan1/10+arctan79/119  

其中，n=4,k=5的公式

π/4 =4arctan1/5-arctan1/239

为Machin公式，一般采用Machin公式进行π的计算。
显然，当tan（narctan1/k）与1最接近时，所得到的计算公式最有效。对于任意给定的n，tan

（narctan1/k）与1最接近，仅当tan（narctan1/x）=1且k=[x]或k=[x]+1，这时，arctan1/x=π/4n，x=

1/(tanπ/4n)。下面的表给出了对于不同的n值，x的取值情况，及tan（narctan1/[x]）和tan

（narctan1/（[x]+1））的值。

n 2 3 4 5 6 7 8 9 10
x 2.4 3.7 5.027 6.3 7.6 8.9 10.1 11.4 12.7
tan（narctan1/（[x]+1）） 0.75 0.9 0.78 0.87 0.92 0.98 0.89 0.93 0.97
tan（narctan1/[x]） 1.33 1.44 1.01 1.08 1.14 1.19 1.02 1.06 1.10

从表中可以看出，当n≤10时，Machin公式最有效。
arctan（x）利用下面的公式进行：
arctan（x）= x-1/3*x^3+1/5*x^5-...+(-1)^(n-1)*(1/(2n-1))*x^(2n-1)+...
在上式中，令x=1，可以得到Libreiz公式，但对于计算π是无效的。
π/4= 1-1/3+1/5-...+(-1)^(n-1)*(1/(2n-1))+...
~~~~~~~~~~~~~~~~
问题：将π展开成arctan1/2^n的一个有穷或无穷和是否可能？
π/4=arctan1/2+arctan1/3
arctan1/3-arctan1/4=arctan1/13
arctan1/13-arctan/16=arctan3/209
arctan3/209-arctan1/128=arctan175/26752
从而
π/4=arctan1/2+arctan1/4+arctan/16+arctan1/128+arctan175/26752
问题是arctan175/26752继续分拆下去能否终止？如果能够终止，我们就找到了一个有效的计算π的办法

，否则，我们可以选择终止在某一项，估计其余项，将误差控制在适当的范围内。