Stability and convergence of the L1 formula on nonuniform time grids are studied for solving linear react ion-sub diffusion equations with the Caputo derivative. A discrete fractional Gronwall ...