diff --git a/README.md b/README.md
index e9b279e057beb329c6f61366e02a40539566f1b6..314ad0f977422481982f02135224a74aba2edc34 100644
--- a/README.md
+++ b/README.md
@@ -69,8 +69,9 @@ The output argument `V` is a matrix with dimensions compatible with `x1vec`and `
 `contour(x1vec, x2vec, V, [0 0])` to plot the bound $`\partial \Omega`$.
 
 ## Plotting during Simulation
-With the input argument `plot_progress = 1`, a surface plot of $`V`$ and a plot of $`\partial \Omega`$ will be updated every 10th iteration of the simulation.
+With the input argument `plot_progress = 1`, a surface plot of $`V`$ and a plot of $`\partial \Omega`$ will be updated every 10th iteration of the simulation. An example is shown below:
 
+![](images/plot_progress_example.png)
 
 ## Printing during Simulation
 With the input argument `print_progress = 1`, information of every 10th iteration will be printed to the command window during the simulation. An example is shown below: