package require Tcl 8.5 variable ::pi [expr acos(0)] proc ::tcl::mathfunc::nrand {} { expr {sqrt(-2*log(rand())) * cos(2*$::pi*rand())} } set mean 1.0 set stddev 0.5 for {set i 0} {$i < 1000} {incr i} { lappend result [expr {$mean + $stddev*nrand()}] }