add support_density.py