RosettaCodeData/Task/Random-numbers/Zkl/random-numbers-1.zkl

6 lines
276 B
Text
Raw Permalink Normal View History

2017-09-23 10:01:46 +02:00
fcn mkRand(mean,sd){ //normally distributed random w/mean & standard deviation
pi:=(0.0).pi; // using the BoxMuller transform
rz1:=fcn{1.0-(0.0).random(1)} // from [0,1) to (0,1]
return('wrap(){((-2.0*rz1().log()).sqrt() * (2.0*pi*rz1()).cos())*sd + mean })
}