This book discusses numerical methods for solving time-fractional evolution equations. The approach is based on first discretizing in the spatial variables by the Galerkin finite element method, using piecewise linear trial functions, and then applying suitable time stepping schemes, of the type either convolution quadrature or finite difference. The main concern is on stability and error analysis of approximate solutions, efficient implementation and qualitative properties, under various regularity assumptions on the problem data, using tools from semigroup theory and Laplace transform. The...