This paper proposes a novel numerical approach for handling fractional boundary value problems. Such an approach is established on the basis of two numerical formulas; the fractional central formula for approximating the Caputo differentiator of order α and the fractional central formula for approximating the Caputo differentiator of order 2α, where 0<α≤1. The first formula is recalled here, whereas the second one is derived based on the generalized Taylor theorem. The stability of the proposed approach is investigated in view of some formulated results. In addition, several numerical examples are included to illustrate the efficiency and applicability of our approach.