A novel algorithm is proposed for computing matrix-vector products A^\alpha v, where A is a symmetric positive semidefinite sparse matrix and \alpha > 0. The method can be applied for the efficient implementation of the matrix transformation method to solve space-fractional diffusion problems. The performance of the new algorithm is studied in a comparison with the conventional MATLAB subroutines to compute matrix powers
