import numpy as np
import holoviews as hv
from holoviews import opts
hv.extension('bokeh')

xs = np.linspace(-5, 5, 100)
ys = -(xs-2)**2
ymax = ys.argmax()

curve = hv.Curve((xs,ys))
vspan = hv.VSpan(xs[ymax-5], xs[ymax+5])

curve.opts(color='#D3D3D3') * vspan.opts(color='red')

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