RMPIA: a new algorithm for computing the Lagrange matrix interpolation polynomials