目次         

円周率と黄金比

 

円周率と黄金比を構成する式を示す。両者は10とも密接な関係を持っている。


円周率と10の関係式。

π2 = 10-1/(1*2)3-1/(2*3)3-1/(3*4)3-1/(4*5)3-1/(5*6)3-      ← 既知の式(Ramanujan)

 

π2 = 10/ΠP{1+(-1)^(p^2)(p^2)^(-1)} =10/{(1+2-2)(1-3-2)(1-5-2)(1-7-2)(1-11-2)}

 

次式の1/k^2の符号はkを素因数分解した時、 (-1)^(素数2の個数) である。

π2 = 10(1/12-1/22+1/32+1/42+1/52-1/62+1/72-1/82+1/92-1/102+1/112+ )

 

上記式で素数2を奇数個含むマイナスの項(k=2,6,8,10,)を考える。k q1q2qn と素因数分解してr=Σj=1n (qj-1)とする時、次式のkに対する符号はr4n±1に対応して'±'である。例えばk=6=2*3ではr=1+2=3=4*1-1 '-'k=24=2*2*2*3ではr=1+1+1+2=5=4*1+1 '+'

 π = 10(1/2-1/6-1/8+1/10-1/14+1/18-1/22+1/24+1/26-1/30+ )

 


黄金比と円周率と10の関係式。φπの級数と乗積の各項が殆ど同じ形式である。

 

次の円周率の乗積と黄金比の乗積は似た形式である。

 

(φπ)2=(102/22) { 1+1/(1-10)2+1/(1+10)2+1/(1-20)2+1/(1+20)2+1/(1-30)2+1/(1+30)2+ }

   π2 =(102/2 ) { (0!)2/(2!φ2)+(1!)2/(4!φ4)+(2!)2/(6!φ6)+(3!)2/(8!φ8)+ }


ゼータ関数値で表すと

(φπ)=(10/2) e^{ ζ(2)/(2*102)+ζ(4)/(4*104)+ζ(6)/(6*106)+ζ(8)/(8*108)+ }

 (φπ)2=(102/2 ) { (-1)ζ(0)/100+1ζ(2)/102+3ζ(4)/104+5ζ(6)/106+7ζ(8)/108+ }

 

ベルヌーイ数との関係は

 

フィボナッチ数と円周率と10の関係式。Castellanosの式を利用すると、√10 の小数部分をtとして

   10 = 3+t = 3.16227766

 

φ  = 5/π +π/120+7π3/720000+31π5/3024000000+127π7/12096000000000+

φ2 = 25/π2+1/12+π2/6000+π4/3780000+π6/2700000000+

φとπeの関係式。

φ  = (5/π) e^( π2/600+π4/1800000+π6/2835000000+ )

 

黄金比の次式の分母は1,3,6で始まる数字であるが、それは後ろに続く0の個数と等しい。

φ1+2π/10-π3/3000+π5/6000000 = 1.6180341

 

1,2,3 の数字を使用した黄金比の近似式。

φ(322-1/322)^(1/12)            12=3*2*2

φ(233-1/(322-1/322))/12^2

φ-((12^2/233)^12-233)/12^2


1'+'は黄金比と深い関係にあるが、2'+','-'φ ,π ,e と関係が深い。

有名な      の連続平方根式に対し、 √の個数をnとして

    ←既知の式  

       

         

      ( )の個数はn

これに関連して      (φπ) ≒ √e + 1/ ,   φπ ≒ 2 + e + 1/e   ,   e (1+2/(φπ-2) )2

 4(φπ) 20 + 2-1 = 3/2   ,    φπ ≒ 20 + 22 + 2-4 + 2-6 + 2-8 + 2-10 + 2-12

 


カタラン数 1,1,2,5,14,42,132,… を使用した円周率、黄金比、白銀比の級数。

 

 

黄金比の連分数を乗積化すると、フィボナッチ数による乗積表現となる。

φ=(1+1/12) (1-1/22) (1+1/32) (1-1/52) (1+1/82) (1-1/132)

 

次の乗積は分母のフィボナッチ数が和でも積であっても黄金比となる。


黄金比は5の倍数を除く乗積で表現することもできる。

 

参考:5の倍数を除く円周率の級数は

π=(5)(1/1+1/3+1/7+1/9-1/11-1/13-1/17-1/19+1/21+1/23+1/27+1/29-1/31-1/33-1/37-1/39+)

 

オイラー積は素数のある種の乗積であるが、ζ(2)偶数の似た乗積に等しい。


円周率の近似式で有名な祖冲之による密率。

355/113=3.14159292

数字を逆にした分数は円周率の平方根に近いことが知られている。553/311≒√π

 

1133で構成した式。

11√(333333/1.133)=3.14159252

 

より近い値は

355/113-1/44^4=3.141592653551  = 3-1/44^4+16/113

 

密率の精度を上げた式。

(1-1/11776666)/(113/355)=3.141592653589779

 

ただし上記式はラマヌジャンによる次式と本質的に等しい。

(355/113)(1-0.0003/3533)=3.1415926535897943

01,3,5だけで構成した式で、精度もこちらが少し良い。

 

21,3,5だけで構成した式。数字は小数点第19位まで正しい。

355/(113+1/(11^3*3+5^5*2^5+11*3*355/(52+1/(113*3+5*2^5+11/(5+3/(1+2^(1/2)))))))

=3.141592653589793238478



インドのマーダヴァが発見(後にグレゴリーとライプニッツが再発見)した級数  1-1/3+1/5-1/7+1/9-1/11+=π/4  は収束が遅いことで有名であるが、次の級数は少ない項数で円周率の近似値を得ることができる。[ ]内は 1-1/3+1/5-… の項数。「円周率と級数」を参照。

 

[   1]  4*(1)-2/(2+1/2) = 3.2
[   2
]  4*(1-1/3)+2/(4+1/4) = 3.13725
[   3
]  4*(1-1/3+1/5)-2/(6+1/6) =3.142342
[   4
]  4*(1-1/3+1/5-1/7)+2/(8+1/8) =3.1413919
[  10
]  4*(1-1/3+1/5-1/7+-1/19)+2/(20+1/20) = 3.141590242
[ 100
]  4*(1-1/3+1/5-1/7+-1/199)+2/(200+1/200) = 3.14159265356480
[1000
]  4*(1-1/3+1/5-1/7+-1/1999)+2/(2000+1/2000) = 3.141592653589792988…
最初は1項の追加で正しい数字が1桁増加している。

 

級数と連分数の合計による円周率の近似値。

[   1]  4*(1)-2/(2+12/2) = 3.2
[   2
]  4*(1-1/3)+2/(4+12/(4+22/(4))) = 3.142857
[   3
]  4*(1-1/3+1/5)-2/(6+12/(6+22/(6+32/(6)))) = 3.1416252
[  10
]  4*(1-1/3+-1/19)+2/(20+12/(20+22/(20++102/(20)) = 3.1415926535897937321
[ 100
]  4*(1-1/3+-1/199)+2/(200+12/(200+22/(200++1002/(200)) = 3.14284818194
最後の数値は154桁正しい。


覚えている円周率の続きは三角関数の値で知ることができる。

 

sin(3.14   )=0.00159265291                tan(3.14  )=-0.001592654936
sin(3.14159)=0.00000265358979323534
     2 cos(3.14/2)= 0.001592653421
3+sin(3)+sin(3+sin(3))=3.1415926535721
  
3/2+cos(3/2)+cos(3/2+cos(3/2))=3.1415926535897243
/2

参考:最も単純と考えられる円周率計算のプログラム。

  a=3:for k=1 to 4:a+=sin(a):next:?a        計算結果は100桁正しい

   3.141592653589793238462643383279502884197169399375105820974944592307816406

 

双曲線関数も使用した精度の良い式。

sinh(sin(3.14))=0.00159265358979255

3+sinh(sin(3))+sinh(sin(3+sinh(sin(3))))=3.1415926535897932384626433832279

 

sinとcosによる高精度な式。

3+sin(3)/|cos(3)|^(1/3) = 3.14159392…

sin(3.14)/|cos(3.14)|^(1/3) = 0.001592653589793466…

 

参考: a=3:for k=1 to 4:a+=sin(a)/(-cos(a))^(1/3):next:?a  計算結果は788桁正しい

      1,2,3回では6,31,157桁が正しく、1回につき約5倍の正しい桁が求まる。

 

sinと自然対数を使用した更に精度の良い式。

ln{sin(3.14)+(1+sin2(3.14))}+sin3(3.14)/3 = 0.001592653589793238460322


π を連続する自然数の比の極限で表すと

例えばn=10 のとき

 2102 (1*2*3*…*10)/(11*12*13*…*20)*(1*2*3*…*40)/(41*42*43*…*80*81)= 3.1516…

階乗や二項係数で表すと次の式になる。

 

3の数字が主となる円周率の極限式は

 

          目次