Resolver ecuaciones diferenciales de orden fraccional (es decir, no enteras) de manera precisa, fiable y eficiente es mucho más difícil que en el caso estándar de orden entero; además, la mayoría de las herramientas computacionales no proporcionan funciones integradas para este tipo de problema. En este artículo, revisamos dos de las familias de métodos numéricos más efectivas para problemas de orden fraccional, y discutimos algunos de los principales problemas computacionales como el tratamiento eficiente del término de memoria persistente y la solución de los sistemas no lineales involucrados en los métodos implícitos. Por lo tanto, presentamos un conjunto de rutinas de MATLAB específicamente diseñadas para resolver tres familias de problemas de orden fraccional: ecuaciones diferenciales fraccionarias (EDF) (también para el caso no escalar), sistemas de múltiples órdenes de EDF y EDF de múltiples términos (también para el caso no escalar); se proporcionan algunos ejemplos para ilustrar el uso de las rutinas.