code: def rtruncated_pareto(alpha, m, b, size=None):
      
      u = random_number(size)
      return (-(u * b ** alpha - u * m ** alpha - b ** alpha) /
              (b ** alpha * m ** alpha)) ** (-1. / alpha)
docs: Random bounded Pareto variates.