'Compute the n-th non-central moment.'

def _munp(self, n):

    integrals = (((self._hbins[1:] ** (n + 1)) - (self._hbins[:(-1)] ** (n + 1))) / (n + 1))
    return np.sum((self._hpdf[1:(-1)] * integrals))
