J阶梅欣公式及其算法说课.doc

J阶梅钦公式及其简式的算法分析与应用 吴敏金 在诸多的反正切和式应用于计算圆周率pi中，梅钦公式[1]占有特别重要的位置。本文将在作者[2]的基础上，进一步将梅钦公式加以拓广，提出j阶梅钦公式及其简式，它包含了日本人在1981年将圆周率pi计算到2*10^6位所采用的公式（即1阶梅钦公式），导出其实现的算法，并给出基于万进制数的计算圆周率pi的应用结果。 反正切和式的回顾 所谓反正切和式[2]为关系式 arctan 1/a=arctan 1/b+arctan 1/c……………(1) 其中a,b,c为正实数。在应用于计算圆周率pi时，要求a,b,c取正整数（或有理分数）。简记为 [a=b+c]………………………………………………(2) 可以证明[2]： a=(bc-1)/(b+c)……………………………………(3) 1+a^2=(b-a)(c-a)…………………………………(4) 反正切和式有多种形式，如： [a=b-c]表示arctan 1/a=arctan 1/b-arctan 1/c，此时， a=(bc+1)/(c-b)……………………………………(5) 。又如,[a=n*b]表示arctan 1/a=n*arctan 1/b。（注意：在[ ]中，n,b不可交换。符号“*”不是算术乘“· ”，是n个反正切的和）。 例如,梅钦公式:[1=4*5-239]，而日本人用的则是：