diff --git a/src/elasticity/assemble/IElasticity.cpp b/src/elasticity/assemble/IElasticity.cpp
index 8341f4fa7a4cf49366461437c92b89bf72655343..15e2c91bd6d66c13f9e18200455fe3acf9905d39 100644
--- a/src/elasticity/assemble/IElasticity.cpp
+++ b/src/elasticity/assemble/IElasticity.cpp
@@ -83,13 +83,17 @@ void IElasticity::Finalize(Vectors &vecs) {
   GetInvariant(vecs[0], vecs[2]);
 }
 
-void IElasticity::CalculateAndPrintMinGammaF() const{
+void IElasticity::CalculateAndPrintMinMaxGammaF() const{
   if (gammaEvaluation.empty()) THROW("gammaEvaluation is empty, should not be used!")
+
+  double max_gamma = gammaEvaluation[0][0];
   double min_gamma = gammaEvaluation[0][0];
   for (VectorField v : gammaEvaluation){
     min_gamma = std::min(min_gamma, v[0]);
+    max_gamma = std::max(max_gamma, v[0]);
   }
   vout(0) << "minimal value of gamma_f is " << min_gamma << endl ;
+  vout(0) << "maximal value of gamma_f is " << max_gamma << endl ;
 }
 
 void IElasticity::PointEvaluationGamma(const std::string &varname, const vector<Point> &point)  {