added error bars

plots/densityprofile2.py

@@ -29,6 +29,8 @@ rhoHernquist = M / (2*np.pi) * a / (bins[1:]-dr/2) / (bins[1:]-dr/2 + a)**3
 print("Number of bins:",nbins)
 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')
 # plot density profile of data