- ベストアンサー
円周率の大桁数の計算
円周率 3.1415~延々と続きますが、この莫大な桁の数字はどうやって計算するのでしょうか? ・計算機ならではの特別なアルゴリズムがあるのか? ・手計算で無理やり算出する方法はあるのか? 正多角形を内接、外接させる方法は桁数が小さい場合に有効ですが、何十桁も求めるには不向きと思います。世に出ている数値はどうやって計算したのでしょうか?素人考えですが、今日は40~50桁数、明日は51~60桁などという方法があるのでしょうか?
- over_the_galaxy
- お礼率87% (424/485)
- 数学・算数
- 回答数1
- ありがとう数1
- みんなの回答 (1)
- 専門家の回答
質問者が選んだベストアンサー
いろいろあるみたいですよ。arctan x を変形してグレゴリ級数で計算するというのが ポピュラーみたいですね。 「円周率の歴史」 http://ja.wikipedia.org/wiki/%E5%86%86%E5%91%A8%E7%8E%87%E3%81%AE%E6%AD%B4%E5%8F%B2 私はガウスの公式で三万桁求めたことがあります。当時(486 33MHz)で20分くらいでした。 いまなら秒殺でしょうね(^^;
関連するQ&A
- 円周率の計算はなんのため?
円周率の計算はなんのため? 最近、自作PCで円周率を5兆桁まで計算したというのが話題になりました。 円周率の話題が出るたびに思いますが、何のために計算するのですか? 学者やマニアがおもしろがってやってるだけですか? どうも、円周率の桁数はどうでもよく、計算機の性能を自慢したいだけのような気がしますが・・・どうなのでしょうか。
- ベストアンサー
- 数学・算数
- アルキメデスが円周率を計算したやり方は?
Blue Backs「パソコンで挑む円周率」で教えられたのですが、世界で最初に円周率を計算により求めたのはアルキメデスとのことです。彼は円に内接・外接する正96角形の周の長さから円周率の近似値を計算し、3.14までは正確に求めたとのことです。 大変ためになる情報ですが、残念ながら私には正96角形の周の長さを求めるやり方が分かりません。アルキメデスは三角関数を知っていたのですか? 三角関数を知っているとしても、それを計算できたのでしょうか。 たぶん簡単なやり方があるのでしょうが、どなたか親切な方、教えてください。
- ベストアンサー
- 数学・算数
- 円周(率)の計算方法
文系人間です。 高校の数学は、計算問題は解けるけど…といったレベルでした。 円周率は、「円周÷直径」と習いました。 しかし、実生活において円周や直径を実測しようと思ったら、 定規や巻き尺を使ってせいぜい10分の1ミリまでが限界です。 それで正確な(小数点以下何兆ケタの)円周率が算出できるとは思えません。 ということは円周も計算式で求めなければならないということです。 で、円周は、「直径×円周率」…これでは堂々巡りですね。 円周を円周率を使わずに求める方法ってあるんでしょうか。 ある値の近似値ってその値そのものですか? 小数点以下何兆ケタの円周率を算出する公式は調べれば出てくると思いますが、 その公式が意味するところは、一般人にも分かるように説明できるのでしょうか。 また、そのような公式は、円周と関係あるのでしょうか、ないのでしょうか。 円周と関係ないとしたら、どうして定義から離れたところで、円周率が算出できるのでしょうか。 円周率に対するもやもやした気持ちを言葉にすること自体が難しいのですが、 あえて質問にするとしたらこんな形です。 推察するに聞きたいことはそうじゃないだろ、というご意見でも結構です。 よろしくお願いいたします。
- ベストアンサー
- 数学・算数
- NASAが使う円周率って???
お世話になります。 単純にあくまでも、好奇心だけの質問なのですが、 NASAやJAXAなんかの宇宙規模の技術計算をする場合、 円周率は、何桁ぐらいのものを使っているのでしょうか? 私も仕事等で円周率を使う場合がありますが、この場合は、 PCのdouble型の有効数字10桁少々のもので十分 事足りています。 宇宙開発に限らず、実用で使用している円周率の桁数について 情報をお持ちの方がいらっしゃいましたら、よろしくお願いします。
- 締切済み
- 数学・算数
- 表計算(エクセル)で、円周率の近似値を求めようとすると・・・
コンピュータのカテゴリに書き込もうかとも迷いました。 場違いでしたらすみません。 表計算ソフト(エクセル)で、遊びで円周率の近似をやってみることにしました。半径=0.5の円に内接する正多角形の周の長さを求めるやり方です。半径=0.5にしたのは、直径1の円の円周率は、周の長さをそのまま円周率とすることができるからという理由からです。 半径0.5の円に内接する正多角形の隣り合う2点ABと円の中心Oとを結んで出来る二等辺三角形OABの辺ABの長さは、 =√(0.5^2+0.5^2-2*0.5*0.5*cos∠AOB)(余弦定理) =√(0.5^2*2*(1-cos∠AOB)) エクセルでの具体的な計算の仕方 (1) A1セルに「=3」 B1セルに「=SQRT(0.5^2*2*(1-COS(RADIANS(360/A1))))*A1」 ※これでA2セルには、半径0.5の円に内接する正三角形の周の長さが表示されます。 (2) A2セルに「=A1+1」 B2セルに「=SQRT(0.5^2*2*(1-COS(RADIANS(360/A2))))*A2」 (3) A2、B2を選択して、下方向へオートフィルします。 オートフィルを続ければ続けるほど、正n角形のnが増大するので、3.14にB列に表示される数値は、”下の行に行くほどどんどん円周率πに近づく”はずです。 なのに、正4316角形と正4317角形(セルB4315とセルB4315)では、 正4316角形の周の長さ=3.14159237622779 正4317角形の周の長さ=3.14159237622464 となっており、正4316角形の周の長さよりも正4315角形の周の長さのほうが長いことになっています。 正∞角形の周の長さ÷直径=円周率というのは、数学の教科書にも載っているようなことなので、”下の行に行くほどどんどん円周率πに近づく”という考え方自体は間違っていないと思うのですが・・・ コンピュータの限界とか、そういう問題でしょうか?
- ベストアンサー
- 数学・算数
お礼
ありがとうございます。 ガウスの式、調べました。この方法なら確かに出来そうです。 扱える数値の最小値は決まってるでしょうから、コツが必要ですね。