diff --git a/dune/tectonic/nicefunction.hh b/dune/tectonic/nicefunction.hh
index 81625152d6426a2e666be05d7c8fd6f731818e46..fdcdc0e4b4f115b9b475a0c4a3ecdc1a63ff0323 100644
--- a/dune/tectonic/nicefunction.hh
+++ b/dune/tectonic/nicefunction.hh
@@ -51,18 +51,26 @@ class RuinaFunction : public NiceFunction {
     y *= coefficient;
   }
 
+  /*
+    (leaving some terms aside): with s > rho
+
+    1/eta d/dx [ a * (s log s - s) + mu s ] where s = eta x
+    = 1/eta [ a * (log (eta x) * eta) + eta mu ]
+    = a * log(eta x) + mu
+  */
   double virtual leftDifferential(double s) const {
     if (eta * s <= rho)
       return 0;
 
-    return coefficient * normalStress * (a * std::log(s) + mu);
+    return coefficient * normalStress * (a * std::log(eta * s) + mu);
   }
 
+  /* see above */
   double virtual rightDifferential(double s) const {
     if (eta * s <= rho)
       return 0;
 
-    return coefficient * normalStress * (a * std::log(s) + mu);
+    return coefficient * normalStress * (a * std::log(eta * s) + mu);
   }
 
 private: