docs: def dist_points(bin_edges, d):
      """Return an array of values according to a distribution
  
      Points are calculated at the center of each bin
      """
      bc = bin_centers(bin_edges)
      if d is not None:
          d = DISTS[d['type']](d, bc)
      return d, bc