Series Solution of the Pantograph Equation and Its Properties

In this paper, we discuss the classical pantograph equation and its generalizations to include fractional order and the higher order case. The special functions are obtained from the series solution of these equations. We study different properties of these special functions and establish the relation with other functions. Further, we discuss some contiguous relations for these special functions.


Introduction
Ordinary differential equations (ODE) are widely used by researchers to model various natural systems.However, it is observed that such equations cannot model the actual behavior of the system.Since the ordinary derivative is a local operator, it cannot model the memory and hereditary properties in real-life phenomena.Such phenomena can be modeled in a more accurate way by introducing some nonlocal component, e.g., delay in it.
The characteristic equation of delay differential equations (DDE) is a transcendental equation in contrast with the polynomial in the case of ODE.Hence, the DDEs are difficult to analyze as compared with ODEs.
Various special functions viz.exponential, sine, cosine, hypergeometric, Mittag-Leffler and gamma are obtained from ODEs [1].If the equations have variable coefficients, then we may get the Legendre polynomial, the Laguerre polynomial, Bessel functions, and so on [2].However, there are vary few papers that are devoted to the special functions arising in DDEs [3].This motivates us to work on the special functions emerging from the solution of DDE with proportional delay.We analyze different properties of such special functions and present the relationship with other functions.

Definition 2 ([5]
).The Riemann-Liouville integral of order µ, µ > 0 is given by: (2) Definition 3 ([5]).The Caputo fractional derivative of f is defined as: Definition 4 ([6]).Two special functions are said to be contiguous if their parameters differ by integers.The relations made by contiguous functions are said to be contiguous function relations.

Pantograph Equation
The pantograph is a current collection device, which is used in electric trains.The mathematical model of the pantograph is discussed by various researchers [8][9][10][11].
The differential equation: with proportional delay modeling these phenomena is discussed by Ockendon and Tayler in [12].Equation ( 4) is called the pantograph equation.Kato and McLeod [13] showed that the problem (4) is well-posed if q < 1.Further, the authors discussed the asymptotic properties of this equation.The coefficients in the power series solution are obtained by using a recurrence relation.
In [14], Fox et al. showed that the solution of ( 4) is given by a power series: Iserles [15] considered a generalized pantograph equation y (t) = Ay(t) + By(qt) + Cy (qt), y(0) = y 0 , where q ∈ (0, 1).The condition for the well-posedness of the problem is given in terms of A, B and C. The solution of the problem is expressed in the form of the Dirichlet series.Further, the advanced pantograph equation y (n) (t) = ∑ l j=0 ∑ m−1 k=0 a j,k y (k) (α j t), t ≥ 0, where a j,k ∈ C and α j > 1 for all j = 0, 1, . . ., l, is also analyzed by Derfel and Iserles [16].In [17], Patade and Bhalekar discussed the pantograph equation with incommensurate delay.Various properties of the series solution obtained are discussed.
In this paper, we write the solution (5) in the form of a special function and study its properties.We discuss the generalization of (4) to fractional order and the higher order case, as well.

Special Function Generated from the Pantograph Equation
We write solution (5) in the form of following special function: The notation R is used in memory of Ramanujan [18].
Corollary 1.The power series ( 6) is absolutely convergent for all t, and hence, it is uniformly convergent on any compact interval on R.
Theorem 4. For m ∈ N, we have: Proof.We prove the result by using the induction hypothesis on m.From Theorem 3, the result is true for m = 1.Suppose, Consider, Using Theorem 3, we have: aq r R(a, b, q, q r t) + bq r R(a, b, q, q r+1 t) This completes the proof.
Theorem 6.If q ∈ (0, 1), then we can express the integral of e −t R in the following series containing q−binomial coefficients: Proof.Consider: Theorem 7. If q ∈ (0, 1), then we have the following relation: Proof.Consider: Theorem 8.The function R shows the following relationship with the lower incomplete gamma function γ: Proof.We have: By using the property [20]: Theorem 9.The following integral gives the relation between the function R and the upper incomplete gamma function Γ: By using the property [20]:

Properties, Relations and Identities of R(a, b, q, t)
We state the following properties of R(a, b, q, t) without proof.

Relation to Other Functions
where B(t, n, 1) and I t (n, 1) are the incomplete beta function and regularized incomplete beta function [4].

Generalizations to Include the Fractional Order Derivative
Consider the fractional delay differential equation with proportional delay: where 0 < α ≤ 1, q ∈ (0, 1), a ∈ R and b ∈ R.
The solution of ( 8) is: We denote the series in (9) by: R α (a, b, q, t) Theorem 10.If q ∈ (0, 1), then the power series: is convergent for all finite values of t.

Generalizations to the Higher Order Case
Consider the system of delay differential equation with proportional delay: where q ∈ (0, where A ∈ M n×n and B ∈ M n×n is convergent for t ∈ R. 6.3.Properties and Relations of R(A, B, q, tI) We get the following analogous properties and relations as in Section 5 under the condition that A, B, A ± I and B ± I are invertible.

Theorem 14 .
1), A = a ij n×n , B = b ij n×n and Y = (12) is:R(A, B, q, tI) = I + If (A + Bq n ) is invertible for each n, then the power series: R(A, B, q, tI) = I +