“…Also, they can be differentiated and integrated easily and can approximate a great variety of functions to any accuracy desired. For these reasons, for the numerical solution of problem (1)-(2), many methods have been developed which uses polynomial basis, such as spline functions ( [1,13,18,28,29]), Bernstein polynomials ( [3]), Hermite polynomials ( [30]), classical Taylor polynomials ( [19,31]), two-point Taylor polynomials ( [9,26]), Lagrange polynomials ( [8]), Lidstone polynomials ( [10]). In [33] the author uses Chebyshev polynomials, along with fast Fourier transform and baricentric interpolation, in order to construct a powerful tool, Chebfun ( [12]), for the solution of differential equations.…”