fix PairCorrelationFunction_3D()

- where() does not scale at all. Replaced with fancy indexing approach
- input Nx3 matrix rather than individual x,y,z axis arrays
- function returns the used subset of the input matrix
- example script example_3D.py