RosettaCodeData/Task/Random-numbers/Phix/random-numbers.phix
2026-02-01 16:33:20 -08:00

12 lines
301 B
Text

with javascript_semantics
requires("1.0.6") -- std_dev
function random_normal()
return sqrt(-2*log(rnd())) * cos(2*PI*rnd())
end function
sequence s = repeat(0,1000)
for i=1 to length(s) do
s[i] = 1 + 0.5 * random_normal()
end for
printf(1,"mean: %g, std dev: %g\n",{average(s),std_dev(s)})