diff --git a/graph_tools/base.py b/graph_tools/base.py
index 7d5f30cf96e8c29ba88123191a9777f0519925a0..a6e6088263a5a63291e23d1e6aeb8896a54d7252 100644
--- a/graph_tools/base.py
+++ b/graph_tools/base.py
@@ -238,7 +238,7 @@ def _init_():
     def graph(self, n):
       return self._graph_cache.get(n, lambda: self._make_graph(self.gadget(n)))
 
-    def stabilize(self, parameter, start_at=None, jordan=True):
+    def stabilize(self, parameter, start_at=None, jordan=False):
       from sage.all import matrix, QQ, binomial, SR
 
       if start_at is None: start_at = self.sequence_start
diff --git a/graph_tools/tests.py b/graph_tools/tests.py
index 98aa5b6f544858cf207c34f0729b86fcee995c62..f00843a67d5d7b4716123e2fecfff619113af157 100644
--- a/graph_tools/tests.py
+++ b/graph_tools/tests.py
@@ -20,12 +20,12 @@
 True
 
 
->>> Necklace.stabilize(CircuitDoubleCover)
+>>> Necklace.stabilize(CircuitDoubleCover, jordan=True)
 Loop 1 done - 1 boundaries (0 new)
 {'variables': [((1, 2), (1, 2), None)], 'step_matrix': [4], 'initial_vector': [2], \
 'finalize': [1], 'simplified': ([1], [4], [2]), 'formula': 1/2*4^k, 'formula_if': k - 1 >= 0}
 
->>> f = CyclicLadder().stabilize(CircuitDoubleCover, start_at=1)['formula']; f
+>>> f = CyclicLadder().stabilize(CircuitDoubleCover, start_at=1, jordan=True)['formula']; f
 Loop 1 done - 6 boundaries (6 new)
 Loop 2 done - 13 boundaries (6 new)
 Loop 3 done - 15 boundaries (2 new)
@@ -33,7 +33,7 @@ Loop 4 done - 15 boundaries (0 new)
 1/2*(-1)^k*(k - 1) + 1/6*4^(k - 1) + 3*2^(k - 1) + 1/3*(-2)^(k - 2) - 5/2*k - 3/2
 >>> assert all( f(k=k) == cdc_count(CyclicLadder().graph(k)) for k in range(3, 10) )
 
->>> f = CyclicLadder(crossed=True).stabilize(CircuitDoubleCover, start_at=3)['formula']; f
+>>> f = CyclicLadder(crossed=True).stabilize(CircuitDoubleCover, start_at=3, jordan=True)['formula']; f
 Loop 1 done - 12 boundaries (4 new)
 Loop 2 done - 14 boundaries (0 new)
 -1/2*(-1)^k*k + 1/6*4^(k - 1) + 3*2^(k - 1) + 1/3*(-2)^(k - 2) - 5/2*k