Why use the different polynomial for newton forward and lagrange
Answers
Hiii
Frankly, Lagrange interpolation is mostly just useful for theory. Actually computing with it requires huge numbers and catastrophic cancellations. In floating point arithmetic this is very bad. It does have some small advantages: for instance, the Lagrange approach amounts to diagonalizing the problem of finding the coefficients, so it takes only linear time to find the coefficients. This is good if you need to use the same set of points repeatedly. But all of these advantages do not make up for the problems associated with trying to actually evaluate a Lagrange interpolating polynomial.
With Newton interpolation, you get the coefficients reasonably fast (quadratic time), the evaluation is much more stable (roughly because there is usually a single dominant term for a given x), the evaluation can be done quickly and straightforwardly using Horner's method, and adding an additional node just amounts to adding a single additional term. It is also fairly easy to see how to interpolate derivatives using the Newton framework.
Here are two differences:
Lagrange's form is more efficient when you have to interpolate several data sets on the same data points.
Newton's form is more efficient when you have to interpolate data incrementally.
Lagrange method is mostly a theoretical tool used for proving theorems. Not only it is not very efficient when a new point is added (which requires computing the polynomial again, from scratch), it is also numerically unstable.
Therefore, Newton's method is usually used. However, there is a variation of the Lagrange interpolation, which is numerically stable and computationally efficient! Unfortunately, this method is not very known...
I am attaching a link to a paper from SIAM Review called "Barycentric Lagrange interpolation", which is not difficult to read. I hope you will find it interesting.