diff --git a/src/one-body-sample.cc b/src/one-body-sample.cc
index cfe1c64f4759ab48e5964739480e360c0dc90f10..faf33e228024fede19d29b77576785bfd17720d9 100644
--- a/src/one-body-sample.cc
+++ b/src/one-body-sample.cc
@@ -7,6 +7,8 @@
 #include <exception>
 #include <iostream>
 
+#include <dune/common/shared_ptr.hh>
+
 #include <dune/common/exceptions.hh>
 #include <dune/common/stdstreams.hh>
 
@@ -26,6 +28,7 @@
 #include <dune/fufem/boundarypatch.hh>
 #include <dune/fufem/functions/constantfunction.hh>
 #include <dune/fufem/functionspacebases/p1nodalbasis.hh>
+#include <dune/grid/io/file/vtk/vtkwriter.hh>
 
 #include <dune/tnnmg/iterationsteps/genericnonlineargs.hh>
 
@@ -34,6 +37,8 @@
 #include <dune/solvers/norms/energynorm.hh>
 #include <dune/solvers/solvers/loopsolver.hh>
 
+#include <dune/vtkgridfunction.hh> // currently homeless
+
 #include <dune/tectonic/globalnonlinearity.hh>
 #include <dune/tectonic/myconvexproblem.hh>
 #include <dune/tectonic/myblockproblem.hh>
@@ -237,6 +242,13 @@ int main() {
       }
     }
 
+    Dune::shared_ptr<Dune::VTKBasisGridFunction<P1Basis, VectorType>> ptr(
+        new Dune::VTKBasisGridFunction<P1Basis, VectorType>(p1Basis, u1,
+                                                            "displacement"));
+    Dune::VTKWriter<GridType::LeafGridView> writer(leafView);
+    writer.addVertexData(ptr);
+    writer.write("weeee");
+
     VectorType diff = u2;
     diff -= u1;
     std::cout << "Infinity norm of the difference of the two solutions: "