先把误差定死,再去切八亿边形
他要的精度是:拿这个 π 去算一个直径六十万倍地球的圆,周长错不过一根马鬃。倒推下来,多边形得切到 3 × 2²⁸ 边。
阿基米德用正 96 边形把 π 夹在 223/71 与 22/7 之间,此后一千六百多年,凡是靠多边形算 π 的人走的都是同一条路:边数翻倍,每翻一倍就要开一次平方(喀拉拉的马德哈瓦是例外,他走的是无穷级数)。1424 年 7 月的《论圆周》把这条路走到了尽头,而新的东西不在算法里,在次序上——他先写下要多准,再回头算这需要多少条边。他给的标准是一句可以验收的话:取一个直径为六十万个地球直径的圆,算出来的周长与真值相差不得超过一根马鬃。倒推并不难:边数每翻一倍,内接与外切两条周长的夹缝就缩到四分之一,也就是多定下 log₁₀4 = 0.602 位小数;要十六位,就得从正六边形起翻二十八次,于是 n = 3 × 2²⁸ = 805306368。他取半径 60,用六十进制一路算到十位,内接与外切各算一遍,最后报出 2π = 6;16,59,28,1,34,51,46,14,50,并亲手把它换算成十进制的 6.2831853071795865。今天复算得到的数对得上:这一档的上下界前十六位小数完全相同,夹缝宽 2.4 × 10⁻¹⁷;他给的值与真值差 7.9 × 10⁻¹⁸,末位那个 50 正是逐位算下去的 49、55 四舍五入的结果。那句马鬃的话也验得了——按今天的地球直径,那个圆的周长误差是 0.03 毫米,而一根马鬃在 0.1 毫米上下。此前的纪录是马德哈瓦的十一位(约 1400),此后要等一百七十二年,才由荷兰的范·科伊伦以二十位超过。
拖动翻倍次数:看内接与外切两串数字一位位对齐;底下那条尺按他的验收标准回算,这一档到底够不够细
他的做法与阿基米德同源:内接多边形给下界,外切多边形给上界,边数每翻一倍,两条线就往中间收一次。 新的是次序——他先定下要多准,再倒推该切多少边。翻一倍只买 0.602 位小数,要十六位就得从正六边形起翻二十八次,于是 n = 3 × 2²⁸ = 805306368。这一档的上下界前十六位小数完全相同, 而他报出的 2π = 6;16,59,28,1,34,51,46,14,50(六十进制)换成十进制是 6.2831853071795865, 与真值差 7.9 × 10⁻¹⁸。底下那条尺按今天的地球直径回算他的验收标准:3 × 2²⁸ 正是第一个达标的档位——上一档那个圆还要差 0.37 毫米,到这一档只差 0.09 毫米。 拖到最左边看看起点:正六边形只给出 3 到 3.4641,一位小数都还没有。
