code: def rexponweib(alpha, k, loc=0, scale=1, size=None):
      
  
      q = np.random.uniform(size=size)
      r = flib.exponweib_ppf(q, alpha, k)
      return loc + r * scale
docs: Random exponentiated Weibull variates.