diff --git a/cosserat-continuum.cc b/cosserat-continuum.cc
index d4b68675f3d2d2116e6a3bfa33a6511acf3e26fe..1eb2d2650e68a631880c369d9266488d5aced397 100644
--- a/cosserat-continuum.cc
+++ b/cosserat-continuum.cc
@@ -201,6 +201,9 @@ int main (int argc, char *argv[]) try
     //   Main homotopy loop
     ////////////////////////////////////////////////////////
 
+    // Output initial iterate (of homotopy loop)
+    CosseratVTKWriter<GridType>::write(*grid,x, resultPath + "cosserat_homotopy_0");
+
     for (int i=0; i<numHomotopySteps; i++) {
 
         double homotopyParameter = (i+1)*(1.0/numHomotopySteps);
@@ -274,6 +277,11 @@ int main (int argc, char *argv[]) try
 
         x = solver.getSol();
 
+        // Output result of each homotopy step
+        std::stringstream iAsAscii;
+        iAsAscii << i+1;
+        CosseratVTKWriter<GridType>::write(*grid,x, resultPath + "cosserat_homotopy_" + iAsAscii.str());
+
     }
 
     // //////////////////////////////