code: def volume_to_distance_with_errors(vol, vol_err):
      
      dist = (vol * 3.0/4.0/numpy.pi) ** (1.0/3.0)
      ehigh = ((vol + vol_err) * 3.0/4.0/numpy.pi) ** (1.0/3.0) - dist
      delta = numpy.where(vol >= vol_err, vol - vol_err, 0)
      elow = dist - (delta * 3.0/4.0/numpy.pi) ** (1.0/3.0)
      return dist, ehigh, elow
docs: Return the distance and standard deviation upper and lower bounds
  
      Parameters
      ----------
      vol: float
      vol_err: float
  
      Returns
      -------
      dist: float
      ehigh: float
      elow: float