High-accuracy numerical integration methods for fractional order derivatives and integrals computations

In this paper the authors present highly accurate and remarkably efficient computational methods for fractional order derivatives and integrals applying Riemann-Liouville and Caputo formulae: the Gauss-Jacobi Quadrature with adopted weight function, the Double Exponential Formula, applying two arbitrary precision and exact rounding mathematical libraries (GNU GMP and GNU MPFR). Example fractional order derivatives and integrals of some elementary functions are calculated. Resulting accuracy is compared with accuracy achieved by applying widely known methods of numerical integration. Finally, presented methods are applied to solve Abel’s Integral equation (in Appendix).