Radial basis functions are a very accurate means of solving interpolation and partial differential equations problems. The global radial basis functions collocation technique produces ill-conditioning matrices when using multiquadrics, making the choice of the shape parameter a crucial issue. The use of local numerical schemes, such as finite differences produces much better conditioned matrices. However, finite difference schemes are limited to special grids. For scattered points, a combi...