This study outlines a modified implicit finite difference method for approximating the local stable manifold near a hyperbolic equilibrium point for a nonlinear systems of fractional differential equations. The fractional deriva- tive is described in the Caputo sense of order (0 < 1) which is approxi- mated based on the modified trapezoidal quadrature rule of order O ( △ t 2
