- /* 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);