On the numerical computation of the matrix exponential

Let us consider the initial-value problem of dimension m: $$ \frac{d\tau}{d}y(\tau) = f(\tau, Y(\tau)), y(0) = y_0, \tau \geq 0, (1.1) $$ Where $ = (f_1, f_2, \cdots, f_m) and y = (y_1, y_2, \cdots, y_m)$.