Este documento investiga la solución a ecuaciones diferenciales fraccionarias unidimensionales con dos tipos de operadores de derivadas fraccionarias de órdenes en el rango de . Se construyen dos esquemas linealizados del método numérico. Las EDF consideradas son transformadas de manera equivalente por la integral de Riemann-Liouville en sus problemas diferenciales parciales integro para reducir el requisito de suavidad en el tiempo. Se discute rigurosamente el análisis de estabilidad y convergencia. Finalmente, se describen experimentos numéricos que confirman el análisis teórico obtenido.