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 Laplace-Padé 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 after-treatment to the approximate solutions with the Laplace-Padé (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 Laplace-Padé 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
(2.1)A(u)-f(r)=0,wherer∈Ω
with the boundary condition given by
(2.2)B(u,∂u∂η)=0,wherer∈Γ,
where A is a general differential operator, f(r) is a known analytic function, B is a boundary operator, Γ is the boundary of the domain Ω, and ∂u/∂η denotes differentiation along the normal drawn outwards from Ω [49]. In most cases, the A operator can be split into two operators, namely L and N, which represent the linear and the nonlinear operators, respectively. Hence, (2.1) can be rewritten as
(2.3)L(u)+N(u)-f(r)=0.
Now, a homotopy formulation can be given as
(2.4)H(v,p)=(1-p)[L(v)-L(u0)]+p(L(v)+N(v)-f(r))=0,p∈[0,1],
where u0 is the trial function (initial approximation) for (2.3) that satisfies the boundary conditions, and p is known as the perturbation homotopy parameter.
The equation above exhibits specific behaviors at the limit values p=0 and p=1, as given in
(2.5)H(v,0)=L(v)-L(u0)=0,H(v,1)=L(v)+N(v)-f(r)=0.
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 p(2.6)v=p0v0+p1v1+p2v2+⋯.
In the limit, when p→1, the approximate solution for (2.1) is give as
(2.7)u=limp→1v=v0+v1+v2+⋯,
where v0,v1,v2,… 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
(2.8)H(v,p)=(1-Q)[L(v)-L(u0)]+Q(L(v)+N(v)-f(r))=0,p∈[0,1],q∈[0,1],(2.9)Q=ap+(1-a)q,0<a<1,
where p and q are the homotopy parameters, and a 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:
(2.10)v=p0v0+p1v1+p2v2+⋯q0w0+q1w1+q2w2+⋯,
where v0,v1,v2,… and w1,w2,… are unknown functions to be determined by the RBHPM method. In addition, w0 is an arbitrary trial function, which is chosen in order to improve the RBHPM convergence in the same way as the trial function u0 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 p considered in the formulation. On the other side, the order for the RBHPM method is given as [i,k], where i and k are the highest power of p and q employed in the numerator and denominator of (2.10), respectively.
The limit of (2.10), when p→1 and q→1, provides an approximate solution for (2.3) given as
(2.11)u=limp→1q→1v=v0+v1+v2+⋯w0+w1+w2+⋯.
The limit above exists in the case that both limits
(2.12)limp→1(∑i=0∞vi),limq→1(∑i=0∞wi),where∑i=0∞wi≠0
exist.
3. Convergence of RBHPM Method
In order to analyse the convergence of the RBHPM method, (2.8) is rewritten as
(3.1)L(v)=L(u0)+Q[f(r)-N(v)-L(u0)]=0.
After applying the inverse operator, L-1, on both sides of (3.1), we obtain
(3.2)v=u0+Q[L-1f(r)-L-1N(v)-u0].
By assuming that (see (2.10))
(3.3)v=∑i=0∞pivi∑i=0∞qiwi,
and after substituting (3.3) in the right-hand side of (3.2) in the following form:
(3.4)v=u0+Q[L-1f(r)-(L-1N)[∑i=0∞pivi∑i=0∞qiwi]-u0].
The exact solution of (2.3) is obtained in the limit when p→1 and q→1 in (3.4), which results in
(3.5)u=limp→1q→1(QL-1f(r)-Q(L-1N)[∑i=0∞pivi∑i=0∞qiwi]+u0-Qu0),=L-1f(r)-[∑i=0∞(L-1N)(viβ)],β=∑i=0∞wi,
where
(3.6)limp→1q→1Q=limp→1q→1(ap+(1-a)q)=1.
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 N.
Theorem 3.1 (Sufficient Condition for Convergence).
Suppose that X and Y are Banach spaces and N:X→Y is a contractive nonlinear mapping, that is,
(3.7)∀w,w*∈X;∥N(w)-N(w*)∥≤γ∥w-w*∥;0<γ<1.
Then, according to Banach Fixed Point Theorem, N has a unique fixed point u, that is, N(u)=u. Assume that the sequence generated by the RBHPM method can be written as
(3.8)Wn=N(Wn-1),Wn-1=∑i=0n-1(viβ),n=1,2,3…,
and suppose that W0=v0/β∈Br(u), where Br(u)={w*∈X|∥w*-u∥<r}, then, under these conditions:
Wn∈Br(u),
limn→∞Wn=u.
Proof.
(i) By inductive approach, for n=1 we have
(3.9)∥W1-u∥=∥N(W0)-N(u)∥≤γ∥w0-u∥.
Assuming that ∥Wn-1-u∥≤γn-1∥w0-u∥, as induction hypothesis, then
(3.10)∥Wn-u∥=∥N(Wn-1)-N(u)∥≤γ∥Wn-1-u∥≤γn∥w0-u∥.
Using (i), we have
(3.11)∥Wn-u∥≤γn∥w0-u∥≤γnr<r⇒Wn∈Br(u).
(ii) Because of ∥Wn-u∥≤γn∥w0-u∥ and limn→∞γn=0, limn→∞∥Wn-u∥=0, that is,
(3.12)limn→∞Wn=u.
4. Padé Approximants
A rational approximation to f(x) on [a,b] is the quotient of two polynomials PN(x) and QM(x) of degrees N and M, respectively. We use the notation RN,M(x) to denote this quotient. The RN,M(x) Padé approximations [51, 52] to a function f(x) are given by
(4.1)RN,M(x)=PN(x)QM(x),fora≤x≤b.
The method of Padé requires f(x) and its derivative to be continuous at x=0. The polynomials used in (4.1) are expressed as
(4.2)PN(x)=p0+p1x+p2x2+⋯+pN(x),QM(x)=q0+q1x+q2x2+⋯+qM(x).
The polynomials in (4.2) are constructed so that f(x) and RN,M(x) agree at x=0 as well as their derivatives up to N+M agree at x=0. A special case occurs for Q0(x)=1, wherein the approximation in (4.1) becomes the Maclaurin expansion for f(x). For a fixed value of N+M the error is smallest when PN(x) and QM(x) have the same degree or when PN(x) has degree one higher than QM(x).
Notice that the constant coefficient of QM is q0-1. This is permissible, because it can be noted that 0 and RN,M(x) are not changed when both PN(x) and QM(x) are divided by the same constant. Hence the rational function RN,M(x) has N+M+1 unknown coefficients. Assume that f(x) is analytic and has the Maclaurin expansion
(4.3)f(x)=a0+a1x+a2x2+⋯+akxk+⋯.
And from the difference f(x)QM(x)-PN(x)=Z(x)(4.4)⌈∑i=0∞aixi⌉⌈∑i=0Mqixi⌉-⌈∑i=0Npixi⌉=⌈∑i=N+M+1∞cixi⌉.
The lower index j=N+M+1 in the summation on the right side of (4.4) is chosen because the first N+M derivatives of f(X) and RN,M(x) should agree at x=0.
When the left side of (4.4) is multiplied out and the coefficients of the powers of xi are set equal to zero for k=0,1,2,…,N+M, the result is a system of N+M+1 linear equations:
(4.5)a0-p0=0,q1a0+a1-p1=0,q2a0+q1a1+a2-p2=0,q3a0+q2a1+q1a2+a2-p3=0,qMaN-M|qM-1aN-M-1|aN-pN=0,(4.6)qMaN-M+1+qM-1aN-M+2+⋯+q1aN+aN+2=0,qMaN-M+2+qM-1aN-M+3+⋯+q1aN+1+aN+3=0,⋮qMaN+qM-1aN+1+⋯+q1aN+M+1+aN+M=0.
Notice that the sum of the subscripts of the terms of each product is the same in each equation, and it increases consecutively from 0 to N+M. The M equations in (4.6) involve only the unknowns q1,q2,…,qM and must be firstly solved. Then the equations in (4.5) are used successively to find p1,p2,…,pN [51].
5. Laplace-Padé 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.
First, Laplace transformation is applied to power series.
Next, s is substituted by 1/x in the resulting equation.
After that, we convert the transformed series into a meromorphic function by forming its Padé approximant of order [N/M]. N and M 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.
Then, x is substituted by 1/s.
Finally, by using the inverse Laplace s 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 Laplace-Padé 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 Laplace-Padé rational biparameter homotopy perturbation method (LPRBHPM).
6. Case Study 1: Riccati Nonlinear Differential Equation
Consider the Riccati nonlinear differential equation [7, 47]
(6.1)y′(x)-y(x)2+1=0,y(0)=0,
having the exact solution given as
(6.2)y(x)=-tanh(x).
6.1. Solution Calculated by RBHPM
We establish the following homotopy equation:
(6.3)(1-(ap+(1-a)q))(L(v)-L(u0))+(ap+(1-a)q)(v′-v2+1)=0,
where we have chosen a=0.25. Besides, the linear operator L is given as
(6.4)L(v)=v′+v+1,
and the trial function is
(6.5)u0=-1+exp(-x).
We suppose that the solution of (6.3) is of order [2,2], which is expressed as follows
(6.6)v=v0+v1p+v2p2w0+w1q+w2q2,
where the trial w-function is chosen as w0=1.
After substituting (6.6) into (6.3), regrouping, and equating the terms having the following powers: p0q0, p1, p2, q1, and q2, it can be solved for v0, v1, v2, w1, and w2. In order to fulfil the initial conditions of (6.1), it follows that v0(0)=0, v1(0)=0, v2(0)=0, w1(0)=0, and w2(0)=0.
The results are recast in the following system of differential equations:
(6.7)p0q0:w0v0′+v0w0-v0w0′+w02=0,v0(0)=0,p1:w0v1′+v1w0-v1w0′-av0w0-av02=0,v1(0)=0,p2:w0v2′-2av0v1+v2w0-v2w0′-av1w0=0,v2(0)=0,q1:-v0w1′-v0w0-v02+2w0w1+v0′w1+av02+av0w0+v0w1=0,w1(0)=0,q2:-v0w2′+2w0w2+av0w1+w12+v0w2+v0′w2-v0w1=0,w2(0)=0.
By substituting (6.8) into (6.6), and calculating the limits when p→1 and q→1, we obtain
(6.9)y(x)=limp→1q→1v=v0+v1+v2w0+w1+w2.
6.2. Solution by HPM
We establish the following homotopy equation:
(6.10)(1-p)(L(v)-L(u0))+p(v′-v2+1)=0,
where the linear operator L and the trial function u0 are (6.4) and (6.5), respectively.
Substituting (2.6) into (6.10), regrouping, and equating the terms with identical powers of p, it can be solved for v0, v1, v2, and so on (in order to fulfil initial conditions from v(0)=y(0)=0, it follows that v0(0)=0, v1(0)=0, v2(0)=0 and so on).
The result is recast in the following system of differential equations:
(6.11)p0:v0′+v0+1=0,v0(0)=0,p1:v1′+v1-v02-v0=0,v1(0)=0,p2:v2′+v2-2v0v1-v1=0,v2(0)=0,p3:v3′+v3-v12-v2-2v0v2=0,v3(0)=0,p4:v4′+v4-2v1v2-v3-2v0v3=0,v4(0)=0,⋮.
By solving (6.11), we obtain
(6.12)v0=-1+exp(-x),v1=-(x+exp(-x)-1)exp(-x),v2=12((4x-2)exp(-x)+2exp(-2x)+x2-2x)exp(-x),v3=-16(-12xexp(-x)+x3-3x2-6exp(-2x)+12x2exp(-x)+18xexp(-2x)+6exp(-3x)3x2)exp(-x),v4=124(-4x3-48x2exp(-x)-72xexp(-2x)-24exp(-3x)+x4+96xexp(-3x)+108x2exp(-2x)+32x3exp(-x)+24exp(-4x))exp(-x),⋮.
By using (6.12) and (2.6) and calculating the limit when p→1, we obtain the second and four order approximations
(6.13)y(x)=limp→1(∑i=02vipi)=v0+v1+v2,(6.14)y(x)=limp→1(∑i=04vipi)=v0+v1+v2+v3+v4,
respectively.
7. Case Study 2: Van Der Pol Oscillator
Consider the Van der Pol Oscillator problem [8, 48] without external forcing
(7.1)u′′+u'+u+u2u'=0,u(0)=0,u′(0)=1,
7.1. Solution by the RBHPM Method
We establish the following homotopy equation:
(7.2)(1-(ap+(1-a)q))(L(v)-L(u0))+(ap+(1-a)q)(v′′+v′+v+v2v′)=0,
where the linear operator L is
(7.3)L(v)=v′′,
and the trial function is
(7.4)u0=t.
We assume that the solution for (7.2) is order [2,2], which is expressed as follows:
(7.5)v=v0+v1p+v2p2w0+w1q+w2q2,
where the trial w-function is chosen as w0=1.
Substituting (7.5) into (7.2), regrouping, and equating the terms having the following powers: p0q0, p1, p2, q1, and q2, it can be solved for v0, v1, v2, w1, and w2. In order to fulfil the initial conditions of (7.1), it follows that v0(0)=0, v0′(0)=1 and the rest are vi(0)=0, vi′(0)=0, wi(0)=0 and wi′(0)=0, (i=1,2,3,…). Then, considering that w0=1, we establish the following system:
(7.6)p0:v0′′=0,v0(0)=0,v0′(0)=1,p1:v1′′+av02v0′+av0+av0′=0,v1(0)=0,v1′(0)=0,p2:v2′′+av02v1′+av1′+2av0v1v0′+av1=0,v2(0)=0,v2′(0)=0,q1:-v0w1′′+3v0′′w1+v0′+v0-2v0′w1′-av02v0′-av0-av0′+v02v0′=0,w1(0)=0,w1′(0)=0,q2:-v0w2′′+v02v0′w1+av03w1′-v03w1′-av02v0′w1+av0w1′-3av0w1-4v0′w1′w1-2v0w1′′w1-3av0′w1-2v0′w2′+2v0(w1′)2+3v0w1+3v0′w1+3v0′′w12-v0w1′+3v0′′w2=0,w2(0)=0,w2′(0)=0.
By solving (7.6), we obtain
(7.7)v0=t,v1=-112at2(t2+2t+6),v2=12520a2t3(30t4+77t3+315t2+210t+420),w1=112(1-a)t3+16(1-a)t2-12at+k1,w2=-51008(-1+a)2t6-1360(-1+a)2t5+15040(-28a+98-70a2)t4+112(-1+a)(a-16+k1)t3+15040(420a2+(-840-840k1)a+840k1)t2-12k1(a-3)t+k2,
where k1 and k2 are integration constants.
Substituting (7.7) into (7.5), and calculating the limits when p→1 and q→1, we obtain
(7.8)u(t)=limp→1q→1v=v0+v1+v2w0+w1+w2.
We select the parameters as a=27/100, k1=22/53, and k2=-17/39 by using the numerical procedure reported in [35–37].
In order to guarantee the validity of the approximate solution (7.8) for large t, the quotient of series solutions is transformed by the Laplace-Padé after-treatment. 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 1/t is written in place of s in the equation. Then, the Padé approximant [3/3] is applied and 1/s is written in place of t. Finally, by using the inverse Laplace s transformation, we obtain the modified approximate solution for numerator
(7.9)un(t)=-0.0566603453624exp(-1.71426194764t)+0.0566603453622exp(-0.0402536974457t)cos(0.736244931229t)+1.22941439189exp(-0.0402536974457t)sin(0.736244931229t).
Now, we repeat the same process to the denominator by changing the order of the Padé approximant to [3/4], to obtain
(7.10)ud(t)=0.003937675067exp(-0.535832417112t)cos(2.71442186582t)+0.003048633725exp(-0.535832417112t)sin(2.71442186582t)+0.228229226307exp(-0.192147218554t)+0.747030002351exp(0.628210432031t).
Then, the approximate solution calculated by LPRBHPM is
(7.11)u(t)=un(t)ud(t).
7.2. Solution by HPM
We establish the following homotopy equation:
(7.12)(1-p)(L(v)-L(u0))+p(v′′+v'+v+v2v')=0,
where the linear operator L and the trial function u0 are (7.3) and (7.4), respectively.
Substituting (2.6) into (7.12), regrouping, and equating the terms with identical powers of p. In order to fulfil initial conditions of (7.1), it follows that v0(0)=0, v0′(0)=1, and the rest are vi(0)=0 and vi′(0)=0 (i=1,2,3,…). Then, we establish the following system:
(7.13)p0:v0′′=0,v0(0)=0,v0′(0)=1,p1:v1′′+v0+v02v0′+v0′=0,v1(0)=0,v1′(0)=0,p2:v2′′+v02v1′+v1+2v0v1v0′+v1′=0,v2(0)=0,v2′(0)=0,p3:v3′′+v02v2′+2v0v1v1′+v12v0′+v2′+v2+2v0v2v0′=0,v3(0)=0,v3′(0)=0,p4:v4′′+2v1v2v0′+v12v1′+2v0v3v0′+v3′+v02v3′+2v0v1v2′+2v0v2v1′+v3=0,v4(0)=0,v4′(0)=0.
By solving (7.13), we obtain
(7.14)v0=t,v1=-112t4-16t3-12t2,v2=184t7+11360t6+18t5+112t4+16t3,v3=-1910080t10-6710080t9-531680t8-31720t7-67720t6-140t5-124t4,v4=7372358720t13+18891330560t12+83471108800t11+14039907200t10+13921362880t9+11336t8+25504t7+1180t6+1120t5.
By substituting the solutions from (7.14) into (2.6) and calculating the limit when p→1, we can obtain the second order approximation
(7.15)u(t)=limp→1(∑i=02vipi)=t-12t2+184t7+11360t6+18t5.
In the same way, we obtain the four order approximation
(7.16)u(t)=limp→1(∑i=04vipi)=t-12t2+311680t7-41720t6+13120t5+12329907200t10+11509362880t9+1840t8-124t4+7372358720t13+18891330560t12+83471108800t11.
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 [2,2] 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).
Comparison between exact solution (6.2) for (6.1) to the results of approximations (6.9), (6.13), and (6.14).
Homotopy (6.3)
Homotopy (6.10)
Homotopy (6.10)
Order
[2,2]
2
4
x
Exact (6.2)
RBHPM (6.9)
HPM (6.13)
HPM (6.14)
0.00
0.000000000000
0.000000000000
0.000000000000
0.000000000000
0.50
-0.462117157260
-0.462720433177
-0.462402288846
-0.462119211588
1.50
-0.905148253645
-0.918326771288
-0.911212665075
-0.90475203875
2.00
-0.964027580076
-0.975747936207
-0.960889970046
-0.965901380964
2.50
-0.986614298151
-0.991622046961
-0.968972449805
-0.992586779171
3.50
-0.998177897611
-0.989972317008
-0.961440997373
-1.00443324553
4.00
-0.999329299739
-0.98793701472
-0.961349802243
-1.0000486249
4.50
-0.999753210848
-0.987368801973
-0.96441914623
-0.993370392728
5.00
-0.999909204263
-0.988053048524
-0.969315733164
-0.986865576296
5.50
-0.999966597156
-0.98949557545
-0.974818141376
-0.981957532724
6.00
-0.999987711651
-0.991249069742
-0.980108525233
-0.979148093068
6.50
-0.999995479351
-0.993004075085
-0.984752811149
-0.97832221949
7.00
-0.999998336944
-0.994584980578
-0.988591496328
-0.979053617564
A.A.E.
8.15 (−3)
1.84 (−2)
8.93 (−3)
Exact solution (6.2) (diagonal cross) for (6.1) and its approximate solutions (6.9) (solid line), (6.13) (solid circles), and (6.14) (solid diamonds).
On the other hand, Figure 2 and Table 2 show a comparison between the Fehlberg fourth-fifth order Runge-Kutta method with degree four interpolant (RKF45) [61, 62] solution (built-in 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 [2,2] (see (7.11)) yields the smallest A.A.E. of all solutions, followed by (7.8) obtained by the RBHPM method of order [2,2], 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.
Comparison numerical solution (RKF45) for (7.1) and its first order approximate solutions given by RBHPM (7.8), LPRBHPM (7.11), HPM (7.15), and HPM (7.16).
Order
[2,2]
[2,2]
2
4
t
RKF45 for (7.1)
RBHPM (7.8)
LPRBHPM (7.11)
HPM (7.15)
HPM (7.16)
0.00
0.000000000000
0.000000000000
−2.04249011 (−13)
0.000000000000
0.000000000000
0.10
0.094996749389
0.0963305722762
0.0963305680033
0.095001281746
0.0949968616125
0.20
0.179964338949
0.181413957989
0.181413692649
0.180042107937
0.179964612584
0.30
0.254886117960
0.255755918036
0.255752978888
0.255328628571
0.254889069792
0.40
0.319824133934
0.319864038194
0.319847953401
0.321424660317
0.319850515572
0.50
0.374962243912
0.374261968751
0.374202151797
0.379476686508
0.3751195795
0.60
0.420620714410
0.419505461842
0.419331298484
0.431478857143
0.421336349604
0.70
0.457249830610
0.456198011306
0.455769940137
0.480583988889
0.459927483684
0.80
0.485408608033
0.485004250862
0.484075378004
0.531466565079
0.49405125938
0.90
0.505735525317
0.506659764874
0.504828319804
0.590743735714
0.53057647436
1.00
0.518917385708
0.521976574206
0.518630236646
0.66746031746
0.583931084348
1.10
0.525660073662
0.531844221386
0.526097880422
0.773643793651
0.683143422756
1.20
0.526664476070
0.537227026265
0.527855737463
0.924935314286
0.884097382101
1.30
0.522608206722
0.539158647966
0.524527343795
1.14130269603
1.2900013214
1.40
0.514133103782
0.538735520529
0.516726377800
1.44784142222
2.08441286596
1.50
0.501837942107
0.537111006964
0.505048313285
1.87566964286
3.58297069111
1.60
0.486275223293
0.535492254682
0.490063208382
2.4629231746
6.31238329453
1.70
0.467950980699
0.535141785803
0.472309972142
3.25585650079
11.128360908
1.80
0.447326548414
0.537385900716
0.452292231487
4.31005577143
19.3882241197
1.90
0.424821677896
0.54363212026
0.430475743238
5.69176980317
33.1990855168
2.00
0.400817960941
0.555398274173
0.407287171967
7.47936507937
55.769016169
2.10
0.375662771298
0.574356631824
0.383113985298
9.76491075
91.8967515218
2.20
0.349672525187
0.602397906631
0.358305196653
12.6558996317
148.645576503
2.30
0.323136129463
0.641722425466
0.333172699524
16.2771112079
236.259417439
2.40
0.296317411714
0.694969880802
0.307992974314
20.7726226286
369.394267712
2.50
0.269457558120
0.765406030734
0.28300899685
26.3079737103
568.756347393
2.60
0.242776634418
0.857196643743
0.258432227399
33.0724919365
863.260364544
2.70
0.216475028642
0.975820118585
0.234444603854
41.2817834571
1292.84749061
2.80
0.190734203918
1.12870912808
0.211200498849
51.1803960889
1912.13383476
2.90
0.165717437177
1.32628670185
0.188828626412
63.0446603151
2795.09702573
3.00
0.141570064350
1.58371504978
0.167433899647
77.1857142857
4041.05178634
A.A.E.
0.231
0.00641
12.1
401
RKF45 solution for Van Der Pol Oscillator problem (7.1) (diagonal cross) and its approximate solutions obtained by RBHPM (7.8) (asterisk), LPRBHPM (7.11) (solid line), HPM (7.15) (solid circles), and HPM (7.16) (empty squares).
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. Laplace-Padé after-treatment [46] was successfully applied, and it yields better accuracy and it allows for larger ranges of the domain. Additionally, the trial functions w0 and u0 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 L and N 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 CB-2010-01 no. 157024. The authors would like to thank Roberto Castaneda-Sheissa, Rogelio-Alejandro Callejas-Molina, and Roberto Ruiz-Gomez for their contribution to this project.
BiazarJ.AminikhahH.Study of convergence of homotopy perturbation method for systems of partial differential equations20095811-122221223010.1016/j.camwa.2009.03.0302557352ZBL1189.65246BiazarJ.GhazviniH.Convergence of the homotopy perturbation method for partial differential equations20091052633264010.1016/j.nonrwa.2008.07.0022523226ZBL1173.35395AğırsevenD.ÖzişT.An analytical study for Fisher type equations by using homotopy perturbation method201060360260910.1016/j.camwa.2010.05.0062665660ZBL1201.65187ÖzişT.AğırsevenD.He's homotopy perturbation method for solving heat-like and wave-like equations with variable coefficients2008372385944595010.1016/j.physleta.2008.07.0602458895ZBL1223.35294YildirimA.ÖzişT.Solutions of singular IVPs of Lane-Emden type by homotopy perturbation method20073691-270762-s2.0-3454824906410.1016/j.physleta.2007.04.072ÖzişT.YildirimA.A note on He's homotopy perturbation method for van der Pol oscillator with very strong nonlinearity20073439899912-s2.0-3424715484210.1016/j.chaos.2006.04.013Vázquez-LealH.Rational homotopy perturbation method201220121410.1155/2012/490342490342BarariA.OmidvarM.GhotbiA. R.GanjiD. D.Application of homotopy perturbation method and variational iteration method to nonlinear oscillator differential equations200810421611712-s2.0-5294913262410.1007/s10440-008-9248-9HeJ. H.Comparison of homotopy perturbation method and homotopy analysis method2004156252753910.1016/j.amc.2003.08.0082087530ZBL1062.65074HeJ. H.An elementary introduction to the homotopy perturbation method200957341041210.1016/j.camwa.2008.06.0032488613ZBL1165.65374HeJ. H.Homotopy perturbation technique19991783-42572622-s2.0-0032672778HeJ. H.The homotopy perturbation method for nonlinear oscillators with discontinuities200415112872922-s2.0-124228758710.1016/S0096-3003(03)00341-2HeJ. H.Homotopy perturbation method: a new nonlinear analytical technique2003135173792-s2.0-003744057910.1016/S0096-3003(01)00312-5HeJ. H.Coupling method of a homotopy technique and a perturbation technique for non-linear problems200035137432-s2.0-003370238410.1016/S0020-7462(98)00085-7HeJ. H.Application of homotopy perturbation method to nonlinear wave equations20052636957002-s2.0-1884442601610.1016/j.chaos.2005.03.006HeJ. H.Homotopy perturbation method for solving boundary value problems20063501-287882-s2.0-3064446035710.1016/j.physleta.2005.10.005KhanY.Vazquez-LealH.WuQ.An efficient iterated method for mathematical biology modelNeural Computing and Applications. In press10.1007/s00521-012-0952-zBiazarJ.EslamiM.A new homotopy perturbation method for solving systems of partial differential equations20116212252342-s2.0-7995952136510.1016/j.camwa.2011.04.070BiazarJ.GhanbariB.The homotopy perturbation method for solving neutral functional-differential equations with proportional delays201224133372-s2.0-7795609933410.1016/j.jksus.2010.07.026GanjiD. D.TariH.JooybariM. B.Variational iteration method and homotopy perturbation method for nonlinear evolution equations2007547-8101810272-s2.0-3474881559110.1016/j.camwa.2006.12.070SheikholeslamiM.AshorynejadH. R.GanjiD. D.YildirimA.Homotopy perturbation method for three-dimensional
problem of condensation film on inclined rotating disk2012193437442GanjiD. D.SadighiA.Application of homotopy-perturbation and variational iteration methods to nonlinear heat transfer and porous media equations2007207124342-s2.0-3425063670110.1016/j.cam.2006.07.030GanjiD. D.The application of He's homotopy perturbation method to nonlinear equations arising in heat transfer20063554-53373412-s2.0-3364689348110.1016/j.physleta.2006.02.056RafeiM.GanjiD. D.DanialiH.Solution of the epidemic model by homotopy perturbation method20071872105610622-s2.0-3424726652810.1016/j.amc.2006.09.019GanjiD. D.RajabiA.Assessment of homotopy-perturbation and perturbation methods in heat radiation equations20063333914002-s2.0-3314447846210.1016/j.icheatmasstransfer.2005.11.001FarazN.KhanY.Analytical solution of electrically conducted rotating flow of a second grade fluid over a
shrinking surface2011234221226KhanY.WuQ.FarazN.YildirimA.The effects of variable viscosity and thermal conductivity on a thin film flow over a shrinking/stretching sheet20116111339133992-s2.0-7995634415110.1016/j.camwa.2011.04.053MadaniM.FathizadehM.KhanY.YildirimA.On the coupling of the homotopy perturbation method and Laplace transformation2011539-10193719452-s2.0-7995195161210.1016/j.mcm.2011.01.023AminikhahH.The combined laplace transform and new homotopy perturbation methods for stiff systems of odes201236836383644CompeanF. I.OlveraD.CampaF. J.Lopez de LacalleL. N.Elias-ZunigaA.RodriguezC. A.Characterization and
stability analysis of a multivariable milling tool by the enhanced multistage homotopy perturbation method2012572733El-SayedA. M. A.ElsaidA.El-KallaI. L.HammadD.A homotopy perturbation technique for solving partial differential
equations of fractional order in finite domains20122181783298340KhanY.Vazquez-LealH.FarazN.An auxiliary parameter method using adomian polynomials and laplace transformation for nonlinear differential equationsApplied Mathematical Modelling. In press10.1016/j.apm.2012.06.026KhanY.Vazquez-LealH.Hernandez-MartinezL.Removal of noise oscillation term appearing in the nonlinear equation solution20122012910.1155/2012/387365387365Mohyud-DinS. T.YildirimA.IncM.Coupling of homotopy perturbation and modified Lindstedt-Poincaré methods for traveling wave solutions of the nonlinear Klein-Gordon equation20122421871912-s2.0-7795623883510.1016/j.jksus.2010.08.002Vazquez-LealH.Castaneda-SheissaR.Filobello-NinoU.Sarmiento-ReyesA.Sanchez OreaJ.High accurate simple
approximation of normal distribution related integrals201220122210.1155/2012/124029124029Vázquez-LealH.Filobello-NiñoU.Castañeda-SheissaR.Hernández-MartínezL.Sarmiento-ReyesA.Modified HPMS inspired by homotopy continuation methods201220121910.1155/2012/309123309123Filobello-NinoU.Vazquez-LealH.Castaneda-SheissaR.YildirimA.Hernandez-MartinezL.Pereyra-DiazD.Perez-SesmaA.Hoyos-ReyesC.An approximate solution of blasius equation by using HPM method201255059BiazarJ.BadpeimaF.AzimiF.Application of the homotopy perturbation method to Zakharov-Kuznetsov equations20095811-12239123942-s2.0-7035078618310.1016/j.camwa.2009.03.102BiazarJ.GhazviniH.Homotopy perturbation method for solving hyperbolic partial differential equations20085624534582-s2.0-4434919263010.1016/j.camwa.2007.10.032BiazarJ.GhazviniH.Exact solutions for non-linear Schrödinger equations by He's homotopy perturbation method20073661-279842-s2.0-3424786217510.1016/j.physleta.2007.01.060BiazarJ.GhazviniH.EslamiM.He's homotopy perturbation method for systems of integro-differential equations2009393125312582-s2.0-6294918578810.1016/j.chaos.2007.06.001BiazarJ.EslamiM.AminikhahH.Application of homotopy perturbation method for systems of Volterra integral equations of the first kind2009425302030262-s2.0-6765122718210.1016/j.chaos.2009.04.016BiazarJ.GhazviniH.He's homotopy perturbation method for solving systems of Volterra integral equations of the second kind20093927707772-s2.0-6274909099010.1016/j.chaos.2007.01.108BiazarJ.GhazviniH.Numerical solution for special non-linear Fredholm integral equation by HPM200819526816872-s2.0-3744903030110.1016/j.amc.2007.05.015Filobello-NinoU.Vazquez-LealH.KhanY.Castaneda-SheissaR.YildirimA.Hernandez-MartinezL.Sanchez-OreaJ.Castaneda-SheissaR.Rabago BernalF.Hpm applied to solve nonlinear circuits: a study case2012685–8843314344JiaoY. C.YamamotoY.DangC.HaoY.An aftertreatment technique for improving the accuracy of Adomian's decomposition method2002436-77837982-s2.0-003649900310.1016/S0898-1221(01)00321-2El-TawilM. A.BahnasawiA. A.Abdel-NabyA.Solving Riccati differential equation using Adomian's decomposition method200415725035142-s2.0-444437471110.1016/j.amc.2003.08.049BehiryS. H.HashishH.El-KallaI. L.ElsaidA.A new algorithm for the decomposition solution of nonlinear differential equations20075444594662-s2.0-3425087343410.1016/j.camwa.2006.12.027WangY.-G.LinW.-H.LiuN.A homotopy perturbation-based method for large deflection of a cantilever beam under a terminal follower force201213197201HeJ. H.An elementary introduction to recently developed asymptotic methods and nanomechanics in textile engineering20082221348735782-s2.0-5094911226610.1142/S0217979208048668MerdanM.GökdoǧanA.YildirimA.On the numerical solution of the model for HIV infection of CD4+ T cells20116211181232-s2.0-7995950632110.1016/j.camwa.2011.04.058BakerG. A.1975London, UKAcademic ExpressSweilamN. H.KhaderM. M.Exact solutions of some coupled nonlinear partial differential equations using the homotopy perturbation method20095811-12213421412-s2.0-7035057474510.1016/j.camwa.2009.03.059MomaniS.ErjaeeG. H.AlnasrM. H.The modified homotopy perturbation method for solving strongly nonlinear oscillators20095811-12220922202-s2.0-7035057116510.1016/j.camwa.2009.03.082KhanY.FarazN.Application of modified Laplace decomposition method for solving boundary layer equation20112311151192-s2.0-7995185836910.1016/j.jksus.2010.06.018BahugunaD.UjlayanA.PandeyD. N.A comparative study of numerical methods for solving an integro-differential equation2009579148514932-s2.0-6304912329110.1016/j.camwa.2008.10.097MomaniS.ErtürkV. S.Solutions of non-linear oscillators by the modified differential transform method20085548338422-s2.0-3764902044810.1016/j.camwa.2007.05.009GökdoĝanA.MerdanM.YildirimA.The modified algorithm for the differential transform method to solution of Genesio systems201217145512-s2.0-7996102758210.1016/j.cnsns.2011.03.039TsaiP. Y.ChenC. K.An approximate analytic solution of the nonlinear Riccati differential equation201034710185018622-s2.0-7865003420210.1016/j.jfranklin.2010.10.005EbaidA. E.A reliable aftertreatment for improving the differential transformation method and its application to nonlinear oscillators with fractional nonlinearities20111615285362-s2.0-7795592259510.1016/j.cnsns.2010.03.012EnrightW. H.JacksonK. R.NorsettS. P.ThomsenP. G.Interpolants for runge-kutta formulas19861231932182-s2.0-002278412010.1145/7921.7923FehlbergE.Klassische runge-kutta-formeln vierter und niedrigerer ordnung mit schrittweiten-kontrolle und ihre anwendung
auf waermeleitungsprobleme197066171ReifK.WeinzierlK.ZellA.UnbehauenR.A homotopy approach for nonlinear control synthesis1998439131113182-s2.0-0032163761AslamN.SunolA. K.Homotopy continuation based prediction of azeotropy in binary and multicomponent mixtures through equations of state200469232023262-s2.0-294255124310.1039/b315579jKublerF.SchmeddersK.Computing equilibria in stochastic finance economies2000151-21451722-s2.0-0007041864MelvilleR. C.TrajkovicL.FangS. C.WatsonL. T.Artificial parameter homotopy methods for the DC operating point problem19931268618772-s2.0-002760855510.1109/43.229761YamamuraK.SekiguchiT.InoueY.A fixed-point homotopy method for solving modified nodal equations19994666546652-s2.0-0032661141GeogheganR.LagariasJ. C.MelvilleR. C.Threading homotopies and DC operating points of nonlinear circuits1999911591782-s2.0-0032245109LeeJ.ChiangH. D.Constructive homotopy methods for finding all or multiple dc operating points of nonlinear circuits and systems200148151662-s2.0-003508527510.1109/81.903186UshidaA.YamagamiY.NishioY.KinouchiI.InoueY.An efficient algorithm for finding multiple DC solutions based on the SPICE-oriented Newton homotopy method20022133373482-s2.0-003649448110.1109/43.986427Vázquez-LealH.Hernández-MartínezL.Sarmiento-ReyesA.Double-bounded homotopy for analysing nonlinear resistive circuitsProceedings of the IEEE International Symposium on Circuits and Systems (ISCAS '05)May 2005320332062-s2.0-6764911155310.1109/ISCAS.2005.1465309Van BarelM.IkramovK. D.ChesnokovA. A.A continuation method for solving symmetric Toeplitz systems20084812212621392-s2.0-5974909249510.1134/S0965542508120026Vázquez-LealH.Hernández-MartínezL.Sarmiento-ReyesA.Castañeda-SheissaR.Numerical continuation scheme for tracing the double bounded homotopy for analysing nonlinear circuitsProceedings of the International Conference on Communications, Circuits and SystemsMay 2005112211262-s2.0-29844434602TrajkovicL.MelvilleR. C.FangS. C.Passivity and no-gain properties establish global convergence of a homotopy method for DC operating pointsProceedings of the IEEE International Symposium on Circuits and SystemsMay 19909149172-s2.0-0025665877GreenM. M.Efficient continuation method for use in globally convergent DC circuit simulationProceedings of the 3rd URSI International Symposium on Signals, Systems and Electronics (ISSSE '95)October 19954975002-s2.0-0029484722GoldgeisserL. B.GreenM. M.A method for automatically finding multiple operating points in nonlinear circuits20055247767842-s2.0-1814440203810.1109/TCSI.2005.844359WolfD. M.SandersS. R.Multiparameter homotopy methods for finding dc operating points of nonlinear circuits1996431082483810.1109/81.5389891416903ChuaL. O.UshidaA.A switching-parameter algorithm for finding multiple solutions of nonlinear resistive circuits1976432152392-s2.0-0016973951ChaoK. S.LiuD. K.PanC. T.A systematic search method for obtaining multiple solutions of simultaneous nonlinear
equations19752297487522-s2.0-0016548852BraninF. H.Jr.Widely convergent method for finding multiple solutions of simultaneous nonlinear equations1973159277927802-s2.0-0015585399VandenbergheL.VandewalleJ.Variable dimension algorithms for solving resistive circuits19901854434742-s2.0-0025488541KojimaM.YamamotoY.Variable dimension algorithms: basic theory, interpretations and extensions of some existing methods198224217721510.1007/BF01585103674630ZBL0509.90070YamamotoY.A variable dimension fixed point algorithm and the orientation of simplices198430330131210.1007/BF02591935769234ZBL0575.65043Vázquez-LealH.Hernández-MartínezL.Sarmiento-ReyesA.Murphy-ArteagaR. S.Improving multi-parameter homotopy via symbolic analysis techniques for circuit simulation2Proceedings of the European Conference on Circuit Theory and Design2003402405RoychowdhuryJ. S.MelvilleR. C.Homotopy techniques for obtaining a DC solution of large-scale MOS circuitsProceedings of the 33rd Annual Design Automation ConferenceJune 19962862912-s2.0-0029710304RoychowdhuryJ.MelvilleR.Delivering global DC convergence for large mixed-signal circuits via homotopy/continuation methods200625166782-s2.0-3134445191610.1109/TCAD.2005.852461TrajkovicL.MelvilleR. C.FangS. C.Finding DC operating points of transistor circuits using homotopy methodsProceedings of the IEEE International Symposium on Circuits and SystemsJune 19917587612-s2.0-0026308616TrajkovicL.WillsonA. N.Jr.Theory of dc operating points of transistor networks19924642282412-s2.0-0026898498Vazquez-LealH.Castaneda-SheissaR.YıldırımA.KhanY.Sarmiento-ReyesA.Jimenez-FernandezV.Herrera-MayA. L.Filobello-NinoU.Rabago-BernalF.Hoyos-ReyesC.Biparameter homotopybased direct current simulation of multistable circuits201223137150InoueY.KusanobuS.YamamuraK.TakahashiT.Newton-fixed-point homotopy method for finding dc operating points of nonlinear circuits1Proceedings of the International Technical Conference on Circuits/Systems, Computer and Communications (ITC-CSCC '01)July 2001Tokushima, Japan370373TrajkovićL.MathisW.Parameter embedding methods for finding dc operating points: formulation and implementationProceedings of the NOLTADicember 199511591164TrajkovićL.MelvilleR. C.FangS. C.Passivity and no-gain properties establish global convergence of a homotopy method for DC operating pointsProceedings of the IEEE International Symposium on Circuits and SystemsMay 19909149172-s2.0-0025665877Vazquez-LealH.Hernandez-MartinezL.Sarmiento-ReyesA.Castaneda-SheissaR.Gallardo-Del-AngelA.Homotopy method with a formal stop criterion applied to circuit simulation201182118081815