Solving BVPs using two-point Taylor formula by a symbolic software

For the numerical solution of boundary value problems a global method, based on two-point Taylor formula is proposed. A Mathematica package to compute approximate solution of BVPs is presented.Numerical examples provide favorable comparisons with other existing methods, especially with respect to accuracy.