算術幾何平均の Python コード

「拡張された算術幾何平均」の数値実験コードを Essays 15 に含めようとしたのだが，MathJax ではなぜかうまく表示されない。仕方ないので参考までに画像として添付する；将来 Essays を pdf 化する際には埋め込む予定です。