density_histogram function

Histogram density estimator