Las ecuaciones diferenciales de orden fraccional se consideran más desafiantes de calcular en comparación con las ecuaciones diferenciales de orden entero debido a sus propiedades arbitrarias. Este estudio propone un método de múltiples pasos para resolver ecuaciones diferenciales fraccionarias. El método se deriva basado en el concepto de un esquema numérico de Adam-Bashforth de tercer orden implementando interpolación de Lagrange para el caso fraccional, donde las derivadas fraccionarias están definidas en el sentido de Caputo. Además, el estudio incluye una discusión sobre el análisis de estabilidad y convergencia del método. Se proporcionan varios ejemplos numéricos para validar la confiabilidad y eficiencia del método propuesto. Los ejemplos en este estudio cubren la resolución de ecuaciones diferenciales fraccionarias lineales y no lineales para el caso tanto de orden único como de orden superior, donde denota el orden de las derivadas fraccionarias de . La comparación en términos de precisión entre el método propuesto y otros métodos existentes demuestra que el método propuesto ofrece un rendimiento competitivo en comparación con los métodos existentes.