Research Article  Open Access
Hector VazquezLeal, Arturo SarmientoReyes, Yasir Khan, Uriel FilobelloNino, Alejandro DiazSanchez, "Rational Biparameter Homotopy Perturbation Method and LaplacePadé Coupled Version", Journal of Applied Mathematics, vol. 2012, Article ID 923975, 21 pages, 2012. https://doi.org/10.1155/2012/923975
Rational Biparameter Homotopy Perturbation Method and LaplacePadé Coupled Version
Abstract
The fact that most of the physical phenomena are modelled by nonlinear differential equations underlines the importance of having reliable methods for solving them. This work presents the rational biparameter homotopy perturbation method (RBHPM) as a novel tool with the potential to find approximate solutions for nonlinear differential equations. The method generates the solutions in the form of a quotient of two power series of different homotopy parameters. Besides, in order to improve accuracy, we propose the LaplacePadé rational biparameter homotopy perturbation method (LPRBHPM), when the solution is expressed as the quotient of two truncated power series. The usage of the method is illustrated with two case studies. On one side, a Ricatti nonlinear differential equation is solved and a comparison with the homotopy perturbation method (HPM) is presented. On the other side, a nonforced Van der Pol Oscillator is analysed and we compare results obtained with RBHPM, LPRBHPM, and HPM in order to conclude that the LPRBHPM and RBHPM methods generate the most accurate approximated solutions.
1. Introduction
Solving nonlinear differential equations is an important issue in sciences because many physical phenomena are modelled using such classes of equations. One of the most powerful methods to approximately solve nonlinear differential equations is the homotopy perturbation method (HPM) [1–45]. The HPM method is based on the use of a power series, which transforms the original nonlinear differential equation into a series of linear differential equations. In this work, we propose a generalization of the aforementioned concept by using a quotient of two power series of different homotopy parameters, which will be denominated as the rational biparameter homotopy perturbation method (RBHPM). In the same fashion, like HPM, the use of this quotient transforms the nonlinear differential equation into a series of linear differential equations. The generated solutions are expressed as the quotient of two truncated power series and they constitute the approximate solutions. Besides, we propose an aftertreatment to the approximate solutions with the LaplacePadé (LP) transform [46] in order to improve the accuracy of the solutions. This coupled method will be denominated as the LPRBHPM. In addition, the method is applied to two case studies, a Ricatti nonlinear differential equation [47] and a Van der Pol Oscillator [8, 48], without external forcing.
This paper is organized as follows. In Section 2, the basic idea of the RBHPM method is given. Section 3 presents a convergence analysis for the proposed method. In Section 4, the basic concept of Padé approximants is explained. In Section 5, the coupling of the RBHPM method with the LaplacePadé transform is recast. In Section 6, the first case study, a Riccati nonlinear differential equation, is solved by using the proposed method and HPM. In Section 7, the second case study, a nonlinear oscillator problem, is treated. Section 8 is devoted to discuss the resulting solutions and some numeric analysis issues. Finally, a brief conclusion is given in Section 9.
2. Basic Concept of RBHPM
The RBHPM and HPM methods share common foundations. For both methods, it can be considered that a nonlinear differential equation can be expressed as with the boundary condition given by where is a general differential operator, is a known analytic function, is a boundary operator, is the boundary of the domain , and denotes differentiation along the normal drawn outwards from [49]. In most cases, the operator can be split into two operators, namely and , which represent the linear and the nonlinear operators, respectively. Hence, (2.1) can be rewritten as
Now, a homotopy formulation can be given as where is the trial function (initial approximation) for (2.3) that satisfies the boundary conditions, and is known as the perturbation homotopy parameter.
The equation above exhibits specific behaviors at the limit values and , as given in
For the HPM method [11–14], we assume, without loss of generality, that the solution for (2.4) can be expressed as a power series of
In the limit, when , the approximate solution for (2.1) is give as where are unknown functions to be determined by the HPM method. The series in (2.7) is convergent in most cases [1, 2, 11, 14, 36].
For the RBHPM method, the homotopy in (2.4) can be rewritten as where and are the homotopy parameters, and is a weight factor that affects them in complementary proportions.
Now, we assume that the solution for (2.8) can be written as the quotient of the power series of both homotopy parameters: where and are unknown functions to be determined by the RBHPM method. In addition, is an arbitrary trial function, which is chosen in order to improve the RBHPM convergence in the same way as the trial function does for the HPM method [50].
On one side, the order of the approximation for the HPM method is determined by the highest power of considered in the formulation. On the other side, the order for the RBHPM method is given as , where and are the highest power of and employed in the numerator and denominator of (2.10), respectively.
The limit of (2.10), when and , provides an approximate solution for (2.3) given as
The limit above exists in the case that both limits exist.
3. Convergence of RBHPM Method
In order to analyse the convergence of the RBHPM method, (2.8) is rewritten as
After applying the inverse operator, , on both sides of (3.1), we obtain
By assuming that (see (2.10)) and after substituting (3.3) in the righthand side of (3.2) in the following form:
The exact solution of (2.3) is obtained in the limit when and in (3.4), which results in where
In order to study the convergence of the RBHPM method, we use the Banach Theorem as reported in [1, 2, 7, 36]. This theorem relates the solution of (2.3) to the fixed point problem of the nonlinear operator .
Theorem 3.1 (Sufficient Condition for Convergence). Suppose that and are Banach spaces and is a contractive nonlinear mapping, that is,
Then, according to Banach Fixed Point Theorem, has a unique fixed point , that is, . Assume that the sequence generated by the RBHPM method can be written as and suppose that , where , then, under these conditions:(i),(ii).
Proof. (i) By inductive approach, for we have Assuming that , as induction hypothesis, then Using (i), we have (ii) Because of and , , that is,
4. Padé Approximants
A rational approximation to on is the quotient of two polynomials and of degrees and , respectively. We use the notation to denote this quotient. The Padé approximations [51, 52] to a function are given by
The method of Padé requires and its derivative to be continuous at . The polynomials used in (4.1) are expressed as
The polynomials in (4.2) are constructed so that and agree at as well as their derivatives up to agree at . A special case occurs for , wherein the approximation in (4.1) becomes the Maclaurin expansion for . For a fixed value of the error is smallest when and have the same degree or when has degree one higher than .
Notice that the constant coefficient of is . This is permissible, because it can be noted that and are not changed when both and are divided by the same constant. Hence the rational function has unknown coefficients. Assume that is analytic and has the Maclaurin expansion
And from the difference
The lower index in the summation on the right side of (4.4) is chosen because the first derivatives of and should agree at .
When the left side of (4.4) is multiplied out and the coefficients of the powers of are set equal to zero for , the result is a system of linear equations:
Notice that the sum of the subscripts of the terms of each product is the same in each equation, and it increases consecutively from to . The equations in (4.6) involve only the unknowns and must be firstly solved. Then the equations in (4.5) are used successively to find [51].
5. LaplacePadé Transform and RBHPM Method Coupling
The coupling of Laplace transform and Padé approximant [46] is used in order to recover part of the lost information due to the truncated power series [51, 53–60]. The process can be recast as follows.(1)First, Laplace transformation is applied to power series. (2)Next, is substituted by in the resulting equation. (3)After that, we convert the transformed series into a meromorphic function by forming its Padé approximant of order . and are arbitrarily chosen, but they should be of smaller value than the order of the power series. In this step, the Padé approximant extends the domain of the truncated series solution to obtain better accuracy and convergence. (4)Then, is substituted by . (5)Finally, by using the inverse Laplace transformation, we obtain the modified approximate solution.
In order to improve the approximate rational solutions generated by the RBHPM method, we propose to apply the LaplacePadé method, separately, to the denominator and the numerator in (2.11), only when they are expressed as power series. We will denominate to this process as the LaplacePadé rational biparameter homotopy perturbation method (LPRBHPM).
6. Case Study 1: Riccati Nonlinear Differential Equation
Consider the Riccati nonlinear differential equation [7, 47] having the exact solution given as
6.1. Solution Calculated by RBHPM
We establish the following homotopy equation: where we have chosen . Besides, the linear operator is given as and the trial function is
We suppose that the solution of (6.3) is of order , which is expressed as follows where the trial function is chosen as .
After substituting (6.6) into (6.3), regrouping, and equating the terms having the following powers: , , , , and , it can be solved for , , , , and . In order to fulfil the initial conditions of (6.1), it follows that , , , , and .
The results are recast in the following system of differential equations:
Solving (6.7) yields
By substituting (6.8) into (6.6), and calculating the limits when and , we obtain
6.2. Solution by HPM
We establish the following homotopy equation: where the linear operator and the trial function are (6.4) and (6.5), respectively.
Substituting (2.6) into (6.10), regrouping, and equating the terms with identical powers of , it can be solved for , , , and so on (in order to fulfil initial conditions from , it follows that , , and so on).
The result is recast in the following system of differential equations:
By solving (6.11), we obtain
By using (6.12) and (2.6) and calculating the limit when , we obtain the second and four order approximations respectively.
7. Case Study 2: Van Der Pol Oscillator
Consider the Van der Pol Oscillator problem [8, 48] without external forcing
7.1. Solution by the RBHPM Method
We establish the following homotopy equation: where the linear operator is and the trial function is
We assume that the solution for (7.2) is order , which is expressed as follows: where the trial function is chosen as .
Substituting (7.5) into (7.2), regrouping, and equating the terms having the following powers: , , , , and , it can be solved for , , , , and . In order to fulfil the initial conditions of (7.1), it follows that , and the rest are , , and , (). Then, considering that , we establish the following system:
By solving (7.6), we obtain where and are integration constants.
Substituting (7.7) into (7.5), and calculating the limits when and , we obtain
We select the parameters as , , and by using the numerical procedure reported in [35–37].
In order to guarantee the validity of the approximate solution (7.8) for large , the quotient of series solutions is transformed by the LaplacePadé aftertreatment. The procedure is applied separately to numerator and denominator of the expression in (7.8). First, Laplace transformation is applied to numerator of (7.8) and then is written in place of in the equation. Then, the Padé approximant is applied and is written in place of . Finally, by using the inverse Laplace transformation, we obtain the modified approximate solution for numerator Now, we repeat the same process to the denominator by changing the order of the Padé approximant to , to obtain
Then, the approximate solution calculated by LPRBHPM is
7.2. Solution by HPM
We establish the following homotopy equation: where the linear operator and the trial function are (7.3) and (7.4), respectively.
Substituting (2.6) into (7.12), regrouping, and equating the terms with identical powers of . In order to fulfil initial conditions of (7.1), it follows that , , and the rest are and (). Then, we establish the following system:
By solving (7.13), we obtain
By substituting the solutions from (7.14) into (2.6) and calculating the limit when , we can obtain the second order approximation
In the same way, we obtain the four order approximation
8. Numerical Simulation and Discussion
On one hand, Figure 1 and Table 1 show a comparison between the exact solution (6.2) for the Riccati nonlinear differential equation (6.1) and the analytic approximations (6.9), (6.13), and (6.14). The RBHPM method of order gives the smallest average absolute error (A.A.E.) of all solutions, followed by the solutions obtained with HPM of order 4 and 2. It is remarkable to observe that the RBHPM method generates an accurate rational approximation that successfully replicates the asymptotic behaviour of (6.2).
On the other hand, Figure 2 and Table 2 show a comparison between the Fehlberg fourthfifth order RungeKutta method with degree four interpolant (RKF45) [61, 62] solution (builtin function of Maple software) for the Van der Pol Oscillator (7.1) and the analytic approximations (7.8), (7.11), (7.15), and (7.16). The LPRBHPM method of order (see (7.11)) yields the smallest A.A.E. of all solutions, followed by (7.8) obtained by the RBHPM method of order , and finally by approximations calculated by HPM of order 4 and 2. In this case study, both LPRBHPM and RBHPM generate more accurate solutions than HPM.

The coupling of Laplace and Padé with RBHPM was a key factor for improving the accuracy in the Van der Pol problem. Both numerator and denominator of the rational expression (7.8) were considered as truncated power series. LaplacePadé aftertreatment [46] was successfully applied, and it yields better accuracy and it allows for larger ranges of the domain. Additionally, the trial functions and play an important role in the behaviour of the RBHPM and LPRBHPM methods; therefore, future research must be done in order to understand what kinds of trial functions produce better results; in this context, the operators and will be a key aspect to consider. Furthermore, it is important to point out that the RBHPM and LPRBHPM methods do not resort to linearization, a perturbation parameter, or assumptions of weak nonlinearity, and it clearly results that the generated solution has a general character and it is more realistic compared to the method of simplifying the physical problems.
Finally, the homotopy in (2.4) can be replaced by homotopy schemes such as those reported in the literature [63–93]. This future line of research can lead us to improve the performance of RBHPM or LPRBHPM methods.
9. Conclusions
This paper presented the RBHPM and LPRBHPM methods as a novel tool with high potential to solve nonlinear differential equations. Furthermore, a comparison between the results of applying the proposed methods and HPM was given. For the Riccati nonlinear asymptotic problem, a comparison between the RBHPM and HPM methods was presented, showing how the RBHPM method generates highly accurate approximate solutions, similar to the results obtained when using the HPM method. Additionally, a Van der Pol Oscillator problem was tackled with the proposed methods and compared to solutions obtained by HPM. It resulted that the RBHPM and LPRBHPM methods produced the most accurate approximated solutions. Because RBHPM, LPRBHPM, and HPM are closely related, it is possible that differential equations solved by HPM can be solved by using RBHPM or LPRBHPM. Besides, further research should be done to apply the proposed methods to the calculation of approximate solutions of nonlinear partial differential equations, nonlinear fractional equations, and boundary value problems, among others. An important remark is that the RBHPM and LPRBHPM methods do not resort to any kind of linearization procedure or perturbation parameter. Thereupon, these methods promise to become important mathematical tools, useful for scientist and engineers working in the area of mathematical modelling and computer simulation.
Acknowledgments
The authors gratefully acknowledge the financial support of the National Council for Science and Technology of Mexico (CONACyT) through Grant CB201001 no. 157024. The authors would like to thank Roberto CastanedaSheissa, RogelioAlejandro CallejasMolina, and Roberto RuizGomez for their contribution to this project.
References
 J. Biazar and H. Aminikhah, “Study of convergence of homotopy perturbation method for systems of partial differential equations,” Computers & Mathematics with Applications, vol. 58, no. 1112, pp. 2221–2230, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Biazar and H. Ghazvini, “Convergence of the homotopy perturbation method for partial differential equations,” Nonlinear Analysis. Real World Applications, vol. 10, no. 5, pp. 2633–2640, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. Ağırseven and T. Öziş, “An analytical study for Fisher type equations by using homotopy perturbation method,” Computers & Mathematics with Applications, vol. 60, no. 3, pp. 602–609, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 T. Öziş and D. Ağırseven, “He's homotopy perturbation method for solving heatlike and wavelike equations with variable coefficients,” Physics Letters A, vol. 372, no. 38, pp. 5944–5950, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 A. Yildirim and T. Öziş, “Solutions of singular IVPs of LaneEmden type by homotopy perturbation method,” Physics Letters A, vol. 369, no. 12, pp. 70–76, 2007. View at: Publisher Site  Google Scholar
 T. Öziş and A. Yildirim, “A note on He's homotopy perturbation method for van der Pol oscillator with very strong nonlinearity,” Chaos, Solitons and Fractals, vol. 34, no. 3, pp. 989–991, 2007. View at: Publisher Site  Google Scholar
 H. VázquezLeal, “Rational homotopy perturbation method,” Journal of Applied Mathematics, vol. 2012, Article ID 490342, 14 pages, 2012. View at: Publisher Site  Google Scholar
 A. Barari, M. Omidvar, A. R. Ghotbi, and D. D. Ganji, “Application of homotopy perturbation method and variational iteration method to nonlinear oscillator differential equations,” Acta Applicandae Mathematicae, vol. 104, no. 2, pp. 161–171, 2008. View at: Publisher Site  Google Scholar
 J. H. He, “Comparison of homotopy perturbation method and homotopy analysis method,” Applied Mathematics and Computation, vol. 156, no. 2, pp. 527–539, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. He, “An elementary introduction to the homotopy perturbation method,” Computers & Mathematics with Applications, vol. 57, no. 3, pp. 410–412, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 34, pp. 257–262, 1999. View at: Google Scholar
 J. H. He, “The homotopy perturbation method for nonlinear oscillators with discontinuities,” Applied Mathematics and Computation, vol. 151, no. 1, pp. 287–292, 2004. View at: Publisher Site  Google Scholar
 J. H. He, “Homotopy perturbation method: a new nonlinear analytical technique,” Applied Mathematics and Computation, vol. 135, no. 1, pp. 73–79, 2003. View at: Publisher Site  Google Scholar
 J. H. He, “Coupling method of a homotopy technique and a perturbation technique for nonlinear problems,” International Journal of NonLinear Mechanics, vol. 35, no. 1, pp. 37–43, 2000. View at: Publisher Site  Google Scholar
 J. H. He, “Application of homotopy perturbation method to nonlinear wave equations,” Chaos, Solitons and Fractals, vol. 26, no. 3, pp. 695–700, 2005. View at: Publisher Site  Google Scholar
 J. H. He, “Homotopy perturbation method for solving boundary value problems,” Physics Letters A, vol. 350, no. 12, pp. 87–88, 2006. View at: Publisher Site  Google Scholar
 Y. Khan, H. VazquezLeal, and Q. Wu, “An efficient iterated method for mathematical biology model,” Neural Computing and Applications. In press. View at: Publisher Site  Google Scholar
 J. Biazar and M. Eslami, “A new homotopy perturbation method for solving systems of partial differential equations,” Computers and Mathematics with Applications, vol. 62, no. 1, pp. 225–234, 2011. View at: Publisher Site  Google Scholar
 J. Biazar and B. Ghanbari, “The homotopy perturbation method for solving neutral functionaldifferential equations with proportional delays,” Journal of King Saud University, vol. 24, no. 1, pp. 33–37, 2012. View at: Publisher Site  Google Scholar
 D. D. Ganji, H. Tari, and M. B. Jooybari, “Variational iteration method and homotopy perturbation method for nonlinear evolution equations,” Computers and Mathematics with Applications, vol. 54, no. 78, pp. 1018–1027, 2007. View at: Publisher Site  Google Scholar
 M. Sheikholeslami, H. R. Ashorynejad, D. D. Ganji, and A. Yildirim, “Homotopy perturbation method for threedimensional problem of condensation film on inclined rotating disk,” Scientia Iranica, vol. 19, no. 3, pp. 437–442, 2012. View at: Google Scholar
 D. D. Ganji and A. Sadighi, “Application of homotopyperturbation and variational iteration methods to nonlinear heat transfer and porous media equations,” Journal of Computational and Applied Mathematics, vol. 207, no. 1, pp. 24–34, 2007. View at: Publisher Site  Google Scholar
 D. D. Ganji, “The application of He's homotopy perturbation method to nonlinear equations arising in heat transfer,” Physics Letters A, vol. 355, no. 45, pp. 337–341, 2006. View at: Publisher Site  Google Scholar
 M. Rafei, D. D. Ganji, and H. Daniali, “Solution of the epidemic model by homotopy perturbation method,” Applied Mathematics and Computation, vol. 187, no. 2, pp. 1056–1062, 2007. View at: Publisher Site  Google Scholar
 D. D. Ganji and A. Rajabi, “Assessment of homotopyperturbation and perturbation methods in heat radiation equations,” International Communications in Heat and Mass Transfer, vol. 33, no. 3, pp. 391–400, 2006. View at: Publisher Site  Google Scholar
 N. Faraz and Y. Khan, “Analytical solution of electrically conducted rotating flow of a second grade fluid over a shrinking surface,” Ain Shams Engineering Journal, vol. 2, no. 34, pp. 221–226, 2011. View at: Google Scholar
 Y. Khan, Q. Wu, N. Faraz, and A. Yildirim, “The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet,” Computers and Mathematics with Applications, vol. 61, no. 11, pp. 3391–3399, 2011. View at: Publisher Site  Google Scholar
 M. Madani, M. Fathizadeh, Y. Khan, and A. Yildirim, “On the coupling of the homotopy perturbation method and Laplace transformation,” Mathematical and Computer Modelling, vol. 53, no. 910, pp. 1937–1945, 2011. View at: Publisher Site  Google Scholar
 H. Aminikhah, “The combined laplace transform and new homotopy perturbation methods for stiff systems of odes,” Applied Mathematical Modelling, vol. 36, no. 8, pp. 3638–3644, 2012. View at: Google Scholar
 F. I. Compean, D. Olvera, F. J. Campa, L. N. Lopez de Lacalle, A. EliasZuniga, and C. A. Rodriguez, “Characterization and stability analysis of a multivariable milling tool by the enhanced multistage homotopy perturbation method,” International Journal of Machine Tools and Manufacture, vol. 57, pp. 27–33, 2012. View at: Google Scholar
 A. M. A. ElSayed, A. Elsaid, I. L. ElKalla, and D. Hammad, “A homotopy perturbation technique for solving partial differential equations of fractional order in finite domains,” Applied Mathematics and Computation, vol. 218, no. 17, pp. 8329–8340, 2012. View at: Google Scholar
 Y. Khan, H. VazquezLeal, and N. Faraz, “An auxiliary parameter method using adomian polynomials and laplace transformation for nonlinear differential equations,” Applied Mathematical Modelling. In press. View at: Publisher Site  Google Scholar
 Y. Khan, H. VazquezLeal, and L. HernandezMartinez, “Removal of noise oscillation term appearing in the nonlinear equation solution,” Journal of Applied Mathematics, vol. 2012, Article ID 387365, 9 pages, 2012. View at: Publisher Site  Google Scholar
 S. T. MohyudDin, A. Yildirim, and M. Inc, “Coupling of homotopy perturbation and modified LindstedtPoincaré methods for traveling wave solutions of the nonlinear KleinGordon equation,” Journal of King Saud University, vol. 24, no. 2, pp. 187–191, 2012. View at: Publisher Site  Google Scholar
 H. VazquezLeal, R. CastanedaSheissa, U. FilobelloNino, A. SarmientoReyes, and J. Sanchez Orea, “High accurate simple approximation of normal distribution related integrals,” Mathematical Problems in Engineering, vol. 2012, Article ID 124029, 22 pages, 2012. View at: Publisher Site  Google Scholar
 H. VázquezLeal, U. FilobelloNiño, R. CastañedaSheissa, L. HernándezMartínez, and A. SarmientoReyes, “Modified HPMS inspired by homotopy continuation methods,” Mathematical Problems in Engineering, vol. 2012, Article ID 309123, 19 pages, 2012. View at: Publisher Site  Google Scholar
 U. FilobelloNino, H. VazquezLeal, R. CastanedaSheissa et al., “An approximate solution of blasius equation by using HPM method,” Asian Journal of Mathematics and Statistics, vol. 5, pp. 50–59, 2012. View at: Google Scholar
 J. Biazar, F. Badpeima, and F. Azimi, “Application of the homotopy perturbation method to ZakharovKuznetsov equations,” Computers and Mathematics with Applications, vol. 58, no. 1112, pp. 2391–2394, 2009. View at: Publisher Site  Google Scholar
 J. Biazar and H. Ghazvini, “Homotopy perturbation method for solving hyperbolic partial differential equations,” Computers and Mathematics with Applications, vol. 56, no. 2, pp. 453–458, 2008. View at: Publisher Site  Google Scholar
 J. Biazar and H. Ghazvini, “Exact solutions for nonlinear Schrödinger equations by He's homotopy perturbation method,” Physics Letters A, vol. 366, no. 12, pp. 79–84, 2007. View at: Publisher Site  Google Scholar
 J. Biazar, H. Ghazvini, and M. Eslami, “He's homotopy perturbation method for systems of integrodifferential equations,” Chaos, Solitons and Fractals, vol. 39, no. 3, pp. 1253–1258, 2009. View at: Publisher Site  Google Scholar
 J. Biazar, M. Eslami, and H. Aminikhah, “Application of homotopy perturbation method for systems of Volterra integral equations of the first kind,” Chaos, Solitons and Fractals, vol. 42, no. 5, pp. 3020–3026, 2009. View at: Publisher Site  Google Scholar
 J. Biazar and H. Ghazvini, “He's homotopy perturbation method for solving systems of Volterra integral equations of the second kind,” Chaos, Solitons and Fractals, vol. 39, no. 2, pp. 770–777, 2009. View at: Publisher Site  Google Scholar
 J. Biazar and H. Ghazvini, “Numerical solution for special nonlinear Fredholm integral equation by HPM,” Applied Mathematics and Computation, vol. 195, no. 2, pp. 681–687, 2008. View at: Publisher Site  Google Scholar
 U. FilobelloNino, H. VazquezLeal, Y. Khan et al., “Hpm applied to solve nonlinear circuits: a study case,” Applie