quant_bounds function

Calculate confidence bands for lower and upper quartile