Important applications in fractional calculus require the numerical computation of the Mittag-Leffler function with matrix arguments. This topic presents sensitive issues that need to be properly addressed. Furthermore, ad hoc numerical strategies are required when large matrices are involved. Numerical methods will be presented, along with numerical tests to demonstrate their effectiveness.