In this paper, the non-standard finite difference method (in short NSFD) is implemented to give numerical solutions for linear and nonlinear systems of differential equations of fractional order. The Grunwald-Letnikov ...