diff --git a/cpp/full_code/ornstein_uhlenbeck_process.hpp b/cpp/full_code/ornstein_uhlenbeck_process.hpp
index 2d264121d3f95613625b22de00eccafbdcd7cf2b..750c3ccc358970047e55038de435036a48796874 100644
--- a/cpp/full_code/ornstein_uhlenbeck_process.hpp
+++ b/cpp/full_code/ornstein_uhlenbeck_process.hpp
@@ -52,10 +52,8 @@ class ornstein_uhlenbeck_process{
 
         inline double gamma(double kabs)
         {
-            return this->kolmogorov_constant*sqrt(this->kolmogorov_constant/
+            return pow(kabs,2./3.)*this->kolmogorov_constant*sqrt(this->kolmogorov_constant/
                            this->ou_energy_amplitude);
-            //return 1./10.*pow(kabs,2./3.)*this->kolmogorov_constant*sqrt(this->kolmogorov_constant/
-            //               this->ou_energy_amplitude); //TODO take out 10
         }