Evolution equations containing fractional derivatives can provide suitable mathemati- cal models for describing important physical phenomena. In this paper, we propose a fast and accurate method for numerical solutions ...