slightly changed error bars, probably incorrect

plots/densityprofile.py

@@ -37,7 +37,7 @@ rho /= 4*np.pi*r**2*dr
 print("Number of bins:",len(rho))
 
 # plot poissonian error bars
-plt.errorbar(r, rhoHernquist, yerr=np.sqrt(rhoHernquist), fmt='grey', label='Hernquist with Poisson error')
+plt.errorbar(r, rhoHernquist, yerr=np.sqrt(rhoHernquist), color='grey', fmt='.', label='Hernquist with Poisson error', markersize=0)
 # plot Hernquist profile
 plt.loglog(r, rhoHernquist, 'r.-', label='Hernquist Density Profile')
 # plot density profile of data