diff --git a/figures/SchueckerSchmidt2017/utils.py b/figures/SchueckerSchmidt2017/utils.py
index ef8e03affc1f325ed928c8e5fe9466d5f79fe10e..2fde31d21d4352bad7b9398c2d310a4d5f539d83 100644
--- a/figures/SchueckerSchmidt2017/utils.py
+++ b/figures/SchueckerSchmidt2017/utils.py
@@ -113,11 +113,11 @@ def compute_iteration(max_iter, fac_nu_ext_5E_list, theory_params, M_base, threa
         r = iteration_results(fac_nu_ext_5E_list, M_base, threads=threads)
         results = np.vstack((results, r))
         par_list = np.append(par_list, fac_nu_ext_5E_list)
-        i = np.argmax(np.diff(np.mean(r, axis=1)[:, -1]))
+        j = np.argmax(np.diff(np.mean(r, axis=1)[:, -1]))
         i += 1
-        fac_nu_ext_5E_list = np.arange(fac_nu_ext_5E_list[i],
+        fac_nu_ext_5E_list = np.arange(fac_nu_ext_5E_list[j],
                                        # to ensure that the array includes the last value, we add a small epsilon
-                                       fac_nu_ext_5E_list[i+1] + 1.e-10,
+                                       fac_nu_ext_5E_list[j+1] + 1.e-10,
                                        10**(-(i+2.)))
 
     ind = np.argsort(par_list)