Abstract
In this article, we develop a powerful method for the numerical solution of boundary value problems (BVPs) of fractional order differential equations (FDEs). Omitting the discretization of data by using Bernstein polynomials, we construct the required scheme. With the help of this scheme we convert the concerned FDEs to algebraic equations whose solutions led us to the numerical solution of the considered problem. Numerical examples are provided to illustrate our main results. Also comparison of the results with the exact solutions and other method like (Haar wavelets) is provided to justify the efficiency of the proposed scheme.