Всегда по другому делал: cos(pi*Rnd2)*sqrt(-ln(1-Rnd1))
В принципе -- то же самое. НО: без "неймановских обрезаний" (так это называлось во времена моего студенчества)

А без него (обрезания) -- может и "на математике" не так уж и медленно получиться ...


------------ Дoбавленo в 12.42:
Комментарий к doRandomNorm -- странный.
И использование Min/Max -- тоже странное... Логично (для этих свойств) было бы что-о типа: Центр-Сигма/Центр+Сигма