González Rodríguez, PedroMoscoso, MiguelKindelan Segura, Manuel2016-09-142017-10-152015-10-15Journal of Computational Physics, 2015, 299, pp. 307–319.0021-9991https://hdl.handle.net/10016/23560In this paper we describe a numerical algorithm to compute the Laurent expansion of the inverse of singularly perturbed matrices. The algorithm is based on the resolvent formalism used in complex analysis to study the spectrum of matrices. The input of the algorithm are the matrix coefficients of the power series expansion of the perturbed matrix. The matrix coefficients of the Laurent expansion of the inverse are computed using recursive analytical formulae. We show that the computational complexity of the proposed algorithm grows algebraically with the size of the matrix, but exponentially with the order of the singularity. We apply this algorithm to several matrices that arise in applications. We make special emphasis to interpolation problems with radial basis functions.13application/pdfeng© Elsevier 2015Atribución-NoComercial-SinDerivadas 3.0 EspañaLaurent seriesInverseRadial basis functionsInterpolationresearch articleBiología y Biomedicina10.1016/j.jcp.2015.07.006open access307319Journal of computational physics299AR/0000016824