def E_sq(z, om0):
    """
    A function giving Hubble's law for flat cosmology

    :param z: redshift value
    :param om0: matter density
    :return: A value for the Hubble parameter
    """
    return om0 * (1 + z) ** 3 + 1 - om0