code: def _areadist(ax, v, xr, c, bins=100, by=None, alpha=1, label=None):
      
      y, x = np.histogram(v[~np.isnan(v)], bins)
      x = x[:-1]
  
      if by is None:
          by = np.zeros((bins,))
  
      ax.fill_between(x, y, by, facecolor=c, alpha=alpha, label=label)
      return y
docs: Plot the histogram distribution but as an area plot