'Critical density in grams per cubic cm at redshift ``z``.
Parameters
z : array-like
Input redshifts.
Returns
rho : `~astropy.units.Quantity`
Critical density in g/cm^3 at each input redshift.'

def critical_density(self, z):

    return (self._critical_density0 * (self.efunc(z) ** 2))
