;π

1#HltCr 2#EndCr %PI %D0 500#OtCnt #MDtRc[0,0,50,8,LZ,*L,1,T] ;小数点以下から揃える *L ⇒ #L
500 n10/8+2=Im50 1280zM50>_E_ 1701m75 9900m76 VZ

; �納 (-1)^n * 48/( (2n+1) * 18^2n+1 )]
N10/[18la*2] +1/2=I+3=m30
7100m55-M50+1=m10 48E+7m22 mM10 8500m56 18m8500 9900m57 V/
324m1402 1m9
1m1401 18m3000 3000m51 5700m56 4400m57m42m46 9900m41m43m45m47
:L0 2mi1401 1402m52m54 M51m53V* 1401m52m54 M56m53V* V/ M9/2=Ff Z=_G1_V-G2 :L1 V+
:L2 M30 LM9I_G0_G3

; �納 (-1)^n * 32/( (2n+1) * 57^2n+1 )]
:L3 N10/[57la*2] +1/2=I+2=m30
7100m55-M50+1=m10 32E+7m22 mM10 8500m56 57m8500 3000m57 V/ 9900m41m43 3000m42 V+

1701m75 5700m76 VZ 3249m1402 1m9
1m1401 57m3000 3000m51 5700m56 4400m57m42m46 9900m41m43m45m47
:L4 2mi1401 1402m52m54 M51m53V* 1401m52m54 M56m53V* V/M9/2=Ff Z=_G5_V-G6 :L5 V+
:L6 M30 LM9I_G4_G7

; -�納 (-1)^n * 20/( (2n+1) * 239^2n+1 )]
:L7 N10/[239la*2] +1/2=I+1=m30
7100m55-M50+1=m10 20E+7m22 mM10 8500m56 239m8500 3000m57 V/ 9900m45m47 3000m46 V-

1701m75 5700m76 VZ 57121m1402 1m9
1m1401 239m3000 3000m51 5700m56 4400m57m42m46 9900m41m43m45m47
:L8 2mi1401 1402m52m54 M51m53V* 1401m52m54 M56m53V* V/M9/2=Ff Z=_G9_V+G10 :L9 V-
:L10 M30 LM9I_G8_
:L90 9901-M50=PM50*M #MDtRc[MN,0,80,D,RS,L,0,Z] 1jN10 2j;桁
0jRL 0j+0,1,2 0jRL+ 0J" p>>@" 1.5T%D1, E :E


 5000桁　[WindowsXP　Core2Duo_2.5GHz_11分　Pentium4_2.5GHz_40分]　 [Windows95　200MHz_03h22m]

10000桁　[WindowsXP　Core2Duo_3.5GHz_65分　Pentium4_2.5GHz_316分]　[Windows95　200MHz_26h46m]

桁数=8*M50　n=M30　　　処理時間≒基準桁の処理時間×桁数比^3


ガウスの公式
π＝48*Arctan(1/18) ＋ 32*Arctan(1/57) − 20*Arctan(1/239)
　＝[0--N:Lm]�納(-1)^n/(2n+1) * [(48/18^2n+1)+(32/57^2n+1)−(20/239^2n+1)]]
π＝3.1415 92653 58979 32384 62643 38327 95028 84197 17

マーチンの公式
π＝16*Arctan(1/5) − 4*Arctan(1/239)
　＝[0--N:Lm]�納 (-1)^n/(2n+1) * [(16/5^2n+1)−(4/239^2n+1)] ]


πの書籍
円周率πの不思議　　　(堀場芳数著)[講談社]　　πの歴史，πの展開公式，他
パソコンで挑む円周率　(大野栄一著)[講談社]　　πの計算式，BASICプログラム，πの値(3万桁)，他
πのはなし 　 　　　　(金田康正著)[東京図書]　スーパーコンピューターによるπの計算，πの値(10万桁)，他

スーパーπ (東京大学金田研究室 フリーソフト)で、πの値が得られます。　(Ver 1.1　3355万桁)



VPi

[ 2007年, 2008年, 2009年 ]


前のページ ＡＴＨＰトップページ