- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

International Journal of Differential Equations

Volume 2012 (2012), Article ID 521750, 38 pages

http://dx.doi.org/10.1155/2012/521750

## Generalized Multiparameters Fractional Variational Calculus

Mechanical Engineering and Energy Processes, Southern Illinois University, Carbondale, IL 62901, USA

Received 2 April 2012; Accepted 8 August 2012

Academic Editor: Fawang Liu

Copyright © 2012 Om Prakash Agrawal. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper builds upon our recent paper on generalized fractional variational calculus (FVC). Here, we briefly review some of the fractional derivatives (FDs) that we considered in the past to develop FVC. We first introduce new one parameter generalized fractional derivatives (GFDs) which depend on two functions, and show that many of the one-parameter FDs considered in the past are special cases of the proposed GFDs. We develop several parts of FVC in terms of one parameter GFDs. We point out how many other parts could be developed using the properties of the one-parameter GFDs. Subsequently, we introduce two new two- and three-parameter GFDs. We introduce some of their properties, and discuss how they can be used to develop FVC. In addition, we indicate how these formulations could be used in various fields, and how the generalizations presented here can be further extended.

#### 1. Introduction

For over a century, many researchers have been in search for a fundamental law that can be used to describe the behavior of the nature. One law that comes very close to it is the universal law of extremum which states that the nature always behaves in a way such that some quantity is an extremum. A catenary takes a shape so that the total potential energy is minimum, light travels from a point to another so that the travel time is minimum, a particle in a flow takes a path of least resistance, and even in social settings, we behave so that our conflict with others within our conviction is minimum. Related laws, principles, and theories have been developed in almost every field of science, engineering, mathematics, biology, economics and social science. For example, applications of such laws, principles, and theories in continuum mechanics, classical and quantum mechanics, relativistic quantum mechanics, and electromagnetics could be found in [1–5] and many other textbooks, monographs, and papers. Opponent of the universal principle may argue that nature behaves in its own way, and the extremum principles are our creations where we design a functional that is extremum for the nature’s trajectory. Whatever may be the reality, the underlying theories have advanced our understanding of the nature tremendously.

The field that deals with the mathematical theories of the extremum principles is known as the variational calculus. Excellent books have been written in this field, see for example [6, 7]. These books provide not only the foundations for theoretical work in the field, but they have also been a basis for many numerical techniques (see, [8, 9]). However, the traditional variational calculus subject has one major drawback; it deals with functionals containing integer-order derivative terms only. Recent progress in last two decades have demonstrated that many phenomena in various fields of science, mathematics, engineering, bioengineering, and economics are more accurately described using fractional derivatives. As a result, many books, monographs, and papers have been written recently on this subject (see, e.g. [10–21]). We assume that many fractional models would follow the universal law of extremum. If this is true, then it is very likely that a variational calculus that deals with fractional derivatives would be necessary. In other words, we need fractional variational calculus.

The subject of fractional variational calculus was initiated by Riewe [22, 23] in 1996. Riewe was interested in developing a variational formulation for a linear damper. He observed that a quadratic term of type in a functional leads to a second-order derivative term of type in the resulting differential equation. Here is the derivative operator. Therefore, he argued that a first-order derivative term of type in a differential equation will come from a quadratic term of type in the functional. Here is a half-order derivative operator (Precise definitions of fractional derivative would be discussed later in the paper.) Using this hypothesis, he proceeded to develop a variational formulation in terms of fractional derivatives. Subsequently, he developed fractional Lagrangian, fractional Hamiltonion, and fractional mechanics.

Klimek [24, 25] and Agrawal [26] brought this subject to the main stream and initiated the field of fractional variational calculus. These authors identified the key integration-by-parts formulas for fractional derivatives, and showed that using these formulas a fractional variational formulation can be obtained in the same way as it is done for integer variational formulation. Note that variational calculus has been applied to an extensively large number of problems, theories, and formulations most of which could be reexamined in the light of fractional variational calculus. Thus, the above work has opened significant opportunities for many new research.

Recently, the field of fractional variational calculus has indeed grown very rapidly. A citation search of [26] and some related papers suggests that in the last 10 years over 300 papers have been published which are directly related to fractional variational calculus; here we cite a few of them [27–36]. These papers further (1) develop fractional variational calculus and Fractional Euler-Lagrange Equations in terms of fractional derivatives not considered earlier, (2) derive trasnversality conditions for fractional problems, (3) propose new fractional Lagrangians and fractional Hamiltonians and fractional mechanics, and (4) develop applications of fractional derivatives in optimal control and fractional inverse problems. Recently, Klimek has written a book dedicated to fractional variational calculus and analytical techniques to solve problems resulting from fractional variational formulations [37]. Reference [38] presents a two- and a three-parameter generalizations of fractional variational calculus. It is demonstrated that by setting these parameters to different values we obtain several different fractional variational formulations presented previously. These citations are clear indications that significant progress has been made in the area of fractional variational calculus. However, compared to the progress that has already been made in ordinary variational calculus, this progress in fractional variational calculus is very small.

In this paper, we provide further generalization of fractional variational calculus. Specifically we introduce Hadamard type, Erdélyi-Kober type fractional integral and fractional derivatives, and fractional integrals and derivatives of a function (say ) with respect to another function (say ). We also introduce some new fractional integrals and derivatives Caputo, Hilfer, and Riesz types. This leads to two-, three-, and four-parameter generalized fractional derivatives of a function with respect to another function. We develop integration-by-parts formulas and Euler-Lagrange equations in terms of these new two-, three- and four-parameters fractional derivatives. It is demonstrated that by taking different values for different parameters and different function , we obtain many old and many new fractional derivatives and fractional variational formulations. We also introduce fractional Lagrangians and Hamiltonians in terms of these derivatives, and develop a more general fractional mechanics. Since, function could be selected from a large set, it provides a large number of fractional variational formulations for modeling purpose. Finally, we discuss how this work can be extended further.

At this point, we would like to emphasize that a comprehensive treatment and an excellent review of many generalized fractional operators proposed in the field could be found in [13, 39] (In this regards, please also see [19] and many references cited in [13, 39]). Generalizations of integral operators with specific weights could be found in [12]. In contrast to these references, some of the fractional operators proposed here are more general. In addition, our focus here is to develop some of the theories for generalized fractional variational calculus in terms of these operators, and provide an outline for other formulations in the field. We consider here functions dependent on one parameter. However, the theories developed here could easily be extended to field variable and distributed order systems ([21]).

#### 2. Preliminaries

In this paper, we will introduce several general multiparameter fractional integrals and derivatives, and show that many specific integrals and derivatives can be obtained from these general derivatives. For ease in the discussion to follow and to make this paper self-contained, we first introduce several symbols and notations, and provide some preliminaries. A large part of these symbols and definitions could be found in [12, 19, 38]. We shall denote the order of the fractional integrals and derivatives as . In [12, 19], in general, is taken as a complex number, and restrictions are imposed on it as necessary. Same approach can be taken here. However, in our discussion to follow, we shall implicitly consider that is a positive real number. We shall assume that and are real such that , and consider the domain of the functions and operators as , although in some cases, the domain of the functions and operators may not contain some isolated points of . Further, we shall consider to be a finite domain. But, in many cases, could be and could be , and some special cases could be derived by setting . However, these would be left as an exercise. Exception to these would be noted as necessary. We shall further assume that our functions are “sufficiently good” so that the operations considered on them are valid.

In this paper, we consider several fractional integrals and derivatives. We begin with the Riemann-Liouville and the Caputo fractional integrals and derivatives.

##### 2.1. Riemann-Liouville and Caputo Fractional Integrals and Derivatives

Many fractional derivatives are defined using the left/forward and the right/backward Riemann-Liouville fractional integrals (RLFIs). These integrals are defined as follows.

*Left/Forward Riemann-Liouville fractional integral of order *

* Right/Backward Riemann-Liouville fractional integral of order *
where is the Gamma function, and is the order of integration. Here , however, we shall restrict our attention to . The fractional integral operators and satisfy the semigroup property, namely,

We now consider the *reflection operator * which is defined as
Operator satisfies the following identities,
An advantage of operator is that one needs to examine the properties of only, and obtain the properties of by using the properties of and .

Using (2.1) and (2.2), the left/forward and the right/backward Riemann-Liouville fractional derivatives (RLFDs) of order , (), an integer, are defined as follows.

*Left/Forward Riemann-Liouville fractional derivative of order *

*Right/Backward Riemann-Liouville fractional derivative of order *
where represents the ordinary differential operator. Operators and are the left inverse of operators and , that is, they satisfy the following identities
where is the identity operator. However, in general

The left/forward and the right/backward Caputo fractional derivatives (CFDs) of order are defined as

*Left/Forward Caputo fractional derivative of order *

*Right/Backward Caputo fractional derivative of order *
Note that the ordinary derivative operators are applied in the RLFDs after the fractional integrals whereas in the CFDs before the fractional integrals. Therefore, the differentiability required of by the CFDs is higher than those by the RLFDs.

Operators , , , , and satisfy the following identities: The RLFDs and the CFDs are related by the following formulas [12, 19]: where represents the th derivative of with respect to . Equation (2.13) can be obtained using the identity and the semigroup property of the integral operators, and (2.14) can be obtained using the property of the operator.

The integral operators and and the derivative operators , and satisfy the following fractional integration by parts formulas [12] Similarly, operators and satisfy and The conditions under which (2.16) and (2.17) are valid can be found in [12]. It will be implicitly assumed that these conditions are satisfied. In a more general setting, operators , , , and satisfy the following fractional integration by parts formula, Here, and must be interpreted as and , respectively. Equations (2.18) and (2.19) can be obtained by considering the Riemann-Liouville and the Caputo fractional derivatives in terms of fractional integrals and ordinary derivatives and (2.16). Conditions under which (2.18) and (2.19) are true can be found in [37]. Using the properties of functions and , these equalities can be further specialized. Equations (2.17) to (2.19) have played key roles in developing fractional variational calculus.

##### 2.2. Hadamard Fractional Integrals and Derivatives

The left/forward and the right/backward Hadamard fractional integrals (HFIs) of order are defined as follows [12, 19].

*Left/Forward Hadamard fractional integral of order *

*Right/Backward Hadamard fractional integral of order *
Like the RLFIs, the HFIs also satisfy the semigroup properties, that is,

The left/forward and the right/backward Hadamard fractional derivatives (HFDs) of order are defined as follows [12, 19].

*Left/Forward Hadamard fractional derivative of order *

*Right/Backward Hadamard fractional derivative of order *
where, as before, . The left and the right Hadamard fractional derivative operators (HFDOs) and are the left inverse of the left and the right Hadamard Fractional Integral Operator (HFIOs) and , respectively, that is, they satisfy the following identities
Note that the left and the right HFDOs are obtained by applying and to the left- and the right HFIOs, respectively. We shall define another set of fractional derivatives by changing the order of the derivative and the integral operators in (2.23) and (2.24). Due to their similarity with Caputo derivatives, we call them Hadamard-Caputo derivatives. Thus, the left/forward and the right/backward Hadamard-Caputo fractional derivatives (HCFDs) of order are defined as follows:

*Left/Forward Hadamard-Caputo fractional derivative of order *

*Right/Backward Hadamard-Caputo fractional derivative of order *

The HFDOs and are related to the HCFDOs and by the following relations where is the time derivative operator. The proofs of (2.28) are the same as those for (2.13) and (2.14).

We can develop integration by parts formula for Hadamard operators also. It can be demonstrated that the HFIOs and satisfy the following integration by parts formula: One can prove (2.29) using the definitions of the HFIOs and the Dirichlet formula.

Similar to (2.17) to (2.19), the HFDOs and the HCFDOs satisfy the following identity and in a more general setting, the HFDOs and the HCFDOs satisfy the following fractional integration by parts formula

As stated earlier, one can take or and , and include some weight functions in the integral to obtain some other types of HFIs and HFDs. However, these will be considered later.

Equations (2.20) to (2.32) provide sufficient number of formulas to develop Euler-Lagrange formulations in terms of Hadamard fractional derivatives. For the time being, we proceed to define the Erdélyi-Kober type fractional integrals and derivatives, and develop some of their properties pertinent to fractional variational calculus.

##### 2.3. Erdélyi-Kober Fractional Integrals and Derivatives

The left/forward and right/backward Erdélyi-Kober fractional integrals (EKFIs) are defined as [12, 19] as follows:

*Left/Forward Erdélyi-Kober fractional integral of order *
and *Right/Backward Erdélyi-Kober fractional integral of order *

As pointed out earlier, by setting to or and to , one can obtain several other types of EKFIs. Indeed, many of such integrals are defined and discussed in [12, 19]. Many of the formulations discussed here can directly be applied to these other EKFIs.

These integrals satisfy the following semigroup properties:

The Erdélyi-Kober fractional derivatives (EKFDs) of order corresponding to the EKFIs are defined as.

*Left/Forward Erdélyi-Kober fractional derivative of order *
and *Right/Backward Erdélyi-Kober fractional derivative of order *
For sufficiently good function , the Erdélyi-Kober fractional differential operators (EKFDOs) and are left inverse of the Erdélyi-Kober fractional integral operators (EKFIOs) and , respectively, that is, they satisfy the following identities,
The EKFIs and the EKFDs satisfy the following integration by parts formulas:
Equation (2.39) could be found in [12, 19], and (2.40) could be derived using the definitions of EKFIs and EKFDs and the Dirichlet formula.

Equations (2.33) to (2.40) are sufficient to develop some fractional variational formulations in terms of EKFIs and EKFDs. Note that in the EKFDs, the EKFIOs are applied first and some derivative operators are applied next. One can also take some derivative operators first and the EKFIOs next to define some Caputo type EKFIs and EKFDs. Indeed, one can show that the following relation is valid
where
is the left/forward Erdélyi-Kober-Caputo type *fractional derivative of order *. Following above discussion, a right/backward Erdélyi-Kober-Caputo type Fractional Derivative (EKC-FD) of order and related integration by parts formula could be developed. However, note that Erdélyi-Kober operators may lead to some nonstandard cases. For example, in (2.35), one cannot simply interchange the first two operator, as it can be done in (2.3). Therefore, the EKC-FDs and related identities and their generalizations will be considered in the future. Note that references [12, 19] also introduce operators and defined as and . These operators link the Erdélyi-Kober operators to the Riemann-Liouville operators, which could simplify some of the formulations. Properties of these formulations in the context of fractional variational calculus will be examined in the future.

##### 2.4. Modified Erdélyi-Kober Fractional Integrals and Derivatives

It was pointed out above that the Erdélyi-Kober operators defined by (2.33), (2.34), (2.36), and (2.37) lead to some nonstandard cases. For example, in (2.35) the first two operators cannot be interchanged. Furthermore, in order for the relations in (2.35) to be valid, some of the parameters of the operators must be related (e.g., notice the presence of in the first two operators in (2.35)). To overcome this difficulty, we define modified Erdélyi-Kober fractional integrals and derivatives in the following way.

*Left/Forward modified Erdélyi-Kober fractional integral (MEKFI) of order *
and *right/backward modified Erdélyi-Kober fractional integral (MEKFI) of order *
It should be pointed out that like in the case of EKFIs, one can obtain several other types of MEKFIs by setting to or and to . Such substitutions may lead to specialized cases, and may result in simplified formulations. These cases will be considered in the future. In any case, the formulation here will also be applicable to these specialized MEKFIs.

These integrals satisfy the following semigroup properties The proof of these identities follows the same steps as those for the proof of (2.35). Note that in the modified definition, the operators and commute.

We define the modified EKFDs of order corresponding to the modified EKFIs as follows.

*Left/Forward modified Erdélyi-Kober fractional derivative of order *
and *right/backward modified Erdélyi-Kober fractional derivative of order *
Note that in the case of the modified EKFIOs and the modified EKFDOs, the value of the parameters and are not modified.

For sufficiently good function , the modified EKFDOs and are left inverse of the modified EKFIOs and , respectively, that is, they satisfy the following identities: The modified EKFIs and the modified EKFDs satisfy the following integration by parts formulas The method of proof for (2.49) and (2.50) are the same as that for (2.39) and (2.40).

Equation (2.50) can be used to obtain a fractional variational formulation in terms of modified EKFDs. However, such formulations limit the terminal conditions and obscure the posibility of other terminal and transversality conditions. Further, note that in (2.46) and (2.47), the integral operators are applied first and the derivative operators are applied next. We can interchange the operation, and these leads to Caputo type derivatives. Accordingly, we define the modified Caputo type Erdélyi-Kober fractional derivatives as follows.

*Left/Forward modified Caputo type Erdélyi-Kober fractional derivative (MCEKFD) of order *
and *right/backward modified Caputo type Erdélyi-Kober fractional derivative (MCEKFD) of order *
The left MEKFD and the left MCEKFD are related in the following way:
A similar relation exits between the right MEKFD and the right MCEKFD. The MEKFDs and MCEKFDs satisfy the following integration by parts formula:
A similar relationship can be obtained relating the right MEKFD and the left MCEKFD.

##### 2.5. Weighted/Scaled Fractional Integrals and Fractional Derivatives of a Function with Respect to Another Function

In this section we define the left/forward and the right/backward fractional integrals and fractional derivatives of a function with respect to another function and weight/scale , and investigate some of their properties. We assume that is an increasing positive monotone function on having a continuous derivative on . We further assume that function is “sufficiently good.”

We now define the left/forward and the right/backward weighted/scaled fractional integrals of a function with respect to another function as follows.

*Left/Forward weighted/scaled fractional integral of order ** of a function ** with respect to another function ** and weight *

*Right/Backward weighted fractional integral of order ** of a function ** with respect to another function ** and weight *
We shall denote them as the left/forward and the right/backward generalized fractional integrals (GFIs). It should be pointed out that these integrals also contain the scaling/weight function and the function with respect to which the function is integrated. In our case, functions and will often remain the same, and therefore, it is not necessary to explicitly mention of the presence of these functions in the definitions of the integrals.

Like the GFIs, we define the left/forward and the right/backward weighted/scaled fractional derivatives of a function with respect to another function and weight/scale as follows.

*Left/Forward weighted fractional derivative of order ** of a function ** with respect to another function ** and weight *

*Right/Backward weighted fractional derivative of order ** of a function ** with respect to another function ** and weight *
where , and . These derivatives contain one parameter only, namely the order of the derivative. (In reality, functions and would also introduce additional parameters. For the time being, we shall keep these functions the same, and therefore, calling these derivatives as one parameter derivatives is justified.) In later sections, we shall define derivatives containing many more parameters. Furthermore, note that in these derivatives, the GFI Operators (GFIOs) are applied first and some derivative operators are applied next, and accordingly, they are like Riemann-Liouville fractional derivatives. Later, we shall consider Caputo type derivatives in which the derivative operators would be applied first, and the GFIOs next. A name that describes all these would be very long. For brevity, we denote these derivatives as One-Parameter Type-1 generalized fractional derivatives (1PT1GFDs). We choose not to call these derivatives as the Riemann-Liouville fractional derivatives for the reason given below. We further define
Here, and are like the left/forward and the right/backward operators and , respectively, except that they also contain functions and . These are new integer order operators. variational calculus in terms of these operators will be considered somewhere else. In this paper, we shall focus on the role of these operators in formulating generalized fractional variational calculus. Using (2.59), (2.57), and (2.58) can be written as
Here subscript “1” is added to indicate that these are Type-1 fractional derivatives. In the same faashion, we introduce One-parameter Type-2 generalized fractional derivatives (1PT2GFDs) as
Note that in these fractional derivatives, the differential operators are applied first, and the GFIOs are applied next, accordingly they are Caputo type fractional derivatives.

Before we proceed further, we would like to note that the fractional integrals and derivatives of a function with respect to another function defined in [12, 19] does not consider the weight/scale function. Further, by taking in (2.55) to (2.58), we obtain the fractional integrals and derivatives , , , and defined in [12, 19]. However, in our case, need not be 1. Thus, the fractional integrals and derivatives defined here are more general than those given in [12, 19]. Also note that in [19], the authors define a substitution operator and its inverse such that . One of the major advantages of these operators is that they link the fractional integrals and derivatives , , , and to the Riemann-Liouville fractional integrals and derivatives defined by (2.1), (2.2), (2.6), and (2.7). Accordingly, the properties of , , , and could be obtained from the operators defined in (2.1), (2.2), (2.6), and (2.7). Operators similar to and could be defined for the present case. However, in this paper, we shall use the definitions directly to obtain the properties of the operators defined in (2.55) to (2.62).

We now list several properties of operators , , , , , and . Operators and satisfy the following semigroup properties: Operators and are left inverse of the operators and , that is, they satisfy the following relations: Operators and are related by the following formula One can develop a similar relation relating to . The above relation shows that like in the case of Riemann-Liouville and Caputo derivatives, if and all its derivatives of order upto are 0 at , then Type-1 and Type-2 derivatives and are the same, that is, Similarly, the following relationship holds: Operators , , , , , and satisfy the following integration by parts formulas: Equations (2.69) to (2.72) would play a key role in deriving a generalized fractional variational formulation. Equation (2.70) will implicitly account for some boundary conditions. In contrast, (2.71) and (2.72) would provide either natural or the boundary conditions. Equations (2.69) to (2.72) would play a key role in developing adjoint equations. However, this issue will be considered elsewhere.

References [12, 19] obtain some special cases of fractional integrals and derivatives of a function with respect to another function by setting or and . Similarly, one can obtain special cases of (2.55) to (2.58) by setting or and . In particular, by setting and we obtain generalized Liouville type fractional integrals and derivatives on the half-axis and on the whole axis .

The generalized fractional integrals and derivatives defined in (2.55) to (2.62) encompass many of the fractional integrals and derivatives defined earlier. For example, for and , the generalized fractional integrals and derivatives reduce to Riemann-Liouville and Caputo fractional integrals and derivatives, that is, we get For and , the generalized fractional integrals and derivatives reduce to Hadamard type fractional integrals and derivatives, that is, we get For and , the generalized fractional integrals and derivatives reduce to Modified Erdélyi-Kober type fractional integrals and derivatives, that is, we get It could also be verified that for equal to , and the semigroup, the left inverse, and the integration by parts type identities for the generalized fractional integrals and derivatives, namely (2.63) to (2.72), reduce to those for Riemann-Liouville, Hadamard, and Modified Erdélyi-Kober fractional integrals and derivatives. It should be pointed out that functions are not limited to , and . In fact this choice is significantly large, and for each choice, one would obtain a different set of fractional integrals and derivatives. The rational for not calling and as the Riemann-Liouville fractional derivatives should now be clear; for different and , they lead to different fractional integrals and derivatives. In literature, these fractional integrals and derivatives have been called by different names. Thus, to avoid confusion, we prefer to call and as the 1PT1GFDs, and and as the 1PT2GFDs.

#### 3. Fractional Variational Formulation in Terms of One-Parameter Generalized Fractional Derivatives

In this section, we present several fractional variational formulations in terms of one parameter generalized fractional derivatives. We shall consider formulations in terms of one variable and one fractional derivative term, specified and unspecified terminal conditions, one variable and multiple fractional derivative terms with different order of derivatives, multiple variables and multiple derivative terms but the same order of derivatives, geometric constraints, and parametric constrains to name a few. The approach presented here will also be applicable to multivariables and multiple fractional derivative terms with different order of derivatives, free end points, free end-point constraints, multidimensions, and many other formulations can also be considered. As a matter of fact almost all variational formulations can be recast in terms of generalized fractional derivatives, and these will be considered in the future. The derivations of almost all formulations follow the same pattern, and for this reason, a fractional variational formulation would be given in detail for a simple fractional variational problem only. For other fractional variational problems, the final Euler-Lagrange equation will be given but the details would be omitted.

##### 3.1. A Simple Fractional Variational Formulation

In this subsection, we develop an Euler-Lagrange formulation for a simple fractional variational formulation. The functional considered in this case may contain the left and the right integrals and derivatives, and the derivatives could be of Type-1, Type-2, or both. In the simplest fractional variational problem considered here, we take only one fractional derivative term, namely the term . The approach for functional containing other fractional integrals and derivatives would be the same. Accordingly, the fractional variational problem is defined as follows: among all functions which are continuously differentiable on find the function for which the functional is an extremum. For simplicity, we assume for the time being that . In case the functional is given in terms of , one can use (2.66) to write the functional in terms of .

At this stage, two points can be made which are the same as those made in [38] for another set of fractional derivatives. We repeat them here for completeness and to indicate that those points equally apply here. First, note that fractional integration introduces some degree of continuity, and therefore the differentiability requirements of could be relaxed. For certain class of problems and in several numerical schemes, the differentiability of would be required on only a finite subset of , and the derivative of could be discontinuous at a finite set of points. Further, we have implicitly assumed that function has continuous first and second partial derivatives with respect to all its arguments. In all formulation to follow, these conditions will be implicitly assumed. In many cases, mathematical operations performed will determine the class of functions being considered. Second, we have not specified the terminal conditions yet, because in fractional variational formulations, the forms of the necessary conditions are tied to the specified terminal conditions. This will be demonstrated shortly, and subsequently the boundary conditions will be specified.

To derive the necessary conditions, we define and substitute it in (3.1) to obtain in terms of and . Here is the desired solution, is an arbitrary function consistent with the boundary conditions, and is a small real number. Thus, for the specified , would be a function of only. We differentiate with respect to , and set and the resulting equation to 0 to obtain Using integration by parts formula given by (2.72), we obtain Here, . Typically, in variational formulations the boundary terms suggest the geometric and the natural boundary conditions. Thus, (3.4) suggests that for this case, and would be the appropriate geometric boundary conditions. These boundary conditions are the same as those considered in ordinary variational calculus. However, this leads to fractional natural boundary conditions. It can be verified that if (2.71) is used for integration by parts formula, then one obtains fractional geometric boundary conditions and regular natural boundary conditions.

Let us assume that and are specified. In this case, , and using the fundamental lemma of variational calculus, (3.4) leads to the following Euler-Lagrange equation: It is the necessary condition for the extremum. Let us now assume that is not specified. In this case, (3.4) leads to A similar condition is obtained if is not specified. Equation (3.6) and its variations are known as the natural boundary conditions.

As a quick modification of the above problem, assume that the functional contains an additional term . Clearly, would be unknown, otherwise would be a constant which could be removed from the functional without altering the problem. Term would lead to an additional term in (3.4) of type and the natural boundary conditions as In the case of greater than 1, (3.4) is replaced with where . Equation (3.9) leads to the same Euler-Lagrange equation as that given by (3.5). It further suggests that the geometric boundary conditions should be specified or the natural boundary conditions should be considered. Here , and . The above discussion assumes that all , , and are independent. In case they are not independent, then the natural boundary conditions are modified accordingly.

It should be pointed out that and are not ordinary differential operators. Therefore, derivatives , , are not the same as the ordinary derivatives . To demonstrate this, take , and . In this case, we have , that is, is a weighted/scaled combination of and , and . Therefore, when writing the geometric and/or natural boundary conditions, one must exercise caution.

In the above functional, we considered only the left fractional derivative. If the functional contains the right fractional derivative also, then the functional is written as where and . For simplicity, we shall assume that . It can be shown that for this functional the above approach leads to the following Euler-Lagrange equation and the following condition at the boundary points Note that in this case, we have left the natural boundary conditions, and together because we cannot treat both and as independent and they are not equal either. As a result, we cannot even write where , , . However, for , , , , and are linearly related to and . These linear relations would be necessary to separate the natural boundary conditions. Terms and could still be taken as the geometric boundary conditions.

In the formulation above, we have considered Type-2 fractional derivatives in the functional. Using (2.66), the functional can be written in terms of Type-1 fractional derivatives. However, Type-2 fractional derivatives were considered for two reasons. First, Type-2 fractional derivatives lead to geometric boundary conditions in terms of the desired function and its ordinary derivatives at the boundary points. In contrast, Type-1 fractional derivatives lead to fractional geometric boundary conditions. Many engineers and scientists avoid fractional geometric boundary conditions with remarks that these conditions are nonphysical. For this reason, Type-2 fractional derivatives would be more appealing. Second, treatment of the geometric and the natural boundary conditions is quite involved due to the presence of and . Fractional geometric boundary conditions and the associated variational formulations will be considered elsewhere.

##### 3.2. Multiorder and Multiterm Fractional Variational Formulation

We now consider several variations of the above formulation. As a first variation, assume that the functional contains left Type-2 fractional derivatives where the order of the derivatives , satisfy the following condition, . For this case, the functional is given as For this functional, following the above approach, we obtain the Euler-Lagrange equation as This is a straight forward generalization of (3.6). Further, assume that both and are independent. In this case, at point either must be specified (geometric boundary condition), or the identity (natural boundary condition) must be satisfied. A similar condition applies at point .

Assume now that the functional is of the following type: such that , and , . In this case, the Euler-Lagrange equation is given as and the boundary terms must satisfy the identity,