In this paper we develop an iterative numerical method based on Bernstein splines for solving two-point boundary problems associated to differential equations of fractional order . The convergence of the method is proved by providing the error estimate and it is tested on a numerical example.