An algorithm is described which computes the exponential function of a matrix multiplying a vector by applying the Arnoldi algorithm to the matrix and computing the exponential function of the resulting small Hessenberg matrix. Error estimates based on the theory for Lagrange interpolation are given and a numerical test is reported.