def _munp(self, n, c):
        k = arange(0, n+1)
        vals = 1.0/c**n * np.sum(
            comb(n, k) * (-1)**k * special.gamma(c*k + 1),
            axis=0)
        return where(c*n > -1, vals, inf)
dataset/ETHPy150Open scipy/scipy/scipy/stats/_continuous_distns.py/genextreme_gen._munp/original