code: def rho0_c(self, c):
          
          return 200./3*self.rhoc*c**3/(np.log(1.+c)-c/(1.+c))
docs: computes density normalization as a function of concentration parameter
          :return: density normalization in h^2/Mpc^3 (comoving)