Radial Distribution Function for a Hard Sphere Gas in the Percus-Yevick Approximation - Is There a Quick Way to Evaluate it Numerically?

The question is as in the title. 


I know that for a one-component hard-sphere gas, there is an analytic solution in the $\vec{k}$-space (PRA, 5, 939). I know that to obtain the correct solution in $r$-space I just need to do something like a Fast Fourier Transform. The problem is that in doing so I can only obtain values for the radial function on a grid, and problem arises as to the magnitude of errors on the grid values, the errors in interpolating between the grid values, and so on.


I have to apologize for the fact that I rarely do fast fourier transforms, and even less estimate their errors.


Is there a standard way (e.g fitting formula) to do this? Thanks.


