diff --git a/plots/densityprofile.png b/plots/densityprofile.png
index aee049b8c2267b68684a020a504d73641ff59fe6..0dc0a8f3a8552a1da01f0b449b6b2de498d97813 100644
Binary files a/plots/densityprofile.png and b/plots/densityprofile.png differ
diff --git a/plots/densityprofile.py b/plots/densityprofile.py
index 1bf376648c37bd45a01d28f2da37a234e1cb83b6..a9870dffe0d2c8b9cacb298dc08d08cc21e086b1 100644
--- a/plots/densityprofile.py
+++ b/plots/densityprofile.py
@@ -45,10 +45,10 @@ print("Bin size:",dr)
 
 # plot poissonian error bars
 plt.errorbar(bins[1:], rhoHernquist, yerr=np.sqrt(rhoHernquist), color='grey', fmt='.', label='Hernquist with Poisson error', markersize=0)
-# plot Hernquist profile
-plt.loglog(bins[1:]-dr/2, rhoHernquist, 'r-', label='Hernquist Density Profile', linewidth=1)
 # plot density profile of data 
 plt.loglog(bins[1:]-dr/2, hist, 'b-', label='Density Profile from data', linewidth=1)
+# plot Hernquist profile
+plt.loglog(bins[1:]-dr/2, rhoHernquist, 'r-', label='Hernquist Density Profile', linewidth=1)
 
 plt.xlabel(r'Radius $r$')
 plt.ylabel(r'Density $\rho(r)$')