The subject of this work is the study of Birkhoff’s polynomial interpolation. Given a list \( Z_n=[z_{0},...,z_{n}] \) with \( (n+1) \) distinct nodes, of \( \mathbb {K}=\mathbb {R} \) or \( \mathbb {C} \) , we will seek to study the questions of existence and uniqueness of a polynomial P such that P and a number of its derivatives take, in these nodes, given values. Recently, Messaoudi et al. (Numer. Algorithms 80, 253–278 2019) presented a new algorithm for computing the Hermite interpolation polynomial called the Generalized Recursive Polynomial Interpolation Algorithm (GRPIA). In this paper, we will give a new formulation of the Birkhoff polynomial interpolation problem and derive a new algorithm, called the Recursive Hermite-Birkhoff Polynomial Interpolation Algorithm (RHBPIA), to solve the Birkhoff interpolation problem and generalize the GRPIA. A new existing result will be established. The numerical stability of this algorithm will also be studied and some examples will be given.