diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index efe421eb1c565f77d20bf2a120cb103709ba1d9d..cf4ef4416d09e248d4b7710d2f5cfd0a181fd71d 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -74,16 +74,16 @@ model = Ruina
 [boundary.friction.state]
 evolve = true
 # log(alpha(0))
-initial = 0
+initial = -2.42037
 
 [boundary.friction.ruina]
 # "For rocks, typical values of A and B range from 0.005 to 0.015"
 # -- Ronaldo I. Borja and Craig D. Foster:
 # Continuum mathematical modeling of slip weakening in geological systems
-a = 0.015
-b = 0.005
+a = 0.0075
+b = 0.015
 # If L is too large, reaching a steady state takes too long
-L = 1e-4
+L = 1e-3
 
 # The second-order a and b terms in [..] are very important in
 # determining whether sliding is stable or unstable and (in appli-
diff --git a/src/one-body-sample.py b/src/one-body-sample.py
index 81d09658a8eebbe6a01479ab4805917c7627e848..3ab9b22d2ae28a6bc1d72616b751d671e4d2da1a 100644
--- a/src/one-body-sample.py
+++ b/src/one-body-sample.py
@@ -14,9 +14,9 @@ class neumannCondition:
 class dirichletCondition:
     def __call__(self, x):
         # return 0
-        fst = 5e-4
-        snd = 2e-4
-        trd = 5e-4
+        fst = 3e-4
+        snd = 5e-4
+        trd = 3e-4
         if x < 1.0/3:
             return fst
         elif x < 2.0/3: