diff --git a/include/util/math.hpp b/include/util/math.hpp
index 1d68856..a69b575 100644
--- a/include/util/math.hpp
+++ b/include/util/math.hpp
@@ -103,8 +103,9 @@ void reservoir_sampling_A_ExpJ(
 
 /**
  * @brief Power law distribution for random numbers.
- * A continuous random distribution on the range [eps, infty) with equal
- * distribution
+ * A continuous random distribution on the range [eps, infty)
+ * with p(x) ~ x^(-gamma)
+ * Including normalization the PDF is
  * p(x) = (1-gamma)/(1-eps^(1-gamma)) * x^(-gamma)
  */
 template<typename ScalarT = double>