Interpolation based on Radial Basis Functions (RBF) is very often used for scattered scalar data interpolation in n-dimensional space in general. RBFs are used for surface reconstruction of 3D objects, reconstruction of corrupted images etc. As there is no explicit order in data sets, computations are quite time consuming that leads to limitation of usability even for static data sets. Generally the complexity of computation of RBF interpolation for N points is of O(N3 ) or O(k N2 ), k is a number of iterations if iterative methods are used, which is prohibitive for real applications. The inverse matrix can also be computed by the Strassen algorithm based on matrix block notation with O(N2.807 ) complexity. Even worst situation occurs when interpolation has to be made over non-constant data sets, as the whole set of equations for determining RBFs has to be recomputed. This situation is typical for applications in which some points are becoming invalid and new points are acquired. In t...