Numerical solution of fractional oscillator equation

We focus on a numerical scheme applied for a fractional oscillator equation in a finite time interval. This type of equation includes a complex form of left- and right-sided fractional derivatives. Its analytical solution is represented by a series of left and right fractional integrals and therefore is difficult in practical calculations. Here we elaborated two numerical schemes being dependent on a fractional order of the equation. The results of numerical calculations are compared with analytical solutions. Then we illustrate convergence and stability of our schemes.