code: def K(self):
          
          K1 = np.float_power(pi, 0.5 * self.n_min * self.n_min)
          K1 /= (
              np.float_power(2, 0.5 * self.n_min * self._n_max)
              * self._mgamma(0.5 * self._n_max, self.n_min)
              * self._mgamma(0.5 * self.n_min, self.n_min)
          )
  
          K2 = np.float_power(
              2, self.alpha * self.size + 0.5 * self.size * (self.size + 1)
          )
          for i in xrange(self.size):
              K2 *= gamma(self.alpha + i + 1)
          return K1 * K2
docs: Normalizing constant for wishart CDF.