School of Mathematics and Statistics, Carleton University, Ottawa, ON, Canada.

Abstract: A new algorithm for computing integrals involving the matrix exponential is given. The method employs diagonal Padé approximation with scaling and squaring ...