Created pair distribution function from parts of function from Anders and Jørgen