The Fractional Complex Step Method

It is well known that the complex step method is a tool that calculates derivatives by imposing a complex step in a strict sense. We extended the method by employing the fractional calculus differential operator in this paper. The fractional calculus can be taken in the sense of the Caputo operator, Riemann-Liouville operator, and so forth. Furthermore, we derived several approximations for computing the fractional order derivatives. Stability of the generalized fractional complex step approximations is demonstrated for an analytic test function.


Introduction
The concept of derivative is one of the most important concepts in science and engineering.It can be described from two equally valid points of view: the geometrical point of view and the physical one.From the geometrical point of view, the derivative can be seen as the tangent line to a function in a certain evaluation point.From the physical point of view, the derivative can be seen as a measure of the rate of change of the function in this point.Further method of computing the derivative of a function comes from its expansion in a Taylor series.Most naturally, derivatives of real functions are evaluated using real numbers, but the less intuitive idea of using an imaginary number in real functions differentiation has been shown capable of overcoming the term cancellation inherent to the ordinary FD method, as well as reducing the associated approximation error.The utilize of complex variables in numerical differentiation was imposed by Lyness and Moler [1], describing a method for calculating the derivatives of any analytic function.Lai and Crassidis [2] used the complex representation of the Taylor series to avoid using the real part for computing the second derivative.Cerviño and Bewley [3] extended the method with an application to pseudospectral simulation codes.Kim et al. [4,5] employed the complex step perturbation in nonlinear robust performance analysis.Recently, the complex set method has been applied by many authors [6][7][8][9][10].
Fractional calculus (real and complex) is a rapidly growing subject of interest for physicists and mathematicians.The reason for this is that problems may be discussed in a much more stringent and elegant way than using traditional methods.Fractional differential equations have emerged as a new branch of applied mathematics which has been used for many mathematical models in science and engineering.In fact, fractional differential equations are considered as an alternative model to nonlinear differential equations.Varieties of them play important roles and tools not only in mathematics but also in physics, dynamical systems, control systems, and engineering to create the mathematical modeling of many physical phenomena.Furthermore, they are employed in social science such as food supplement, climate, and economics.Several different derivatives were introduced: Riemann Liouville, Hadamard, Grunwald Letnikov, Riesz, and Erdelyi-Kober operators and Caputo [11][12][13][14][15][16][17].It is well known that the physical interpretation of the fractional derivative is an open problem today.There is no formal interpretation of the physical meaning of the fractional derivative.Since the appearance of the idea of differentiation and integration of arbitrary order, there was not any acceptable geometric and physical interpretation of these operations for more than 300 year.Recently, the physical interpretation is an open problem.In [18], it is shown that geometric interpretation of fractional integration is "Shadows on the walls" and its Physical interpretation is "Shadows of the past." In this work, we extended the complex step method by employing the fractional calculus differential operator.Furthermore, we derived several approximations for computing the fractional order derivatives.Stability of the generalized fractional complex step approximations is demonstrated for an analytic test function.Moreover, examples are illustrated.

Fractional Calculus
The concept of the fractional calculus (i.e., calculus of integrals and derivatives of any arbitrary real or complex order) was performed over 300 years ago.Abel in 1823 studied the generalized tautochrone problem and for the first time applied fractional calculus techniques in a physical problem.Later Liouville considered fractional calculus to problems in potential theory.Since that time, the fractional calculus has haggard the attention of many researchers in all areas of sciences.

The Fractional Complex Step Method
The local fractional Taylor formula has been generalized by many authors [19][20][21][22].This expansion takes the following formula: where    is the Riemann-Liouville differential operator and .
As a generalization of the above approximate method, we let ( +  +  √ V), the corresponding fractional Taylor series expansion, become where  and V are real numbers related to the real and imaginary differential steps.It is clear that when  = 0, (9)  17)-( 20):  1/5  2 = 1.8 *  4/5 * .

Conclusion
We extended the complex step method by employing the fractional calculus differential operator (the fractional complex step).The approximation is provided for  ∈ (0, 1), and in the same manner, we can consider the approximation for all  ∈ ( − 1, ).This derivative concept imposes improving many different approximations for the fractional derivatives of any complex valued analytic function utilizing its real and imaginary parts.We provided different approximations for the operator   .Moreover,  2 is approximated in (8).This work can be applied in physics and computer sciences such as image processing and signal processing.