Need help with understanding Thin Plate Splines as Radial Basis function for Meshfree simulation

I want to try applying TPS as the RBF for the simulation of groundwater flow and transport. but the i'm getting error because of zeros in diagonal elements in the global matrix. Please help with the shape function