A new algorithm based on compressed Legendre polynomials for solving boundary value problems

In this paper, we discuss a novel numerical algorithm for solving boundary value problems. We introduce an orthonormal basis generated from compressed Legendre polynomials. This basis can avoid Runge phenomenon caused by high-order polynomial approximation. Based on the new basis, a numerical algorithm of two-point boundary value problems is established. The convergence and stability of the method are proved. The whole analysis is also applicable to higher order equations or equations with more complex boundary conditions. Four numerical examples are tested to illustrate the accuracy and efficiency of the algorithm. The results show that our algorithm have higher accuracy for solving linear and nonlinear problems.