2.3 参考: 色々な $\arctan$ 公式を比較する

http://nalab.mind.meiji.ac.jp/~mk/syori2-2011/jouhousyori2-2011-06/node18.htmlで紹介した公式を比較してみるのも面白いです。

こういうのは実際に試してみて初めて分かることが多く、 挑戦できる腕力のある人は是非やってみることをお勧めします。

いくつかアドバイスをします。

(1)

たくさんの $\arctan$ を計算するので、 $\arctan$ の計算を副プログラム (関数またはサブルーチン) にすることを お勧めします。 次のサンプル・プログラムは、 piarctan.bas を外部副関数を使って計算するように書き換えたものです。

REM piarctan2.bas --- �ޡ������������쥴�꡼���饤�ץ˥åĵ���ǦФ�׻�
REM arctan�ε���η׻��ϳ������ؿ���Ǥ���뤳�Ȥˤ�����
DECLARE EXTERNAL FUNCTION arctan
INPUT x
INPUT n
LET  s=arctan(x,n)
PRINT "arctan(x)��";s
print "����4��";4*s
PRINT USING "�ФȤκ�=-%.###^^^^^^":4*s-PI
END

EXTERNAL FUNCTION arctan(x,n)
REM arctan x �ε������ n ��ޤǷ׻�
LET  f=-x*x
LET  t=x
LET  s=0
for j=1 to n
   LET  a=t/(2*j-1)
   LET  s=s+a
   LET  t=f*t
NEXT j
LET  arctan=s
END FUNCTION

(2)

級数の部分和 $s_n(x)$ は、 $n$ をどこまで大きく取れば十分な精度を持つのか？ 一般にはそれほど簡単ではありませんが、 今考えている $\arctan$ の級数の場合は、 各項の絶対値が単調に減少する交代級数になるので、 次の Leibniz の定理によって簡単に見通しが立ちます ($s_n$ の誤差は、次に足す項の絶対値 $a_{n+1}$ 以下である)。

交代級数に関する Leibniz の定理

数列 $\{a_n\}_{n\ge 1}$ は単調減少 ( $a_1\ge a_2\ge \cdots$ ) で、 $\dsp\lim_{n\to\infty}a_n=0$ を満たすとき、級数

$$s=\sum_{n=1}^\infty (-1)^{n-1}a_n =a_1-a_2+a_3-a_4+\cdots$$

は収束する。第 $n$ 部分和を $s_n$ とすると、

$$s_{2n}\le s\le s_{2n-1}$$   $$\mbox{($n\in\N$)}$$$$,$$

$$\left\vert s-s_n\right\vert\le a_{n+1}$$   $$\mbox{($n\in\N$)}$$

が成り立つ。

次に掲げる結果は、 100桁程度の精度を得るために、 級数の和を何項加える必要があるか概算し、 実際にそれを実行するとどれくらいので精度になるかを計算したものです。

大体どれくらいか…

���㡼��
�� 205 ��ޤDzä����׵����٤�ã�����ޤ�����
 9.81E-101

�ޥ��� 2�����
�� 71 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 22 ��ޤDzä����׵����٤�ã�����ޤ�����
 1.20E-101

Gauss 3�����
�� 40 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 29 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 22 ��ޤDzä����׵����٤�ã�����ޤ�����
 1.23E-102

Gauss 4�����
�� 32 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 29 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 22 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 21 ��ޤDzä����׵����٤�ã�����ޤ�����
 1.09E-103

Gauss 9�����
�� 14 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 13 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 13 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 12 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 12 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 11 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 11 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 10 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 10 ��ޤDzä����׵����٤�ã�����ޤ�����
 7.45E-106

������ͺ��4�����
�� 30 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 29 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 22 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 11 ��ޤDzä����׵����٤�ã�����ޤ�����
 7.29E-105

Stormer��4�����
�� 29 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 22 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 18 ��ޤDzä����׵����٤�ã�����ޤ�����
�� 13 ��ޤDzä����׵����٤�ã�����ޤ�����
 7.73E-104

この結果をどう分析すべきか… 書きたくてむずむずするけれど、 誰かにやって欲しいので書きません。

(一言だけ： Gauss は色々大理論を作った偉い人だけど、 こういう計算も大好きだったのが良く分かります。)