First draft of dgp_density.py