diff --git a/dune/microstructure/energies/discretekirchhoffbendingenergyprestrained.hh b/dune/microstructure/energies/discretekirchhoffbendingenergyprestrained.hh
index a512bbd20397f209da7875ca2b8917871623fdc2..def8d5b91d5cca4811bb5d01a7e8e5775ed875b4 100644
--- a/dune/microstructure/energies/discretekirchhoffbendingenergyprestrained.hh
+++ b/dune/microstructure/energies/discretekirchhoffbendingenergyprestrained.hh
@@ -592,19 +592,28 @@ namespace Dune::GFE
        */
       RT harmonicEnergy = 0;
 
-      // Gauss-Quadrature:
+      // Gauss-Legendre-Quadrature (default):
       const auto &quadRule = QuadratureRules<double, gridDim>::rule(lagrangeLFE_.type(), quadOrder);
+      //TEST Gauß-Jacobi-Quadrature formula: 
+      // const auto &quadRule = QuadratureRules<double, gridDim>::rule(lagrangeLFE_.type(), quadOrder, QuadratureType::GaussJacobi_1_0);
+
+
       // Trapezoidal-rule:
       // std::vector<Dune::QuadraturePoint<double,2>> quadRule = { {{0.0,0.0}, 1.0/6.0}, {{1.0,0.0}, 1.0/6.0}, {{0.0,1.0}, 1.0/6.0} };
       // Trapezoidal-rule + evaluation on edge centers (only for testing purposes):
       // std::vector<Dune::QuadraturePoint<double,2>> quadRule = { {{0.0,0.0}, 1.0/6.0}, {{0.5,0.0}, 1.0/6.0}, {{1.0,0.0}, 1.0/6.0}, {{0.0,0.5}, 1.0/6.0}, {{0.5,0.5}, 1.0/6.0}, {{0.0,1.0}, 1.0/6.0} };
       // int count = 0;
-      // std::cout << "Number of quadrature points: " << quadRule.size() << std::endl;
+
+
+      std::cout << "Number of quadrature points: " << quadRule.size() << std::endl;
+
+
       for (auto&& quadPoint : quadRule)
       {
         const auto& quadPos = quadPoint.position();
         const auto integrationElement = geometry.integrationElement(quadPos);
 
+        std::cout << "quadPos:" << quadPos << std::endl;
 
         // exit(0);