Solve a different problem

src/one-body-sample.parset

 # -*- mode:conf -*-
-timesteps = 500
+timesteps = 1000
 
 verbose = false
 printCoefficient = false
@@ -60,7 +60,7 @@ acceptFactor = 1.0
 requiredResidual = 1e-12
 
 [boundary.friction]
-normalstress = 0.1 # laursen depends a lot more on this
+normalstress = 10000000
 # "Most dry materials in combination have friction coefficient values
 #  between 0.3 and 0.6"
 # -- http://en.wikipedia.org/wiki/Friction#Coefficient_of_friction
@@ -84,7 +84,7 @@ model = Dieterich # Ruina
 a = 0.0075
 b = 0.015
 # If L is too large, reaching a steady state takes too long
-L = 1e-3 # 1e-3 for the Dirichlet case, 1e-5 - 1e-7 for the Neumann case
+L = 1e-8 # 1e-3 for the Dirichlet case, 1e-5 - 1e-7 for the Neumann case
 
 # The second-order a and b terms in [..] are very important in
 # determining whether sliding is stable or unstable and (in appli-

src/one-body-sample.py

@@ -13,6 +13,7 @@ class neumannCondition:
 
 class dirichletCondition:
     def __call__(self, x):
+        return 3e-3
         # return 0
         fst = 3e-4
         snd = 5e-4