'Return the density parameter for dark energy at redshift ``z``.
Parameters
z : array-like
Input redshifts.
Returns
Ode : ndarray, or float if input scalar
The density of non-relativistic matter relative to the critical
density at each redshift.'

def Ode(self, z):

    if isiterable(z):
        z = np.asarray(z)
        if (self._Ode0 == 0):
            return np.zeros(np.asanyarray(z).shape, dtype=np.float)
    elif (self._Ode0 == 0):
        return 0.0
    return ((self._Ode0 * self.de_density_scale(z)) * (self.inv_efunc(z) ** 2))
