Critical density

临界密度

 from astropy import constants as c, units as u
 rhocrit0=(3*1e4*(u.km/u.Mpc/u.s)**2/8/np.pi/c.G).to(u.Msun / u.Mpc**3).value
 #Msun h^2/Mpc^3