diff --git a/cpp/full_code/ornstein_uhlenbeck_process.hpp b/cpp/full_code/ornstein_uhlenbeck_process.hpp
index d122677b2e887e45a0061e88162f3dfdd5b97d38..2d264121d3f95613625b22de00eccafbdcd7cf2b 100644
--- a/cpp/full_code/ornstein_uhlenbeck_process.hpp
+++ b/cpp/full_code/ornstein_uhlenbeck_process.hpp
@@ -52,10 +52,10 @@ class ornstein_uhlenbeck_process{
 
         inline double gamma(double kabs)
         {
-            //return 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
+            return 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
         }