#### Abstract

The solution of the fractional hyperbolic partial differential equation is obtained by means of the variational iteration method. Our numerical results are compared with those obtained by the modified Gauss elimination method. Our results reveal that the technique introduced here is very effective, convenient, and quite accurate to one-dimensional fractional hyperbolic partial differential equations. Application of variational iteration technique to this problem has shown the rapid convergence of the sequence constructed by this method to the exact solution.

#### 1. Introduction

It is known that various problems in fluid mechanics (dynamics, elasticity) and other areas of physics lead to fractional partial differential equations. Methods of solutions of problems for fractional differential equations have been studied extensively by many researchers (see [1–11]).

The variational iteration method (VIM), which was proposed by He (see, e.g., [12–21]), was successfully applied to autonomous ordinary and partial differential equations and other fields. He [15] was the first research who applied the VIM to fractional differential equations. Odibat and Momani [22] extended the application of this method to provide approximate solutions for initial value problems of nonlinear partial differential equations of fractional order. VIM [23–25] is relatively a new approach to provide an analytical approximation to linear and nonlinear problems which is particularly valuable tools for scientists and applied mathematicians. Yulita et al. [26] used the VIM to obtain analytical solutions of fractional heat- and wave-like equations with variable coefficients. In the Ashyralyev et al. [27], the mixed boundary value problem for the multidimensional fractional hyperbolic equation is considered. The first order of accuracy in and the second order of accuracy in space variables for the approximate solution of problem were presented. The stability estimates for the solution of this difference scheme and its first- and second-order difference derivatives were established. A procedure of modified Gauss elimination method [28] was used for solving this difference scheme in the case of one-dimensional fractional hyperbolic partial differential equations.

In this paper, we apply variational iteration method to fractional hyperbolic partial differential equations and then we compare the results with those obtained using modified Gauss elimination method [27].

#### 2. Definitions

*Definition 2.1. *A reel function , , is said to be in the space , , if there exists a real number , such that , where , and it is said to be in the space if , .

*Definition 2.2. *If and , then
where is called the Riemann-Liouville fractional integral operator of order .

*Definition 2.3. *For , we let
which is called the Riemann-Liouville fractional derivative operator of order .

#### 3. Variational teration Method

In the present paper, the mixed boundary value problem for the multidimensional fractional hyperbolic equation is considered. Here is the unit open cube in the -dimensional Euclidean space with boundary , , , , and are given smooth functions and .

The correction functional for (3.1) can be approximately expressed as follows: where is a general Lagrangian multiplier [29] and is considered as a restricted variation as a restricted variation [21], that is, , and is its initial approximation. Using the above correction functional stationary and noticing that , we obtain From the above relation for any , we get the Euler-Lagrange equation: with the following natural boundary conditions: Therefore, the Lagrange multiplier can be identified as follows: Substituting the identified Lagrange multiplier into (3.2), following variational iteration formula can be obtained: In this case, let an initial approximation be. Then approximate solution takes the form .

##### 3.1. The Difference Scheme

The discretization of problem (3.1) is carried out in two steps. In the first step, let us define the grid space We introduce the Banach space of the grid functions defined on equipped with the norm To the differential operator generated by problem (3.1), we assign the difference operator by the formula acting in the space of grid functions satisfying the conditions for all It is known that is a self-adjoint positive definite operator in . With the help of we arrive at the initial boundary value problem for an finite system of ordinary fractional differential equations.

In the second step, applying the first order of approximation formula for and using the first order of accuracy stable difference scheme for hyperbolic equations (see [30]), one can present the first order of accuracy difference scheme: for the approximate solution of problem (3.1). Here

##### 3.2. Example 1

For the numerical result, the mixed problem for solving the one-dimensional fractional hyperbolic partial differential equation is considered.

According to the formula (3.7), the iteration formula for (3.13) is given by Now we start with an initial approximation:

Using the above iteration formula (3.14), we can obtain the other components as and so on; in the same manner the rest of the components of the iteration formula (3.14) can be obtained using the Maple package.

##### 3.3. Example 2

We consider one-dimensional fractional hyperbolic partial differential equation as follows: The iteration formula for (3.17) is given by When an initial approximation is , we have the other components as and so on. For (3.18), the rest of the components of the iteration formula can be obtained using the Maple 10 package.

#### 4. Conclusions

Variational iteration method is a powerful and efficient technique in finding exact and approximate solutions for one-dimensional fractional hyperbolic partial differential equations. The solution procedure is very simple by means of variational theory, and only a few steps lead to highly accurate solutions which are valid for the whole solution domain. The results of applying variational iteration method are exactly the same as those obtained by modified Gauss elimination method [27]. All Computations are performed by Maple 10 package program.

Figure 1 shows the exact solution of (3.13). Figure 2 shows difference scheme solution of (3.13). Figure 3 shows approximate solution by VIM for (3.13). Figure 4 shows the exact solution of (3.17). Figure 5 shows approximate solution by VIM for (3.17).