interval_hdr function

Highest Density Regions