code: def redshift_from_comoving_volume(vc, **kwargs):
      r
      cosmology = get_cosmology(**kwargs)
      return z_at_value(cosmology.comoving_volume, vc, units.Mpc**3)
docs: r"""Returns the redshift from the given comoving volume.
  
      Parameters
      ----------
      vc : float
          The comoving volume, in units of cubed Mpc.
      \**kwargs :
          All other keyword args are passed to :py:func:`get_cosmology` to
          select a cosmology. If none provided, will use
          :py:attr:`DEFAULT_COSMOLOGY`.
  
      Returns
      -------
      float :
          The redshift at the given comoving volume.