Las ecuaciones diferenciales fraccionarias multitermino (MTFDEs) representan hoy en día una herramienta ampliamente utilizada para modelar muchos procesos importantes, particularmente para sistemas multirrate. Su solución numérica es, por lo tanto, un tema convincente que merece una gran atención, no solo por las dificultades para aplicar métodos de propósito general para ecuaciones diferenciales fraccionarias (FDEs) a este caso. En este documento, primero transformamos las MTFDEs en sistemas equivalentes de FDEs, como lo hicieron Diethelm y Ford; de esta manera, la solución puede expresarse en términos de funciones de Mittag-Leffler (ML) evaluadas en argumentos de matriz. Luego proponemos calcularlo recurriendo al enfoque de matriz propuesto por Garrappa y Popolizio. Se presentan varios tests numéricos que muestran claramente que este enfoque de matriz es muy preciso y rápido, también en comparación con otros métodos numéricos.