diff --git a/ohmpi/sequence.py b/ohmpi/sequence.py
index a5066cd4bd932a8cd8e1385ba7114023d22b372d..683bc0d08114cd01412dfc306e7eb233dd776cd6 100644
--- a/ohmpi/sequence.py
+++ b/ohmpi/sequence.py
@@ -191,7 +191,7 @@ def multigrad(nelec, a, n, s):
     if isinstance(n, int):
         n = np.arange(n) + 1
     if isinstance(s, int):
-        n = np.arange(s) + 1
+        s = np.arange(s) + 1
     elec = np.arange(nelec) + 1
     abmn = []
     for aa in a: