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

CLS
FOR I=1 TO F
　Z[I]=FLOOR(B/5)
NEXT

N=F
C=0
FOR I=1 TO H　　　　　　　// 外側ループは結果格納配列数でよい
　LOCATE 1,1:PRINT N
　T=0
　FOR J=N-1 TO 1 STEP -1
　　D=2*J-1
　　W=T*J+Z[J]*B
　　T=FLOOR(W/D)
　　Z[J]=W-T*D
　NEXT

　W=FLOOR(T/B)
　A0[I]=C+W
　C=T-W*B

  N=N-E
NEXT

PRINT "PI=";　　　　　　// 計算結果表示
FOR I=1 TO H
　PRINT A0[I];
NEXT

PRINT "END"