Added routines to compute pair correlation function

Handles 2D and 3D distributions of particles
Includes a demo "driver" script