Spectral method based on new operational matrices for system of differential equations with Jumarie’s modified Riemann-Liouville fractional derivative

In this paper, we present a new operational matrix by Bernstein polynomials and apply it for solving the nonlinear system of differential equations with Jumarie’s modified Riemann-Liouville fractional derivative. This method reduces the problems to a system of algebraic equations. The results obtained are with high accuracy. Also, the solutions approach to classical solutions as the order of the fractional derivatives approach to 1.