diff --git a/conmech/solvers/optimization/optimization.py b/conmech/solvers/optimization/optimization.py
index b6ae519a..1e5c4eae 100644
--- a/conmech/solvers/optimization/optimization.py
+++ b/conmech/solvers/optimization/optimization.py
@@ -124,9 +124,11 @@ def _solve_impl(
                 "qsmlm",
             ):
                 # pylint: disable=import-outside-toplevel,import-error)
-                from kosopt import qsmlm
+                from kosopt import qsmlmi
 
-                solution = qsmlm.minimize(self.loss, solution, args=args, maxiter=maxiter)
+                solution = qsmlmi.minimize(
+                    self.loss, solution, args=args, maxiter=maxiter
+                )
             else:
                 result = scipy.optimize.minimize(
                     self.loss,
diff --git a/examples/Bagirov_Bartman_Ochal_2023.py b/examples/Bagirov_Bartman_Ochal_2023.py
index 14f4ba9f..a09bf002 100644
--- a/examples/Bagirov_Bartman_Ochal_2023.py
+++ b/examples/Bagirov_Bartman_Ochal_2023.py
@@ -103,13 +103,15 @@ def main(config: Config):
     mesh_descr = RectangleMeshDescription(
         initial_position=None,
         max_element_perimeter=0.25 * 10 * mm,
-        scale=[2.5 * 10 * mm, 10 * mm],
+        scale=[8 * 10 * mm, 10 * mm],
     )
 
     setup = StaticSetup(mesh_descr=mesh_descr)
 
-    for method in ("Powell", "BFGS", "CG", "qsm"):
-        for force in np.arange(25e3 * kN, 26e3 * kN + 1, 1e3 * kN) * surface:
+    for method in ("Powell", "BFGS", "CG", "qsm")[3:]:
+        for force in (
+            np.asarray([23e3 * kN, 26.2e3 * kN, 27e3 * kN, 30e3 * kN]) * surface
+        ):
 
             def outer_forces(x, t=None):
                 if x[1] >= 0.0099: