多数桁のπ(円周率)の値を計算するFORTRAN及びCのソースプログラムを掲載する。本プログラムは内容を理解してもらうことを目的とし、計算時間及び使用メモリ量の効率化は目的としていない。本格的なプログラムは金田研究室から他のプログラムと合わせて公開する予定。4バイト整数に10進5桁を入れて計算。演算量は桁数の2乗に比例する。)と桁数を画面上から応答する。画面上に利用したArctanとπの計算時間を表示する。πの値はAtanOrg?.txtに出力する。?は1〜5の利用公式の番号である。倍精度実数に10進5桁を入れて計算。演算量は桁数nのn・log2に比例する。多数桁乗算は複素数のFMT(高速剰余法)を使用している。桁数を画面上から応答する。計算桁数は2のべき乗の5倍である。現在の利用公式はウス・ルジャンドル公式であるが、利用公式の応答要。画面上にAGMの各stepとπの計算時間を表示する。πの値はAGM-2nd?.txtに出力する。Errorが出力される。これは、FMTによる多数桁乗算時に丸め誤差のため実数を整数にする過程で発生した最大誤差である。これが0.1以上の値となる場合は1要素に詰める桁数を10進5桁から4桁にする必要がある。倍精度実数に10進5桁を入れて計算。演算量は桁数nのn・log2〜3に比例する。多数桁乗算は複素数のFMT(高速剰余法)を使用している。)と桁数を画面上から応答する。画面上に利用したArctanとπの計算時間を表示する。πの値はAtanDRM?.txtに出力する。?は1〜5の利用公式の番号である。Errorが出力される。乗算時に丸め誤差のため実数を整数にする過程で発生した最大誤差である。これが0.1以上の値となる場合は1要素に詰める桁数を10進5桁から4桁にする必要がある。倍精度実数に10進5桁を入れて計算。演算量は桁数nのn・log2〜3に比例する。多数桁乗算は複素数のFMT(高速剰余法)を使用している。)と桁数を画面上から応答する。画面上に利用したDRMの各step単位とπの計算時間を表示する。πの値はRamnDRM?.txtに出力する。?は1〜3の利用公式の番号である。Errorが出力される。乗算時に丸め誤差のため実数を整数にする過程で発生した最大誤差である。これが0.1以上の値となる場合は1要素に詰める桁数を10進5桁から4桁にする必要がある。倍精度実数に10進5桁を入れて計算。演算量は桁数nのn・log2〜3に比例する。多数桁乗算は複素数のFMT(高速剰余法)を使用している。原理プログラムDRM法によるAtan公式計算プログラムとマヌジャン公式計算プログラムを載せる。3〜62進数までの任意の基底による円周率の多数桁計算を行う。任意基底への変換は小数点以下に対してで、最初の値は常に3と表示する。結果の表示方法は下記の2種類がある。英小文字(26進数)や英大小文字(52進数)で表示するときに都合が良い。倍精度実数に指定した基底数(3〜62)を3〜11桁入れて計算。倍精度実数に入れる桁数は自動調整する。演算量は桁数nのn・log2〜3に比例する。多数桁乗算は複素数のFMT(高速剰余法)を使用している。)と桁数を画面上から応答する。画面上に利用したArctanとπの計算時間を表示する。πの値はAtanDRM?.txtに出力する。?は1〜5の利用公式の番号である。倍精度実数に指定した基底数(3〜62)を3〜11桁入れて計算。倍精度実数に入れる桁数は自動調整する。演算量は桁数nのn・log2〜3に比例する。多数桁乗算は複素数のFMT(高速剰余法)を使用している。)と桁数を画面上から応答する。画面上に利用したDRMの各step単位とπの計算時間を表示する。πの値はRamnDRM?.txtに出力する。?は1〜3の利用公式の番号である。xdrmrmnj.txt。