Code to compute the correlation function from the P(k)