equation in a bounded domain. A fully discrete scheme combining a Petrov–Galerkin
spectral method for the spatial discretization and L 1-approximation for the Caputo temporal
derivative is proposed. Stability and convergence of the fully discrete scheme are
rigourously established. Numerical results are presented to confirm the theoretical results.