Hard Sphere Monte Carlo Simulations for Pair Correlation Functions