In this talk, we develop efficient and accurate algorithms for evaluating both  and  where () is the function related to the exponential,  is an   matrix and  is an  dimensional vector. Such matrix functions play a key role in a class of numerical methods well-known as exponential integrators. The algorithms use the modified scaling and squaring procedure combined with truncated Taylor series. A quasi-backward error analysis is presented to find the optimal value of the scaling and the degree of the Taylor approximation. Some useful techniques are employed for reducing the computational cost. Numerical comparisons with state-of-the-art algorithms show that the algorithms perform well in both accuracy and efficiency.