compute_density function

Compute density of data.