diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc
index 30be33a4f3970a16bb6227787c0407e79e9cc8b9..c611292c37379d3c68169c4ad0ea6f06cee5f66f 100644
--- a/cosserat-continuum.cc
+++ b/cosserat-continuum.cc
@@ -1,5 +1,7 @@
 #include <config.h>
 
+#define SECOND_ORDER
+
 #include <fenv.h>
 
 // Includes for the ADOL-C automatic differentiation library
@@ -25,6 +27,7 @@
 #include <dune/fufem/functiontools/boundarydofs.hh>
 #include <dune/fufem/functiontools/basisinterpolator.hh>
 #include <dune/fufem/functionspacebases/p1nodalbasis.hh>
+#include <dune/fufem/functionspacebases/p2nodalbasis.hh>
 
 #include <dune/solvers/solvers/iterativesolver.hh>
 #include <dune/solvers/norms/energynorm.hh>
@@ -191,7 +194,11 @@ int main (int argc, char *argv[]) try
     typedef GridType::LeafGridView GridView;
     GridView gridView = grid->leafGridView();
 
+#ifdef SECOND_ORDER
+    typedef P2NodalBasis<GridView,double> FEBasis;
+#else
     typedef P1NodalBasis<GridView,double> FEBasis;
+#endif
     FEBasis feBasis(gridView);
 
     // /////////////////////////////////////////