Numerical method for fractional Bagley-Torvik equation

In this paper, we solve the fractional Bagley-Torvik equation by using discrete cubic spline and a fourth order approxi-mation based on weighted shifted Grunwald-Letnikov difference operator. By matrix analysis, the numerical scheme is proved to be uniquely solvable and third order accurate. An example is presented to verify the efficiency of the numerical scheme and to compare with other methods in the literature.In this paper, we solve the fractional Bagley-Torvik equation by using discrete cubic spline and a fourth order approxi-mation based on weighted shifted Grunwald-Letnikov difference operator. By matrix analysis, the numerical scheme is proved to be uniquely solvable and third order accurate. An example is presented to verify the efficiency of the numerical scheme and to compare with other methods in the literature.