--- a/pymedcouple
+++ b/pymedcouple
@@ -70,7 +70,7 @@
         else:
             return trial
 
-def medcouple_1d(X, eps1 = 2**-52, eps2 = 2**-1022, debug=False ):
+def medcouple_1d(X, eps1 = 2**-52, eps2 = 2**-1022):
     """Calculates the medcouple robust measure of skewness.
 
     Parameters
@@ -129,7 +129,7 @@
     n_plus = len(Zplus)
     n_minus = len(Zminus)
 
-    def h_kern(i, j, debug=False):
+    def h_kern(i, j):
         """Kernel function h for the medcouple, closing over the values of
         Zplus and Zminus just defined above.
 
@@ -208,7 +208,7 @@
     A = []
     for (i, (l, r)) in enumerate(izip(L, R)):
         for j in xrange(l, r + 1):
-            A.append(h_kern(i, j, debug=debug))
+            A.append(h_kern(i, j))
     A.sort()
     A.reverse()