39 lines
942 B
Plaintext
39 lines
942 B
Plaintext
module statistics1 {
|
|
open "out1.txt" for output as #f
|
|
sample(100)
|
|
sample(1000)
|
|
sample(10000)
|
|
sample(100000)
|
|
sample(1000000)
|
|
close #f
|
|
Win "NotePad", dir$+"out1.txt"
|
|
Sub sample(n)
|
|
local samp(1 to n) as double<<rnd
|
|
|
|
' calculate mean, standard deviation
|
|
local sum, sumSq, i, max, mean, bins, b
|
|
for i = 1 to n
|
|
sum += samp(i)
|
|
sumSq += samp(i)^2
|
|
next i
|
|
print #f, n; " Samples used."
|
|
mean = sum / n
|
|
print #f, "Mean = "; mean
|
|
|
|
print #f, "Std Dev = "; (sumSq /n -mean^2)^0.5
|
|
|
|
'------- Show histogram
|
|
bins = 10
|
|
dim bins(bins)
|
|
for i = 1 to n
|
|
bins(int(bins * samp(i)))++
|
|
next i
|
|
max=bins()#max()
|
|
for b = 0 to bins -1
|
|
print #f, b;" ";string$("*",bins(b)/max*60);" ";bins(b)/n*100;"%"
|
|
next b
|
|
print #f
|
|
end sub
|
|
}
|
|
statistics1
|