A Systematic Approach to Delay Functions

: We present a systematic introduction to a class of functions that provide fundamental solutions for autonomous linear integer-order and fractional-order delay differential equations. These functions, referred to as delay functions, are deﬁned through power series or fractional power series, with delays incorporated into their series representations. Using this approach, we have deﬁned delay exponential functions, delay trigonometric functions and delay fractional Mittag-Lefﬂer functions, among others. We obtained Laplace transforms of the delay functions and demonstrated how they can be employed in ﬁnding solutions to delay differential equations. Our results, which extend and unify previous work, offer a consistent framework for deﬁning and using delay functions.


Introduction
Special functions play a vital role in simplifying the analysis of many problems in mathematical physics and applied mathematics [1], especially in the solution of differential equations.For example, the exponential function and the Mittag-Leffler function [2,3] are fundamental in composing solutions for linear integer-order and linear fractionalorder differential equations, respectively.Given the importance and prevalence of delay differential equations in mathematical modeling across a diverse array of disciplines, including biology [4,5], economics [4], engineering [6], pharmacokinetics [7], physics [8], population dynamics [9], and traffic modeling [4], it is natural to seek similar special functions for linear integer-order and linear fractional-order delay differential equations.
Indeed, there has been some progress in this direction, such as the Lambert W function [10] and delayed exponential functions for linear first-order delay differential equations [11,12], delay sine and cosine functions for linear second-order delay differential equations [13], and delayed Mittag-Leffler functions for fractional-order delay differential equations [14][15][16][17].The introduction of such functions has proliferated in recent years; to date, this has been conducted in a non-uniform way, with varying notations and different types of series representations.
Certain special delay functions have been previously introduced for particular delay equations.In this work, we provide a systematic approach for obtaining representations to delay functions that provide solutions to delay differential equations.This provides a uniform representation of delay functions and it enables the identification of a large class of new delay functions that generalize special functions.First, we show how we can use the known series expansion solution of a given differential equation to create a corresponding delay function series expansion solution for a corresponding delay differential equation.In this way, we introduce a large class of previously unknown delay functions that provide solutions to delay differential equations.We then introduce a new solution method for delay differential equations by showing how solutions for arbitrary autonomous linear delay differential equations can be constructed from trial delay function series expansions.We also present the first table of Laplace transforms for delay functions, facilitating the use of Laplace transform methods in solving linear delay differential equations.Finally, we have generalized our approach using fractional power series and generalized power series to derive solutions for linear fractional-order delay differential equations and certain non-autonomous delay differential equations.This generalization is a further novel aspect of our systematic approach to delay functions as solutions of associated delay differential equations.
We hope that our approach will highlight the importance of delay functions, including those defined previously, and make them more accessible to a broader mathematics community.In general terms, delay functions are represented as truncated power series, truncated fractional power series, or truncated generalized power series, each incorporating a delay parameter σ, with the variable in the nth term of the series delayed by nσ.If the delay function is taken to be a function of a real variable, the series expansion can be truncated using a Heaviside function, which also facilitates the calculation of Laplace transform properties.
In Section 2, we define delay functions through truncated power series with delays incorporated in the series representation.Explicit definitions are given for delay exponential functions, delay trigonometric functions, delay hyperbolic trig functions, and delay hypergeometric functions.We calculated the derivatives of these functions and explored some of their basic properties.In Section 3, we define fractional delay functions through truncated fractional power series with delays incorporated in the series representation.Delay Mittag-Leffler functions and generalized delay Mittag-Leffler functions have been defined through this approach.In Section 4, we provide Laplace transform results for delay functions and fractional delay functions.In Section 5, we show how standard delay functions arise as solutions of autonomous linear integer-order delay differential equations, and in Section 6, we show how fractional delay functions arise as solutions of autonomous linear fractional-order delay differential equations.In Section 7, we introduce generalized power series delay functions and illustrate their application to a non-autonomous delay differential equation with a periodic coefficient.We provide a brief summary in Section 8.

Standard Delay Functions
Here, we define standard delay functions as truncated power series with the nth term delayed by nσ.The Heaviside function, which facilitates the analysis of truncated power series, is introduced first.Definition 1.The Heaviside function of a real variable x is defined as See, for example, Lighthill [18], for further details on Heaviside and related functions.

Definition 2.
A delay function of a variable x with delay parameter σ is defined as The Heaviside function truncates the series expansion; thus, we can write This explicit truncated form is preferred for defining delay functions over complex variables, x, σ ∈ C, with x σ ∈ R. Note that if σ ∈ R >0 , then the domain for the delay function is restricted to x ∈ R ≥0 , and conversely, if σ ∈ R <0 , then the domain for the delay function is restricted to x ∈ R ≤0 .
The definition of the delay function allows for a natural correspondence between a function and a delay function.This correspondence is noteworthy because some properties and relations of the original function have analogs in the delay function.In particular, this will allow us to correlate the solutions of certain ordinary differential equations with those of delay differential equations.In brief, for a real-valued function, f (x), defined by a power series we can define a corresponding delay function which is convergent on the domain of f (x).The convergence is assured by the truncation of the series in the definition of the delay function.The introduction of delay functions through infinite series expansions without the Heaviside function would typically lead to non-convergence.For example, if f (x) is an analytic function with a n = f (n) (0) n! , then ∑ +∞ n=0 a n (x − nσ) n would fail the nth term test at x = 0, σ = 1; viz lim n→+∞ n n n! → +∞.
Example 1.If we consider the exponential function, then we can define a corresponding delay function If we consider the exponential function, then we can define a corresponding delay function Note that with these definitions, df(−x; −σ) = dg(x; σ).
To avoid possible ambiguity, we define the delay exponential function below.

Definition 3.
The delay exponential function is defined as A delayed exponential function, similar to the above, has been previously introduced using a slightly different notation in [11,12].It immediately follows from Definition 3 that which identifies the parameter λ as a common scale factor for x and σ.Alternative representations of this function are also possible.For example, if t > 0 and 0 < τ < 1 e , then we can use Laplace transform methods, as shown below, to write where W k (z) denotes the kth branch of the Lambert W function [10], which satisfies Figure 1 shows plots of the delay exponential functions dexp(x; σ) and dexp(−x; −σ) for three different delay values.For positive x and σ, the behavior is similar to the standard exponential function, with the growth slower for larger delays.For negative x and σ, the behavior is similar to the standard exponential function for very small delays, but for larger delays, the delay exponential function exhibits oscillatory behavior.
The following lemma shows that the derivative of the delay exponential function is proportional to the function itself but with x replaced by x − σ.Lemma 1.The derivative of the delay exponential function is given by Proof.We first note that the nth-generalized derivative of the Heaviside function is given by where δ (n−1) (x) denotes the (n − 1)th derivative of the Dirac delta generalized function δ(x).Hence, by differentiating Equation (11), term-by-term, gives = λdexp(λ(x − σ); λσ), (19) where the only contribution from the Dirac delta sum is when n = 0 but this term vanishes on the domain x σ ∈ R\{0}.We next focus on delay trigonometric functions and begin with the definition of the delay cosine and delay sine functions.Definition 4. The delay cosine function is defined as Definition 5.The delay sine function is defined as The delay sine and cosine functions were previously introduced by [13] in connection with second-order delay differential equations.
It is evident from Definitions 4 and 5 that dcos(−x; −σ) = dcos(x; σ) and dsin(−x; −σ) = −dsin(x; σ), which mirrors the even and odd nature of the cosine and sine functions, respectively.Plots of the delay cosine and delay sine functions are shown in Figure 2 for three different delay values.These delay functions are oscillatory but not periodic for finite delays.Also, the amplitude of the oscillations increases as the delay increases.The next lemma shows that the derivatives of the delay cosine and delay sine functions are analogous to the derivatives of their non-delay counterparts.
Lemma 2. The derivatives of the delay cosine and delay sine functions are given by Proof.We verify Equation ( 22) by differentiating Equation ( 20) term-by-term, with x and σ scaled by the parameter λ, to give where we note that the only contribution from the Dirac delta sum is when n = 0 but this term vanishes on the domain x σ ∈ R\{0}.Similarly, Equation ( 23) follows by differentiating Equation ( 21) term-by-term, with x and σ scaled by the parameter λ, to give = λdcos(λ(x − σ); λσ), (29) where we note that the Dirac delta sum vanishes for all n ≥ 0.
The following proposition shows that Euler's formula e ix = cos(x) + i sin(x) can be extended to delay functions.Proposition 1. Euler's formula for delay functions is given by the equation Proof.Starting with Equation (11), where λ = i, we observe that dexp(ix; iσ) Using Proposition 1, we can obtain other identities involving delay trigonometric functions.An example of such an identity is given in the following corollary.
We can also formulate the delay cosine and delay sine functions in terms of delay exponential functions in a similar manner as their non-delay counterparts.
The delay cosine and delay sine functions both have an analog delay hyperbolic function, which we define next.Definition 6.The delay hyperbolic cosine function is defined as Definition 7. The delay hyperbolic sine function is defined as It is clear from the above definitions that dcosh(−x; −σ) = dcosh(x; σ) and dsinh(−x; −σ) = −dsinh(x; σ), mirroring the even and odd nature of the hyperbolic cosine and hyperbolic sine functions, respectively.
Plots of the delay hyperbolic cosine and delay hyperbolic sine functions are shown in Figure 3 for three different delay values.The derivatives of the delay hyperbolic cosine and delay hyperbolic sine functions closely resemble the derivatives of their non-delay counterparts.Lemma 3. The derivative of the delay hyperbolic cosine and delay hyperbolic sine functions are given by Proof.The method to verify Equations ( 45) and ( 46) is almost identical to that presented in Lemma 2 but with the initial removal of the (−1) n term from Equations ( 24) and ( 27).
The next theorem shows that the hyperbolic identity exp(x) = cosh(x) + sinh(x) can be extended to delay functions.Proposition 2. The delay hyperbolic cosine and delay hyperbolic sine functions satisfy the equation Proof.Starting with Equation ( 10), we observe that Using Proposition 2, we can obtain other identities involving delay hyperbolic trig functions.An example of such an identity is stated in the next corollary.

Corollary 3. A relation between the delay hyperbolic cosine and delay hyperbolic sine functions is given by
Proof.From Proposition 2, we can write which simplifies to yield Equation (52).
We can also express the delay hyperbolic cosine and delay hyperbolic sine functions in terms of delay exponential functions.
We can use the known series expansions of the inverse hyperbolic trigonometric functions to construct corresponding series expansions for delay inverse hyperbolic trig functions.It is important to note that the delay inverse hyperbolic trig functions are not inverse functions for the delay hyperbolic trig functions.The delay inverse hyperbolic trig functions do however share some of the other properties with their non-delay counterparts.This is highlighted in the following example.

Example 2. Starting with the series expansion for the inverse hyperbolic tangent function
we identify the corresponding delay inverse hyperbolic tangent function The derivative of the inverse hyperbolic tangent function is which we recognize as the series expansion for The derivative of the delay inverse hyperbolic tangent function is where is the corresponding delay function for f (x) defined above.Thus, the derivative of the delay inverse hyperbolic tangent function is the corresponding delay function for the derivative of the standard inverse hyperbolic tangent function.
We now formulate the delay function of the hypergeometric function.We begin with a definition.Definition 8.The delay hypergeometric function is defined as Here, (c) n denotes the rising factorial Pochhammer symbol for n ≥ 0, i.e. (c The next lemma presents the derivative of the delay hypergeometric function.
Lemma 4. The derivative of the delay hypergeometric function is given by Proof.Differentiating Equation (68) term-by-term gives where we note that the only contribution from the Dirac delta sum is when n = 0 but this term vanishes on the domain x σ ∈ R\{0}.

Fractional Delay Functions
Here, we define fractional delay functions based on fractional power series, which have been proven to be useful for solving fractional-order differential equations [19][20][21].We begin with the definition of a fractional power series.Definition 9.The fractional power series of a real function, f (t), is an infinite series of the form We now expand our definition of delay functions by defining fractional delay functions, which may prove useful for solving fractional delay equations.Definition 10.For 0 < α < 1, a fractional delay function of a variable x with delay parameter σ is defined as Arguably, the most fundamental function in fractional calculus is the function e α (x) = E α (−x α ) [3], where is the Mittag-Leffler function and the parameter α ∈ (0, 1).The fractional power series representation of the function e α (x) can be used as a basis to define delay fractional Mittag-Leffler functions.
Definition 11.For 0 < α < 1, the delay fractional Mittag-Leffler function is defined as and The delay Mittag-Leffler function, similar to Equation (77) for x ∈ R + , was introduced as a solution of a linear homogeneous [15] and a linear non-homogeneous [14,17]  Next, we generalize the delay fractional Mittag-Leffler function given by Equation (77) via the inclusion of two additional parameters.Definition 12.For 0 < α < 1 and β, γ ∈ R ≥0 , the generalized delay fractional Mittag-Leffler function is defined as Interestingly, we can express the derivative of the generalized delay fractional Mittag-Leffler function, Equation (79), as a linear combination of three generalized delay fractional Mittag-Leffler functions, but with different parameters and arguments.Lemma 5.For 0 < α < 1 and β, γ ∈ R ≥0 , the derivative of the generalized delay fractional Mittag-Leffler function is given by Proof.Differentiating Equation (79) term-by-term gives where we note that the Dirac delta sum vanishes on the domain x σ ∈ R\{0}, irrespective of the value of γ ∈ R ≥0 .We then simplify the remaining sum in Equation (82) by replacing Substituting the results from Equation (86) into Equation (82), we obtain Equation (80).
The following lemma provides two special cases of Lemma 5.
Lemma 6.For 0 < α < 1, two derivatives of the generalized delay fractional Mittag-Leffler function with different values for β and γ are given by Proof.This follows from an application of Lemma 5 with the specified values of β and γ.

Laplace Transform Results
In this section, we derive results that enable the Laplace transform of both standard and fractional delay functions to be readily determined.These results essentially rely on finding the delay function corresponding to the given power series and then taking the Laplace transform term-by-term of the resulting sum.

Standard Delay Functions
We first consider the Laplace transform of standard delay functions.
Theorem 1.If f (t) permits a power series representation, Equation (4), with Laplace transform then the Laplace transform of the corresponding delay function df(t; τ) is given by Proof.Starting with the power series we take the Laplace transform term-by-term with to obtain From Equation (92), we identify the corresponding standard delay function and we take the Laplace transform term-by-term and use to write = e sτ f (se sτ ) (99) The second result follows in a similar manner.Starting with and proceeding as above, we have The following example demonstrates how the previous theorem can be used to obtain the Laplace transform of a standard delay function.
Example 3. The exponential function, defined by the power series, It then follows from Theorem 1 that the corresponding delay exponential function has Laplace transform (108) Thus, in this example, we have (109) It then follows from Equation (105 (110) Example 4. The representation of the delay exponential function dexp(−t; −τ) in terms of Lambert W functions, Equation ( 12), can now be recovered by finding the inverse Laplace transform of using the complex inversion formula In this formula, γ ∈ R is to the right of all of the singularities of F(z).If 0 < τ < 1 e , then it is known [10] that F(z) has simple poles, at z + e −zτ = 0 given by z * = W k (−τ) τ .It is a simple exercise in complex analysis to show that the Bromwich contour integral is equivalently given as the modified Bromwich contour integral where Γ is a closed path from γ − iL to γ + iL with γ > 0 on the right side plane, completed with a circular arc of radius R on the left-hand side, in the limit L, R → +∞.The modified Bromwich contour integral can then be evaluated using the Cauchy residue theorem, In Table 1, we provide a summary list of some standard delay functions and their Laplace transforms, which can be obtained via Theorem 1.
Table 1.The left-hand column lists some special delay functions with their corresponding Laplace transforms in the right-hand column.

Fractional Delay Functions
Next, we consider the Laplace transform of fractional delay functions.
Theorem 2. If f (t) permits a fractional power series representation (Definition 9) with the Laplace transform then the Laplace transform of the corresponding fractional delay function df α (t; τ) is given by Proof.Starting with the fractional power series we take the Laplace transform term-by-term with From Equation (120), we identify the corresponding fractional delay function and again, we take the Laplace transform term-by-term and use to write The second result follows in a similar manner.Starting with The next example demonstrates how the previous theorem can be used to determine the Laplace transform of a fractional delay function.
Example 5.The Le Roy-type function [22] can be modified to yield the fractional power series which has Laplace transform It then follows from Theorem 2 that the corresponding delay Le Roy-type function, has the Laplace transform (138) Thus, in this example, we have It then follows from Equation ( 133) with λ = 2 that Note that our notation for the delay Le Roy-type function follows the notation used for general fractional delay functions stated in Definition 10.
In Table 1, we present a summary list of some fractional delay functions and their Laplace transforms that can be obtained via Theorem 2.

Extension of the Laplace Transform
Throughout this section, we restrict the time domain of the Laplace transform to nonnegative real numbers.However, since delay functions vanish for all negative real numbers due to the Heaviside function, we can extend the time domain to include all real numbers.Consequently, we can consider more general integral transforms.For instance, the Laplace transform considered previously is equivalent to the bilateral Laplace transform, which we define as If we make the change of variable s = 2πiω with ω ∈ R, then the bilateral Laplace transform is equivalent to the Fourier transform, which we define as Hence, if the Laplace transform of a delay function is known, we can use this variable change to readily obtain its Fourier transform.The next example helps to elucidate this idea.
Example 6. Recall from Equation (108) that the Laplace transform of the delay exponential function is (144) It then follows that the Fourier transform of the delay exponential function is (145)

Integer-Order Delay Differential Equations
Here, we exploit the ability of standard delay functions to solve linear integer-order delay differential equations.As a preliminary to more general results, we first consider examples of first-order differential equations with delay exponential functions in their solutions.
Example 7. The delay differential equation has the solution To show this, we first verify that it is a solution to the equation.From Lemma 1, we have so that, given u(t) in Equation ( 148), we have It remains to show that the solution, Equation ( 148), matches the initial condition.First, we note that hence, dexp(t + τ; τ) = 1 for −τ ≤ t ≤ 0 and dexp(t Clearly, the integral vanishes if t > t + τ, so that In the special case with φ(t) = 1, we recover u(t) = dexp(t + τ; τ).The dynamics of the solution given by Equation ( 148) are illustrated in Figure 5 for various delays.
We can readily construct more complicated equations that involve delay functions in their solutions.This is highlighted by the next example.

Example 8. The delay differential equation
with the initial condition has the solution u(t) = e −at dexp(be aτ t; be aτ τ). (155) To show this, let v(t) = e at u(t), so that Equation ( 153) is consistent with which has solution dexp(be aτ t; be aτ τ), and the result then follows.The dynamics of the solution given by Equation ( 155) are illustrated in Figure 6 for various delays.The following theorem shows the more general utility of delay functions in representing solutions to delay differential equations.Theorem 3. If f (x) permits a power series representation, Equation (4), with the jth derivative ) where df(x; σ) is the corresponding delay function for f (x), and dg(x; σ) is the corresponding delay function for g(x).
Proof.Starting with the power series and differentiating j times, we arrive at Hence, we now have From Equation (159), we identify the corresponding delay function and then Hence, differentiating Equation (163) j times, we obtain From Equation (161), we identify the corresponding delay function and then which matches the first sum on the right-hand side of Equation (164).
Example 9. Consider the logarithmic function f (x) = log(1 + x) with derivative g(x) = 1 1+x .Starting with the Taylor series for f (x) we identify the corresponding delay function The Taylor series for g(x) is with the corresponding delay function From Theorem 3, we have Corollary 5.If dy(x; σ) is a delay function with where dz(x; σ) is a delay function Proof.This follows from an application of the chain rule and Theorem 3 with the identification Example 10.Consider dy(x; σ) = dexp(−x; −σ), then a n = (−1) n n! from which it follows that a n+1 (n + 1) = −a n and dz(x; σ) = −dexp(−x; −σ).Theorem 3 in this case yields Note that in the case of f (x) = x, we have and in the case of f (x) = x 2 2 , we have It is worth comparing these results with the non-delay counterparts.We note that exp(−x) is the solution to the differential equation and from Equation (178), dexp(−x; −σ) is the solution to the delay differential equation We also note that exp − x 2 2 is the solution to the differential equation but from Equation (179), dexp − x 2 2 ; −σ is not the solution to the delay differential equation In the next example, we demonstrate the utility of Theorem 3 in solving first-order delay differential equations.
Example 11.The delay differential equation has a solution dexp(−x; −σ).This is easy to establish via the direct substitution of Definition 3 into the above equation with λ = −1.Alternatively, note that d f (x) dx = f (x) has a solution f (x) = exp(x), and by using Theorem 3 with g(x) = exp(x) and a 0 = 1, we can write If λ = −1 this identifies dexp(−x; −σ) as a solution to the delay differential equation, Equation (184).
The solution can also be represented as a linear combination of Lambert W functions [10,23], and it can be generalized to include more general initial conditions [12] and systems of equations [11,12].
Using Theorem 3, we can solve higher-order delay differential equations, which we consider in the following example.
We observe that f (x) = sin(x) is also a solution of d 2 f (x) dx 2 = − f (x) and by using Theorem 3 with g(x) = − sin(x), a 0 = 0, and a 1 = 1, we can write Interestingly, if λ = 1, this identifies dsin(x; σ) as a solution of a separate delay differential equation The use of delay sine and cosine functions to represent solutions of systems of second-order differential delay equations with Cauchy initial conditions was explored in [13].
Corollary 6.If x = kσ with k ∈ N and k < j, then Equation (158) reduces to d j dx j df(λx; λσ) = λ j dg(λ(x − jσ); λσ), Proof.From Equation (158), we observe that the sum on the right-hand side vanishes for all x ∈ R, except x = kσ with k ∈ N and k < j.
The next example highlights how delay function series expansions can solve specific linear delay differential equations.
We also have We now substitute the delay series expansions from Equations ( 193), (194) Thus, we require a n+2 (n + 1)(n + 2) − a n+1 (n + 1) + a n = 0. (198) If we choose initial conditions that are consistent with a 0 = 1 and a 1 = 1, then the recurrence relation can be solved in closed form, yielding Finally, we have the solution We can generalize the example above by imposing an arbitrary initial condition.
Example 14.The delay differential equation has the solution where y p (x) is the particular solution to the equation with the initial condition y(x) = Θ(x) for x ≤ 0 given by Equation (200), provided that x = kσ for k = 0, 1.To show that this is a solution, let y p (x) = y c (x) + y s (x), where y c (x) and y s (x) are defined as and respectively.It is straightforward to show that so that, given y(x) in Equation (203), we have Next, we verify that the solution, Equation (203), matches the initial condition.We note that y p (x), like dexp(x; σ), satisfies the properties The particular solution y(x) = y p (x) is recovered when φ(x) = Θ(x) for x ≤ 0.
Interestingly, the previous two examples are special cases of a more general theorem.
Proof.Substituting Equation (216) into Equation (215) yields which defines a linear homogeneous mth-order recurrence relation with polynomial coefficients to solve for the sequence, a n .The initial terms of the sequence, a 0 , a 1 , . . ., a m−1 , are defined by the initial conditions for Equation (215).

Fractional-Order Delay Differential Equations
We now demonstrate the ability of fractional delay functions to solve various linear fractional-order delay differential equations.Although we focus on differential equations that involve the Caputo fractional derivative, it is possible to consider equations with alternative fractional derivatives.For readers interested in the history of fractional calculus, we refer them to the preface in the classic text on fractional differential equations by Podlubny [24].
Example 15.For 0 < α < 1, the fractional-order delay differential equation with initial condition has the solution u where C D α t is a Caputo fractional derivative of order α [25], To show this, we use Laplace transform methods, starting with the Laplace transform of the Caputo derivative Consistent with the requirement, u(t − τ) = 0 if t < τ, we introduce a Heaviside function on the right-hand side of the equation and take the Laplace transform using The Laplace transform of Equation (219) then gives which we can solve for û(s) with u(0) = 1 to obtain Recall from Table 1 that this is the Laplace transform of dE − α (−t; −τ).To show this, we can take the Laplace transform term-by-term in the definition of dE − α (−t; −τ) with and then sum over n using the series expansion 1 1+x = ∑ +∞ n=0 (−1) n x n .The solution to this example was essentially given in [26] but using different notations.
We can increase the complexity of the previous example through the addition of a function of time, f (t), on the right-hand side of the equation.
Example 16.For 0 < α < 1, the fractional-order delay differential equation with initial condition has the solution To show this, we can take a Laplace transform of the equation and solve for û(s).We require the Laplace transform results which can be obtained from Table 1, together with the convolution rule for Laplace transforms to carry out the inverse Laplace transforms.
We can further generalize Example 15 by imposing an arbitrary initial condition.
Example 17.For 0 < α < 1, the fractional-order delay differential equation has the solution The study of fractional-order delay differential equations is a recent and rapidly growing field [14][15][16][26][27][28][29][30][31].It is anticipated that special delay functions will facilitate the representation of solutions for these types of problems, more generally, and enable more direct comparisons with their non-delay counterparts.

Generalized Power Series Delay Functions
In the above, we developed definitions of delay functions based on power series expansions and fractional power series expansions.Clearly, it is possible to broaden this class of functions in many ways.Our primary focus above has been to define delay functions that could prove useful for finding solutions to delay differential equations.In the following, we consider delay functions based on an infinite series of powers of a delayed function.We refer to this class of functions as generalized power series delay functions.Our goal here has not been to be exhaustive in the identification of delay functions that can be used in solving delay differential equations, but to provide one of many possible generalizations.Definition 13.A generalized power series delay function of a real-valued function h(x) with delay parameter σ is defined as In previous sections, we showed that delay functions based on standard power series, and delay functions based on fractional power series, provided solutions to autonomous delay differential equations and autonomous fractional delay differential equations.In the following, we will show that we can use delay functions based on generalized power series, Equation (13), to provide solutions to certain delay differential equations with periodic coefficients.
In the corollary below, we show that if the periodic function p(x) = p(x − σ) has a mean of zero, then the infinite sum solution, Equation (238), for Equation (237), can be simplified to a closed-form solution.
Corollary 7. If p(x) is a periodic function with period σ and a mean of zero, then the delay differential equation, Equation (5), has a closed-form solution y(x) = Γ(1 + x σ , q(x)) Γ(1 + x σ ) e q(x) , x σ ∈ R, where q(x), defined by Equation (239), is periodic with period σ, and denotes the upper incomplete gamma function.
Proof.We begin by showing that if p(x) is periodic with period σ and has a mean of zero, then q(x), defined by Equation (239), is periodic with period σ.Note that if p(x) is periodic with period σ, then from Equation (239), we can write q(x) = so that q(x) = q(x − σ) is periodic with period σ, and we can replace q(x − nσ) with q(x) in the solution to the delay differential equation from Theorem 5 to obtain or equivalently which simplifies to give Equation (245) since Note that the above result is simple to establish using Laplace transform methods.
We demonstrate the utility of the previous theorem and corollary in the following final example.The dynamics of the solution given by Equation (258) are illustrated in Figure 7 for various delays.

Summary
In the above, we provided a systematic approach to delay functions, defined as truncated power series or truncated fractional power series, characterized by a delay parameter σ, where the variable in the nth term of the series is delayed by nσ.Here, we concentrated on the properties of the functions over real variables.We provided their Laplace transform properties and demonstrated how the delay functions could be used to represent solutions of autonomous linear delay differential equations and autonomous linear delay fractional differential equations.Starting with the series expansion solution of an autonomous linear differential equation, we identified a corresponding delay series expansion as a solution of a corresponding autonomous linear delay differential equation.We then showed how general delay series expansions could be employed as trial solutions to construct solutions for arbitrary delay linear differential equations with constant coefficients.Finally, we considered more general delay functions based on truncated series expansions with delayed arguments in powers of general functions, h(x), rather than powers of x or powers of x α .This systematic approach to delay functions that we provide here should make them more accessible to a broader mathematics community, and bring further attention to important results that have already been made [11][12][13][14][15][16]26].More generally, we hope that this publication will simplify modeling with delay differential equations and inspire others to extend research on novel delay functions and their properties.
fractionalorder delay differential equation.Note that our notation for the delay fractional Mittag-Leffler function follows the notation used for general fractional delay functions stated in Definition 10.Plots of the delay fractional Mittag-Leffler function dE − α (−x; −σ) with α = 0.8, 0.2 are shown in Figure 4 for three different delay values.The delay fractional Mittag-Leffler function dE − α (−x; −σ) has larger amplitude oscillations for larger delays.The oscillations are smoother for larger values of α.Further study is needed to determine if the oscillations vanish for sufficiently small σ and sufficiently large α 1.

Example 13 .
Consider the delay differential equationd 2 y(x) dx 2 − dy(x − σ) dx + y(x − 2σ) = 0, x σ ∈ R. (191) If x = kσ for k = 0, 1, we can use Corollary 6 to seek a delay function series expansion solution to the by substituting this into the equation.It is straightforward to show that y

Theorem 4 .
The mth-order linear delay differential equation of the form c m d m dx m y(x)

Example 19 .
The delay differential equation dy