in this article we implement an operational matrix of fractional integration for legendre polynomials. we proposed an algorithm to obtain an approximation solution for fractional differential equations, described in riemann-liouville sense, based on shifted legendre polynomials. this method was applied to solve linear multi-order fractional differential equation with initial conditions, and the...