An alternative way to find a polynomial approximation of a function is something of a mix of using the Taylor Series and numerical differentiation. Newton’s Divided Difference method results in a nice, recursive mechanism for calculating coefficients of the polynomial.

Click here for lecture slides based on Kaw, Chapter 5.03.