If a function can be explicitly expressed, then one can easily compute its Caputo derivative by the known methods. If a function cannot be explicitly expressed but it satisfies a differential equation, how to seek Caputo derivative of such a function has not yet been investigated. In this paper, we propose a numerical algorithm for computing the Caputo derivative of a function defined by a classical (integer-order) differential equation. By the properties of Caputo derivative derived in this paper, we can change the original typical differential system into an equivalent Caputo-type differential system. Numerical examples are given to support the derived numerical method.
Financed by the National Centre for Research and Development under grant No. SP/I/1/77065/10 by the strategic scientific research and experimental development program:
SYNAT - “Interdisciplinary System for Interactive Scientific and Scientific-Technical Information”.