hdr function

Compute highest density regions