円周率について語り合おう【π】
lim_[n→∞]n*cos(90-180/n)=π
>>307 >>308
「円周率は3.05より大きいことを証明せよ。(2003)東京大学」
――「 I. Y.」 を名乗るフェイスブックのタイムラインに掲載される、手書きの「過去問」の一節だ。
比較的直近に記されたものらしく、(中略)メモ用紙に、解答までぎっちり書き込まれている。
http://www.j-cast.com/2019/06/17360206.html
円周率ヲタクでも無さそう。 円周率の355/113の近似が精度高いのはπを連分数にしたときに大きな数字がすぐに出るためらしい
なるほど。逆に黄金比は連分数は全部1だから分数近似精度が一番悪い無理数 円周率の正則連分数表示を求めるのってどうやればいいの?
もちろん先に小数を求めてからやるのは無しでπの性質から導く方法が知りたい。
書いてある本とかでも有り難い。 もし>>313へのレスならありがとう。
今ガラケーから見てるんであとでじっくり確認します。 pi^4の連分数展開もでっかい数字が出てくるからラマヌジャンの近似も精度高い
(2143/22)^(1/4)=3.1415926525..
https://www.wolframalpha.com/input/?i=ContinuedFraction%5BPi%5E4%5D (31)^(1/3)=3.1413...
(4930/159)^(1/3)=3.141593... 円周率を11進法で計算していたコンピューターが
1857万桁のところで異変を感知し、その部分を画面に表示し始めた
その表示は0と1のみしか登場せず、ある一定の区間ごとにに折り返され、
0と1によってある図形が浮かび上がった…
0000000011111100000000
0000011110000111100000
0001110000000000111000
0011000000000000001100
0110000000000000000110
1000000000000000000001
1000000000000000000001
0110000000000000000110
0011000000000000001100
0001110000000000111000
0000011110000111100000
0000000011111100000000 ネイピア数
640320^(3/163^0.5pi)=2.718281828459045 Log(640320^3)/163^0.5=3.14159265358979 >>320を少し分かりやすく
ネイピア数
(640320^3)^1/(pi*163^0.5))=2.718281828459045 「おいらの贈り物」 〜人類の至宝 e^(π√163) = 640320^3 + 744 を学ぶ〜 (640320^3)^1/(pi*163^0.5)=2.718281828459045 より高い精度
円周率
Log(640320^3+744)/163^0.5=3.141592653589793238462643383279
ネイピア数
(640320^3+744)^1/(pi*163^0.5)=2.718281828459045235360287471352 π≒2^9≒3.1411
e≒163/(3*4*5)≒2.7166
163(π-e)≒68.9996644963
((2^9)-(163^2/(3*4*5)))≒69.18333... >>321
Log(640320^3)/163^0.5≒3.141592653589793 e^(pi163^0.5)+0.00000000000075≒640320^3-744
の
0.00000000000075の誤差部分はモンストラス・ムーンシャインというやつです。
j(τ)=1/q+744+196884p+21493760q^2+864299970q^3+20245856256q^4
つまり
e^(pi163^0.5)+196884/(e^(pi163^0.5))≒640320^3-744
と評価されます。21493760q^2以降の部分はかなり計算精度が高くないと正確に
求まらないので省略しています。
楕円関数、モジュラー関数、虚二次体、ヘーグナー数、j関数、モンスター群、
モンストラス・ムーンシャインと訳の判らないものがいっぱい出てきます。 芯径10oの感熱紙の芯に糸巻いて切って一周分の長さ測ったらわ。
31.4oぐらいになるはず。 訂正
e^(pi163^0.5)+0.00000000000075≒640320^3-744
の
0.00000000000075の誤差部分はモンストラス・ムーンシャインというやつです。
j(τ)=1/q+744+196884p+21493760q^2+864299970q^3+20245856256q^4
つまり
e^(pi163^0.5)-196884/(-e^(pi163^0.5))≒640320^3-744
と評価されます。21493760q^2以降の部分を追加するとこうなります。
e^(pi163^0.5)-196884/(-e^(pi163^0.5))-21493760/(-e^(pi163^0.5))^2-864299970/(-e^(pi163^0.5))^3
-20245856256/(-e^(pi163^0.5))^4+333202640600/(-e^(pi163^0.5))^5 e^(pi163^0.5)-196884/(-e^(pi163^0.5))≒640320^3+744=262537412640768744
e^(pi67^0.5)-196884/(-e^(pi67^0.5))≒5280^3+744=147197952744
e^(pi43^0.5)-196884/(-e^(pi43^0.5))≒960^3+744=884736744
e^(pi19^0.5)-196884/(-e^(pi19^0.5))≒96^3+744=885480
e^(pi11^0.5)-196884/(-e^(pi11^0.5))≒32^3+744=33512
e^(pi7^0.5)-196884/(-e^(pi7^0.5))-21493760/(-e^(pi7^0.5))^2-864299970/(-e^(pi7^0.5))^3 -20245856256/(-e^(pi7^0.5))^4-333202640600/(-e^(pi7^0.5))^5-4252023300096/(-e^(pi7^0.5))^6≒15^3+744=4119
e^(pi3^0.5)-196884/(-e^(pi3^0.5))-21493760/(-e^(pi3^0.5))^2-864299970/(-e^(pi3^0.5))^3 -20245856256/(-e^(pi3^0.5))^4-333202640600/(-e^(pi3^0.5))^5-4252023300096/(-e^(pi3^0.5))^6≒0^3+744=744
j(τ)=1/q+744+196884q+21493760q^2+864299970q^3+20245856256q^4+333202640600q^5+4252023300096q^6+44656994071935q^7+401490886656000q^8+3176440229784420q^9+22567393309593600q^10+…
j関数
http://mathworld.wolfram.com/j-Function.html 再訂正
e^(pi163^0.5)+0.00000000000075≒640320^3+744
の
0.00000000000075の誤差部分はモンストラス・ムーンシャインというやつです。
j(τ)=1/q+744+196884p+21493760q^2+864299970q^3+20245856256q^4
つまり
e^(pi163^0.5)-196884/(-e^(pi163^0.5))≒640320^3-744
と評価されます。21493760q^2以降の部分を追加するとこうなります。
e^(pi163^0.5)-196884/(-e^(pi163^0.5))-21493760/(-e^(pi163^0.5))^2-864299970/(-e^(pi163^0.5))^3
-20245856256/(-e^(pi163^0.5))^4+333202640600/(-e^(pi163^0.5))^5 再再訂正
e^(pi163^0.5)+0.00000000000075≒640320^3+744
の
0.00000000000075の誤差部分はモンストラス・ムーンシャインというやつです。
j(τ)=1/q+744+196884p+21493760q^2+864299970q^3+20245856256q^4
つまり
e^(pi163^0.5)-196884/(-e^(pi163^0.5))≒640320^3+744
と評価されます。21493760q^2以降の部分を追加するとこうなります。
e^(pi163^0.5)-196884/(-e^(pi163^0.5))-21493760/(-e^(pi163^0.5))^2-864299970/(-e^(pi163^0.5))^3
-20245856256/(-e^(pi163^0.5))^4+333202640600/(-e^(pi163^0.5))^5 "分母が小さいにも関わらず考えている数にかなり近い"
有理数を作れるかが勝負なのです
314159265/100000000=3.14159265
355/113≒3.14159292 118132人目の素数さん2019/06/22(土) 06:55:51.56ID:mnSGZhQY
π - e = 69/163
円周率スレ【π】 - 326 >>304 >>305
バーゼル問題について
藤田岳彦: 数学セミナー, 51(3), p.30-36 (2012/Mar)
「リーマン・ゼータ関数の特殊値を確率論で求める」 314159265/99999999を約分してみる。
314159265/99999999
=104719755/33333333
=34906585/11111111 >>258
> 円周率の小数点以下762桁目から”9”が6回連続で現れるファインマン・ポイント。
これ”確率”で考えるとこの早い段階で見つかるのはレアやな 1415
ふうL@Fu_L12345654321
学コン1傑いただきました!
とても嬉しいです!
https://pbs.twimg.com/media/D-IuUuqVUAALnAB.jpg
https://twitter.com/Fu_L12345654321/status/1144528199654633477
https://twitter.com/5chan_nel (5ch newer account) >>240 >>243 >>316 より
π^4 = 2143/22
= 100 - 3(19/22)
= 100 - 3(20/23) + 9/(22・23)
≒ 100 - 20(3/23) + (3/23)^2
= (10 - 3/23)^2,
π^2 ≒ 10 - 3/23 = 9.869565 2020/03/14 15:09:26.5359
(公財)日本数学検定協会(数検)が「数学の日」制定(1997)
日本パイ協会 の「パイの日」
http://www7a.biglobe.ne.jp/~pienohi/index.htm
A.アインシュタイン (1879/03/14〜1955/04/18) なぜ2πが円一周なのか
なぜπでは半円にしかならないのか
それはおっぱいは2つで一つだからである
一つでは不完全だからである >>243
下の近似式はモジュラー関数に基づく公式
1/π = {(2√2)/(99^2)}Σ[n=0,∞] (4n)!(1103+26390n)/{(4^n)(99^n) n!}^4
の初項から。 Emma Haruka Iwao is a Japanese computer scientist
https://en.wikipedia.org/wiki/Emma_Haruka_Iwao
31.4兆桁出したGoogleの人抜かれちゃったのね。 https:/twitter.com/nami_twun2
ニホンザルヒトモドキを焼き殺せ
https://twitter.com/5chan_nel (5ch newer account) コンピューターの表示環境にもよるが、
n(エヌ)とπ(パイ)が似ていてまぎらわしい。
なんとかしてほしい。 Windows XPのパソコンで2ちゃんねるをしていた時代は、
パイがエヌみたいな形に表示されることはなかったのに。
誰だよ、「エヌにそっくりなパイ」を考案した奴は。 3月の宿題で(1)のみ正解の数弱@shukudai_sujaku
昨年度の大学への数学(大数)での勝率は、
学コンBコースが 1/1 = 100% ,
宿題が 3/10 = 30% でした!
宿題の勝率が低すぎると思うので、
これからは一層精進していきたいです!
https://twitter.com/shukudai_sujaku
https://twitter.com/5chan_nel (5ch newer account) 3月の宿題で(1)のみ正解の数弱@shukudai_sujaku
昨年度の大学への数学(大数)での勝率は、
学コンBコースが 1/1 = 100% ,
宿題が 3/10 = 30% でした!
宿題の勝率が低すぎると思うので、
これからは一層精進していきたいです!
https://twitter.com/shukudai_sujaku
https://twitter.com/5chan_nel (5ch newer account) 円周率で今計算されている数値で連続するn個の10^n個の順列すべてがあるっていう
条件を満たすnの最大値ってどれくらいでしょうか? >>312
黄金比は φ = √(π/1.2) だから同じだよね >>337
π - e = 69/163
π = 2・37・173/(163・25) = 12802/(163・25),
e = 11・19・53/(163・25) = 11077/(163・25) から。 1からNまでの自然数から、無作為に二つを選ぶ、
二つが互いに素となる確率は
N→∞ とするとバーゼル問題の逆数です。 (π/2 - 1)^8 + (4/3)^8 = 10
より
π = 2(1 + [10 - (4/3)^8]^{1/8}),
だから無理数。(8次の代数的数?) π = 3 + 0.1√2 = √(9 +0.6√2 +0.02),
とおく。
√2 の 近似分数を「ペル方程式」を使って求める。
7^2 - 2・5^2 = -1 より
√2 ≒ 7/5,
p = 3 + 0.1√2 ≒ 3 + 7/50 = 3.14
10^2 - 2・7^2 = 2 より
√2 ≒ 10/7,
q = 3 + 0.1√2 ≒ 3 + 1/7 = 3.142857
{p,p,p,p, q,q,q,q,q} の相加平均、調和平均より
π' = (4p+5q)/9 = 3 + 223/1575 = 3.1415873
π” = 9/(4/p + 5/q) = 3 + 1401/9895 = 3.14158666
17^2 - 2・12^2 = 1 より
√2 ≒ 17/12,
π = 3 + 0.1√2 ≒ 3 + 17/120 = 3.1416667
π = √(9 +0.6√2 +0.02)
≒ √(9 +0.85 +0.02) = √(9.87) = 3.1416556 π^2 + (1/π)^2 = (π - 1/π)^2 + 2 ≒ 10,
π - 1/π ≒ 2√2,
これを改良して
π - 1/π + 1/(2π^4) = 2√2,
∴ π = 3.141603 π^2 + (1/π)^2 = (π + 1/π)^2 - 2 ≒ 10,
π + 1/π ≒ 2√3,
これを改良して
π + 1/π + 1/((√6)π^4) = 2√3,
∴ π = 3.1416016
また
π = √3 + √2 - (√3 + √2)/(4(√3)π^4),
1/π = √3 - √2 + (√3 - √2)/(4(√3)π^4), √3 = 1 + (1/2) 1.1^4 = 1.73205
√2 = 1 + 0.8^4
π = √3 + √2 = 2 + (1/2) 1.1^4 + 0.8^4 = 3.14165 >>346
π^2 = (3 + 14/99)^2 = 9 + 28/33 + 1/50
= 10 - 5/33 + 16/(23・33) = 10 - 3/23
= 9.86956522
π^2 = (3 + 14/99)^2 = 9 + 28/33 + 1/50
= 10 - 5/33 + 23/(33^2) = 10 - 142/(33^2)
= 9.86960514 >>363
φ = √(π/a) より
π = aφ^2,
π - a - √(aπ) = 0,
0 = (π-a)^2 - aπ
= π^2 - 3aπ + aa
= (π - 3a/2)^2 - 5aa/4,
π = aφ^2 = (3+√5)a/2 = 1.8 + √1.8 >>373
下から2行目
π = (√3 + √2) {1 - 1/(4√3・π^4)},
を4乗して
π^4 = (√3 + √2)^4 {1 - 1/(√3・π^4)},
これを解いて
{π/(√3 + √2)}^4 = (1 + √(1 - 196/√3 + 80√2))/2 = 0.994072927
∴ π = (0.994072927)^(1/4) * (√3 + √2)
= 0.998514926 (√3 + √2)
= 3.14159194 a = 0.00727079154 に対して
π = √(3-a) + √(2-a) = 3.141591246
1/π = √(3-a) - √(2-a) = 0.318310029
π + 1/π = 2√(3-a) = 3.459901275
π - 1/π = 2√(2-a) = 2.823281217
(a = 1 - exp(-α),
α = 0.00729735257 は Sommerfeld の微細構造定数) >>368 と π ≒ (20/9)√2 から
π = (160/9)√2 - 22 = 3.14157444 π ≒ (20/9)√2 = 3.1427
π ≒ (64/27)(√2 - 4/45) = 3.14151
から
π = (64/29){(77/72)√2 - 4/45} = 3.14159216 >>374
p = √3 + √2 = 3.14626437
π = p - (√2) /p^5 - √(2/3) /p^8 - … = 3.14159223 (3 + √5)(√7 + √11) = 31.2194 > π^3,
(3 + √11)(√5 + √7) = 30.8366 < π^3,
π = [(3 + √5)(√5 + √7)(√7 + √11)(√11 + 3) - (√7)/2]^{1/6}
= 3.141587 π = [(3 + √5)(√5 + √7)(√7 + √11)(√11 + 3)]^{1/6} - 8/(571√385)
= 3.1415926518 tan(1) < π/2.
(略証)
1 = π/3 - δ,
δ = 0.04719755
加法公式で
tan(1) = tan(π/3 - δ)
= {tan(π/3) - tanδ}/{1 + tan(π/3)tanδ}
= (√3 - tanδ)/{1 + (√3)tanδ}
< 3/(√3 + 4 tanδ)
< 3/(√3 + 4 δ)
< 3/(1.732 + 4・0.047)
= 3 / 1.92
= 25/16
= 1.5625
< π/2, >>381
√2+√3を使う規則性のあるπの公式
π = (2/3)/(a/p^2 - (1/2)^3*(a+28)/p^10 + ((1*3)/(2*4))^3*(a+56)/p^18 - ((1*3*5)/(2*4*6))^3*(a+84)/p^26 + ...),
p=√2+√3, a=7-2√6 単位円に内接する正8角形の半周の長さを求める。
(1,0) と (0,1) の間にある頂点を (x,y) とすると
x=y>0, xx+yy=1 より (x,y)=(1/√2, 1/√2) となる。
一辺の長さは
L = √{(1/2) + (1-1/√2)^2} = √(2-√2),
ところで
99^2 - 2・70^2 = 1 より √2 < 99/70,
4・70・41 - 107^2 = 31 より 41/70 > (107/140)^2,
したがって
π > 4L
= 4√(2-√2)
> 4√(2 - 99/70)
= 4√(41/70)
> 4(107/140)
= 107/35
= 3 + 2/35
> 3.05 正解っすね。他の解き方としては三角関数を用いたり、ルート2の近似値を直接用いて解く方法もあったり。
これは有名なので知ってるかもしれないんですが、円周率が3.14であった最後の年に東大で出された受験問題。
円周率を3とした文科省への抗議とされてる有名な問題。 抗議したの誰だろうね。
変な人だね。
初めは概略だけにして、徐々に詳しく…というのが常道だと思うけどな。 行列が終わる年にも、まさかの大問2つまるごと行列だったりしてるから、東大は不満を試験にぶつけてくるw 今知ったんだけど行列なくなったんか
高校レベルの行列って数式の表記の話でしかないし、削る程のものでもないと思うんだけどな 行列が消えて複素数平面が15年ぶりくらいに復活。
複素数平面はかなりえげつないので、行列のままなら受験生がかなりラクなんやけどね。 単位円に内接する8角形の半周の長さを求める。
(1,0) (12/17, 12/17) (0,1) (-12/17, 12/17) (-1,0)
に頂点があるとする。
一辺の長さは
L = √{(12/17)^2 + (1-12/17)^2} = 13/17,
一方、
2・(12/17)^2 = 288/289 < 1,
∴ 頂点および辺は、円周上または内側にある。
凸な折れ線は、外側を通る曲線より短いから
π > 4 L
= 4・(13/17)
= 52/17
= 3 + 1/17
= 3.0588… 単位円に内接する12角形の半周の長さを求める。
A(1,0) B(41/48, 1/2) C(1/2, 41/48) D(0,1)
に頂点があるとする。
辺の長さは
AB = CD = (1/48)√(7^2 + 24^2) = 25/48,
BC = (41/48 - 1/2)√2 = 17/(24√2) > 1/2,
ここで
17^2 - 2・12^2 = 1, 17 > 12√2 を使った。
凸な折れ線は、外側を通る曲線より短いから
π > 4AB + 2BC
> 25/12 + 1
= 3 + 1/12
= 3.0833…
もう秋田? Ln(640320^3+744+196884/(-(640320^3+744))+21493760/(-(640320^3+744))^2)/sqrt(163)
=3.14159265358979323846264338327950288419716939937510582097494
https://hebi.5ch.net/test/read.cgi/news4vip/1621516129/1 >>398の式は微妙に間違っているらしい
正しくは
ln(640320^3+744+196884/(-640320)^3+(744*196884+21493760)/(-640320)^(3*2))/sqrt(163) >>330
>>332-335
結局q^2の所で使う係数を間違えていたから精度が上がらなかったのでしょう。
ムーンシャインとかj-不変量の係数をそのまま突っ込んでも駄目で項を増やしていくと
複雑な計算が必要になってくるようです。
A178449 - OEIS:
https://oeis.org/A178449
↑に係数が乗ってるけど、どうやって計算しているかはよく分からん >>399
とりあえずごちゃごちゃしているのでシンプル?に
ln(640320^3+744-196884/640320^3+167975456/640320^6)/163^0.5
更に精度を上げると
ln(640320^3+744-196884/640320^3+167975456/640320^6-180592706130/640320^9)/163^0.5 さらにもっと精度を上げると
pi-ln(640320^3+744-196884/640320^3+167975456/640320^6-180592706130/640320^9+217940004309743/640320^12)/163^0.5
q^3以上計算はかなり重たくなってきます Log[((3+√5)(√5+√7)(√7+√11)(√11+3)/(4√2))^12-24]/√385
=3.1415926535897932384626433832795028841971693993751058...
https://rio2016.5ch.net/test/read.cgi/math/1613029064/92
Conjecture:
π = Log[k - 24 - 276/k - 8672/k^2 - 344658/k^3 -...]/√385,
k = ((3+√5)(√5+√7)(√7+√11)(√11+3)/(4√2))^12 >>351
1年以上経っても変わっていないな。
誰か新たに計算しているのかな? >>402
なんか精度の挙動がおかしいと思たら
https://oeis.org/A178449
は近似で求めた係数で√163にしか使えない(4次以上が誤り)
q-展開からきちんと求めた係数は以下の通り
[1,744,-196884,167975456,-180592706130,
217940004309744,-282054965806724344,
382591095354251539392,-536797252082856840544683,
772598111838972001258770120,
-1134346327935015067651297762308,
1692324738742597705005194275401888,
-2558136060792026773012451913035887538,
3909566534059719280565543662082528637552,
-6030806348626044568366137322595811547663800,
9377648421379464305085605549750143357652168640,
-14683413510495912973021347501907744913788055440950]
この修正した係数をa[n]とすると
Log[640320^3 + 744 + Σ[n=1,15] a[n]/640320^(3n)]/√163
で誤差が10^(-244)以下
Log[t + 744 + Σ[n=1,15] a[n]/t^n]/√1435,
t = (108 (2+√5)^10 (9559+2212√5+1315√41+425√205)^2 - 12)^3
で誤差が10^(-827)以下になる ちょっと強引なπの1000桁近似をしてみた
近似1(j-invariant, ε<10^(-1055)):
Log[t+744-196884/t+167975456/t^2-180592706130/t^3+217940004309744/t^4-282054965806724344/t^5+382591095354251539392/t^6-536797252082856840544683/t^7+772598111838972001258770120/t^8]/√6307,
t = (27 (4+√17)^16 (1272659166+396488754√17+175423977√53+54313779√901+(564772430+112064468√17+79752371√53+14865745√901)√((27+4√53)/7))^2 - 12)^3
近似2(eta function, ε<10^(-1109)):
Log[u+24-276/u+8672/u^2-344658/u^3+15390480/u^4-737293560/u^5+37026698304/u^6-1923581395371/u^7+102518730258488/u^8-5573961072647172/u^9+307952836032412512/u^10-17239165406937117618/u^11+975709822658417655696/u^12]/√3502,
u = (2 (a+√(a^2-1))^2 (b+√(b^2-1))^2 (c+√(c^2-1)) (d+√(d^2-1)))^6,
a = (23+4√34)/2, b = (19√2+7√17)/2, c = 429+304√2, d = (627+442√2)/2 ただのメモ
e^π-π+1/(1111+(1/(11+((1/√2)))))=20.000000000001214