The two-parameter Mittag–Leffler function Eα,β is of fundamental importance in fractional calculus, and it appears frequently in the solutions of fractional differential and integral equations. However, the expense of calculating this function often prompts efforts to devise accurate approximations that are more cost-effective. When α>1, the monotonicity property is largely lost, resulting in the emergence of roots and oscillations. As a result, current rational approximants constructed mainly for α∈(0,1) often fail to capture this oscillatory behavior. In this paper, we develop computationally efficient rational approximants for Eα,β(−t), t≥0, with α∈(1,2). This process involves decomposing the Mittag–Leffler function with real roots into a weighted root-free Mittag–Leffler function and a polynomial. This provides approximants valid over extended intervals. These approximants are then extended to the matrix Mittag–Leffler function, and different implementation strategies are discussed, including using partial fraction decomposition. Numerical experiments are conducted to illustrate the performance of the proposed approximants.