Random numbers: Difference between revisions

Content added Content deleted
(add RPL)
Line 785: Line 785:


<syntaxhighlight lang="text">
<syntaxhighlight lang="text">
numfmt 5 0
e = 2.7182818284590452354
for i = 1 to 1000
for i = 1 to 1000
a[] &= 1 + 0.5 * sqrt (-2 * logn randomf) * cos (360 * randomf)
a[] &= 1 + 0.5 * sqrt (-2 * log10 randomf / log10 e) * cos (360 * randomf)
.
.
print a[]
for v in a[]
avg += v
.
avg = avg / len a[]
print "Average: " & avg
for v in a[]
s += pow (v - avg) 2
.
s = sqrt (s / len a[])
print "Std deviation: " & s
</syntaxhighlight>
</syntaxhighlight>