Abstract
This paper proposes a new method of numerical calculation for fractional derivatives. Conventional methods using difference approximation have some defects such as rounding errors and computation time. In this paper, fractional derivatives are expressed as the combination of ordinary derivatives and integrals which are calculated using automatic differentiation and double exponential formula, respectively. Numerical examples show that the proposed method gives high-precision results with practical computational cost.
