Fractional Definite Integral

This paper proposes the definition of fractional definite integral and analyses the corresponding fundamental theorem of fractional calculus. In this context, we studied the relevant properties of the fractional derivatives that lead to such a definition. Finally, integrals on R2 and R3 are also proposed.


Introduction
Fractional Calculus has evolved considerably during the last 30 years and has become popular in many scientific and technical areas [1][2][3][4].The progress in applications vis-a-vis theoretical developments motivated the re-evaluation of past formulations.The concepts of fractional derivative (FD) and fractional integral (FI) assume various forms not always equivalent and also not compatible with each other.Notwithstanding this development a singular situation exists, since there is no definition of "fractional definite integral" [5,6].Even more strangely, it seems that no-one has the goal to propose such a adefinition [7,8].In fact, the term FI is used for an indefinite integral, which, as we will see, stands for a primitivation operation.
A review of the literature reveals two indirect attempts to formulate the "fractional generalization of the fundamental theorem of calculus" [9,10].However, in [9] several formulations of the theorem were described, one for each derivative; the approach followed in [10] was more directed into application but also did not address the definite integral definition.None of those formulations introduced the notion of a definite integral that will be the topic of this paper.Starting from a revision of classic results, an approach based on a generalisation of the Barrow formula is used to propose the "definite fractional integral" and, from it, to formulate the "fractional fundamental theorem of calculus".These developments allowed the definition of double and triple integrals on rectangular spaces.
The paper is organised as follows.In Section 2, we introduce the "definite fractional integral" and the corresponding "fractional fundamental theorem of calculus".For this purpose we study the admissibility of fractional derivatives and we find those that are suitable for achieving the proposed integral.In Section 3 the definite integrals on R, R 2 , and R 3 are presented.Finally, in Section 4 the conclusions are drawn.

On the One-Sided Integer Order Derivatives and Their Inverses
Let f be a function defined on a closed interval [c, d] ∈ R where f is continuous.Then, f is said to have a lefthand derivative at any point of the closed interval if exists as a finite value or ±∞.The righthand derivative is defined similarly by the expression If both derivatives exist at x = x 0 and are equal, then we say that f is differentiable at that point.In the following we work preferably with the left derivative and will omit the subscript "l" unless it is needed to clarify the formulae.Proceeding as in [11] it is possible to obtain an operator, J, that is the inverse of the derivative: where (x − c)/h is the integer part of (x − c)/h.The summation in (3) is the so-called Riemann sum.
Using (1) in ( 3) and vice-versa it is straightforward to show that This means that J is simultaneously the left and right inverse of D.
As the derivative of a constant is zero, there are many right inverses, P, of the derivative having the form P f (x) = D −1 f (x) + C, C ∈ R. In fact, DP f (x) = f (x), but PD f (x) = C.These inverses are called primitives of f (x).The particular primitive P = J = D −1 can be called "proper primitive" [12].Here, it will be named anti-derivative and when needed a subscript will be inserted to indicate that it is a left operator D l .
On the other hand, if we represent J by D −1 , the Formulaes (2) and (3) can be joined in one: where (a a) 0 = 1, is the Pochhammer symbol for the rising factorial.In a similar way, the right hand derivative and inverse can be formulated as This unified formulation opens a new perspective into the generalisation for fractional derivatives.

Order 1 Integral
There are several equivalent definitions of integral of a function in [a, b], a ≥ c, b ≤ d (in particular c can be −∞ and d can be +∞), but there is no doubt about the meaning of b a f (x)dx and its computation.Moreover, b a f (x)dx enjoys relevant properties such as follows.
Let f (x) as above and F(x) be a function defined, for all x ∈ [a, b], by differentiable on the open interval (a, b), and for all x ∈ (a, b).This constitutes the fundamental theorem of integral calculus and has as important consequence the well-known Barrow formula These results show that the function F(x) is the anti-derivative and can be expressed by ( 4) If the upper limit is variable relation (10) gives establishing the connection and simultaneously the difference between the anti-derivative and the integral.
These expressions were presented in terms of the left derivative and anti-derivative.The corresponding right operators are readily obtained.However, in an integer order formulation, they give the same result in all the points where f (x) is a continuous function.

Definite Fractional Integral
According to the above discussion, it is not clear how to obtain a generalisation of the definite integral.However, expression (10) suggests a shortcut for the solution: introduce it using a generalization of the Barrow formula.
Let f (x), x ∈ R, such that exist their left and right derivatives of any order.Denote the left and right anti-derivatives by f Definition 1.We define α-order fractional integral (FI) of f (x) over the interval (a, b) through the fractional Barrow formula In coherence with the integer order result, this FI must lead to a fractional formulation of the fundamental theorem of integral calculus.Making variable the upper limit, b = x ∈ R, and putting f (x) = D α lr g(x) in ( 13) it comes but it is expected to obtain This means that we should have On the other hand, according to (8), the right hand side in the equation, should be equal to f (x).This implies that and that the derivative of a constant is zero These results show that not all the FD can be used to generalise the notion of FI.For now, let us assume that the derivatives to be adopted verify (15), (17), and (18).
In this case, ( 12) can be readily generalised to This definite integral makes sense if suitable FD are used.Therefore, the admissible derivatives have to agree with relations (15) to (19).In the next section, the most important derivatives are analysed to establish their admissibility.

Which Fractional Derivative?
In a previous paper we answered to the question "What is a fractional derivative?" [13] using two criteria as guidelines for deciding if a given operator is a FD.The wide sense criterion is as follows.
An operator is considered as a FD under 1 P criterion if it enjoys the properties: The index law holds for negative orders 1 P5 Generalised Leibniz rule The index law property can be modified to include positive orders.This led to the strict sense criterion.Therefore, criterion 2 P keeps four conditions and 1 P4 is modified to: 2 P4 The index law holds for any real α and β.The difference between the two criteria lies in the validity of the index law for positive values only, or for any real order.
Let us consider relation (20) again and the case α > 0 with β = −α.We conclude that for each derivative operator, D α , there is an anti-derivative, D −α .As in the classic framework we can say that, given a function, it can have an infinite number of primitives, but only one of them does not depend on arbitrary constants.

The Riemann-Liouville and Caputo Derivatives
For the RL and C derivatives the negative order operator is given by the Riemann-Liouville fractional integral (RL-FI) that can be introduced as follows.Let [a, b] (−∞ < a < b < ∞) be an interval on the real axis R. The left-sided RL-FI, I α a+ f , of order α ∈ R + is given by [6] RL The left Riemann-Liouville fractional derivative (RL-FD), D α a+ f , of order α ∈ R + 0 is defined as The C-FD, C D α a+ f (x), of order α ∈ R + on [a, b] is introduced via the above RL-FD as [5]: With the RL-FI formula above introduced, we verify [5,13] that:

•
The RL-FD is an inverse operator of the left RL-FI • The C-FD is also an inverse operator of the left RL-FI Therefore, the RL-FI is the right inverse of the RL and C derivatives.The RL-FI is not a left inverse, because we have [5]: In order to study the admissiblity of RL and C derivatives we compare their properties with those required in the previous sub-section, namely conditions (15) to (18).
The RL derivative is not acceptable because it does not give the result zero for the derivative of a constant.On the other hand, (26) is not suitable also since it introduces terms that do not verify (19).Concerning the C derivative, it is acceptable from the point of view of (25) and gives zero for the derivative of a constant, but ( 27) is only acceptable if α ≤ 1.Also, it creates difficulties when trying to use (19), since we must have that is true if, and only if the following equation is valid However, in general, this relation is not valid.

Grünwald-Letnikov Derivatives
In a recent paper [14], we considered the problem of selecting derivatives suitable for generalizing standard laws of physics and results while maintaining fully compatibility with them.It was shown that the Grünwald-Letnikov (GL) derivative verifies the above strict sense criterion and is fully compatible with the classic laws and tools of Physics.Here it is important to study the GL derivative from the point of view of the fractional integral introduced in Section 2.3.
Let f (x) be a function defined on R. We introduce the left (forward) FD by means of the fractional incremental ratio where h > 0. Relation (28) expresses the called left GL-FD [6,15,16].Without intending to explore existence problems (see [16]) we can say that f (x) must decrease to zero as x goes to −∞.It is interesting to remark that (28) generalises (5) and this corresponds to making α = ±1.
The GL derivative enjoys relevant characteristics, namely the index law (20) [15,16].This means that given an FD of order α > 0, there is an FD, of negative order, that we will call "anti-derivative" and verifying It can be shown [14] that the FD of the constant function is identically null.The corresponding righthand (backward) GL derivative is defined by Contrarily to the above situation, here a suitable function, f (x), must decrease to zero when x → +∞.
When dealing with space derivative we do not need to impose causality, because we can move in all the directions in the 3D space.Therefore, the factor e −iπα can be removed.This allows us to write a unified formula joining the lefthand and righthand derivatives Attending to (29) and to the fact that the GL derivative of a constant is zero, we arrive to the conclusion that the GL derivative is suitable for computing the definite integral introduced in Section 2.3.

Liouville Derivatives
The name "Liouville derivative" is frequently attached to the RL derivative when the integration domain is R [6].This derivative has the inconvenience of being divergent when the function is constant.
Another alternative is given by the so called Liouville-Caputo derivative [1,17,18] that is also similar to the C derivative but is also defined on R.This derivative poses strong requirements to the integrand function, since this must have the nth order derivative with n ≥ α.A third alternative was proposed by Liouville in the 19th century and can be obtained from the system interpretation of the GL derivative.Let h(x), x ∈ R be the impulse response of the GL derivative [15] where u(x) is the Heaviside unit step.In a linear system perspective, when the input is a given function, f (x), the output is a α-order derivative that generalises the well known result f (x) = ∞ −∞ f (x − y)δ (y)dy,and represents an integral formulation of the left side FD.For α > 0, the above integral is singular.However, it can be regularised, leading to what we will call the left Liouville (L) derivative, defined as follows: where N is the integer part of α, so that α − This derivative verifies also relations (29).The equivalence of GL and L derivatives is valid for at least functions with Laplace transform, although we can assume it for a broader class of functions and use them interchangeably.Joining the left and right Liouville derivatives, it becomes

Definite Fractional Integrals
Using the Liouville anti-derivative expression we can write From the standard (integer order) Barrow formula b a f (x)dx = f (b) − f (a) we obtain the expression Using (37) it comes because the GL and Liouville derivatives yield the result 0 for the derivative of the constant function.This formula reduces the α-order integral to one integer order.
3.1.Integrals in R 2 and R 3 The above result can be enlarged for double and triple integrals.This will be done next for rectangular domains, meaning that in R 2 the domain is obtained by joining rectangles so that the limiting line is made of segments parallel to one of the coordinated axes.Similarly in R 3 the domain is made of parallelepipeds with a boundary constituted by rectangles parallel to one of the coordinated planes.
Let us assume that the function f is dependent of a variable, x 1 , and a parameter, x 2 , that is kept fixed, f (x 1 , x 2 ).Consider the parametric integral and similarly, fixing x 1 , another one This leads to the following definition Definition 3. The FI on a rectangular region (a 1 , b 1 ) × (a 2 , b 2 ) is given by where each integration is performed by the fractional Barrow formula.
If the orders of the FD are equal (i.e., α 1 = α 2 = α) we put x α 1 x α 2 = S α to get then The integral in R 3 is obtained in a similar way.where each integration is performed by means of the fractional Barrow formula.

Conclusions
In this paper, a definition of fractional definite integral was introduced and the corresponding fundamental theorem of fractional calculus formulated.For this purpose, the mostly used fractional derivatives were studied and a compatibility with the definite fractional integral was saught.It was shown that the Grünwald-Letnikov and Liouville derivatives were suitable for such a definition.Finally, integrals on R 2 and R 3 were also proposed.