def sample(self):
        icdf = self._inverse_cdf_expression()
        return floor(icdf(random.uniform(0, 1)))
A random realization from the distribution
dataset/ETHPy150Open sympy/sympy/sympy/stats/drv.py/SingleDiscreteDistribution.sample