diff --git a/gtsam/nonlinear/NonlinearFactorGraph.cpp b/gtsam/nonlinear/NonlinearFactorGraph.cpp
index 0b9ae35d5..3d51ef3ac 100644
--- a/gtsam/nonlinear/NonlinearFactorGraph.cpp
+++ b/gtsam/nonlinear/NonlinearFactorGraph.cpp
@@ -42,6 +42,30 @@ void NonlinearFactorGraph::print(const std::string& str, const KeyFormatter& key
   }
 }
 
+/* ************************************************************************* */
+void NonlinearFactorGraph::saveGraph(std::ostream &stm, const KeyFormatter& keyFormatter) const {
+  stm << "graph {\n";
+
+  // Create nodes for each variable in the graph
+  BOOST_FOREACH(Key key, this->keys()) {
+    // Label the node with the label from the KeyFormatter
+    stm << "  var" << key << "[label=\"" << keyFormatter(key) << "\"];\n"; }
+  stm << "\n";
+
+  // Create factors and variable connections
+  for(size_t i = 0; i < this->size(); ++i) {
+    // Make each factor a dot
+    stm << "  factor" << i << "[label=\"\", shape=point];\n";
+
+    // Make factor-variable connections
+    if(this->at(i)) {
+      BOOST_FOREACH(Key key, *this->at(i)) {
+        stm << "  var" << key << "--" << "factor" << i << ";\n"; } }
+  }
+
+  stm << "}\n";
+}
+
 /* ************************************************************************* */
 double NonlinearFactorGraph::error(const Values& c) const {
   gttic(NonlinearFactorGraph_error);
diff --git a/gtsam/nonlinear/NonlinearFactorGraph.h b/gtsam/nonlinear/NonlinearFactorGraph.h
index 3b5f38eb3..945189a41 100644
--- a/gtsam/nonlinear/NonlinearFactorGraph.h
+++ b/gtsam/nonlinear/NonlinearFactorGraph.h
@@ -46,6 +46,9 @@ namespace gtsam {
     /** print just calls base class */
     void print(const std::string& str = "NonlinearFactorGraph: ", const KeyFormatter& keyFormatter = DefaultKeyFormatter) const;
 
+    /** Write the graph in GraphViz format for visualization */
+    void saveGraph(std::ostream& stm, const KeyFormatter& keyFormatter = DefaultKeyFormatter) const;
+
     /** return keys as an ordered set - ordering is by key value */
     FastSet<Key> keys() const;