23 lines
447 B
Plaintext
23 lines
447 B
Plaintext
10 ' Random numbers
|
|
20 randomize timer
|
|
30 dim r(999)
|
|
40 sum = 0
|
|
50 for i = 0 to 999
|
|
60 r(i) = 1+randomnormal()/2
|
|
70 sum = sum+r(i)
|
|
80 next
|
|
90 mean = sum/1000
|
|
100 sum = 0
|
|
110 for i = 0 to 999
|
|
120 sum = sum+(r(i)-mean)^2
|
|
130 next
|
|
140 sd = sqr(sum/1000)
|
|
150 print "Mean is ";mean
|
|
160 print "Standard Deviation is ";sd
|
|
170 print
|
|
180 end
|
|
|
|
500 sub randomnormal()
|
|
510 randomnormal = cos(2*pi*rnd(1))*sqr(-2*log(rnd(1)))
|
|
520 end sub
|