“…An automatic quadrature method based on the Chebyshev polynomials was presented for approximating the Caputo derivative in [26]. Some other computational schemes, such as the L1, L2 and L2C schemes, etc., are also introduced [8,11,13,14,16,18,19,21,23,25,28,29,31].…”