density.py