This paper is devoted to the numerical simulation for solving a special class of fractional differential equations. Based on the Grünwald–Letnikov definition of a fractional derivative, a numerical scheme for the approximation of the solution is discussed. By using this scheme, we solve the fractional Van der Pol equation. The results obtained here compare well with the analytical solutions and this shows that the numerical scheme is stable.
