OPTION BASE 0
PRINT
PRINT "平均値が５０，標準偏差が２０ の正規様分布をしている ０ 〜 ９９ の整数が多数集まった母集合があります。"
PRINT "この中から ｎ 個の標本（数）を無作為に抽出し、 この標本の平均値を求めます。"
PRINT "それを１００回繰り返し、 ｎ 個の標本の平均値が１００個集まった集合を作ります。"
PRINT "そして、その集合の平均値と標準偏差を求めます。"
PRINT
RANDOMIZE
DEF nrandom=RND+RND+RND+RND+RND+RND+RND+RND+RND+RND+RND+RND-6
LET width=20
INPUT PROMPT "標本数？ （ 2 〜 1000 ）": n
DIM Q(100,n),R(100,n),U(100),S(100)
FOR i=1 TO 100
   LET U(i)=0
   LET S(i)=0
NEXT i
LET X=0
LET Y=0
PRINT "1回目の結果:"
PRINT "　　";
FOR i=1 TO 100
   FOR j=1 TO n
      LET Q(i,j)=INT(nrandom*20+50)
      IF Q(i,j)<0 THEN LET Q(i,j)=0
      IF Q(i,j)>99 THEN LET Q(i,j)=99
      IF i=1 THEN PRINT Q(i,j);" ";
      LET U(i) = U(i) + Q(i,j)/n
   NEXT j
   FOR j=1 TO n
      LET S(i) = S(i) + (Q(i,j)-U(i))^2/n
   NEXT j
   IF i=1 THEN
      PRINT
      PRINT "　　　平均値；"; ROUND(U(i)*100)/100  
      PRINT "　　標準偏差："; ROUND(SQR(S(i))*100)/100
      PRINT
   END IF
NEXT i
PRINT "標本の平均の集合（１００個）の要素たち："
PRINT "　";
FOR i=1 TO 100
   PRINT ROUND(U(i)*100)/100;"　";
NEXT i
FOR i=1 TO 100
   LET X = X + U(i)/100
NEXT i
FOR i=1 TO 100
   LET Y = Y + (U(i)-X)^2/100
NEXT i
PRINT
PRINT
PRINT "標本平均の集合の平均値：　";ROUND(X*100)/100
PRINT "標本平均の集合の平均値（理論値）＝ 50"
PRINT
PRINT "標本平均の集合の標準偏差：";ROUND(SQR(Y)*100)/100
PRINT "標本平均の集合の標準偏差（理論値）：";ROUND(20/SQR(n)*100)/100
PRINT
END