plot_density function

Plot density