INPUT "DIGITS=";A　　　// 計算する桁数
INPUT "BASE=";K　　　　　// 変数あたりの桁数
B=POW(10,K)
F=FLOOR(A/LOG(2,10)+1)　　// 計算項数
H=FLOOR((A+K-1)/K+1)　　//　結果格納配列数
DIM　Z[H+1]　　　　　　　　//　配列確保

CLS
FOR N=F TO 1 STEP -1
　LOCATE 1,1:PRINT N
　C=0
　FOR I=H TO 1 STEP -1
　　W=Z[I]*N+C
　　C=INT(W/B)
　　Z[I]=W-C*B
　NEXT
　G=N*2+1
　FOR I=1 TO H
　　T=Z[I]+C*B
　　W=FLOOR(T/G)
　　C=T-W*G
　　Z[I]=W
　NEXT
　Z[1]=Z[1]+2
NEXT

CLS           　　　// 計算結果表示
PRINT "PI=3.";
FOR I=2 TO H
　PRINT Z[I];
NEXT

PRINT "END"