A polynomial approximation for arbitrary functions

We describe an expansion of Legendre polynomials, analogous to the Taylor expansion, for approximating arbitrary functions. We show that the polynomial coefficients in the Legendre expansion, and thus, the whole series, converge to zero much more rapidly compared to those in the Taylor expansion of the same order. Furthermore, using numerical analysis with a sixth-order polynomial expansion, we demonstrate that the Legendre polynomial approximation yields an error at least an order of magnitude smaller than that of the analogous Taylor series approximation. This strongly suggests that Legendre expansions, instead of Taylor expansions, should be used when global accuracy is important.