code: def simulate(
      N=100,
      k=10000,
  ):
      truth = gen_prob_dist(N)
  
      area_ratio = truth
      accept, alias = create_alias_table(area_ratio)
  
      ans = np.zeros(N)
      for _ in range(k):
          i = alias_sample(accept, alias)
          ans[i] += 1
      return ans / np.sum(ans), truth
  
imports: import matplotlib.pyplot
  import numpy
  from ge.alias import alias_sample, create_alias_table
label_desc: Use image node_id 2 for calling a global function with example usage: simulate(N, k) and returns: unknown,  truth