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