Abstract This paper outlines a reliable strategy to approximate the local stable manifold near a hyperbolic equilibrium point for nonlinear systems of differential equations of fractional order. Furthermore, the local behavior of these systems near a hyperbolic equilibrium point is investigated based on the fractional Hartman-Grobman theorem. The fractional derivative is described in the Caputo sense. The solution existence, uniqueness, stability and convergence of the proposed scheme is discussed. Finally, the validity and applicability of our approach is examined with the use of a solvable model method.