Solving Heat and Wave-Like Equations Using He ’ s Polynomials

We use He’s polynomials which are calculated form homotopy perturbation method HPM for solving heat and wave-like equations. The proposed iterative scheme finds the solution without any discretization, linearization, or restrictive assumptions. Several examples are given to verify the reliability and efficiency of the method. The fact that suggested technique solves nonlinear problems without using Adomian’s polynomials is a clear advantage of this algorithm over the decomposition method.


Introduction
The heat and wave-like models are the integral part of applied sciences and arise in various physical phenomena.Several techniques including spectral, characteristic, modified variational iteration and Adomian's decomposition have been used for solving these problems; see 1-3 and the references therein.Most of these techniques encounter a considerable size of difficulty.He 4-11 developed and formulated homotopy perturbation method HPM by merging the standard homotopy and perturbation.The homotopy perturbation method HPM proved to be compatible with the versatile nature of the physical problems and has been applied to a wide class of functional equations; see 1, 4-19 and the references therein.In this technique, the solution is given in an infinite series usually converging to an accurate solution; see 1, 4-19 and the references therein.It is worth mentioning that HPM is applied without any discretization, restrictive assumption or transformation and is free from round off errors.The HPM is applied for all the nonlinear terms in the problem without discretizing either by finite difference or by spline techniques at the nodes and involves laborious calculations coupled with a strong possibility of the illconditioned resultant equations which are a complicated problem to solve.Moreover, unlike the method of separation of variables that requires initial and boundary conditions, the homotopy perturbation method HPM provides an analytical solution by using the initial conditions only.The fact that HPM solves nonlinear problems without using Adomian's polynomials is a clear advantage of this technique over the decomposition method.It is worth mentioning that 12, 13 introduced He's polynomials by splitting the nonlinear term and proved their compatibility with Adomian's polynomials coupled with the conclusion that He's polynomials are easier to calculate, are more user friendly, and are independent of the complexities arising in calculating the so-called Adomian's polynomials.It is to be highlighted that He's polynomials are calculated from homotopy perturbation method HPM .Inspired and motivated by the ongoing research in this area, we use He's polynomials for solving heat and wave-like equations.It is worth mentioning that Noor and Mohyud-Din 20 introduced a homotopy approach which involves an additional term and consequently leads towards laborious and redundant calculations, whereas the approach used in Section 2 is more precise and easier to implement.Moreover, it reduces the huge unnecessary calculation arising in 20 .Several examples are given to verify the reliability and efficiency of the algorithm.

Homotopy Perturbation Method and He's Polynomials
To explain the homotopy perturbation method, we consider a general equation of the type where L is any integral or differential operator.We define a convex homotopy H u, p by where F u is a functional operator with known solutions v 0 , which can be obtained easily.It is clear that, for This shows that H u, p continuously traces an implicitly defined curve from a starting point H v 0 , 0 to a solution function H f, 1 .The embedding parameter monotonically increases from zero to unit as the trivial problem F u 0 is continuously deforms the original problem L u 0. The embedding parameter p ∈ 0, 1 can be considered as an expanding parameter 1, 2, 4-19 .The homotopy perturbation method uses the homotopy parameter p as an expanding parameter 4-11 to obtain if p → 1, then 2.5 corresponds to 2.2 and becomes the approximate solution of the form It is well known that series 2.6 is convergent for most of the cases and also the rate of convergence is dependent on L u ; see 4-11 .We assume that 2.6 has a unique solution.
The comparisons of like powers of p give solutions of various orders.In sum, according to 12, 13 , He's HPM considers the solution, u x , of the homotopy equation in a series of p as follows: and the method considers the nonlinear term N u as where H n 's are the so-called He's polynomials 12, 13 , which can be calculated by using the formula , n 0, 1, 2, . . . .2.9

Numerical Applications
In this section, we use He's polynomials which are calculated from homotopy perturbation method HPM for solving heat and wave-like equations.
Example 3.1 1, 2 .Consider the one-dimensional initial boundary value problem which describes the heat-like models and initial conditions u x, 0 x 2 .

3.3
Apply the convex homotopy Compare the coefficient of like powers of p

3.10
Apply the convex homotopy method

3.11
Compare the coefficient of like powers of p p 0 : u 0 x, y, t y

3.16
and the initial conditions u x, y, z, 0 0.

3.17
Apply the convex homotopy method

3.18
Compare the coefficient of like powers of p p 0 : u 0 x, y, z, t x 4 y 4 z 4 t,

3.23
Apply the convex homotopy method

3.25
where p i s are He's polynomials.The series solution is given by

3.30
Apply the convex homotopy method Compare the coefficient of like powers of p p 0 : u 0 x, t x x 2 t,

3.37
Apply the convex homotopy method

3.38
Compare the coefficient of like powers of p p 0 : u 0 x, y, t x 4 y 4 t, p

3.44
Apply the convex homotopy method

3.45
Compare the coefficient of like powers of p and proceede as before, the series solution is given by

Conclusion
In this paper, we use He's polynomials which are calculated from homotopy perturbation method HPM for solving heat and wave-like equations.The method is applied in a direct way without using linearization, transformation, discretization, or restrictive assumptions.
It may be concluded that the proposed scheme is very powerful and efficient in finding the analytical solutions for a wide class of boundary value problems.The method gives more realistic series solutions that converge very rapidly in physical problems.It is also observed that He's polynomials are compatible with Adomian's Polynomials but are easier to calculate and are more user friendly.