數(shù)學(xué)家在計算任何東西的近似值時,都會給出一個區(qū)間,像左手和右手一樣,把心愛的人擁抱在懷中。如果只給一側(cè)的值,那么就會像楊過一樣,也是頂尖的高手。
阿基米德不是楊過。所以,他計算出圓周率小于 22/7 以后,他還要計算,圓周率大于多少。劉徽也不是楊過,他計算出圓周率大于3.141024以后,還要計算圓周率小于多少。盈朒二數(shù)都給出,才算完美。
阿基米德計算朒數(shù)的時候,也采用計算內(nèi)接多邊形邊長的方法。因為,內(nèi)接多邊形,周長一定小于圓的周長。

只是要計算邊長。如圖,AC是圓的一條直徑,O是圓心。AD是內(nèi)接N邊形的一條邊,AE是內(nèi)接2N邊形的一條邊。這一次,計算邊長,阿基米德使用的是直徑所在的三角形,而非半徑。直徑對的圓周角是直角,這一點(diǎn)可以帶來便捷。邊AD對的圓周角為ACD,圓周角是圓心角的一半。EC平分角ACD。
知道這些關(guān)系,按照前面的方法。在三角形ADC中,使用角平分線的性質(zhì),可以算出三角形MDC三邊的比例。而三角形AEC同三角形MDC是相似的,這一點(diǎn)很重要。也就是說,對應(yīng)邊的比例是一樣的。
在得到AEC三邊的比例以后,繼續(xù)平分角ACE,繼續(xù)計算比例。注意,阿基米德計算的是比例,不是長度。所以有這樣便捷的方法。
這次計算的起點(diǎn)是角ACD為30度,AD還是六邊形的一條邊。始于3的算術(shù)平方根小于 1351/780 。也就是DC/DA小于1351/780。同樣,使用計算半角余切和余割的方法,一直計算最終得到直徑與內(nèi)接正96邊形一邊的比例。該值小于(2017又1/4)比66,化成小數(shù)是 30.5643939。

在計算內(nèi)接多邊形的時候,更關(guān)注余割值。
圓周率大于(96 × 66) :2017.25 = 6366 比 2017.25,最后這個比例,阿基米德把它縮小一點(diǎn)點(diǎn),寫作 223/71 ,也就是 3又 10/71。
到此,阿基米德的工作完成了,他證明了圓周率大于 223/71,小于22/7 。在此基礎(chǔ)上,再接再厲,他給出了圓同其外切正方形面積之比近似值 11/14。 但在他的著作中,先給出了圓與其 外切正方形的面積比為 11/14 這個結(jié)論,明顯,他用的是西方人慣用的手法,先說重要的事情。11/14為何重要呢?因為,從11/14也可以逆推出圓周率,設(shè)圓的半徑為R,其面積為 PI*R*R ;外切正方形面積為 (2R)*(2R)= 4*R*R,兩者的比例是 PI:4,也就是說 PI/4 = 11/14,從而,PI =44/14,約分以后,正好就是 22/7 。可見,應(yīng)用中,阿基米德更喜歡這個分?jǐn)?shù)。
理論上講,阿基米德可以使用精度更高的3的平方根開始計算,平分到更小的角度,從而得到更加精密的數(shù)值。但那個年代,22/7 的精度用在生活中足夠了,而且很方便使用。每當(dāng)覺得自己做的不夠好的時候,想想“美國人,在1897年準(zhǔn)備立法,用3.2做圓周率”這件事。笑一笑,就過去了。
劉徽在得到 3.141024以后,也不滿足于僅僅使用3.14,他指出,這個數(shù)字還是偏小,“圓率猶為微小”。那么,要找到一個比圓率大一點(diǎn)點(diǎn)的數(shù)字,才行。怎樣尋找呢?依然從面積入手。

如圖,設(shè)BC是正N邊形的一條邊,BE是正2N邊形的一條邊。圖中,綠色的是圓弧CEB。正2N邊形面積比正N邊形大,大多少呢?在這個局部看,大出的部分是紫色三角形BCE。那么,整體上看,大了N個這樣的三角形。
假如把紫色三角形復(fù)制一份,切成兩半,貼在2N邊形的外面,那么,2N邊形的面積加上這份復(fù)制,就可以比圓大一點(diǎn)點(diǎn)。
不得不嘆服,這個處理太有技巧了!那么,好,192邊形的面積計算過了,是3.141024。
96邊形的面積呢?用48邊形邊長計算,

正48邊形邊長是 0.130806,那么正96邊形面積是 48*0.130806/2 = 3.139344
192邊形比96變形大的數(shù)值是 3.141024 - 3.139344 =0.00168
把這個數(shù)值加到192邊形面積上,就可以比圓大了。
3.141024+0.00168 = 3.142704
到此,劉徽得到圓周率的范圍是 3.141024 至 3.142704。劉徽的工作已經(jīng)圓滿完成了。盈朒二數(shù)都有了。
從這個范圍可以發(fā)現(xiàn),3.14已經(jīng)確定了,無論計算多少邊形,都不會再變了,再變,也是3.14后面的部分。3.14比“周三徑一”已經(jīng)是很偉大的突破了。因此劉徽對 3.14作為圓周率近似值,已經(jīng)比較滿意了。
但還不是很滿意。他做了更深遠(yuǎn)的計算。一方面,驗證割圓術(shù)的正確性;一方面,驗證3.14不會再改變;一方面,為了獲得更高的精度。這就是,為什么他的文章中前前后后的引用,很多人不知道他在做什么。
這部分計算,分成兩個部分。一是,直接使用前面的數(shù)據(jù),在正192邊形的面積上直接做“消息”,就是根據(jù)192邊形的面積,以及前面各多邊形的面積,對圓面積做更精確的估計;另一方面,繼續(xù)增加邊的數(shù)量,一直計算到正1536邊形的邊長,得到正3072邊形的面積。通過消息,獲得估計值 3.1416,通過增加邊數(shù),他驗證了自己的估值,并且說“若此者,蓋盡其纖微矣”,他猜測,3.1416可以把圓最細(xì)小的部分包括進(jìn)來,然后,通過計算到3072邊形,發(fā)現(xiàn)面積仍然沒有超出3.1416,于是,他對新的估值 3.1416 更滿意。

計算多邊形的方法是一樣的。那么,如何估值呢?從最初的幾個數(shù)字,可以發(fā)現(xiàn),隨著邊數(shù)量的增加,面積增加的幅度越來越小,每次增量大約是前一次的1/4。192邊形面積是3.141024,比96邊形大0.00168。那么可以估計,384邊形會比192邊形大 0.00168/4,再往下一級,大0.00168/16...
于是,極限的狀態(tài)應(yīng)該是 3.141024 + 0.00168(1/4+ 1/16 +1/64 + 1/256...) = 3.141024 + 0.00168/3
=? 3.1415842,畢竟,后面的增長比0.25稍大,因此,對3.14之后兩位進(jìn)行“入”的操作是合理的,估值就可以是 3.1416。這是我對劉徽估值方法的猜測。因為古人喜歡用比例,如果看到每次增長都變?yōu)樯弦淮蔚脑鲩L 1/4這樣的事實,不可能不顧及。
計算,估值,驗證。再計算,再估值,再驗證。在中國,在劉徽之前,沒有人計算過圓周率。所以,劉徽會如此謹(jǐn)慎的前行。
后來的祖沖之有福了,依然用這個方法,向下計算,直接獲得7位小數(shù)??上У氖牵潜窘凶觥毒Y術(shù)》的書在哪里?期待著有一天,考古工作者從古墓里挖出一本來。