一个有趣的计算圆周率方法.pdfVIP

數學傳播 33 卷 1期, pp. 57-64 一個有趣的計算圓周率方法 一一雙反正切項尤拉型公式 黃見利 本文謹獻給作者的恩師 Dar´ıo Castellanos 教授 (生於 1937年 12 月4 日; 逝於 1995年 11 月23 日)。 作者永遠懷念他的教誨。 一. 前言 在 1671 年, 蘇格蘭數學家 James Gregory 發現了下列的反正切型級數 ∞ n 2n+1 tan−1 x = (−1) x , −1 ≤ x ≤ 1. 2n + 1 n=0 當 x = 1 時, 另一位偉大的數學家, 也是微積分的創造者之一, Gottfried Wilhelm Leibniz π 1 1 1 1 1 在三年後的 1674 年發現了所謂的 Leibniz 級數 = 1 − + − + − + · · · “不 4 3 5 7 9 11 過, 一位反正切型級數的發現者, 而且還著手於圓周率超越性方面的研究, 會忽略了在他的級數 中替代 x = 1 這麼如此明顯的舉手之勞, 還真是令人不可思議的一件事。 較可能的情況似乎是, 他不認為替代 x = 1 在他的級數中是重要的, 因為它的收斂速度實在太慢了, 因而無法使用在 圓周率的計算上。”[1] 如果他真的如此考慮過, 那真是絕對正確且嚴謹周密的想法! 現在就讓我 們來瞧一瞧這條級數的收斂速度到底有多慢。 首先, 我們知道 Leibniz 級數是一種所謂的交替級數 (alternating series)。 因此, 我們可 以利用下列著名的定理來估計誤差: 交替級數估計定理 (Alternating Series Estimation Theorem) ∞ n 假設 s = (−1)n−1an 為一交替級數之和, 而 sn = (−1)k−1ak 為其部份項之和, n=1 k=1 且滿足下列兩條件: 57 58 數學傳播 33 卷 1期 民 98 年 3 月 (i) 0 ≤ an+1 ≤ an (ii) lim an = 0 n→∞ =⇒ |Rn | = |s − sn | ≤ an+1 . 由這個定理 , 我們得到了一個令人驚異的事實。 假設我們現在想得到德國數學家 Ludolph van Ceulen 在 1615年死後才發表的 35位圓周率小數 3.14159265358979323846264338327950288。 由於第 36 位小數是 4, 所以 |Rn | ≤ 25 ·10−37 1 不會影響到第 35 位小數的實際值。 所以 |Rn | ≤ 25 · 10−37 = 4 · 1035 , 亦即利用 Leibniz 級 4 4 4 4 4 4 數計算 π ≈ 4 − 3 + 5 − 7 + 9 − 11 + · · · + 16 · 1035 − 1 時, 絕對有 35 位圓周率小數 的正確值。 不過, 話說回來, 僅僅因為想求得 35 位圓周率小數位數, 我們利用 Leibniz 級數竟 然需計算 8 · 1035 項才能得到!!! 我們所要下的苦工比 Ludolph van Ceulen 當年好不到那裏 去, 甚至可能更糟! 因此, 在 1983 年之前, 幾乎所有圓