next up previous
Next: 3 研究課題2 Up: 情報処理2 第5回 十進BASIC (2) Previous: 1.5 レポート課題5B

2 研究課題1 -- アルキメデスの円周率計算

これはレポートを提出するかどうか任意(余裕がある人向けの「挑戦課題」)。 締め切りはこの講義の最終回まで。 提出方法は、 syori2@math.meiji.ac.jp (@はASCIIの@) に電子メールを送ること。

半径 $ 1/2$ の円の円周 (の長さ) は、円周率 $ \pi$ に等しい ($ r=1/2$ なので $ 2\pi r=2\pi\cdot 1/2=\pi$ )。 有名なシラクサのアルキメデス (BC287?〜BC212) は、 内接正 $ n$ 角形の周の長さ $ p_n$ と、 外接正 $ n$ 角形の周の長さ $ q_n$ を考え、

$\displaystyle p_n<\pi<q_n,\quad \lim_{n\to\infty}p_n=\lim_{n\to\infty}q_n=\pi
$

をもとにして、 実際に $ n=6,12,24,48,96$ に対して $ p_n$ , $ q_n$ を計算して (正六角形から始めて、辺の数を倍にしていった)、 $ \pi$ の評価を得た (「円の測定」という短い論文の中で、 $ 3+\frac{10}{71}<\pi<3+\frac{1}{7}$ を証明している)。 アルキメデスは $ \{p_n\}$ , $ \{q_n\}$ について成り立つ漸化式 (?)

$\displaystyle q_{2n}=\frac{2 p_n q_n}{p_n+q_n},\quad
p_{2n}=\sqrt{p_n q_{2n}}
$

を用いて計算したという5 $ n=3\cdot 2^k$ の場合、 すなわち内接・外接正 $ 3\cdot 2^k$ 角形の周長をそれぞれ $ P_k$ , $ Q_k$ とおく:

$\displaystyle P_k:=p_{3\cdot 2^k},\quad
Q_k:=q_{3\cdot 2^k}.
$

このとき $ \{P_k\}$ , $ \{Q_k\}$ について成り立つ漸化式を導き、

$\displaystyle P_1=p_6=3,\quad Q_1=q_6=2\sqrt{3}
$

と合わせて用いて $ \{P_k\}$ , $ \{Q_k\}$ を計算するプログラムを作成し、 以下の (1), (2) に答えよ。
(1)
正96角形 ($ k=5$ の場合) を利用すると、 $ \pi$ のどのような評価が得られるか。
(2)
ネットで検索すると、 内接正多角形の周長で円周率の近似値を求めた歴史上の話が色々見つかる (円周率マニアは多く、歴史のうんちく蘊蓄を語ってくれます)。 それらの話をいくつか選んで確認せよ (正 $ n$ 角形の周の長さを用いて小数点以下○○桁の値を求めたとある場合、 本当にそれができるか -- ときどきウソが書いてある)。
(3)
$ P_kQ_k$ $ 1:2$ または $ 2:1$ の比の内分点は、 $ \pi$ のより高精度な近似値となる。そのことを確めよ。 (どうしてそうなるか説明できたら素晴らしい。)

以下はヒント。


next up previous
Next: 3 研究課題2 Up: 情報処理2 第5回 十進BASIC (2) Previous: 1.5 レポート課題5B
桂田 祐史
2013-05-22