27 lines
611 B
Plaintext
27 lines
611 B
Plaintext
# Generates normally distributed random numbers with mean 0 and standard deviation 1
|
|
function randomNormal()
|
|
return cos(2.0 * pi * rand) * sqr(-2.0 * log(rand))
|
|
end function
|
|
|
|
dim r(1000)
|
|
sum = 0.0
|
|
# Generate 1000 normally distributed random numbers
|
|
# with mean 1 and standard deviation 0.5
|
|
# and calculate their sum
|
|
for i = 0 to 999
|
|
r[i] = 1.0 + randomNormal() / 2.0
|
|
sum += r[i]
|
|
next i
|
|
|
|
mean = sum / 1000.0
|
|
sum = 0.0
|
|
# Now calculate their standard deviation
|
|
for i = 0 to 999
|
|
sum += (r[i] - mean) ^ 2.0
|
|
next i
|
|
sd = sqr(sum/1000.0)
|
|
|
|
print "Mean is "; mean
|
|
print "Standard Deviation is "; sd
|
|
end
|