Recursive algorithms for computing Birkhoff interpolation polynomials
Abstract
As a generalization of Hermite interpolation problem, Birkhoff interpolation is an important subject in numerical approximation. This paper generalizes the existing Generalized Recursive Polynomial Interpolation Algorithm (GRPIA) that is used to compute the Hermite interpolation polynomial. Based on the theory of the Schur complement and the Sylvester identity, the proposed recursive algorithms are applicable to a broader class of Birkhoff interpolation problems, where each interpolation condition is given by the composition of an evaluation functional and a differential polynomial. The approach incorporates a judgment condition to ensure the problem's well-posedness and computes a lower-degree Newton-type interpolation basis (which is also a strongly proper interpolation basis) along with the corresponding interpolation polynomial. Following the numerical examples, we analyze and compare the computational process and complexity of the proposed algorithm against traditional interpolation methods based on Gaussian elimination, and thus demonstrate that the proposed recursive approach reduces both computational cost and storage space requirements.
Cite
@article{arxiv.2511.09014,
title = {Recursive algorithms for computing Birkhoff interpolation polynomials},
author = {Xue Jiang and Yuanhe Li and Zhe Li},
journal= {arXiv preprint arXiv:2511.09014},
year = {2026}
}