diff --git a/binarycpython/utils/distribution_functions.py b/binarycpython/utils/distribution_functions.py
index 63459057a0b48695ac0b95df0e4ce767769dd721..0ef0c158a9aa9972495f99d5f4087cbdfdd0f425 100644
--- a/binarycpython/utils/distribution_functions.py
+++ b/binarycpython/utils/distribution_functions.py
@@ -394,11 +394,10 @@ def gaussian(
         # normalize over given range
         # TODO: add loading into global var
         normalisation = gaussian_normalizing_const(mean, sigma, gmin, gmax)
-        prob = normalisation * gaussian_func(x, mean, sigma)
+        prob = gaussian_func(x, mean, sigma)/normalisation
 
     return prob
 
-
 #####
 # Mass distributions
 #####