DECLARE EXTERNAL FUNCTION RNDNormal
DECLARE NUMERIC m,s,i
RANDOMIZE
LET m=100
LET s=40
LET C=0
LET D=0
PRINT 
PRINT "　平均 １００   標準偏差 ４０  の正規分布に従う母集団より"
PRINT "抽出した標本から、母集団の平均を推定します。"
PRINT
INPUT PROMPT "標本の大きさ？ （ 20 〜 100万 ）":n
PRINT
DIM A(n)
FOR i=1 TO n
   LET A(i) = RNDNormal(m,s)
   LET C=C+A(i)
NEXT i
FOR i=1 TO n
   LET D=D+(A(i)-C/n)^2
NEXT i
FOR i=1 TO 20
   PRINT "　";A(i)
NEXT i
LET x1 =  ROUND ( C/n - 1.96 * SQR(D/n) / SQR(n), 4 )
LET x2 =  ROUND ( C/n + 1.96 * SQR(D/n) / SQR(n), 4 )
PRINT
PRINT "　　　以上、最初に抽出された２０個の標本の内容でした。"
PRINT
PRINT "抽出された標本について"
PRINT "　　平均 ： ";ROUND(C/n,4)
PRINT "    標準偏差 ： ";ROUND(SQR(D/n),4)
PRINT
PRINT "　したがって、母集団の平均値は、信頼度 95% で、"
PRINT x1;" から ";x2;" の範囲内にあると推定されます。"
PRINT
END
MERGE "random.lib"