def METHOD_NAME(n, k, m, vol):
    result = [[] for i in xrange(n)]
    beta = (vol/n*k)**(-1.0/m)
    while (min([len(x) for x in result]) == 0):
        i = gen.randint(0, n - 1)
        result[i].append((gen.weibullvariate(beta, m)))
    return result