def sample(self, n=1):
        return self.dist().rvs(size=n)[:, np.newaxis]