+ /* Include cases which are just on the boundary */
+ if (y == gsl_histogram_max (gslh))
+ {
+ double lower, upper;
+ gsl_histogram_get_range (gslh, gsl_histogram_bins (gslh)-1, &lower, &upper);
+ gsl_histogram_accumulate (gslh, lower + (upper - lower)/2.0, c);
+ }
+ else
+ gsl_histogram_accumulate (hist->gsl_hist, y, c);