diff --git a/2bmgnonlincontact.cc b/2bmgnonlincontact.cc
index 482943cd2d4435fdd699a48e39a8d0bc3890e66b..26debacef27c384154820a16a8729aeff83b2aff 100644
--- a/2bmgnonlincontact.cc
+++ b/2bmgnonlincontact.cc
@@ -27,7 +27,7 @@
 #include <dune/fufem/utilities/dirichletbcassembler.hh>
 #include <dune/fufem/improvegrid.hh>
 
-#include <dune/functions/functionspacebases/pqknodalbasis.hh>
+#include <dune/functions/functionspacebases/lagrangebasis.hh>
 
 #include <dune/solvers/norms/energynorm.hh>
 
@@ -237,7 +237,7 @@ int main (int argc, char *argv[]) try
     }
 
     // Create the materials
-    using P1Basis =  DuneFunctionsBasis<typename Dune::Functions::PQkNodalBasis<typename GridType::LeafGridView, 1> >;
+    using P1Basis =  DuneFunctionsBasis<typename Dune::Functions::LagrangeBasis<typename GridType::LeafGridView, 1> >;
     vector<std::shared_ptr<P1Basis> > p1Basis(nGrids);
 
     // using MaterialType = GeomExactStVenantMaterial<P1Basis>;
@@ -351,7 +351,7 @@ int main (int argc, char *argv[]) try
     // make contact problem for P2 functions
 
     // Create the materials
-    //using P2Basis =  DuneFunctionsBasis<typename Dune::Functions::PQkNodalBasis<typename GridType::LeafGridView, 2> >;
+    //using P2Basis =  DuneFunctionsBasis<typename Dune::Functions::LagrangeBasis<typename GridType::LeafGridView, 2> >;
     using P2Basis =  P2NodalBasis<typename GridType::LeafGridView,field_type>;
     vector<std::shared_ptr<P2Basis> > p2Bases(nGrids);
     for (size_t j=0; j<nGrids; j++)