diff --git a/libs/random/include/psemek/random/normal.hpp b/libs/random/include/psemek/random/normal.hpp index 65d8626a..cb761fb1 100644 --- a/libs/random/include/psemek/random/normal.hpp +++ b/libs/random/include/psemek/random/normal.hpp @@ -64,8 +64,8 @@ namespace psemek::random T const r = std::sqrt(- 2 * std::log(u1)); T const th = 2 * geom::pi * u2; - T const z1 = r * std::cos(th); - T const z2 = r * std::sin(th); + T const z1 = mean_ + stddev_ * r * std::cos(th); + T const z2 = mean_ + stddev_ * r * std::sin(th); cached_ = z2; return z1;