[SciPy-user] Numerical stability and Legendre polynomials

Yaroslav Bulatov yaroslavvb at gmail.com
Wed Sep 7 12:28:05 CDT 2005

I'm trying to approximate a function using an expansion over few dozen
of Legendre polynomials, but I'm running into a large numerical error.

For instance the following
evaluates to 2e-4, instead of 0

So my question is -- is such numerical error normal? Is it realistic
to try to approximate a function using 30 Legendre polynomials?

Yaroslav Bulatov
bulatov at cs.oregonstate.edu
Dearborn 102
Oregon State University
Corvallis, OR

More information about the SciPy-user mailing list