diff --git a/doc/source/code/figure/closed_system.py.diff b/doc/source/code/figure/closed_system.py.diff
index d70b1cc7e649a6b901afc56ba520322e1bd2c00b..1eb0080c95e479086fd7861db1399650763345cb 100644
--- a/doc/source/code/figure/closed_system.py.diff
+++ b/doc/source/code/figure/closed_system.py.diff
@@ -71,7 +71,7 @@
          ham_mat = syst.hamiltonian_submatrix(args=[B], sparse=True)
  
          # we only calculate the 15 lowest eigenvalues
-         ev = sla.eigsh(ham_mat, k=15, which='SM', return_eigenvectors=False)
+         ev = sla.eigsh(ham_mat, k=15, sigma=0, return_eigenvectors=False)
  
          energies.append(ev)