This paper sheds the light on the vital role that fractional ordinary differential equations(FrODEs) play in the mathematical modeling and in real life, particularly in the physical conditions. Furthermore, if the problem is handled directly by using numerical method, it is a far more powerful and efficient numerical method in terms of computational time, number of function evaluations, and precision. In this paper, we concentrate on the derivation of the direct numerical methods for solving fifth-order FrODEs in one, two, and three stages. Additionally, it is important to note that the RKM-numerical methods with two- and three-stages for solving fifth-order ODEs are convenient, for solving class's fifth-order FrODEs. Numerical examples have been analyzed to demonstrate the efficacy of the new methods in comparison to the analytical method. Therefore, the numerical compression is carried out to confirm the efficiency and precision of the modified numerical methods. Significantly, the study demonstrates that the numerical outcomes of the proposed derived and modified numerical applied methods proved to be brilliant. Finally, based on the findings of the study, it could be said that the numerical outcomes of the test-problems using proposed and modified methods agree well with the analytical solutions. Hence, we can conclude that the proposed numerical methods that are derived or modified in the analytic study of this paper are quite efficient.
