グレゴリー級数での円周率計算

π/4=atn(1)=1-1/3+1/5-1/7+・・・

　
最初に文字列一致が見られる項の番号を計算
3.14159265358979323846・・・

計算項数	計算結果	計算時間
3 項を計算	3.4666666666666	20050209 23:18:50
19 項を計算	3.1941879092319	20050209 23:18:50
119 項を計算	3.1499958665934	20050209 23:18:50
1688 項を計算	3.1410002365801	20050209 23:18:51
10794 項を計算	3.1415000095284	20050209 23:18:53
136121 項を計算	3.1415999999948	20050209 23:19:20
1530012 項を計算	3.1415920000001	20050209 23:24:27
18660270 項を計算	3.1415926000000	20050210 00:28:25
156001209 項を計算	3.1415926599999	20050210 08:42:42

―――――――――――――――――――――――――――――――――――
!　π=4*atn(1)=4*(1-1/3+1/5-1/7+・・・)       十進ＢＡＳＩＣ

OPTION ARITHMETIC DECIMAL_HIGH
LET CNT=0
FOR I=1 TO 100000000000 STEP 2
LET GGG=GGG+(-1)^((I-1)/2)*1/I

IF left$(STR$(GGG*4),CNT+2)=left$("3.141592653589793",CNT+2) AND CNT=CNT THEN
PRINT (I-1)/2+1;"項目 ",left$(STR$(GGG*4),15);" ";DATE$;" ";TIME$
LET CNT=CNT+1
END IF
NEXT I

PRINT "終了!"
END