diff --git a/assets/example.png b/assets/example.png
new file mode 100644
index 0000000..38eca55
Binary files /dev/null and b/assets/example.png differ
diff --git a/readme.md b/readme.md
index 3d70f6e..e8c6290 100644
--- a/readme.md
+++ b/readme.md
@@ -25,3 +25,9 @@ with this solver is the convergence rate. We have built this as a feasibility st
 made to optimize the solver. Specifically we use a seriese of nested picard iterations to solve the nonlinear system. These do
 seem to converge reliably; however, they do so very slowly. 
 
+
+## Example
+A 2D slice of the solution for a polytropic index of 1.5 and rotating just below the critical rotation rate is shown below. The color shows
+the effective potential. This visualization was made with GLVis.
+
+![Example solution](./assets/example.png)