Refactor and test find_tps_coefficients
- Create a function that computes pairwise rho(r). This function is used as a part of constructing the linear system of equations.