Variational Iteration Method for Solving Fractional Integro-Differential Equations with Conformable Differointegration

: Multidimensional integro-differential equations are obtained when the unknown function of several independent variable and/or its derivatives appear under an integral sign. When the differentiation or integration operators or both are of fractional order, the integral equation in this case is called a multidimensional fractional integro-differential equation. Such equations are difﬁcult to solve analytically; therefore, as the main objective of this paper, an approximate method—which is the variational iteration method—will be used to solve this type of equation with conformable fractional-order derivatives and integrals. First, we drive the iterative sequence of approximate solutions using the proposed method, and then, under certain conditions over the kernel of the integro-differential equation, prove its convergence to the exact solution. Two illustrative examples, linear and nonlinear, are given, and their approximated solutions are simulated using computer programs in order to verify from the reliability and applicability of the proposed method.


Introduction
Recently, fractional differential equations, fractional calculus and fractional differential integral equations became highly important in several branches of science and engineering because many mathematical models are used to formulate different phenomena, such as mechanics, physics, chemical kinetics, astronomy, biology, economics, potential theory and electrostatistics, which are modeled using integro-differential equations [1][2][3][4][5][6].In 2016, Shi [7] introduced a formula of mild solutions for impulsive fractional evolution equations.
Many academic researchers continue to be interested in the use of fractional differential equations and/or integral equations, which are based on the development and applications of fractional calculus [8,9].Nonlinear fractional integral equations and integro-differential equations are notoriously difficult to solve analytically.Moreover, accurate solutions for these equations are extremely rare.As a result, various authors have taken an interest in numerically solving these problems, particularly after the big revolution in computer application.Among the techniques used to solve integro-differential equations are the multistep methods [10], Adomian decomposition method (ADM) [11], homotopy perturbation method (HPM) [12][13][14], homotopy analysis method (HAM) [15], variational iteration method (VIM) [16] and so on.
The VIM has been successfully applied to solve many problems in different fields of mathematics and its applications.For example, He was the first researcher to propose the use of VIM to solve linear and nonlinear differential and integral equations [17].In 1998, He used VIM to solve the classical Blasius equation, ref. [18] and in 1999, he provided the approximate solutions for some well-known nonlinear problems [19].In 2000, He used the VIM to solve autonomous ordinary differential equations.Moreover, in 2006, Soliman applied the VIM to solve the kdv-Burger's and Lax's seventh-order kdv equations.In the same year, VIM was applied by Abulwafa and Momani [20] to solve a nonlinear coagulation problem with mass loss and Odibat et al. employed the VIM to solve nonlinear differential equations of fractional order in 2006.In 2006, the VIM was also utilised by Bildiki et al. to solve a variety of problems, including nonlinear partial differential equations, Dehghan and Tateri's Fokker-Plank equation and quadratic Riccati differential equations with constant coefficients.Wang [21] used the VIM to solve integro-differential equations in 2009, while Sweilam used the VIM to solve both linear and nonlinear boundary value problems of the fourth-order integro-differential equations.
In 2009, Wen-Hua Wang used the VIM to solve certain types of fractional integrodifferential equations [21].In 2011, Muhammet and Adin used the VIM to solve the problem of nonlinear fractional integro-differential equations [20].
If an exact solution exists, the VIM provides rapidly convergent consecutive approximations to the precise solution; otherwise, a few approximations might be employed for numerical results.
In this paper, we shall present the VIM used to solve integro-differential equations with conformable fractional order differointegration of the form: where K is given continuous function, 0 here, T α is understood as a conformable fractional derivative of order α, while I β and I γ stands for conformable fractional-order integrals of order β and γ, respectively.

Main Concepts of Factional Calculus
Among the most important definitions of fractional-order derivatives or integrals which will be used next in this paper is the conformable type, which is more simple than other definitions and more stable in comparison with the nonfractional (or integer order) derivatives and integrals.
Definition 1 (Conformable Fractional-Order Derivative [22]).Given a function f : [a, ∞) → R, then the left conformable fractional derivative of order α can be defined as: The right conformable fractional derivative of order α ∈ (0, 1] terminating at b of f , is defined by: Definition 2 (Conformable Fractional-Order Integral [23]).Given a continuous function f : [a, ∞) → R, then the left conformable fractional integral of order α of f is: where the integral is considered as the usual Riemann improper integral and a ≥ 0. On the other hand, in the right case, we have: Among some different properties concerning fractional integrals and derivative which are very useful in applications are the following, where T α and I α refer to fractional-order conformable and integral, respectively: T α (c T α (x p ) = px p−α , for all p ∈ R.

5.
T α (λ) = 0, for all constant functions f (x) = λ.6. T where f is any continuous function in the domain of b α I.

Variational Iteration Method
The essential aspect of the VIM, as previously stated in the literature, is that the solution of a mathematical problem under the linearization assumption is utilized as a starting approximation or trial function for the next successive approximate solution to the problem under certain conditions [2].
Consider the following general nonlinear equation in operator form to demonstrate the VIM's essential concept [3]: and suppose that Equation (1) may be decomposed as: where A is any operator that may be decomposed into linear and nonlinear operators L and N, respectively, and g(x) is any function that is referred to as the nonhomogeneous term.Equation ( 5) may be solved iteratively using the VIM by using the correction functional defined by: where λ is the general Lagrange multiplier that may ideally be discovered using variational theory the nth approximation of the subscript n denotes the solution u, and ũn is considered a restricted variation, i.e., δ ũn = 0 where the δ is the first variation [3].

Applications of the Vim for Multidimensional Integro-Differential Equations of Fractional Order
Consider the fractional integro-differential Equation (1), which may be written as: Multiplying Equation (7) by a general Lagrange multiplier λ, yields to: Now, take I α x to both sides of Equation (8), which give: Then, the correction functional with respect to x will be read as follows: and the problem now is to evaluate λ.The problem of evaluating λ is difficult, since Equation ( 10) consists of functional derivatives and integrals, so to avoid this difficulty, approximate I α x and T a w for 0 < α 1 by the first integral and derivative.Hence: Now, taking the first variation of Equation ( 11) with respect to u n , give: Using integration by parts, = Solving this equation will give λ(w, y) = −1, and substituting in Equation ( 10), we get: , which is a Banach space with a m-th-order con- tinuous partial derivative with respect to x, be the exact and approximate solutions of the integrodifferential equation of fractional order (1).If E n (x, y) = u n (x, y) − u(x, y) and the kernel K satisfy Lipschitz with respect to u, with constant L satisfying L < θ(2θ + 1)(3θ + 1) . . .(nθ + 1)γ(2γ + 1)(3γ + 1) . . .
then, the sequences solutions of approximation {u n } converge to the exact solution u.
Proof.Consider the integro-differential equation of fractional order: T α x u(x, y) = g(x, y) + I β x I γ y K(x, y, s, t, u(s, t)), where u(0, y) = u 0 The approximate solution using the VIM is given by: Since u is the exact solution of the integro-differential of the fractional order, Then, subtracting Equation ( 21) from Equation ( 20), we get: x I α w E n (x, y) = E n (x, y) − E n (0, y) and since E n (0, y) = 0, then If θ = α + β, then: Now, taking the supermum norm for both sides of Equation ( 23) where the K (kernel function) satisfies the Lipschitz condition with constant L, then: Using the conformable definition of integrals in Equation (24) implies: Now, applying mathematical induction over the last inequality: If n = 1, then: Hence, by induction, we have: and upon taking the supremum values of x and y, over [a, b] × [c.d], getting: Hence as as n → ∞, we have E n (x, y) → 0 , i.e., u n (x, y) → u(x, y) as n → ∞.

Illustrative Examples
Two examples of using the VIM to solve linear and nonlinear integro-differential equations with conformable fractional order differo-integration are presented in this section.
(0.00000525264x 8.1 y 3.25 +4.0829e − 11x 5.8 y 2.5 − 7.91503e − 7x 8.9 y 3.25 +0.0379687x 4.3 y 1.75 + 0.00000787187x 6.1 y 5.25 − 0.0000113188x 7.1 y 4.25 + 5.946e − 12x 6.6 y 2.5 + 6.78554e − 11x 4.8 y 3.5 − 9.03963e − 7x 6.9 y 5.25 + 0.00000150668x 7.9 y 4.25 −3.86567e − 8x 8.4 y 6 + 9.63085e−  When the exact solution u e (x, y) = x 2 y 2 .For the simplicity of calculations and in order to use the properties of conformable and integrations e u by the Tayler series After tow terms, we have e u ∼ = 1 + u .By calculus of variation, the initial condition to Now, using the variation iteration method to find the next approximation solution as follows: Carrying out recursively order integrals of order γ to represent y, of order β to represent x and of order α to represent y, we get the final from of u 1 (x, y), which is as follows.

Figure 1 5 y = 1 Figure 1 .
Figure1shows the comparison between the exact and the approximated solution for different values of y.

Figure 2 5 Figure 2 .
Figure 2 compares the exact and the approximated solutions for different values of y.