PRINT
PRINT "ある自然数以下の素数の数の予想"
PRINT
INPUT PROMPT "自然数は？ (10000000以下で） ":n
PRINT
DIM A$(7),B(7)
LET A$(2)="    25   約２５％ )"
LET A$(3)="   168   約１７％ )"
LET A$(4)="  1229   約１２％ )"
LET A$(5)="  9592   約　９％ )"
LET A$(6)=" 78498   約　８％ )"
LET A$(7)="664579   約　７％ )"
LET B(2)=25
LET B(3)=168
LET B(4)=1229
LET B(5)=9592
LET B(6)=78498
LET B(7)=664579
LET y=0
FOR i=2 TO n
   LET y=y+1/LOG(i)
NEXT i
PRINT n;"以下の素数の数はおよそ ";INT(y);" と予想します。"
PRINT
PRINT "もっと粗い予想：";n;"/ LOG (";n;") ＝";INT(n/LOG(n))
PRINT
PRINT "○桁までの素数の数の予想"
FOR j=2 TO 7
   LET z=O
   FOR i=2 TO 10^j
      LET z=z+1/LOG(i)
   NEXT i
   PRINT j;"桁までの素数の数  ";
   PRINT USING "######":INT(z);
   PRINT " （ 実測値 ";A$(j);
   PRINT "  誤差： +";ROUND(100*INT(z)/B(j),2)-100;"%"
NEXT j
PRINT
PRINT "このプログラムは終了しました。"
END