Power Series Solution for Solving Nonlinear Burgers-Type Equations

and Applied Analysis 3 U(x, 0) as a polynomial series applying Taylor’s theorem. Then U (x, 0) = γ 2 + γ 2 (xγ −] + √]2 + 8η 8 − 1 3 (xγ −] + √]2 + 8η 8 ) 3 + 2 15 (xγ −] + √]2 + 8η 8 ) 5


Introduction
Power series solution (PSS) method is an old method that has been limited to solve linear differential equations, both ordinary differential equations (ODE) [1,2] and partial differential equations (PDE) [3,4].Linear PDE have traditionally been solved using the separation of variables method because it permits obtaining a coupled system of ODE easier to solve with the PSS method.Some examples of these are the Legendre polynomials and the spherical harmonics used in Laplace's equations in spherical coordinates or in Bessel's equations in cylindrical coordinates [3,4].It is known that in nonlinear PDE (NLPDE) this procedure is not possible.
In this work we compare the spectral method (SM) with the PSS method solving three versions of nonlinear timedependent Burgers-type equations [5] because we know that the SM is the more accurate numerical method.The SM with collocation points (SMCP) is a numerical technique applied to solve linear and nonlinear differential equations with high accurate approximations to the solution [6].This has been used to solve PDE using polynomial interpolation function with an orthogonal basis such as Fourier, Chebyshev, or Legendre functions [7].The SM has also been very successful to solve any kind of DE problems, including integro-differential problems [8], with Newman boundary values [9], and nonlinear PDE [10].
We use the symbolic computation package Matlab to obtain the algebraic operations for the truncated series approximation.This program helps to do easier the tedious algebraic operations.

Power Series Solution Method
We know that almost the totality of the NLPDE does not have a solution with an analytic expression, that is, a solution in a closed form of known functions.Our goal is to construct a solution using a power series, taking advantage of the capacity of power series to represent any function with algebraic series developing the idea to construct an approximate solution [11][12][13][14][15][16][17].It also has the possibility to approximate a solution, inclusive if an analytic form does not exist, in a similar way like the Taylor's series approximate the functions.The existence of the PSS does not guarantee per se that the represented function has an exact approximation in distant points relative to the central value.However, considering 2 Abstract and Applied Analysis that the PSS needs to satisfy the NLPDE, with initial values condition (IVC) or with boundary values conditions (BVC), therefore we can construct a well posed problem to obtain an accurate solution, constrained with all these limiting conditions [18].Furthermore, the polynomial of the PSS is a smoothed function and this can guarantee the existence of a solution [18].
The PSS method represents a general solution with a series of unknown coefficients.When the PSS polynomial is substituted in the PDE we obtain a recurrence relation for the expansion coefficients.These coefficients should be expressed in function of the coefficients result from IVC or BVC.In this way, we obtain a system of equations depending on these initial value based coefficients.In order to obtain and solve a consistent algebraic system of equations, we also need the same number of coefficients and equations [11].All these conditions, in the beginning, provide a guarantee that the PDE is a well posed problem; that is, existence, uniqueness, and smoothness of the solution are well defined [18].
Finally, the PSS method is a proposal to find a semianalytic solution as an asymptotic approximation (in space and time) of a finite series with minimal error in the expansion of terms of the series.From numerical analysis when a power series ∑      converges on an interval (−, ) to a function , the radius of convergence is .In our work, the radius of convergence is defined by each interval where our error was estimated as we see below.

Numerical Results
First, we consider the nonlinear time-dependent one-dimensional generalized Burgers-Huxley equation [5]: with the initial condition where , ], , and  are real parameters.With  = 1, this equation admits a travelling wave solution.Then (, ) reads with  =  − , where  and  represent the wave number and frequency of the travelling wave, respectively, working as unknown variables.Introducing (3) in (1) we obtain The ansatz for (4) will be a PSS as The respective derivatives and nonlinear terms in (4) result in Substituting the series of ( 6) in (4), we obtain the recurrence relation Solving with Matlab until degree  = 3 of PSS from (4), we obtain the following values for the coefficients: We will use the initial conditions to obtain the unknown coefficients (8).From the initial condition (2), we express Abstract and Applied Analysis 3 (, 0) as a polynomial series applying Taylor's theorem.Then Matching the coefficients of this polynomial with the coefficients (8) of our ansatz, we obtain the next values:  0 = /2,  1 = /2,  2 = 0,  3 = −/6,  4 = 0,  5 = /15,  6 = 0,  7 = −17/630, and so forth.
With  2 and  3 values matched to their respective coefficients in (8), we obtain an algebraic system of 2 equations with two variables.Solving this one, we obtain the value of the unknown variables  and : Then, the complete solution as PSS for the NLDE (1) reads As it usually does when an approximate solution with PSS is obtained, a test of accuracy of the approximation must be performed.In this way, we calculate the absolute difference between exact and approximated solution defined as (, ) = |(, ) −   (, )|, where  is the exact solution obtained from [5], and   is the calculated solution (11), at the point (, ), until the power degree  = 19, respectively.We compute the error, with the parameters values  = 0.001, ] = 1, and  = 1, within the intervals  = [−2000, 2000] and  = [0, 1000].This result is shown in Figure 1.This parameter set was selected because it is the same one used in [5] to do a comparison.The convergence of the power series, () = ∑ ∞ =0     , depends on  and also the coefficient   , and then it is possible to adjust these ones to solve the NLDE and to find a solution that approximates its behavior to any distance and time, at less in the interval where we calculate the solution.We also compare the approximated solution with the spatial part fixed in  = 0.9, with different power degree polynomial ( = 7, 11, 19), relative to the exact result.This comparison is shown in Figure 2. From this figure we note that the improvement is better when the power degree increases.In a similar way, we solved (1)    For our second example, let us consider the nonlinear time-dependent one-dimensional Burgers-type equation [5]: where  = − and  and  are unknown parameters, and replacing ( 13) in (12) then We use the following PSS ansatz for (14): then or moving conveniently the sum index  to , then Substituting the above series into (14), we have This recurrence relation is used in the PSS until the degree  = 5.Using Matlab to solve the recurrence algebraic system, we obtain Expanding the initial condition of ( 12) as a polynomial through the Taylor series, Matching the coefficients of (20) with the coefficients of solution (19), we arrive to , and so forth. (21) The  2 and  3 matched values are used to get a system of 2 equations with 2 variables,  and .Solving this one, we obtain the value of the unknown variables  and : and the complete solution in (15) reads  (12), with ] = 10,  = 0.1, and  = 0.1, for  = 0.9.As in the previous example, we calculate the error (, ) with the parameters values ] = 10,  = 0.1, and  = 0.1 within the intervals  = [−1, 2] and  = [0, 13] until the power degree  = 19.This result is shown in Figure 5.In Figure 6 we make the comparison for a fixed  = 0.9 in the time interval  = [0, 30] with several power degree polynomials ( = 7, 11, 19).
For the third example, we consider the nonlinear timedependent one-dimensional generalized Burgers-Fisher-type equation [5]: with the initial condition where ], , and  are real numbers.When  = 1, this equation admits a travelling wave solution with  =  − .This one is replaced in (24) and we attempt a PSS and so the terms of ( 27) can be written as Substituting ( 29) into (27), we obtain the recurrence relation This recurrence relation is solved with Matlab until the degree  = 4 and the coefficients result The initial condition (25) is expressed like a polynomial series applying Taylor's formula Matching the coefficients of this polynomial with those of polynomial from solution (31), we obtain the value of the coefficients:  0 = 1/2,  1 = −1/2,  2 = 0,  3 = 1/6,  4 = 0,  5 = −2/30,  6 = 0,  7 = 17/630, and so forth.With  2 and  3 values matched to their respective coefficients of (31), we obtain a system of 2 equations with 2 variables to solve and to obtain the values of the unknown parameters  and : Therefore from the solution (28) reads We computed the error with (34) and the exact solution presented in [5] with the parameters values  = 0.01 and Abstract and Applied Analysis ] = 0.01 within the intervals  = [−300, 300] and  = [0, 2200] until the power degree  = 19.This error is shown in Figure 7. Figure 8 shows the trajectories of the approximate solution of problem (24) where  = 1, for  = 0.9, and  = 0.01 and ] = 0.01, in the time interval  = [0, 3000], with the power degree polynomial ( = 7, 11, 19).In Figure 9 We remark that Figures 2, 6, and 8 are shown simple trajectories of the approximate solution due to the fact that they are the solution of this NLDE in the time interval with the axis  fixed (to compare with [5]) where the behavior is smoother.These are nonlinear equations where the error grows in an exponential way, at difference of the behavior of linear equations.We show that by adding more terms to the ansatz series, the error is reduced and the solution is extended to long times and an increased space in comparison with other methods such as the spectral method that is one of the most accurate numerical methods.The PSS method is employed to more complex trajectories as shown in [13], where the method is used to solve chaotic NLDE.

Discussion and Conclusions
In this work we have shown that it is possible to solve nonlinear differential equations with the power series solution method.This method is implemented as a general approximate solution for each nonlinear PDE or ODE, in a similar way to the solution of a Linear DE.We transform a nonlinear DE problem into an algebraic system of equations.
Therefore, with this method, it is possible to obtain a well posed problem when the system of DE is consistent.It means that the system has the same number of variables as equations, and consequently the number of algebraic variables (coefficients of expansion series) obtained with this method is the same as the number of equations.We had shown that the PSS method is a semianalytic technique that leads, in an easier and exact way, to the solution of difficult differential equations with an approximated closed form expression.This is especially useful to solve nonlinear equations and leaves open the possibility to describe in an exact and convergent way the behavior of chaotic dynamical systems [13,15].Then, the solution can be approximated until the necessary degree into the power series.The convergence of the PSS method depends on the number of terms used in the power series.Once we know it, we can determine the domain of the space and time where the solution is valid.Three examples of the nonlinear time-dependent one-dimensional Burgers-Huxley, Burgerstype, and Burgers-Fisher-type equations were solved within this capable method.The errors between the exact analytic solution from [5] and our PSS are compared in each case.We obtained a better approximation in larger space and time intervals (∼10 −3 -10 −7 ) with PSS of order  = 19 but not with the accuracy of [5] in a small space and time interval of the solution.
The advantages of the PSS method over the traditional methods are the same ones that are presented in the work with respect to the spectral method that solves any kind of NLDE.It means that it can solve NLDE with higher accuracy in an interval of time and space than the other traditional methods, adding more terms to the series.
In summary, we have shown that the PSS method is a technique which can be used to solve this kind of nonlinear differential equations [11].In this work, the main goal was to solve the differential partial equation nonlinear Burgers type with dimension 1 (spatial) + 1 (time), doing a transformation of this one to a single nonlinear equation.We compared the PPS solution with that coming from the spectral method.We believe that it is possible to extend the solution to 2 or 3 spatial dimensions plus the time dimension.This possibility will be considered in future works.The procedure provides a much greater accuracy adding more terms into the series.The PSS method opens the possibility to analyze other characteristics of the NLPDE, obtaining better semianalytic approximations that involve less computational efforts.