INTRODUCTION
In this study, we consider the linear and nonlinear wave equation (Inc et al., 2004):
On the finite xinterval [0, π] with Dirichlet boundary condition:
and F (x,t) is a given function of x and t.
This equation describes the propagation of a wave (or disturbance) and it arises in a wide variety of physical problems. some of these problems include a vibrating string, vibrating membrane, longitudinal vibrations of an elastic rod or beam, shallow water waves, acoustic problems for the velocity potential for a fluid flow through which sound can be transmitted, transmission of electric signals along a cable, shock waves, chemical exchange processes in chromatography, sediment transport in rivers and waves in plasmas and both electric and magnetic fields in the absence of charge and dielectric (Debnath, 1997).
We will apply the Variational Iteration Method (VIM) (He, 1999a, 2000; Momani and Abuasad, 2006; Ganji et al., 2007; Sweilam and Khader, 2007; Bildik and Konuralp, 2006) and Homotopy Perturbation Method (HPM) (He, 1999b, 2006; Zhang and He, 2006; Ganji and Sadighi, 2007) for the three types of partial differential wave equations.
Variational Iteration Method (VIM) and HomotopyPerturbation Method (HPM) are the most effective and convenient ones for both linear and nonlinear equations. The VIM is to construct correction functional using general Lagrange multipliers identified optimally via the variational theory and the initial approximation can be freely chosen with unknown constants. The HPM deforms a difficult problem in to a simple problem which can be easily solved using ordinary methods.
Extensive studies have been done regarding this matter. Such as investigation of Helmholtz equation and fifthorder KdV equation or Nonlinear Coupled Systems of reactiondiffusion Equations using Homotopy perturbation method (Ganji and Sadighi, 2006).
BASIC IDEA OF HOMOTOPYPERTURBATION METHOD
Linear and Nonlinear phenomena are of fundamental importance in various fields of science and engineering. Most models of reallife problems are still very difficult to solve. Therefore, approximate analytical solutions such as Homotopy Perturbation Method (HPM) were introduced. This method is the most effective and convenient ones for both linear and nonlinear equations.
Perturbation method is based on assuming a small parameter. The majority of nonlinear problems, especially those having strong nonlinearity, have no small parameters at all and the approximate solutions obtained by the perturbation methods, in most cases, are valid only for small values of the small parameter.
Generally, the perturbation solutions are uniformly valid as long as a scientific system parameter is small. However, we cannot rely fully on the approximations, because there is no criterion on which the small parameter should exists. Thus, it is essential to check the validity of the approximations numerically and/or experimentally. To overcome these difficulties, HPM have been proposed recently.
To explain this method, let us consider the following function:
With the boundary conditions of:
Where, A,B, f (r) and Γ are a general differential operator, a boundary operator, a known analytical function and the boundary of the domain Ω, respectively.
Generally speaking the operator can be divided in to a linear part L and a nonlinear part N (u). Equation 3 can therefore, be written as:
By the homotopy technique, we construct a homotopy v (r, p): Ωx[0,1]→R Which satisfies:
Or
Where, p ε [0,1] is an embedding parameter, while u_{0} is an
initial approximation of Eq. 3, which satisfies the boundary
conditions. Obviously, from Eq. 6 and 7
we will have:
The changing process of p from zero to unity is just that of v (r, p) from u_{0} to u (r). In topology, this is called deformation, while L (v)L (u_{0}) and A (v) f (r) are called homotopy.
According to the HPM, we can first use the embedding parameter p as a small
parameter and assume that the solutions of Eq. 6 and
7 can be written as a power series in p:
Setting p = 1 yields in the approximate solution of Eq. 3 to:
The combination of the perturbation method and the homotopy method is called the HPM, which eliminates the drawbacks of the traditional perturbation methods while keeping all its advantage.
The series (11) is convergent for most cases. However, the convergent rate
depends on the nonlinear operator A (v). Moreover, He (1999b) made the following
suggestions:
• 
The second derivative of N (v) with respect to v must be
small because the parameter may be relatively large, i.e., p→1. 
• 
The norm of must be smaller than one so that the series converges. 
BASIC IDEA OF VARIATIONAL ITERATION METHOD
To clarify the basic ideas of VIM, we consider the following differential equation:
Where:
L 
= 
A linear operator 
N 
= 
A nonlinear operator 
g (t) 
= 
A homogeneous term 
According to VIM, we can write down a correction functional as follows:
Where, λ is a general lagrangian multiplier which can be identified optimally via the variational theory. The subscript n indicates the nth approximation and u_{n} is considered as a restricted variation, i.e., δ " é_{n} = 0.
Example 1: We consider a homogeneous linear wave equation:
With the initial and boundary conditions posed are:
Exact solution of this equation is:
APPLICATION OF HOMOTOPYPERTURBATION METHOD
To solve Eq. 14 by means of HPM, we consider the following process after separating the linear and nonlinear parts of the equation; we apply HomotopyPerturbation to Eq. 6.
A homotopyperturbation method can be constructed as follows:
Substituting Eq. 10 into 17 and rearranging
the resultant equation based on powers of pterms, one has:
With the following conditions:
With the effective initial approximation for v_{0} from the conditions
(21) and solutions of Eq. 1820 may be
written as follows:
In the same manner, the rest of components were obtained using the maple package.
According to the HPM, we can conclude that:
Therefore, substituting the values of v_{0} (x, t), v_{1} (x,
t) and v_{2} (x, t) from Eq. 2224
into 25 yields:
APPLICATION OF VARIATIONAL ITERATION METHOD
Its stationary conditions can be obtained as follows:
We obtain the lagrangian multiplier:
As a result, we obtain the following iteration formula:
Now we start with an arbitrary initial approximation that satisfies the initial condition:
Using the above variational formula (30), we have
Substituting Eq. 31 into 32 and after
simplifications, we have:
In the same way, we obtain u_{2} (x, t) as follows:
And so on. In the same way the rest of the components of the iteration formula
can be obtained.
As it can be seen, using these two approximate methods lead to the exact solution (Fig. 1) using more iteration.
Example 2: We consider the firstorder nonlinear wave equation. The equation of the form:
With the exact solution of:
APPLICATION OF HOMOTOPYPERTURBATION METHOD
A homotopyperturbation method can be constructed as follows:

Fig. 1: 
3D obtained results of u(x,t) by the HPM, VIM and the exact
solution 
Substituting Eq. 10 into 37 and rearranging
the resultant equation based on powers of pterms, one has:
With the following conditions:
With the effective initial approximation for from the conditions (42) and solutions
of Eq. 3841 may be written as follows:
In the same manner, the rest of components were obtained using the maple package.
According to the HPM, we can conclude that:
Therefore, substituting the values of v_{0} (x, t), v_{1} (x,
t), v_{2} (x, t) and v_{3} (x, t) from Eq. 4346
into 47 yields:
APPLICATION OF VARIATIONAL ITERATION METHOD
Its stationary conditions can be obtained as follows:
The lagrangian multiplier can therefore be identified as:
As a result, we obtain the following iteration formula:
Now we start with an arbitrary initial approximation that satisfies the initial condition:
Using the above variational formula (52), we have
Substituting Eq. 53 into 54 and after
simplification, we have:

Fig. 2: 
3D obtained results of u(x,t) by the HPM, VIM and the exact
solution 
In the same way, we obtain u_{2} (x, t), u_{3} (x, t) as follows:
And so on. In the same way the rest of the components of the iteration formula can be obtained.
As it can be seen, using these two approximate methods in solving this equation also lead to the exact solution (Fig. 2) using more iteration.
Example 3: We consider a secondorder nonlinear wave equation:
With the initial and boundary conditions posed are:
Exact solution of this equation is:
APPLICATION OF HOMOTOPYPERTURBATION METHOD
A homotopy can be constructed as follows:
Substituting Eq. 10 into 61 and rearranging
the resultant equation based on powers of pterms, one has:
With the following conditions:
With the effective initial approximation for v_{0} from the conditions
(65) and solutions of Eq. 6264 may be
written as follows:
In the same manner, the rest of components were obtained using the maple package.
According to the HPM, we can conclude that:
Therefore, substituting the values of v_{0} (x, t), v_{1} (x,
t) and v_{2} (x, t) from Eq. 6669
yields:
Application of variational iteration method
Its stationary conditions can be obtained as follows:
The lagrangian multiplier can therefore be identified as:
As a result, we obtain the following iteration formula:
Now we start with an arbitrary initial approximation that satisfies the initial condition:
Using the above variational formula (74), we have:
Substituting Eq. 7576 and after simplification,
we have:
In the same way, we obtain u_{2} (x, t), u_{3} (x, t) as follows:

Fig. 3: 
3D obtained results of u(x,t) by the HPM, VIM and the exact
solution 
And so on. In the same way the rest of the components of the iteration formula can be obtained.
As it can be seen, using these two approximate methods in solving this equation also lead to the exact solution (Fig. 3) using more iteration.
CONCLUSION
The homotopy perturbation method and variational iteration method have been successfully used to study three types of partial differential wave equations. These equations describe the propagation of a wave (disturbance) and it arises in a wide variety of physical problems. The results obtained here were compared with the exact solutions. It can be easily seen that applying these two methods to the mentioned equations were led to the exact solution in all examples, so the results revealed that the homotopy perturbation method and variational iteration method are powerful mathematical tools for solutions of differential equations in terms of accuracy and efficiency.