RosettaCodeData/Task/Random-numbers/Nim/random-numbers.nim

15 lines
375 B
Nim

import math, strutils
const precisn = 5
var rs: TRunningStat
proc normGauss: float {.inline.} = 1 + 0.76 * cos(2*PI*random(1.0)) * sqrt(-2*log10(random(1.0)))
randomize()
for j in 0..5:
for i in 0..1000:
rs.push(normGauss())
echo("mean: ", $formatFloat(rs.mean,ffDecimal,precisn),
" stdDev: ", $formatFloat(rs.standardDeviation(),ffDecimal,precisn))