code: def median_bias(n):
      
      if type(n) is not int or n <= 0:
          raise ValueError('n must be a positive integer')
      if n >= 1000:
          return numpy.log(2)
      ans = 1
      for i in range(1, int((n - 1) / 2 + 1)):
          ans += 1.0 / (2*i + 1) - 1.0 / (2*i)
      return ans
docs: Calculate the bias of the median average PSD computed from `n` segments.
  
      Parameters
      ----------
      n : int
          Number of segments used in PSD estimation.
  
      Returns
      -------
      ans : float
          Calculated bias.
  
      Raises
      ------
      ValueError
          For non-integer or non-positive `n`.
  
      Notes
      -----
      See arXiv:gr-qc/0509116 appendix B for details.