Abstract. This paper takes a look at numerical procedures for computing approximation of the exponential of a matrix of large dimension. Existing approximation methods to evaluate the exponentiation of a matrix will be reviewed paying more attention to Krylov subspace methods and Schur factorization techniques. Some theoretical results on the bounds for the entries of the exponential matrix and some implementation details will be also discussed.