diff --git a/meshpy/triangle.py b/meshpy/triangle.py
index 1e6e9e6..0c0400a 100644
--- a/meshpy/triangle.py
+++ b/meshpy/triangle.py
@@ -61,7 +61,15 @@ def set_facets(self, facets, facet_markers=None):
     def dump(self):
         for name in self._constituents:
             dump_array(name, getattr(self, name))
-
+    
+    def write_vtk(self, filename):
+        import pyvtk
+        vtkelements = pyvtk.VtkData(
+            pyvtk.UnstructuredGrid(
+                self.points,
+                triangle=self.elements),
+            "Mesh")
+        vtkelements.tofile(filename)