Interpolation Polynomials and Linear Algebra
We reconsider the theory of Lagrange interpolation polynomials with multiple interpolation points and apply it to linear algebra. For instance, be a linear operator satisfying a degree polynomial equation . One can see that the evaluation of a meromorphic function at is equal to , where is the degree interpolation polynomial of with the the set of interpolation points equal to the set of roots of the polynomial . In particular, for an matrix, there is a common belief that for computing one has to reduce to its Jordan form. Let be the characteristic polynomial of . Then by the Cayley-Hamilton theorem, . And thus the matrix can be found without reducing to its Jordan form. Computation of the Jordan form for involves many extra computations. In the paper we show that it is not needed. One application is to compute the matrix exponential for a matrix with repeated eigenvalues, thereby solving arbitrary order linear differential equations with constant coefficients.