discint function

Highest probability interval for a discrete distribution