diff --git a/figures/SchueckerSchmidt2017/stabilization.py b/figures/SchueckerSchmidt2017/stabilization.py
index 1e8047aecdd0936b60cc86d2f6df38034cdb541b..dc6f6820c543dccd1415bce6a797f4e18b722b74 100644
--- a/figures/SchueckerSchmidt2017/stabilization.py
+++ b/figures/SchueckerSchmidt2017/stabilization.py
@@ -32,12 +32,12 @@ M_target = MultiAreaModel(network_params_target, theory=True,
                           theory_spec=theory_params)
 
 THREADS = 4
-load_list = [1, 2, 3, 4]
+load_list = []
 
 # This list defines which of the detected minima of the velocity
 # vector is identified as the unstable fixed point. It has to be
 # created manually.
-ind_list = [1, 1, 0, 1]
+ind_list = []
 
 """
 Main loop
diff --git a/figures/SchueckerSchmidt2017/utils.py b/figures/SchueckerSchmidt2017/utils.py
index 556bd576a20b18cccfc7b186be5815a31e7be6c6..ef8e03affc1f325ed928c8e5fe9466d5f79fe10e 100644
--- a/figures/SchueckerSchmidt2017/utils.py
+++ b/figures/SchueckerSchmidt2017/utils.py
@@ -118,7 +118,7 @@ def compute_iteration(max_iter, fac_nu_ext_5E_list, theory_params, M_base, threa
         fac_nu_ext_5E_list = np.arange(fac_nu_ext_5E_list[i],
                                        # to ensure that the array includes the last value, we add a small epsilon
                                        fac_nu_ext_5E_list[i+1] + 1.e-10,
-                                       10**(-(i+2)))
+                                       10**(-(i+2.)))
 
     ind = np.argsort(par_list)
     par_list = par_list[ind]