[SciPy-user] Numerical stability and Legendre polynomials