Statistics: added subroutine prob_range() to determine probability ranges

* Statistics: added subroutine prob_range() to determine probability ranges:
  * returns the lower and upper limit of the narrowest ranges containing a given
    fraction of the provided data points
  * prob_range_sp() is a single-precision wrapper