“…The proposed method is fourth order convergence in space and (2 − )-th order convergence in time. Various numerical schemes based on B-spline collocation methods have been applied to solve a wide variety of problems, see [24][25][26][27][28][29][30][31][32][33][34]. Besides B-spline collocation techniques, there are other methods that can be used to solve classic or fractional differential equations, such as finite difference methods [6,7,14], finite element method [35], boundary element methods [36,37], meshless methods [38,39] etc.…”