Computing the matrix Mittag-Leffler function with applications to fractional calculus