code: def make_quantile_df(data, draw_quantiles):
      
      dens = data['density'].cumsum() / data['density'].sum()
      ecdf = interp1d(dens, data['y'], assume_sorted=True)
      ys = ecdf(draw_quantiles)
  
      # Get the violin bounds for the requested quantiles
      violin_xminvs = interp1d(data['y'], data['xminv'])(ys)
      violin_xmaxvs = interp1d(data['y'], data['xmaxv'])(ys)
  
      data = pd.DataFrame({
          'x': interleave(violin_xminvs, violin_xmaxvs),
          'y': np.repeat(ys, 2),
          'group': np.repeat(np.arange(1, len(ys)+1), 2)})
  
      return data
docs: Return a dataframe with info needed to draw quantile segments