diff --git a/viscoelast.cc b/viscoelast.cc
index 6770420a2dac426ea1f72d96383cfa7b60d10ea6..2bd917464c6a094fe643da8e5b4d53339c4b924b 100644
--- a/viscoelast.cc
+++ b/viscoelast.cc
@@ -29,6 +29,8 @@
 #include <dune/fufem/functiontools/gridfunctionadaptor.hh>
 #include <dune/fufem/utilities/dirichletbcassembler.hh>
 
+#include <dune/functions/functionspacebases/lagrangebasis.hh>
+
 
 /** Solver for a quasi-static (viscoelastic) Kelvin-Voigt material
  *
@@ -152,7 +154,7 @@ int main (int argc, char *argv[]) try {
 
     MatrixType elasticPart,viscousPart;
 
-    using DuneP1Basis  = Dune::Functions::PQkNodalBasis<typename GridType::LeafGridView, 1>;
+    using DuneP1Basis  = Dune::Functions::LagrangeBasis<typename GridType::LeafGridView, 1>;
     using P1Basis = DuneFunctionsBasis<DuneP1Basis>;
     P1Basis p1Basis(grid->leafGridView());
     OperatorAssembler<P1Basis,P1Basis> globalAssembler(p1Basis,p1Basis);