円周率について語り合おう【π】
lim_[n→∞]n*cos(90-180/n)=π >>238 e^{(√163)π/3} = 640320 + 6.04863735049×10^{-10}, e^{(√163)π} = 640320^3 + 744 -7.499274028×10^{-13}, 640320 = (2^6) 3・5・23・29, 「なぜeやπは様々な性質を持つのか?」 053-056 数セミ増刊「数の世界」日本評論社(1982) 高橋秀俊「"163"の不思議」 p.157-161 >>239 >>217 にもある。 ついでに π = (10 - 3/23)^{1/2} = 3.1415864 π = (10 -130/997)^{1/2} = 3.14159336 π = (31 + 1/159)^{1/3} = 3.14159308 π = (97 + 9/22)^{1/4} = 3.14159265 π = (306 + 5/254)^{1/5} = 3.1415926541 π = [(31^2) + 374/(31^2)]^{1/6} = 3.141592645 「πって本当に無理数なの?」 - 171,172 >>18 S = x + 2Σ[n=1,∞] sin(nx)/n = π (0<x<π) = -π (-π<x<0) 奇数項の和 S_o = 2Σ[m=0,∞] sin((2m+1)x)/(2m+1) 偶数項の和(n=0も含む) S_e = x + 2Σ[m=1,∞] sin(2mx)/(2m), は等しい。 S_o = S_e = S/2, ついでに (π/4)|x| + Σ[m=0,∞] cos((2m+1)x)/(2m+1)^2 = ππ/8 (|x|<π) >>48 x^4・(1-x)^4 /(1+xx) = x^6 -4x^5 +5x^4 -4x^2 +4 -4/(1+xx), . ∫[0,1] x^4・(1-x)^4 /(1+xx) dx = [ (1/7)x^7 -(2/3)x^6 +x^5 -(4/3)x^3 +4x -4arctan(x) ](x=0,1) = 22/7 - π, . . . >>50 x^4・(1-x)^8 /(1+xx) = x^10 -8x^9 +27x^8 -48x^7 +43x^6 -8x^5 -15x^4 +16x^2 -16 +16/(1+xx), . ∫[0,1] (1/4)x^4・(1-x)^8 /(1+xx) dx = (1/4)[ (1/11)x^11 -(4/5)x^10 +3x^9 -6x^8 +(43/7)x^7 -(4/3)x^6 -3x^5 +(16/3)x^3 -16x +16arctan(x) ](x=0,1) = π - 2419/770, >>240 π = (97 + 9/22)^{1/4} = 3.1415926525 は大昔に S.Ramanujan が発見してますた。 他にも色々あります。 π= (63/25)[1 + 10/(7+15√5)] = 3.1415926538 π = (99^2)/(2206√2) = 3.141592730 >>238 (上) e^{10π/3} / (2927+1323√5) = 6 - 3.37940×10^{-11}, (e^{10π/3} + 248・e^{-20π/3}) / (2927+1323√5) = 6 - 1.27627×10^{-23}, (e^{2π√190} + 744) / (12+108(1+√2)^12*[154+210√2+144√5+41√10]^2)^3 = 1 - 1.168664×10^{-70} (e^{(2π√190)/3} + 248・e^{-2(2π√190)/3}) / (12+108(1+√2)^12*[154+210√2+144√5+41√10]^2) = 1 - 2.447924×10^{-72} e^{2π√N} の N は 163, 190, 193, 232, 253 … と続くようです。 なぜeやπは様々な性質を持つのか? -061 >>245 e^(2π√N), N=190の場合に対応する πの公式(厳密な式)を作ってみた 1/π = √(760*(1-(12/α)^3)) Σ[n=0,∞] (6n)!*(β+n)/((3n)!(n!)^3*α^(3n)) ここで α=12+108*(1+√2)^12*(154+210√2+144√5+41√10)^2, β=(852020366870471-395185702196000√2+75149192062748√5-84038529457275√10)/16925656058141292 この公式はn=0項目でπと32桁一致し、1項加えるごとに約35桁ずつ増える >>238 (e^{10π/3} + 248・e^{-20π/3})/(2927+1323√5) = 6 - 1.2762708×10^{-23}, (e^{20π/3} + 248・e^{-40π/3})/(2 +9[1201+537√5 +5^{1/4}・(1607+719√5)/2]^2) = 6 - 6.5828772×10^{-51} >>239 e^{(√163)π/3} -248・e^{-2(√163)π/3} = 640320 - 3.8311867×10^{-26} πの公式(N=253): 1/π = 2√(253*(1-(12/a)^3)) Σ[n=0,∞] (6n)!*(b+n)/((3n)!(n!)^3*a^(3n)) ここで a = 12+(27/16)*(3+√11)^10*(470530+117549√11+15*(9181+838√11)√23), b = (135768074392841107-32707446943866240√11+1875*(-78714146177804+16717935852863*√11)/√23)/1818228027142880892 この公式は1項加えるごとに約40桁増える πの公式(N=400): 1/π = 40√(1-(12/A)^3) Σ[n=0,∞] (6n)!*(B+n)/((3n)!(n!)^3*A^(3n)) ここで A = (3/2)*(2+√5)^13*(869800084+703067697√2+430478740√5+265027941√10+5^(1/4)*(42730416+583140762√2+528899760√5+126712674√10)), B = (10/151706805578559)*(1480534452621+37543365334√2-169393026952√5-69939075619√10-2*5^(1/4)*(148486354235+7499787468√2+53215164243√5-1977920683√10)) この公式は1項足すごとに約52桁増加する この性能はChudnovskyの公式よりも桁数/項数比で3.7倍大きい しかし根号が多く含まれるので計算に向くかどうかは不明 >>248 定数aはもう少し簡単な形になって a = 12+27*(9+2√23)*(191469+57730√11+7*(5688+1715√11)√23)^2 この因数分解は一意ではなくもっと簡単に表されるかもしれません >>245 163と190の間に177があります (e^(2π√177)+744)/((23+3√59)^5*(24587023+11657412√3+2634179√59+1851252√177)^3) = 13500 - 6.56*10^(-64) 2から256までの整数のうちe^(2π√N)がほとんど整数または根号を含む有理数になるNは 少なくとも {2から25までの整数,27,28,30から34までの整数, 36,37,39,40,42,43,45,46,48,49,52,55,57,58, 60,63,64,67,70,72,73,78,82,85,88,93,97, 100,102,112,130,133,142,148,163,177,190,193,232,253} 新しいπの公式見つけた(N=760): 1/π = 4√(190*(1-(12/p)^3)) Σ[n=0,∞] (6n)!*(q+n)/((3n)!(n!)^3*p^(3n)) ここで p = 6*(2+√5)^15*(92210562874930+72942820661700√2+46133090462261√5+29159540310024√10+9*(2346786457760+1856794307075√2+1177625294626√5+744442212212√10)√19), q = (38*(168850305099411513534741603514-65584148987153268032666270500√2+7540774496057908735892220904√5-17937625574963862951981217275√10) -(453911230872051456086681856600-253947186186349584616653907250√2-22230644754209419107315374416√5+45689905770524967041148178125 √10)√19)/95410768893023153163385247146614 この公式はQ(√k,√l,√m)上の級数となるタイプで、1項あたり約72桁増加する ちなみに (e^(4π√190)+744)/p^3 = 1 - 6.93*10^(-146) この先は N=772, 793, 862, 928, 1012 と続くが、多くの場合多重根号が付くようです α = (10 -π^2 -1/π^2)/4 = 0.0072686 = 1/137.578 祖沖之が錬金術師だったら 約率 (Ti)/(N) 密率 355/(Nh) と言ったかも。 なお、355番元素は未発見。 >>18 フーリエ級数展開 θ= -2Σ[n=1,∞] (-1)^n sin(nθ)/n, ( |θ| <π) θ=π-1 とおいて 1 を移項すると π = 1 - 2Σ[n=1,∞] (-1)^n sin(n(π-1))/n = 1 + 2Σ[n=1,∞] sin(n)/n ・・・ (1) θ=π-2 とおいて 2 を移項すると π = 2 - 4Σ[n=1,∞] (-1)^n sin(n(π-2))/2n = 2 + 4Σ[n=1,∞] sin(2n)/2n ・・・・ (2) 2・(1) - (2) より π = 4Σ[n=0,∞] sin(2n+1)/(2n+1), >>254 積分と和が一致する例 ∫(0,∞) x^5/(e^(2πx)-1) dx = Σ[n=1,∞] n^5/(e^(2πn)-1) = 1/504 (ラマヌジャン) 積分と和が一致する例のつづき kを正の整数として ∫(0,∞) (sin x/x)^k dx = 1/2 + Σ[n=1,∞] (sin n/n)^k dx aを√-1の整数倍でない数として ∫(-∞,∞) (sin√(x^2+a^2))/√(x^2+a^2) dx = Σ[n=-∞,∞] (sin√(n^2+a^2))/√(n^2+a^2) = π BesselJ(0,a) >>256 すまない、kは1≦k≦6<2πの整数に訂正し和のdxは要らない 等式の続き: ∫(0,∞) sin x/(x(x^2-π^2)) dx = -1/(2π^2) + Σ[n=1,∞] sin n/(n(n^2-π^2)) = -1/π, ∫(0,∞) sin x/(x(x^2-π^2)(x^2-(3π)^2)) dx = 1/(2π^2(3π)^2) + Σ[n=1,∞] sin n/(n(n^2-π^2)(n^2-(3π)^2)) = 1/(9π^3), ∫(0,∞) sin x/(x(x^2-π^2)(x^2-(3π)^2)(x^2-(5π)^2)) dx = -1/(2π^2(3π)^2(5π)^2) + Σ[n=1,∞] sin n/(n(n^2-π^2)(n^2-(3π)^2)(n^2-(5π)^2)) = -1/(225π^5), …… 円周率の小数点以下762桁目から”9”が6回連続で現れるファインマン・ポイント。 100,000,000桁までの円周率で、同じ数字が六回連続する箇所が数回出てくる。 それどころか、“1234567”も数回出てくるし、“1010101”は2回、 “23571113”(素数を6個並べた数字)は2回、“31415926”が1回出てくる。 もし、1グーゴル桁まで計算出来たらな、0が100回並ぶ箇所があるのだろうか? y = (x/π)^2 とおく。 部分分数分解で 1/{(y-1^2)(y-3^2)・・・・(y-(2k-1)^2)} = (-1)^k {1/[(2k-1)!!]^2 + (1/4)^(k-1)・yΣ[j=1,k] (-1)^j /((2j-1)・(k-j)!(k+j-1)![y-(2j-1)^2]) } 1/{x(x^2-π^2)(x^2-(3π)^2)・・・・[x^2-((2k-1)π)^2]} = (-1/π^2)^k {1/(((2k-1)!!)^2・x) + (1/4)^(k-1)・xΣ[j=1,k] (-1)^j /{(2j-1)・(k-j)!(k+j-1)![x^2-((2j-1)π)^2]} ∫[0,∞] sin(x)・x/(x^2 - (Lπ)^2) dx = (π/2)(-1)^L, 特に∫[0,∞] sin(x)/x dx = π/2, 円周率を求める公式ばっかで円周率について語ってはいない。あほ? アホ(亜父)は南を向いて座った。アホ(亜父)とは范増のことである。 ・・・・・ 司馬遷「史記」(項羽本紀「鴻門之会」) 項王、項伯東嚮坐; 亜父南嚮坐 ─ 亜父者,范増也; 沛公北嚮坐; 張良西嚮侍。 初めてこれを読んだとき、4人が大黒柱の4面に背を向けて座るのかとオモタが・・・・ いくら戦国時代でも、やっぱり対面して座るんだろうなぁ。 >>259 部分分数分解で 1/{(y-b_1)(y-b_2)・・・・(y-b_k)} = 1/{Π[j=1,k] (y-b_j)} = Σ[j=1,k] 1/{(y-b_j)Π[i=1,k (但しjを除く)] (b_j-b_i)} * y≒b_j における主要項を考える。(留数) で、8ギガバイトのメモリーのパソコンで円周率10億桁を数分で計算できる昨今だが、 自然数の二乗の逆数の和が なんで 円周率の2乗/6になるの? N88BASICでプログラムを組んで10億回ループさせたが、なかなか収束しない。 今度はBASICコンパイラで100億回ループさせようかな。 >>264 級数を解析的に求める常套手段: ・Σf(n) の計算は g(z)=f(z)π/tan(πz) と置いてg(z)の留数を解析せよ ・Σ(-1)^n f(n) の計算は g(z)=f(z)π/sin(πz) と置いてg(z)の留数を解析せよ が複素関数論の演習問題などでよく知られている Σ1/n^2 のとき g(z)=π/(z^2 tan(πz))と置くと g(z)の極はz=0,±1,±2,±3,...にあり、留数はそれぞれ-π^2/3,1/1^2,1/2^2,1/3^2,... 積分路をN+1/2-(N+1/2)i,N+1/2+(N+1/2)i,-N-1/2+(N+1/2)i,-N-1/2-(N+1/2)i を頂点とする正方形にとると留数定理より (1/(2πi))∫g(z)dz = -π^2/3 + 2Σ[n=1,N] 1/n^2 が成り立ち、N→∞とすると∫|g(z)||dz|=O(1/N)→0 従って Σ[n=1,∞] 1/n^2 = π^2/6 級数を数値的に求める常套手段 --- オイラー・マクローリンの公式: Σ[n=1,∞]f(n) = Σ[n=1,N-1]f(n) + ∫(N,∞)f(x)dx + (1/2)f(N) - Σ[k=1,M] B_{2k} f^(2k-1)(N)/(2k)! - R_{M,N} ただしB_0=1, B_n=-Σ[k=0,n-1](n!/(k!(n+1-k)!))B_kで定義されるベルヌーイ数で この公式の後の和はM→∞で収束するとは限らないが(漸近級数)、 M,Nを適切に大きく取ることで剰余項R_{M,N}を非常に小さくできる Σ[n=1,∞] 1/n^2 のとき M=10, N=10と置いてオイラーの時代に戻った気分で手計算してみると Σ[n=1,9]1/n^2 + ∫(10,∞)1/x^2 dx + (1/2)1/10^2 - Σ[k=1,10] B_{2k} (-(2k)!/10^(2k+1))/(2k)! = 9778141/6350400 + 1/10 + 1/200 - (-1613404548290414275767377/9699690000000000000000000000) = 1.644934066848226436417… となって π^2/6 と20桁一致する YOUTUBEでバーゼル問題を解説している動画があったが、 解説者や、オイラーって本当に天才だな。正直に言う。さっぱりわからん。 https://www.youtube.com/watch?v=liwyGHRr9dk >>267 リーマンゼータ関数と関数等式ってところでおっぱいとリーマンゼータ関数の関係について語っていて面白いですよ。谷間のとことが難しいらしいです。 【超悪質!盗聴盗撮・つきまとい嫌がらせ犯罪者の実名と住所を公開】 @井口・千明(東京都葛飾区青戸6−23−16) ※盗聴盗撮・嫌がらせつきまとい犯罪者のリーダー的存在/犯罪組織の一員で様々な犯罪行為に手を染めている 低学歴で醜いほどの学歴コンプレックスの塊/超変態で食糞愛好家である/醜悪で不気味な顔つきが特徴的である A宇野壽倫(東京都葛飾区青戸6−23−21ハイツニュー青戸202) ※色黒で醜く太っている醜悪黒豚宇野壽倫/低学歴で人間性が醜いだけでなく今後の人生でもう二度と女とセックスをすることができないほど容姿が醜悪である B色川高志(東京都葛飾区青戸6−23−21ハイツニュー青戸103) ※色川高志はyoutubeの視聴回数を勝手に短時間に何百何千時には何万回と増やしたり高評価・低評価の数字を一人でいくつも増やしたり減らしたりなどの youtubeの正常な運営を脅かし信頼性を損なわせるような犯罪的業務妨害行為を行っています ※色川高志は現在、生活保護を不正に受給している犯罪者です/どんどん警察や役所に通報・密告してやってください 【通報先】 ◎葛飾区福祉事務所(西生活課) 〒124−8555 東京都葛飾区立石5−13−1 рO3−3695−1111 C清水(東京都葛飾区青戸6−23−19) ※低学歴脱糞老女:清水婆婆 ☆☆低学歴脱糞老女・清水婆婆は高学歴家系を一方的に憎悪している☆☆ 清水婆婆はコンプレックスの塊でとにかく底意地が悪い/醜悪な形相で嫌がらせを楽しんでいるまさに悪魔のような老婆である D高添・沼田(東京都葛飾区青戸6−26−6) ※犯罪首謀者井口・千明の子分/いつも逆らえずに言いなりになっている金魚のフン/親子孫一族そろって低能 E高橋(東京都葛飾区青戸6−23−23) ※高橋母は夫婦の夜の営み亀甲縛り食い込み緊縛プレイの最中に高橋親父にどさくさに紛れて首を絞められて殺されそうになったことがある F長木義明(東京都葛飾区青戸6−23−20) ※日曜日になると風俗店に行っている Emma Haruka Iwao smashes pi world record with Google help By Zoe Kleinman Technology reporter, BBC News https://www.bbc.com/news/technology-47524760 Google、円周率計算31兆桁達成 世界記録更新 2019年03月14日 19時56分 公開 https://www.itmedia.co.jp/news/articles/1903/14/news148.html Emma Haruka Iwao @Yuryu Neutral Good with Lawful Evil traits / Developer Advocate for Google Cloud Platform / Software engineer, gamer, queer, and feminist / Tweets are my own https://twitter.com/Yuryu https://twitter.com/5chan_nel (5ch newer account) 31兆桁って数字に意味あんの? 良くわからんけど 円周率が直径と円周の比なら 宇宙と同じ大きさの円があったとしてその直径を1としたとき 1/10^31兆とか素粒子以下の長さにならないの? 宇宙の地平線まで行って戻る時に誤差を水素原子程度に抑えるには小数点以下40桁程度わかってれば十分らしい 円周率をアルファベット26文字、空白一文字、区切り,.;:の四種類 合計31進数にしたときに神からのメッセージが現れるかも? ζ(2x)の結果はπを含みます xが大きくなればなるほど計算は大変になりますが正確に値が出せます ζ(20)(Σ[1/x^22,x,1,n)ではnがn=1で六桁、n=2で12桁ほど求まります tfAAKAVSFXE ヒトモドキゴキブリvカスネズミのネトウヨ自殺しろの 355/113で十分な近似値 113分の355、113355でとっても覚えやすいのよ奥さぁん なんだそれだけ? じゃ、なんのために、スパコン使って何万桁も求めてるの? 宇宙の地平線まで行って戻るのに水素原子程度の誤差に収めるには小数点以下40桁程度で十分だそうだから実用面での要請でないのは間違いない >>283 >宇宙の地平線まで行って戻る そなことできるの? 宇宙の地平線なんてわかってないでしょ。どうやって計算するんでしょうかね? お釈迦さまもびっくりだわ。 観測できてるギリギリのことを比喩表現で地平線と呼んだだけなのになんかキリキリしとる鈍い奴が居よる。 >>287 >観測できてるギリギリのこと フーーん。君、詳しいのね。もっと教えて。 138億光年や465億光年ってのは単なる観測限界で果てではないよね それより遠いと光速超えて膨張しているからこっち側に光が届かず観測できないだけ ボイド構造を見たらごく一部なのが誰でも分かるはず https://pi.delivery/ πの31兆桁の値は一応公開してるんだな 無料版だと100桁ずつしか見られんが >>291 の続き γ = 0.5772156649・・・・ をオイラーの定数とする。 (7) e□e□γ□γ = 13.99983・・・・ (8) e□π□e□γ = 20.99962・・・・ (9) e□e□e□π□γ = 15.00035・・・・ (10) e□e□e□π□γ = 28.000040・・・・ (11) e□π□π□π□e□γ = 71.999996・・・・ (12) π□π□γ□π□γ□e = 9.000001・・・・ (13) π□γ□γ□e□π = 4.00019・・・・ >>292 の続き (14) π□π□γ□π□γ□e = 16.000039835677 (15) π□π□π□γ□e□e = 23.000001617056 ------------------------------ e ≒ 19/7, π ≒ 22/7, γ ≒ 4/7 から1次の関係式 e + π + π = 9.0014671356 = (1) 3e - 2γ = 7.000414156 3π + γ = 10.00199363 π - e + γ = 1.00052649 http://rio2016.5ch.net/test/read.cgi/math/1542555999/102-104 >>293 の続き 3e - 2γ = 7.000414156 = 2 [(7)-(9)] + (1) 3π + γ = 10.00199363 = (9) - (7) + (1) π - e + γ = 1.00052649 = (9) - (7) -π + 4e - 3γ = 5.999887665542 = 3 [(7)-(9)] + (1) -3π + 15e - 11γ = 25.0000771522 = 11 [(7)-(9)] + 4・(1) -4π + 19e - 14γ = 30.9999648177 = 14 [(7)-(9)] + 5・(1) http://rio2016.5ch.net/test/read.cgi/math/1542555999/109 >>276 >.>276 Σ[x=1,n] 1/x^20 n=1, 1.0 n=2, 1.000000953674316406 n=3, 1.000000953961113605 n=4, 1.000000953962023100 ・・・・・ n=∞ 1.000000953962033873 = ζ(20) = (7・283・617)/(10・22!) (2π)^20, Σ[x=1,n] 1/x^22 n=1, 1.0 n=2, 1.0000002384185791015 n=3, 1.000000238450445457 n=4, 1.000000238450502300 ・・・・・ n=∞ 1.000000238450502728 = ζ(22) = 22(131・593)/24! (2π)^22, πは中間子です。 種 類 中間子(ゲージボゾン) π± π° --------------------------------------------------------------- 電 荷 ±e 0 クォーク組成 π+ = ud~, π- = u~d, π°= (uu~-dd~)/√2, 質 量 139.5700×10^6 [eV] 134.9764×10^6 [eV] 寿 命 2.603×10^(-8) [s] 8.4×10^(-17) [s] スピン 0 0 (スカラー ボゾン) アイソスピン 1 1 http://rio2016.5ch.net/test/read.cgi/math/1331723086/180 πは(芳香族)有機化合物中の電子軌道です。 ・C-2p軌道(2p_z)が重なって生じる。 ・原子同士の場合、σ結合の60%ぐらいしか安定化しない。(エチレン) ∵ 電子雲が平行なのでσ結合よりも重なりSが小さい。 ・Cの数が多いときは平面状に広がって非局在化する。このため(超高圧でない場合)エネルギー的に有利。 例:グラフェン、グラファイト ・面対称性により、σ軌道や内殻軌道と直交している。 クーロン積分・交換積分など低次の積分は0である。 また、内殻軌道(1s,2s)とのエネルギー差もかなり大きい。 このため孤立性が強く、π電子だけを考慮する近似が可能。(ヒュッケル法) http://rio2016.5ch.net/test/read.cgi/math/1331723086/181 >>48 22/7 = 3.142857 (アルキメデス) >>50 2419/770 = 3.14155844 >>278 355/113 = 3.141592920 (祖沖之) 103993/33102 と 104348/33215 を「平均」すれば 208341/66317 = 3.14159265347 >>298 103993/33102 (3個) と 104348/33215 (5個) を「平均」すれば 833719/265381 = 3.141592653581  ̄ ̄ ̄  ̄ ̄ ̄ >>298 103993/33102 (4個) と 104348/33215 (7個) を「平均」すれば 1146408/364913 = 3.1415926535914 連分数表示 3 + 1・1/{6 + 3・3/[6 + 5・5/(6 + 7・7/(6 + ・・・・))]} >>293 e ≒ 193/71, π ≒ 223/71,γ ≒ 41/71 〔問題〕 10 - 7/48 < 6ζ(2) < (√2 + √3)^2 を示せ。 ただし ζ(2) = 1 + 1/4 + 1/9 + ・・・・ = Σ[k=1,∞] 1/kk である。 〔系〕 3.139134 < √{6ζ(2)} < √2 + √3 = 3.146264 右 ζ(2) = 1 + 1/4 + 1/9 + Σ[k=4,∞] 1/kk < 49/36 + Σ[k=4,∞] 1/(kk-1/4) = 49/36 + Σ[k=4,∞] {1/(k-1/2) - 1/(k+1/2)} = 5/6 + 19/36 + 2/7 = (5 + 205/42) /6, ∴ 6ζ(2) < 5 + 205/42 < 5 + 44/9 < 5 + 2√6 = (√2 + √3)^2, 6 - (22/9)^2 = 2/81 > 0 より √6 > 22/9, 左 ζ(2) = Σ[k=1,∞] 1/kk = 2 - Σ[k=1,∞] {2/(2k-1) -2/(2k+1) -1/kk} = 2 - Σ[k=1,∞] {4/(4kk-1) - 1/kk} = 2 - Σ[k=1,∞] 1/{(4kk-1)kk} = 2 - 1/3 - 1/60 - 1/315 - Σ[k=4,∞] 1/{(4kk-1)kk} > 2 - 89/252 - (1/63)Σ[k=4,∞] 1/kk = 2 - 89/252 - (1/63){ζ(2) - 49/36}, ∴ 6ζ(2) > 10 - 7/48 = 9.854167 〔問題〕 ζ(2) = Σ[k=1,∞] 1/kk = (log 2)^2 + Σ[k=1,∞] 2/(kk・2^k) を示せ。 (バーゼル問題に関連) >>304 マクローリン展開 Σ[k=1,∞] (1/k)x^(k-1) = -(1/x)log(1-x), より Σ[k=1,∞] 1/(kk・2^k) = -∫[0〜1/2] (1/x)log(1-x) dx, Σ[k=1,∞] {1/kk - 1/(kk・2^k)} = -∫[1/2〜1] (1/y)log(1-y) dy, 辺々引く。 ζ(2) - Σ[k=1,∞] 2/(kk・2^k) = -∫[1/2〜1] log(1-y)/y dy + ∫[0〜1/2] (1/x)log(x) dx, = -∫[0〜1/2] log(x)/(1-x) dx + ∫[0〜1/2] (1/x)log(1-x) dx = [ log(x)log(1-x) ](x=0,1/2) = (log 1/2)^2 = (log 2)^2 = 0.4804530139182 http://club.informatix.co.jp/?p=3326 数列総合スレ http://rio2016.5ch.net/test/read.cgi/math/1290234907/203-205 オイラーの贈物スレ http://rio2016.5ch.net/test/read.cgi/math/1417406099/244 ー247 飯森裕次郎は神奈川・川崎市出身!駒澤大学卒で円周率おじさんだった https://koku-byakunews.com/archives/30124 円周率ヲタクは犯罪者になると叩かれる日も近いなw 円周率ヲタって普段から円周率の十六進法でのある桁の値をBBPアルゴリズムで暗算していたとかかな >>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) 「リーマン・ゼータ関数の特殊値を確率論で求める」 read.cgi ver 07.5.1 2024/04/28 Walang Kapalit ★ | Donguri System Team 5ちゃんねる