code: def rhypergeometric(n, m, N, size=None):
      
      if n == 0:
          return np.zeros(size, dtype=int)
      elif n == N:
          out = np.empty(size, dtype=int)
          out.fill(m)
          return out
      return np.random.hypergeometric(n, N - n, m, size)
docs: Returns hypergeometric random variates.