A powerful tool for dealing with high-dimensional fractional-order systems with applications to fractional Emden–Fowler systems

In this study, an approximation solution for a high-dimensional system in terms of the Caputo fractional derivative operator is obtained using the improved modified fractional Euler method, or IMFEM for short. To accomplish this aim, a result that can transform such a system into a double-equation, one-dimensional fractional-order system, is provided theoretically. Some physical applications, including fractional-order systems of equations of Emden–Fowler type, are discussed, and their graphs are plotted using MATLAB to demonstrate the IMFEM schema’s validity.