Neville's interpolation method for tabulated functions is an optimized algorithm for constructing
the unique Lagrange interpolation polynomial for a given set of data points. It provides an error estimate as well.
Coded in C++ and Python.