diff --git a/README.md b/README.md
index 9cf1f721f94854503c4f34ec941f80580247fdc3..a6402da80b197f7b279e23f3c55e0f9e32d25deb 100644
--- a/README.md
+++ b/README.md
@@ -7,6 +7,9 @@ The python routines to reproduce each figure from Ploeckinger et al. ([arXiv](ht
 - numpy
 - unyt
 
+**Additional python packages required for individual plots**
+- astropy (Fig. 4) 
+
 ### Fig. 1: Comparison between Newtonian and softened Jeans mass
 Reproduce Fig. 1 with:
 
@@ -39,8 +42,23 @@ python3 Fig3_runaway_collapse_zone.py
 The plotted resolution parameters can be updated by changing the lines that set `epsarr`, `mBarr`, and `h_min_ratio_arr` (see comments in the python script).  
 
 ### Fig. 4: Runaway collapse zones examples from the literature
+Reproduce Fig. 4 with:
+
+```
+python3 Fig4_runaway_collapse_zone_examples.py 
+```
+
+**Comment:**
+The astropy python package is used to convert the over-density into a physical density for the Springel & Hernquist (2003) effective equation of state. If the relevant lines that add this line to the subplot are commeted out, the astropy package is not necessary. The parameters for some other simulation projects can be found in the python script and can be used for any of the subplots by modifying the call to `make_subplot`. Some examples are commented out under "other options" towards the end of the python script. 
 
-TBC..
+### Fig. 5: 
+Reproduce Fig. 5 with:
+
+```
+python3 Fig5_zones_at_lsmooth.py 
+```
+
+**Comment:**
 
 
 # Authors and acknowledgment