diff --git a/numerical_methods/fast_fourier_transform.cpp b/numerical_methods/fast_fourier_transform.cpp
index 473cdbd2d..198f68f19 100644
--- a/numerical_methods/fast_fourier_transform.cpp
+++ b/numerical_methods/fast_fourier_transform.cpp
@@ -61,7 +61,8 @@ std::complex<double>* FastFourierTransform(std::complex<double>*p,uint64_t  n)
 		y[i]=ye[i]+pow(om,i)*yo[i]; ///Updating the first n/2 elements
 		y[i+n/2]=ye[i]-pow(om,i)*yo[i];///Updating the last n/2 elements
 	}
-
+    free(ye);
+    free(yo);
 	return y;///Returns the list 
 	
 }
@@ -112,6 +113,8 @@ static void test() {
         o2++;
     }
     
+    free(o1);
+    free(o2);
     
 }