2023-12-16 21:33:55 -08:00
|
|
|
numfmt 5 0
|
|
|
|
|
e = 2.7182818284590452354
|
2023-07-01 11:58:00 -04:00
|
|
|
for i = 1 to 1000
|
2023-12-16 21:33:55 -08:00
|
|
|
a[] &= 1 + 0.5 * sqrt (-2 * log10 randomf / log10 e) * cos (360 * randomf)
|
2023-07-01 11:58:00 -04:00
|
|
|
.
|
2023-12-16 21:33:55 -08:00
|
|
|
for v in a[]
|
|
|
|
|
avg += v / len a[]
|
|
|
|
|
.
|
|
|
|
|
print "Average: " & avg
|
|
|
|
|
for v in a[]
|
|
|
|
|
s += pow (v - avg) 2
|
|
|
|
|
.
|
|
|
|
|
s = sqrt (s / len a[])
|
|
|
|
|
print "Std deviation: " & s
|
