This paper proposes an cost-efficient method for the numerical solution of fractional differ- ential equations of multi-order. The suggested algorithm adopts the operational matrix of fractional derivative in term of Caputo sense using modified cubic B-spline functions. The new method leads to a system of algebraic equations that improves the solution accuracy and reduces the computational complexity. Several numerical examples illustrate the applicability of the proposed method.
