A robust scheme for reduction of higher fractional-order systems

The objective of this work is to present a numerical solution to a system of higher fractional-order differential equations with initial value problems. In order to achieve this objective, we develop a novel theoretical result aimed to reduce these higher fractional-order systems to α -fractional systems, where 0<α≤ 1 , and then apply a recent numerical approach called modified fractional Euler method, which is regarded a numerical modification of the fractional Euler Method (FEM). Finally, we will give numerical applications to illustrate our results using MATLAB procedures.