Create histogram_number_density.py