#### Abstract

Using Riemann-Liouville fractional differential operator, a fractional extension of the Lagrange inversion theorem and related formulas are developed. The required basic definitions, lemmas, and theorems in the fractional calculus are presented. A fractional form of Lagrange's expansion for one implicitly defined independent variable is obtained. Then, a fractional version of Lagrange's expansion in more than one unknown function is generalized. For extending the treatment in higher dimensions, some relevant vectors and tensors definitions and notations are presented. A fractional Taylor expansion of a function of -dimensional polyadics is derived. A fractional -dimensional Lagrange inversion theorem is proved.

#### 1. Introduction

The fractional calculus (FC) may be considered as an old and yet novel topic. It dates back to the end of the seventeenth century through the pioneering works of Leibniz, Euler, Lagrange, Abel, Liouville, and many others. In a letter to L'Hospital in 1695, Leibniz raised the possibility of generalizing the operation of differentiation to noninteger orders, and L'Hospital asked what would be the result of half-differentiating . Leibniz replied: It leads to a paradox, from which one day useful consequences will be drawn. The paradoxical aspects are due to the fact that there are several different ways of generalizing the differentiation operator to non-integer powers, leading to inequivalent results.

The fractional calculus (FC) generalizes the ordinary differentiation and integration so as to include any arbitrary real or even complex order instead of being only the positive integers (see, e.g., Samko et al. [1], Kilbas, et al. [2], Magin [3], and Podlubny [4]).

During the second half of the twentieth century till now, FC gained considerable popularity and importance. Many authors have explored the world of FC giving new insight into many areas of scientific research in physics, mechanics, and mathematics. Miller and Ross [5] pointed out that there is hardly a field of science or engineering that has remained untouched by the new concepts of FC.

Fractional derivatives provide an excellent as well as very powerful tool for the description and modeling of many phenomena in nature. There are many applications where the fractional calculus can be widely used, for example, viscoelasticity, electrochemistry, diffusion processes, control theory, heat conduction, electricity, mechanics, chaos and fractals, turbulence, fluid dynamics, stochastic dynamical system, plasma physics and controlled thermonuclear fusion, nonlinear control theory, image processing, nonlinear biological systems, astrophysics, and so forth, see for details [2–12] and the references therein.

In a very good book by Baleanu et al. [13], readers were given the possibility of finding very important mathematical tools for working with fractional models and solving fractional differential equations, such as a generalization of Stirling numbers in the framework of fractional calculus and a set of efficient numerical methods. Moreover, they introduced some applied topics, in particular, fractional variational methods which are used in physics, engineering, or economics. They also discussed the relationship between semi-Markov continuous-time random walks and the space-time fractional diffusion equation, which generalized the usual theory relating random walks to the diffusion equation.

Debbouche and Baleanu [14] introduced a new concept called implicit evolution system to establish the existence results of mild and strong solutions of a class of fractional nonlocal nonlinear integrodifferential system; then they proved the exact null controllability result of a class of fractional evolution nonlocal integrodifferential control systems in Banach space. As an application that illustrates their abstract results, they provided two examples.

Babakhani and Baleanu [15] considered a class of nonlinear fractional order differential equations involving Caputo fractional derivative with lower terminal at 0 in order to study the existence solution satisfying the boundary conditions or satisfying the initial conditions. They derived unique solution under Lipschitz condition. In order to illustrate their results they presented several examples.

Finally and roughly speaking, the fractional calculus may improve the smoothness properties of functions rather than the calculus with integer orders. The development of the FC theory is due to the contributions of many mathematicians such as Euler, Liouville, Riemann, and Letnikov. Several definitions of a fractional derivative have been proposed. These definitions include Riemann-Liouville, Grunwald-Letnikov,Weyl, Caputo, Marchaud, and Riesz fractional derivatives, see Miller and Ross [5] and Riewe [16]. Riemann-Liouville derivative is the most used generalization of the derivatives. It is based on the direct generalization of Cauchy's formula for calculating an -fold or repeated integral, see Oldham and Spanier [17].

In 1770, Lagrange (1736–1813) published his power series solution of the implicit equation. However, his solution used cumbersome series expansions of logarithms. [18, 19]. This expansion was generalized by Bürmann [20–22]. There is a straightforward derivation using complex analysis and contour integration; the complex formal power series version is clearly a consequence of knowing the formula for polynomials; so the theory of analytic functions may be applied. Actually, the machinery from analytic function theory enters only in a formal way in this proof. In 1780, Laplace (1749–1827) published a simpler proof of the theorem, based on the relations between partial derivatives with respect to the variable and the parameter, see [23, 24], Hermite (1822–1901) presented the most straightforward proof of the theorem by using contour integration [25–27].

In mathematical analysis, this series expansions is known as Lagrange inversion theorem, also known as the Lagrange-Bürmann formula, giving the Taylor series expansion of the inverse function. Suppose that where is analytic function at a point and . Then, it is possible to invert or solve the equation for such that on a neighborhood of , where is analytic at the point . This is also called reversion of series. The series expansion of is given by In this work, we will apply the concepts of fractional calculus to obtain a fractional form of the Lagrange expansion and some generalizations.

#### 2. Basic Definitions and Theorems

*Definition 1. *By , we denote the operator that maps a differentiable function onto its integer derivative; that is, ; by , we denote the integer integration operator that maps a function , assumed to be (Riemann) integrable on the compact interval , onto its primitive centered at ; that is, .

*Definition 2. *By and , we denote the -fold iterates of and , respectively. Note that is the left inverse of in a suitable space of functions.

Lemma 3. *Let be Riemann integrable on . Then, for and , one has
*

*Definition 4. *The operator , defined on Lebesgue space , denotes the Riemann-Liouville fractional operator of order . That is,

*Remark 5. *It is evident that , except for the fact that we have extended the domain from Riemann integrable functions to Lebesgue integrable functions (which will not lead to any problems in our development). Moreover, in the case , it is obvious that the integral exists for every because the integrand is the product of an integrable function and the continuous function . One important property of integer-order integral operators, is preserved by this generalization. That is,

*Definition 6. *Let and let , The Riemann-Liouville fractional differential operator of order is defined as such that. Then, . That is,

Lemma 7. *Let and let such that . Then, .*

*Proof. *Since yields .

Thus,

Theorem 8. * Let . Then, for every , .*

*Proof. *For , both operator, are the identity. For , let ; then,

Corollary 9. *Let be analytic in for some , and let .**Then,
**
The binomial coefficients for and are defined as
*

*Proof. * For the first two statements and we use the definition of the Riemann-Liouville integral operator and expand into a power series about . Since , the power series converges in the entire interval of integration and exchanges summation and integration. Then, we use the explicit representation for the fractional integral of the power function:
follows immediately. For the second statement, we proceed in a similar way; but we now expand the power series at and not at . This allows us again to conclude the convergence of the series in the required interval. The analyticity of follows immediately from the second statement.

To prove we use the relation
This allows us to rewrite the statement as
Differentiating times with respect to , we find
The classical version of Leibniz’ formula yields
which yields
By definition, if and . Thus, we may replace the upper limit in the inner sum by without changing the expression. The substitution gives
Using the fact that ,
And the explicit calculation yields
thus, follows directly.

#### 3. Fractional Form of Lagrange's Expansion in One Variable

We can use the standard form of Lagrange's expansion for one implicitly defined independent variable and the Definition 6 to obtain the fractional form of Lagrange's expansion as follows.

Let be a function of and in terms of another function such that Then, for any function for small . If is the identity then This classical result can be obtained using the following integral: Now we are going to introduce a fractional form of the Lagrange inversion formula.

Rewrite the integral (23) in the fractional form as Writing the delta function as an integral, we have The integral over then gives and we have On extracting the first term out of summation, set , and rearranging the terms then gives the result:

#### 4. Generalized Fractional Lagrange's Expansion in One Variable

We can the Lagrange's expansion in more than one unknown function as This classical result can be obtained using the following integral: Now, we are going to introduce a fractional form of the Lagrange inversion formula.

Rewrite the integral (29) in the fractional form as Writing the delta function as an integral, we have The integral over then gives and we have On extracting the first term out of summation, set , and rearranging the terms then gives the result

#### 5. Vector and Tensor Definitions and Notation

For the treatment in higher dimensions, consider the -dimensional space with orthogonal unit base vectors : Let , and the function be -dimensional vectors in this space such that where For any arbitrary differentiable function , we can introduce the following fractional gradient operator as .

*Definition 10. *Let be a domain of . Let is a scalar function that has absolutely continuous derivatives up to order on; then fractional gradient is defined as
where the partial derivatives are taken holding all other components of the argument fixed.

#### 6. The -Dimensional Polyadics (th-Order Tensors)

For arbitrary -dimensional vectors we use an extension of the notion of an -dimensional dyadic (second-order tensor): to define the th-order tensors: We might call and “polyadics,” since the special cases for are known, respectively, as dyadics, triadics, and tetradics [21]. The following defined scalar products then follow quite naturally from (34): and, in general, define the th scalar product: Particular examples of th order tensors to be used are

Theorem 11. *Assume that , and let . Then,
**
holds almost everywhere on . If additionally or , then the identity holds everywhere on .*

*Proof. *We have
Thus, we can write
Using Fubini’s theorem to interchange the order of integration yields
The substitutions yields the new limits of integration as follows: when , and when :
Iterating the Euler Beta integral yields
hold almost everywhere on .

Moreover, by the classical theorems on parameter integrals, if , then also , and therefore , and too. Thus, since these two continuous functions coincide almost everywhere, they must coincide everywhere. Finally, if and , we have, by the result above
almost everywhere. Since is continuous, and once again we may conclude that the two functions on either side of the equality almost everywhere are continuous, thus they must be identical everywhere.

Theorem 12. *Assume that . Moreover let . Then,
*

*Proof. *The key for proof is using the semigroup property of the integral operators, the assumption on , and the definition of the Riemann-Liouville differential operator:
The proof that is quite straightforward
Thus,

Theorem 13. *Assume that . Moreover, let and . And let be the fractional gradient operator. Then,
*

*Proof. *By the assumption on and the successive application of fractional gradient operator, we have :
Thus using the theorem, the results follow directly:

#### 7. Fractional Taylor Expansion of a Function of -Dimensional Polyadics

We have the classical Taylor expansion for the -independent variables as In the light of the above definitions and theorems, we can state the following theorem.

*Definition 14. *Let ; then the fractional Riemann-Liouville multiple integrals and partial derivatives with respect to are;

Theorem 15. *Let and . Assume that is such that , where is the domain of . Then,
*

*Proof. *Because of our assumption about that implies the continuity of , there exists some such that
This is a classical partial differential equation of order for ; its solution is easily seen to be of the form:
Thus, by definition of ,
annihilates every summand in the sum. And due to Theorem 8, we obtain
Next, we apply the operator to (63), finding
Using (64), we obtain
Upon extracting the first term out of summation, we can write

This result can be written conveniently in the Polyadics notation , extending the upper limit of summation to to absorb the remainder, as

#### 8. Fractional -Dimensional Lagrange Expansion

Theorem 16. *Let such that
**
then it is still possible to invert or solve the equation for such that on a neighborhood of using the fractional calculus as follows:
**
where and .*

*Proof. * It is possible to invert or solve the equation for in a neighborhood of . Using the fractional Taylor expansion.
Extracting the first term out of the summations, we obtain the result directly

#### Acknowledgments

The author is deeply indebted and thankful to the deanship of the scientific research and its helpful and distinct team of employees at Taibah University, Al-Madinah Al-Munawwarah, Saudi Arabia This research work was supported by a Grant no. 3017/1434. Also he wishes to express his deep gratitude and indebtedness to the referees of his manuscript for all very fruitful discussions, constructive comments, remarks, and suggestions.