code: def gen_mu(K, delta, c):
      
  
      S = c * log(K/delta) * sqrt(K) 
      tau = gen_tau(S, K, delta)
      rho = gen_rho(K)
      normalizer = sum(rho) + sum(tau)
      return [(rho[d] + tau[d])/normalizer for d in range(K)]
docs: The Robust Soliton Distribution on the degree of 
      transmitted blocks