  We prove a discrete analogue for the composition of the fractional integral
and Caputo derivative. This result is relevant in numerical analysis of
fractional PDEs when one discretizes the Caputo derivative with the so-called
L1 scheme. The proof is based on asymptotic evaluation of the discrete sums
with the use of the Euler-Maclaurin summation formula.
