【１】一般化

　　ｘn＝ｆ（ｘn-1，ｙn-1）

　　ｙn＝ｇ（ｘn-1，ｙn-1）

を繰り返すとき，数列ｘn，ｙnが同じ極限Ｍ（ｘ，ｙ）に収束し，極限関数Ｍ（ｘ，ｙ）が具体的に定められるのはどのような場合でしょうか？

［１］収束性

　答えを先にいうと，ｆ，ｇが比較可能な平均という性質をもつとき，同じ値に収束することが証明されます．

(1)平均の定義

　　ｘ＜ｙならば，ｘ＜ｆ（ｘ，ｙ）＜ｙ

　　λ＞０ならば，ｆ（λｘ，λｙ）＝λｆ（ｘ，ｙ）

(2)比較可能関数の定義

　　２つの関数の間に不等式ｆ（ｘ，ｙ）≦ｇ（ｘ，ｙ）が成り立つ．

この意味でＨ，Ｇ，Ａ，Ｅは平均の条件を満たし，不等式Ｈ＜Ｇ＜Ａ＜Ｅが成り立ちます．

［２］極限関数

　極限関数Ｍ（ｘ，ｙ）は初期値（ｘ0，ｙ0）の如何に関わらず，

　　Ｍ（ｘ，ｙ）＝Ｍ（ｘ0，ｙ0）

したがって

　　Ｍ（ｘ，ｙ）＝Ｍ（ｆ（ｘ，ｙ），ｇ（ｘ，ｙ））

を満たします（不変性）．

（例１）ｆ＝２ｘｙ／（ｘ＋ｙ）＝Ｈ，ｆ＝（ｘ＋ｙ）／２＝Ａ

とおくと，ｘ0ｙ0＝ｘ1ｙ2＝・・・という不変性を示す．この極限関数は

　　Ｍ（ｘ，ｙ）＝（ｘｙ）^(1/2)＝Ｇ

に急速に収束する．→算術調和平均

（例２）ｆ＝Ｇ，ｇ＝Ａならば極限Ｍ（ａ，ｂ）は楕円積分

　　Ｍ（ａ，ｂ）＝１／（2/π∫(0,π/2)ｄφ/√{(acosφ)^2+(bsinφ)^2}）

により表すことができる（ガウス）．→算術幾何平均

　ラグランジュ，ガウス，ルジャンドルは１８世紀に算術幾何平均を熟知していたようです．ガウス・ルジャンドルの算術幾何平均法では，反復のたびに有効数字は２倍になりますから，数値計算アルゴリズムの強力な武器となりえます．

　楕円積分の二重周期は算術幾何平均法を使って計算されます．実際，東京大学の金田康正氏のグループは楕円積分の計算と関係した算術幾何平均法を用いて円周率πの計算の世界記録を樹立しました．その計算量はＯ（ｎｌｏｇｎ）となり，計算能率はarctan（ｘ）の展開公式Ｏ（ｎ^2）よりも格段に優れています．円周率πの計算や巨大な素数の発見はコンピュータシステムの信頼性や処理速度といった性能をテストするのに最適ということです．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

【２】バリエーション

(1)２数ａ0，ｂ0をとり，ａ1＝（ａ0＋ｂ0）／２，ｂ1＝√ａ1ｂ0＝√ｂ0（ａ0＋ｂ0）／２を計算する．次に，ａ2＝（ａ1＋ｂ1）／２，ｂ2＝√ａ2ｂ1とする．ｂnを計算するときａn-1でなく最新のａnを使っていることに注意されたい．

　すると，ａnとｂnは急速に同じ極限Ｍ（ａ，ｂ）に到達する．このとき，極限関数はＭ（ａ，ｂ）＝

（ｂ^2−ａ^2）^(1/2)／ａｒｃｃｏｓ（ａ／ｂ）　　　　（０≦ａ＜ｂのとき）

ａ　　　　　　　　　　　　　　　　　　　　　　　　　（ａ＝ｂのとき）

（ａ^2−ｂ^2）^(1/2)／ａｒｃｃｏｓｈ（ａ／ｂ）　　　（ｂ＜ａのとき）

で与えられる（パッフ）．

(2)０≦ａ0＜ｂ0なる２数ａ0，ｂ0をとり，ａ1＝｛ａ0（ａ0＋ｂ0）／２｝^(1/2)，ｂ1＝｛ｂ0（ａ0＋ｂ0）／２｝^(1/2)を計算する．これを繰り返すとａnとｂnは急速に同じ極限Ｍ（ａ，ｂ）に到達する（カールソン）．

　　Ｍ（ａ，ｂ）＝｛（ｂ^2−ａ^2）／２ｌｏｇ（ａ／ｂ）｝^(1/2)

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

【３】調和幾何平均

「２数ａ0，ｂ0をとり，それらの調和平均ａ1＝２ａ0ｂ0／（ａ0＋ｂ0），幾何平均ｂ1＝√ａ0ｂ0を計算する．次に，ａ1，ｂ1の算術平均と幾何平均を計算し，ａ2＝２ａ1ｂ1／（ａ1＋ｂ1），ｂ2＝√ａ1ｂ1とする．すると，ａnとｂnは急速に同じ極限Ｈ（ａ，ｂ）に到達する．」

　算術幾何平均と調和幾何平均は

　　Ｈ（ａ，ｂ）Ｍ（ａ，ｂ）＝ａｂ

という簡明な関係で結ばれている．

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝