diff --git a/inputs/computeMuGamma.parset b/inputs/computeMuGamma.parset
index 6acf41019269f9686a55716591a4c6c2606d09f6..35fdbeeeab3534d346179fd2809a65fb79308328 100644
--- a/inputs/computeMuGamma.parset
+++ b/inputs/computeMuGamma.parset
@@ -29,7 +29,7 @@
 
 nElements = 32 32
 
-gamma=0.5
+gamma=100
 
 
 #############################################
@@ -40,15 +40,16 @@ write_materialFunctions = true   # VTK mu-functions , lambda-functions
 
 beta=2.0
 
-mu1=10.0
-
+mu1=1.0
 
+lambda1=0.0
 #lambda1=10.0
 
 
 ####  material_implementation("analytical_Example") ? 
 
-material_prestrain_imp= "analytical_Example"
+material_prestrain_imp= "parametrized_Laminate"
+#material_prestrain_imp= "analytical_Example"
 #material_prestrain_imp ="isotropic_bilayer"
 
 
@@ -63,9 +64,9 @@ write_prestrainFunctions = true  # VTK norm of B ,
 
 
 rho1=1.0
-alpha=9.764705882352942
+alpha=2.0
 
-theta=0.15000000000000002
+theta=0.541859296482412
 
 #theta = 0.3    # volume fraction   #default = 1.0/4.0
 #theta = 0.25   # volume fraction
@@ -111,12 +112,3 @@ Solvertype = 1
 
 
 
-)*mu2)
-
-
-
-
-
-
-
-
diff --git a/src/Compute_MuGamma.cc b/src/Compute_MuGamma.cc
index d59db32e028de4796344734f056a71ddec31340a..94ab7e364ef633be51caf940a015d428642f69ec 100644
--- a/src/Compute_MuGamma.cc
+++ b/src/Compute_MuGamma.cc
@@ -744,7 +744,7 @@ int main(int argc, char *argv[])
     using Domain = GridView::Codim<0>::Geometry::GlobalCoordinate;
     
     // ----------- INPUT PARAMETERS -----------------------------
-    std::string imp = parameterSet.get<std::string>("material_prestrain_imp", "analytical_Example");
+    std::string imp = parameterSet.get<std::string>("material_prestrain_imp", "parametrized_Laminate");
     log << "material_prestrain used: "<< imp  << std::endl;
     double gamma   = parameterSet.get<double>("gamma", 1.0);   
     double theta   = parameterSet.get<double>("theta", 1.0/4.0);