Although there have existed some numerical algorithms for the fractional differential equations, developing high-order methods (i.e., with convergence order greater than or equal to 2) is just the beginning. Lubich has ever proposed the high-order schemes when he studied the fractional linear multistep methods, where he constructed the pth order schemes (p = 2, 3, 4, 5, 6) for the αth order Rie...