HPDcrd function

generates the 100 * (1 - alpha)% most probable interval of a distribution of empirical values