diff --git a/src/Makefile.am b/src/Makefile.am
index 22912c6c7c5a3e8eec08ac347e7aa0d1ca1659fd..da531dfbaaa8b1e89a32232657ee9f21f12af376 100644
--- a/src/Makefile.am
+++ b/src/Makefile.am
@@ -18,9 +18,12 @@ cosserat_continuum_LDFLAGS  = $(UG_LDFLAGS) $(IPOPT_LDFLAGS) \
                               $(ADOLC_LDFLAGS) $(PYTHON_LDFLAGS)
 
 harmonicmaps_SOURCES = harmonicmaps.cc
-harmonicmaps_CXXFLAGS = $(UG_CPPFLAGS) $(AMIRAMESH_CPPFLAGS) $(IPOPT_CPPFLAGS) $(PSURFACE_CPPFLAGS)
+harmonicmaps_CXXFLAGS = $(UG_CPPFLAGS) $(AMIRAMESH_CPPFLAGS) $(IPOPT_CPPFLAGS) $(PSURFACE_CPPFLAGS) \
+                        $(ADOLC_CPPFLAGS) $(PYTHON_CPPFLAGS)
 harmonicmaps_LDADD    = $(UG_LDFLAGS) $(AMIRAMESH_LDFLAGS) $(UG_LIBS) $(AMIRAMESH_LIBS) \
-                          $(IPOPT_LDFLAGS) $(IPOPT_LIBS) $(PSURFACE_LDFLAGS) $(PSURFACE_LIBS)
+                        $(IPOPT_LDFLAGS) $(IPOPT_LIBS) $(PSURFACE_LDFLAGS) $(PSURFACE_LIBS) \
+                        $(ADOLC_LIBS) $(PYTHON_LIBS) $(ADOLC_LDFLAGS) $(PYTHON_LDFLAGS)
+
 
 rod3d_SOURCES = rod3d.cc
 rod3d_CXXFLAGS = $(IPOPT_CPPFLAGS)
diff --git a/src/harmonicmaps.cc b/src/harmonicmaps.cc
index 19147ba66d8994bd987065b2fdd0ef5b4cc1dc39..9f089b8bc4dc20856b4045b3ef4b930042ef2843 100644
--- a/src/harmonicmaps.cc
+++ b/src/harmonicmaps.cc
@@ -25,6 +25,7 @@
 
 #include <dune/fufem/boundarypatch.hh>
 #include <dune/fufem/functions/vtkbasisgridfunction.hh>
+#include <dune/fufem/dunepython.hh>
 
 #include <dune/solvers/solvers/iterativesolver.hh>
 #include <dune/solvers/norms/energynorm.hh>
@@ -81,6 +82,17 @@ computeEmbeddedDifference(const std::vector<TargetSpace>& a, const std::vector<T
 int main (int argc, char *argv[]) try
 {
     //feenableexcept(FE_INVALID);
+    // Start Python interpreter
+    Python::start();
+    Python::Reference main = Python::import("__main__");
+    Python::run("import math");
+
+    //feenableexcept(FE_INVALID);
+    Python::runStream()
+        << std::endl << "import sys"
+        << std::endl << "sys.path.append('/home/sander/dune/dune-gfe/src')"
+        << std::endl;
+
 
     typedef std::vector<TargetSpace> SolutionType;