Backward error accurate methods for computing the matrix exponential and its action