Return the Lagrange interpolating polynomial of the data-points (x,w)
Warning: This implementation is numerically unstable; do not expect to be able to use more than about 20 points even if they are chosen optimally.
scipy.interpolate.bisplev
scipy.interpolate.approximate_taylor_polynomial