Exact and Numerical Analysis of the Pantograph Delay Differential Equation via the Homotopy Perturbation Method

: The delay differential equations are of great importance in real-life phenomena. A special type of these equations is the Pantograph delay differential equation. Generally, solving a delay differential equation is a challenge, especially when the complexity of the delay terms increases. In this paper, the homotopy perturbation method is proposed to solve the Pantograph delay differential equation via two different canonical forms; thus, two types of closed-form solutions are determined. The ﬁrst gives the standard power series solution while the second introduces the exponential function solution. It is declared that the current solution agrees with the corresponding ones in the literature in special cases. In addition, the properties of the solution are provided. Furthermore, the results are numerically validated through performing several comparisons with the available exact solutions. Moreover, the calculated residuals tend to zero, even in a huge domain, which reﬂects the high accuracy of the current analysis. The obtained results reveal the effectiveness and efﬁciency of the current analysis which can be further extended to other types of delay equations.


Introduction
The pantograph is a particular device which collects the current in electric trains.The pantograph problem has been extensively studied in the literature.Andrews [1] calculated the behavior of an overhead catenary system for railway electrification.Abbott [2] introduced a numerical method for calculating the dynamic behavior of a trolley wire overhead contact system for electric railways.The Pantograph motion on a nearly uniform railway overhead line and system has been discussed in [3][4][5], respectively.The Pantograph's model has been investigated by Kato and McLeod [6], Iserles [7], Derfel and Iserles [8], Patade and Bhalekar [9], Fox et al. [10], and recently by Alenazy et al. [11].However, this problem still has undiscovered properties as will be declared in this paper.Although a considerable amount of real-world problems were modeled by means of ordinary differential equations (ODEs), it is observed that such ODEs cannot model the actual behavior of natural systems.This is because the ordinary derivative is a local operator which cannot be used to model the memory/hereditary properties in real-life phenomena.Hence, physical phenomena can be accurately modeled by incorporating nonlocal components such as delays.Obtaining an exact/approximate solution of an ODE is regularly easier in contrast with a delay differential equation (DDE).As simple examples, it is easy to find the exact solution of the ODE y (t) − y(t) = 0 by any of the standard methods while finding solutions for the very simple DDEs y (t) − y(−t) = 0, y (t) − y(2t) = 0, y (t) − y t 3 = 0 is not an easy task.Actually, there are no standard methods for solving DDEs and, accordingly, the DDEs are difficult to solve when compared with ODEs.Despite the excellency of the obtained results in the literature for the pantograph delay differential equation (PDDE) [6]: the authors believe that the PDDE (1) still needs an effort to determine its exact solution for all possible real values of parameters a, b, c, and λ.As a special case, Equation (1) reduces to the standard Ambartsumian delay differential equation (ADDE) when a = −1 and b = c = 1 q (q > 1) which has been used for describing the theory of surface brightness in the Milky Way [12].The ADDE has been investigated in the literature using different methodologies such as the Taylor series [13] and the exponential function solution [14][15][16].It was also generalized using different approaches by Khaled et al. [17] and Kumar et al. [18].Moreover, Ebaid and Al-Jeaid [19] obtained the exact periodic solution for Equation (1) for the special case c = −1 such that b > a.In addition, they determined the exact solution when a = b and c = −1.The purpose of this paper is to obtain the closed-form solution of the PDDE for arbitrary values of the proportional delay parameter c.
On the other hand, the homotopy perturbation method (HPM) is a relatively recent method, see [28][29][30][31][32][33][34].It will be shown that the HPM is an effective tool to deal with the current PDDE through utilizing two different canonical forms.The first provides the standard power series solution (PSS) while the second gives the exponential function solution (EFS).In addition, it will be declared that the current PSS agrees with the corresponding PSS in [10] when λ = 1 and also agrees with the EFS in [15,16] when a = −1, b = c = 1 q (q > 1) for the ADDE.Furthermore, the properties of the EFS, in terms of exponential functions, will be provided.Moreover, several comparisons between the EFS and the available exact ones in the literature are to be performed and, accordingly, the advantages of the present analysis over those in previously published papers will be demonstrated.
The paper has the following structure.In Section 2, the PSS is derived, where the existing results in the literature are obtained for particular cases.Section 3 focuses on evaluating the EFS while its characteristics are discussed in Section 4. Section 5 applies the quantum calculus notations on a general-component formula.Section 6 is devoted to analyze and discuss the results which are concluded in Section 7.

The first Canonical Form: PSS
In order to apply the HPM on the PDDE (1), we first rewrite Equation (1) in the following canonical form: where 0 < p ≤ 1 is an embedding parameter which is used to construct the homotopy series solution: On substituting Equation (3) into Equation (2), we have or i.e., Equation ( 6) leads to the following systems of initial value problems (IVPs): and The solution of Equation ( 7) is given by From Equation ( 8), the nth-order component can be obtained as Since Utilizing Equation (11), we obtain . . .
The HPM gives the solution as p → 1 by At λ = 1, the solution (17) becomes which is the same closed-form series solution obtained by Fox et al. [10].In addition, if a = −1 and b = c = 1 q q > 1, then the PDDE (1) becomes the ADDE: Substituting a = −1 and b = c = 1 q into Equation ( 17) gives i.e., which is the same solution obtained in [15,16].

The Second Canonical Form: EPS
In this section, we start with transforming Equation (1) into the following equivalent integral equation: The canonical form of the HPM, in this case, is of the form: or which gives the recurrence scheme: For n = 1, we obtain where y 0 (cτ) = λe acτ .Evaluating the last integral gives Similarly, at n = 2 and n = 3, we can get e ac 2 t − (c + 1)e act + ce at , a = 0, c = ±1, and From above, the calculation of the first few components of the homotopy series is indicated, however, the components of higher-order can be easily derived via any software.It will be shown later in a subsequent section that it is possible to derive a unified formula for the nth-component y n , ∀ n ≥ 1.To do so, at first, we have to list the main characteristics of the preceding components y 1 , y 2 , and y 3 .These characteristics are investigated in the next section and, hence, a unified formula for y n will be possible.

Characteristics of the Components and a Unified Formula for y n (t) 4.1. Characteristics of the Components
From the expressions for the components y 1 (t), y 2 (t), and y 3 (t), it can be noticed that • y 1 (t) involves two terms and can be written as where • y 2 (t) consists of three terms in the form: where • y 3 (t) contains four terms and we can write where • Similarly, y 4 (t) can be evaluated from (26) and is given by where • In view of ( 31), ( 34), (38), and (44), y n (t) takes the form y n (t) = ∑ n m=0 µ n,m e ac m t ∀ n ≥ 1.Moreover, the initial component y 0 (t) = ∑ 0 m=0 µ 0,m e ac m t = µ 0,0 e at where µ 0,0 = λ.Based on the above characteristics, unified formulas for the coefficients µ n,m and, hence, the nth-order component y n (t) will be determined.This is the issue of the next section.

Closed-Form Solution via Quantum Calculus
The objective of this section is to obtain a closed-form solution by means of some properties in the quantum calculus [35].We start with the product: where (α : β) i denotes the Pochhammer symbol.For α = β = c, we have Using the above property, we are able to construct the closed-form solution of Equation ( 1) by means of quantum calculus notations as addressed in the theorem below.
Theorem 2. The solution y(t) of the PDDE (1) is given by Proof.Let us start the proof by rewriting the general component y n (t), from Theorem 1, in the form: The products in the denominator of (64) can be written in view of quantum calculus notations as Substituting ( 65) into (64) and simplifying, then The homotopy series solution is given in terms of the auxiliary parameter p as and, hence, the solution y(t) can be obtained by letting p → 1.Thus, From ( 66) and (68), we obtain i.e., which completes the proofs.

Discussion
In previous sections, two types of closed-form solutions were established.In this discussion, the behaviors of the two types of solutions, mainly the PSS and the EFS, are explored.Before launching, it is better to mention that the PSS in Equation ( 17) reduces to several forms of exact solutions at c = −1, as pointed out in [19].In this case, the PDDE (1) becomes y (t) = ay(t) + by(−t), y(0 For a = b = α, the exact solution of Equation ( 71) reads [20] For b > a, the periodic solution of Equation (71) was obtained in the form [19]: The r-term approximation Θ r (t) of the EFS solution is given from Equation (70) as At the special case c = −1, the approximations Θ r (t) can be obtained by taking the limit of Equation (74) as c → −1: So, we focus in the first part of this discussion on performing comparisons between the approximations Θ r (t) in Equation ( 75) and the exact solutions given by Equation (72) (a = b = α) and Equation (73) (b > a).Figures 1 and 2 show the comparisons between the approximations Θ 5 (t), Θ 6 (t), Θ 7 (t) (Figure 1) and the approximations Θ 15 (t), Θ 20 (t), Θ 25 (t) (Figure 2) with the exact solution (72) when λ = 1, a = b = α = 1.It is seen from these figures that the domain of coincidence between the approximations Θ r (t) and the exact solution increases as the number of terms r increases.Additional comparisons are depicted in Figures 3 and 4 between the approximations Θ 15 (t), Θ 20 (t), Θ 25 (t) and the periodic solution (73) when λ = 1, a = 1, b = √ 2 (Figure 3) and λ = 1, a = 2, b = 3 (Figure 4).It is also observed from these figures that Θ 25 (t) is in full agreement with the exact periodic solution (73) for a complete period indicated by the points in pink and yellow.For c = −1, the accuracy of the approximations Θ r (t) is validated through calculating the residuals RE r (t) defined by The numerical results in Figures 5-10 for the obtained residuals reveal the high accuracy of the present analysis, even in a huge domain.

Conclusions
In this paper, the HPM is applied to solve the PDDE using two different canonical forms.Accordingly, two types of closed-form solutions, the PSS and EFS, were obtained.It was shown that the present PSS reduces to the corresponding ones in the literature [12,15] for the ADDE in special cases.A unified formula for the homotopy components was successfully obtained for the EFS.The results were numerically validated through several comparisons with the available exact solutions.The calculated residuals reflect the high accuracy of the current analysis, especially, residuals were close to zero even in a huge domain of the studied PDDE.The obtained results reveal that the current approach is effective and straightforward and can be further applied to include other DDEs of higher orders [36][37][38][39].