The main purpose of this work is to provide an
optimized version of the Newton divided-difference algorithm
presented by [11] for the polynomial interpolation problem.
This optimized algorithm is used for the computation of the
determinant of a two-variable polynomial matrix.