diff --git a/Framework/CurveFitting/src/Functions/DynamicKuboToyabe.cpp b/Framework/CurveFitting/src/Functions/DynamicKuboToyabe.cpp
index 61c1d0c71376fa5212974c4a0dee0f95e594684c..1bef62eaf54b14fe5b1c1511662f309dc173408f 100644
--- a/Framework/CurveFitting/src/Functions/DynamicKuboToyabe.cpp
+++ b/Framework/CurveFitting/src/Functions/DynamicKuboToyabe.cpp
@@ -128,8 +128,8 @@ double integral(double func(const double, const double, const double),
 
 // f1: function to integrate
 double f1(const double x, const double G, const double w0) {
-	// G = Delta in doc
-	// x = dummy time variable 
+  // G = Delta in doc
+  // x = dummy time variable
   return (exp(-G * G * x * x / 2) * sin(w0 * x));
 }
 
@@ -143,7 +143,7 @@ double ZFKT(const double x, const double G) {
 
 // Static non-zero field Kubo Toyabe relaxation function
 double HKT(const double x, const double G, const double F) {
-	// q = Delta^2 t^2 in doc
+  // q = Delta^2 t^2 in doc
   const double q = G * G * x * x;
   // Muon gyromagnetic ratio * 2 * PI
   const double gm = 2 * M_PI * PhysicalConstants::MuonGyromagneticRatio;
@@ -172,10 +172,10 @@ double HKT(const double x, const double G, const double F) {
       (1 - 2 * r * (1 - exp(-q / 2) * cos(w * x)) + 2 * r * r * w * ig);
 
   if (F > 2 * G) {
-	// longitudinal Gaussian field 
+    // longitudinal Gaussian field
     return ktb;
   } else {
-	 // 
+    //
     const double kz = ZFKT(x, G);
     return kz + F / 2 / G * (ktb - kz);
   }
@@ -228,7 +228,7 @@ double DynamicKuboToyabe::getDKT(double t, double G, double F, double v,
     // Generate dynamic Kubo Toyabe
     for (int k = 0; k < tsmax; k++) {
       double y = gStat[k];
-	  // do integration
+      // do integration
       for (int j = k - 1; j > 0; j--) {
         y = y * (1 - hop) + hop * gDyn[k - j] * gStat[j];
       }