On the Numerical Computation of the Mittag–Leffler Function