Se presenta en este artículo un esquema basado en operadores para la integración numérica de ecuaciones diferenciales fraccionarias. El operador diferencial generalizado se utiliza para construir la solución analítica a la ecuación diferencial ordinaria característica correspondiente en forma de una serie de potencias infinita. La solución numérica aproximada se construye truncando la serie de potencias y cambiando el punto de la expansión. La estrategia de selección de paso de integración adaptativa desarrollada se basa en el error controlado de la aproximación inducida por la truncación. Se utilizan experimentos computacionales para demostrar la eficacia del esquema propuesto.