Radial distribution function interpolation from Numerical Recipies