import numpy as np
import holoviews as hv
from holoviews import opts

hv.extension('matplotlib')

xs = np.random.normal(size=500)
ys = np.random.normal(size=500) * xs
ymean, ystd = ys.mean(), ys.std()

points = hv.Points((xs,ys))
hspan  = hv.HSpan(ymean-ystd, ymean+ystd)

hspan.opts(facecolor='blue') * points.opts(color='#D3D3D3')

(hv.HSpan(1, 3) * hv.HSpan(5, 8)).opts(
    opts.HSpan(apply_ranges=True))