Problem with expand in lagrange interpolation
- Subject: Problem with expand in lagrange interpolation
- From: Daniel Lakeland
- Date: Wed, 13 Jun 2007 11:01:31 -0700
By expanding a high degree polynomial we typically create a
numerically unstable result, since the first term begins to dominate
rapidly, among other things.
The newton divided difference formula is designed to avoid that, I
think it's the same thing as the horner scheme that we get from
"horner" in maxima
http://en.wikipedia.org/wiki/Lagrange_interpolation
http://en.wikipedia.org/wiki/Newton_polynomial
http://en.wikipedia.org/wiki/Horner_scheme
so perhaps instead of "expand" try applying "horner"
--
Daniel Lakeland
dlakelan at street-artists.org
http://www.street-artists.org/~dlakelan