RANDOMIZE
PRINT
PRINT "シュミレーション："
PRINT "    サイコロを１０万回振ったときの出た目の回数"
PRINT
DIM C(6),D(6)
MAT C=ZER
FOR i=1 TO 100000
   LET A=INT(RND*6)+1
   LET C(A)=C(A)+1
NEXT i
FOR i=1 TO 6
   LET D(i)=shisyago(C(i)/1000,4)
   PRINT USING "###":i;
   PRINT " : "; 
   PRINT C(i);"回 出た  → ";
   PRINT " ";D(i);"%"
   PRINT
NEXT i

LET X=AVERAGE(6,D)
LET Y=VARP(6,D)
PRINT "以上の確率の平均：";shisyago(X,4);"%"
PRINT "以上の確率の分散：";shisyago(Y,4);"    参考： 標準正規分布の分散は 1.00"
PRINT

PRINT "     ＊ 1 ÷ 6 ≒";shisyago(100/6,4)
PRINT "       ";shisyago(100/6,4);"× 6 ≒";shisyago(x*6,4)
PRINT
END

EXTERNAL FUNCTION shisyago(n,m)
LET shisyago=INT(n*10^m+0.5)/10^m
END FUNCTION

EXTERNAL FUNCTION AVERAGE(N,A())     !平均
LET  s=A(1)
FOR i=2 TO N
   LET  s=s+A(i)
NEXT i
LET  AVERAGE=s/N
END FUNCTION

EXTERNAL FUNCTION VARP(N,A())     !分散
LET  ma=AVERAGE(N,A)
LET  s=(A(1)-ma)^2
FOR i=2 TO N
   LET  s=s+(A(i)-ma)^2
NEXT i
LET  VARP=s/N
END FUNCTION