diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index dd0b31d8614109152d2bc8f3086ff02b299320ca..b753f0f9b2c7f487fc9fcc8e7e13d1f9fdb265c2 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -452,6 +452,18 @@ int main(int argc, char *argv[]) {
       u4 += u4_diff;
       s4_old = *s4_new;
 
+      if (parset.get<bool>("printEvolution")) {
+        // Print displacement on frictional boundary
+        boost::format const formatter("s[%02d] = %+3e, "
+                                      "%|40t|u[%02d] = %+3e");
+        for (size_t i = 0; i < frictionalNodes.size(); ++i)
+          if (frictionalNodes[i][0]) {
+            std::cout << boost::format(formatter) % i % (*s4_new)[i] % i % u4[i]
+                      << std::endl;
+            break; // Only print the first element
+          }
+      }
+
       // Compute von Mises stress and write everything to a file
       if (parset.get<bool>("writeVTK")) {
         auto const displacement =
diff --git a/src/one-body-sample.parset b/src/one-body-sample.parset
index 9629b741fb5a1486b11f2335c48ffde6d732f73a..b91d6eda2e0485bdf5f891769d4173cdb2b4cf83 100644
--- a/src/one-body-sample.parset
+++ b/src/one-body-sample.parset
@@ -2,6 +2,7 @@
 timesteps = 30
 
 verbose = false
+printEvolution = true
 printFrictionalBoundary = false
 printProgress = false
 printDifference = false