Computing the Matrix Mittag-Leffler Function with Applications to Fractional Calculus