#### Abstract

The variational iteration algorithm combined with the exp-function method is suggested to solve the generalized Benjamin-Bona-Mahony equation (BBM) with variable coefficients. Periodic and soliton solutions are formally derived in a general form. Some particular cases are considered.

#### 1. Introduction

The BBM equation

which describes approximately the unidirectional propagation of long waves in certain nonlinear dispersive systems, has been proposed by Benjamin et al. in 1972  as a more satisfactory model than the KdV equation 

It is easy to see that (1.1) can be derived from the equal width EW-equation :

by means of the change of variable , that is, by replacing with . This last equation is considered as an equally valid and accurate model for the same wave phenomena simulated by (1.1) and (1.2). On the other hand, some researches analyzed the generalized KdV equation with variable coefficients

because this model has important applications in several fields of science .

Motivated by these facts, we will consider here the generalized EW-equation with variable coefficients

Using the solutions of (1.5) we obtain exact solutions to the generalized BBM equation

of order .

#### 2. Exact Solutions to Generalized BBM Equation

##### 2.1. The Variational Iteration Method

Consider the following nonlinear equation:

where and are linear and nonlinear operators, respectively, and is an inhomogeneous term. According to the variational iteration method (VIM) , a functional correction to (2.1) is given by

where is a general Lagrange's multiplier, which can be identified via the variational theory; the subscript denotes the th order approximation and is a restricted variation which means . In this method, we first determine the Lagrange multiplier that will be identified optimally via integration by parts. The successive approximation of the solution will be readily obtained upon using the determined Lagrangian multiplier and any selective function . One of the advantages of the VIM, is the free choice of the initial solution . If we consider a special form to with arbitrary parameters, using the relations

we can obtain a set of algebraic equations in the unknowns given by the parameters that appear in . Solving this system, we have exact solutions to (2.1). To solve (1.5), we construct the following functional equation

where

Taking in (2.4) variation with respect to the independent variable , and noticing that we have

This yields the stationary conditions

Therefore,

Substituting this value into (2.4) we obtain the formula

Using the wave transformation

setting

and performing one integration, (2.9) reduces to

where for sake of simplicity we set the constant of integration equal to zero. With the change of variable

equation (2.12) converts to

Observe that if is a solution to (2.14), then is also a solution to this equation.

##### 2.2. The Exp-Function Method

Recently, He and Wu  have introduced the Exp-function method to solve nonlinear differential equations. In particular, the Exp-function method is an effective method for solving nonlinear equations with high nonlinearity. The method has been used in a satisfactory way by other authors to solve a great variety of nonlinear wave equations . The Exp-function method is very simple and straightforward, and can be briefly revised as follows: Given the nonlinear partial differential equation

it is transformed to ordinary differential equation

by mean of wave transformation . Solutions to (2.16) can then be found using the expression

where , and are positive integers which are unknown to be determined later, and are unknown constants.

After balancing, we substitute (2.17) into (2.16) to obtain an algebraic systems in the variable . Solving the algebraic system we can obtain exact solutions to (2.16) and reversing, solutions to (2.15) in the original variables.

#### 3. Solutions to (2.14) by the Exp-Function Method

Using the Exp-function method, we suppose that solutions to (2.14) can be expressed in the form

We obtain following solutions to (2.14): Some special solutions are obtained if

This choice gives solutions

Solution (3.6) follows from (3.4) with the identifications and

Solution (3.7) follows from (3.4) with the identifications and

#### 4. Particular Cases

##### 4.1. Case 1: Solutions to (2.14) When 𝑝=2

Equation (2.14) takes the form

From (3.2) with :

From (3.3)–(3.7) with :

Other exact solutions are:

##### 4.2. Case 2: Solutions to (2.14) When 𝑝=4

Equation (2.14) takes the form

From (3.2) with :

From (3.3)–(3.7) with :

Other exact solutions are:

It is clear that using (2.13) we obtain solutions to (1.5). Finally, observe that if is a solution of (1.5), then the solutions to the generalized BBM equation (1.6) are obtained as follows:

#### 5. Conclusions

We have considered the generalized EW-equation with variable coefficients and the generalized BBM-equation with variable coefficients. We obtained analytic solutions by using the variational iteration method combined with the exp-function method. With the aid of Mathematica we have derived a lot of different types of solutions for these two models. Combined formal soliton-like solutions as well as kink solutions have been formally derived. The results obtained show that the technique used here can be considered as a powerful method to analyze other types of nonlinear wave equations.

According to , there are alternative iteration alorithms, which might be useful for future work. Furthermore, various modifications of the exp-function method have been appeared in open literature, for example, the double exp-function method [23, 24].

Other methods for solving nonlinear differential equations may be found in .

We think that the results presented in this paper are new in the literature.