java圆周率祖冲之减法次数_祖冲之是如何计算圆周率的（一)

祖冲之计算圆周率，使用的方法依然是刘徽的割圆术。唯一不同的地方是，他可能掌握了独特的开平方办法。那么用割圆术这种方法，每割一次，圆周率大约可以前进几位呢？得到7位小数，需要迭代几次呢？

用数据说话，不再用复杂的微积分计算。

先引用隋书上对祖冲之的记载：

故体有长短，检之以度，则不失毫厘；

物有多少，受之以器，则不失圭撮；

量有轻重，平之以权衡，则不失黍丝；

声有清浊，协之以律吕，则不失宫商；

三光运行，纪以历数，则不差晷刻；

事物糅见，御之以率，则不乖其本。

故幽隐之情，精微之变，可得而综也。

夫所谓率者，有九流焉：

一曰方田，以御田畴界域。

二曰栗米，以御交质变易。

三曰衰分，以御贵贱廪税。

四曰少广，以御积冪方圆。

五曰商功，以御功程积实。

六曰均输，以御远近劳费。

七曰盈肭，以御隐杂互见。

八曰方程，以御错糅正负。

九曰句股，以御高深广远。

皆乘以散之，除以聚之，

齐同以通之，今有以贯之。

则算数之方，尽于斯矣。

古之九数，圆周率三，圆径率一，其术疏舛。

自刘歆、张衡、刘徽、王蕃、皮延宗之徒，各设新率，未臻折衷。

宋末，南徐州从事史祖冲之，更开密法，以圆径一亿为一丈，

圆周盈数三丈一尺四寸一分五厘九毫二秒七忽，

朒数三丈一尺四寸一分五厘九毫二秒六忽，

正数在盈朒二限之间。

密率，圆径一百一十三，圆周三百五十五。

约率，圆径七，周二十二。

又设开差冪，开差立，兼以正圆参之。

指要精密，算氏之最者也。

所著之书，名为《缀术》，学官莫能究其深奥，是故废而不理。

最前面是说数学很重要，日常生活中要用，研究历法、天文、音乐也要用。

然后说九数，就是现在能看到的《九章算术》编排的顺序。刘徽的割圆术就在《九章算术》的注释中出现。具体是第一章《方田》。方田讨论各种形状的土地面积，有长方形的、正方形的、三角形的、梯形的;还有圆形的、凸月形的、弓形、环形的。其中，刘徽对园田计算时的古圆周率不满意，所以开发了割圆术。

下面讲，祖冲之的伟大成就。“朒”读音是“女”第四声，意思是“亏缺，不足”。

割圆术，仍然用刘徽的方法，但数据要计算的更远。

下面是本人用电脑模仿的次序：

置圆径2尺，半径1尺。

割6觚为12觚术曰:

令半径为弦，半面为句，为之求股。

[[正 6 边形:

边长(小弦)=0.1e1

(正 12 边形面积:

0.3e1)]]

置上小弦幂，四而一之,得(正6边形边长平方的四分之一)

0.25e0,余分弃之，即句幂也。

以减弦幂，其余开方除之，得股(边到圆心的距离)

0.86602540378443864676372317075293618347140262690519031402790348972596650845440001854057309337862429e0

以减半径，余(即30.000000度角的正矢)

0.13397459621556135323627682924706381652859737309480968597209651027403349154559998145942690662137571e0,谓之小勾。

觚之半面又谓之小股，为之求小弦。其幂：

0.26794919243112270647255365849412763305719474618961937194419302054806698309119996291885381324275142e0

开方除之，即12觚之一面也。

割12觚为24觚术曰:

令半径为弦，半面为句，为之求股。

[[正 12 边形:

边长(小弦)=0.51763809020504152469779767524809665669813780263986102762800641463011394949760399384473594938849933e0

(正 24 边形面积:

0.3105828541230249148186786051488579940188826815839166165768038487780683696985623963068415696330996e1)]]

置上小弦幂，四而一之,得(正12边形边长平方的四分之一)

0.66987298107780676618138414623531908264298686547404842986048255137016745772799990729713453310687856e-1,余分弃之，即句幂也。

以减弦幂，其余开方除之，得股(边到圆心的距离)

0.96592582628906828674974319972889736763390483900840455040234307631042321397985551634756174185807045e0

以减半径，余(即15.000000度角的正矢)

0.34074173710931713250256800271102632366095160991595449597656923689576786020144483652438258141929548e-1,谓之小勾。

觚之半面又谓之小股，为之求小弦。其幂：

0.68148347421863426500513600542205264732190321983190899195313847379153572040288967304876516283859096e-1

开方除之，即24觚之一面也。

割24觚为48觚术曰:

令半径为弦，半面为句，为之求股。

[[正 24 边形:

边长(小弦)=0.26105238444010318309681245579097802038748140962346450378123389667026565045964725922848073798458425e0

(正 48 边形面积:

0.31326286132812381971617494694917362446497769154815740453748067600431878055157671107417688558150109e1)]]

置上小弦幂，四而一之,得(正24边形边长平方的四分之一)

0.17037086855465856625128400135551316183047580495797724798828461844788393010072241826219129070964774e-1,余分弃之，即句幂也。

以减弦幂，其余开方除之，得股(边到圆心的距离)

0.9914448613738104111445575269285628712777382744481022714587746035289220684050825317632654334532774e0

以减半径，余(即7.500000度角的正矢)

0.85551386261895888554424730714371287222617255518977285412253964710779315949174682367345665467226026e-2,谓之小勾。

觚之半面又谓之小股，为之求小弦。其幂：

0.171102772523791777108849461428742574445234511037