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)})