Eigenfunction Families and Solution Bounds for Multiplicatively Advanced Differential Equations

A family of Schwartz functions W ( t ) are interpreted as eigensolutions of MADEs in the sense that W ( δ ) ( t ) = E W ( q γ t ) where the eigenvalue E ∈ R is independent of the advancing parameter q > 1 . The parameters δ , γ ∈ N are characteristics of the MADE. Some issues, which are related to corresponding q-advanced PDEs, are also explored. In the limit that q → 1 + we show convergence of MADE eigenfunctions to solutions of ODEs, which involve only simple exponentials and trigonometric functions. The limit eigenfunctions ( q = 1 + ) are not Schwartz, thus convergence is only uniform in t ∈ R on compact sets. An asymptotic analysis is provided for MADEs which indicates how to extend solutions in a neighborhood of the origin t = 0 . Finally, an expanded table of Fourier transforms is provided that includes Schwartz solutions to MADEs.


Introduction
The introduction of a relaxing parameter q > 1 in differential equations was found to provide stability properties for their corresponding solutions. This is a phenomenon well-known in numerical analysis where if the Ordinary Differential Equation (ODE) y (t) = f (t, y(t)) , y(t 0 ) = y 0 , is stiff then one can try to use the backward Euler method to obtain the sequence { t n , y n } ∞ n=0 by first considering the algebraic equations t n+1 = t n + ∆t , y n+1 = y n + f (t n+1 , y n+1 ) · ∆t , for small time-steps ∆t > 0. If one can obtain y n+1 explicitly in terms of y n then the iteration scheme often converges much faster, and for longer time intervals, than that provided by the forward Euler method [1], p. 349. That such a principle holds for ODEs as ∆t → 0 + was established through the study of Multiplicatively Advanced Differential Equations (MADEs) as q → 1 + , and will be discussed further in this article. Part of our analysis of stability will require obtaining uniform apriori bounds. This will be achieved in a somewhat general setting, and the consequences will be presented in the form of examples of advanced differential equations.

Solutions of MADEs as Eigenfunctions
In [2] solutions to equations of the form y (t) = ay(qt) + by(t) , y(0) = 1 or 0 (wlog) , (1) were studied for q > 1, a ∈ C, b ∈ R and t ≥ 0. In the case that b = 0, with y(0) = 0, solutions y(t) are referred to as eigenfunctions since y(t) → 0 as t → ∞. Specific asymptotic properties of solutions were obtained in Theorem 10 of [3]. Here we only consider the case that b = 0 and a ∈ R, however the derivatives may be of higher (integer) order than in Equation (1). In addition, we extend solutions of these equations to all t ∈ R so that the eigen equation, referred to as an eigen-MADE, has a solution y(t) ∈ S(R) the Schwartz space of infinitely differentiable functions, with derivatives that decay faster than reciprocal polynomials (as defined in [4] section V. 3). An asymptotic theory near t = 0 can be developed indicating that an extension to t < 0 is quite natural. In this way the special functions that we study are eigenfunctions in L 2 (R), although not in the traditional, local (q = 1) sense. The significance of these functions will be demonstrated by examples, and convergence to familiar functions is obtained on compact subsets of R, as q → 1 + .
In the current work we continue by focusing on global solutions of a MADE on R. In particular, we discuss several techniques for starting with a given global solution to an original MADE and then generating solutions of new related MADEs. This theme will be developed as follows: In Section 2, a known MADE solution first introduced in [25], namely q Cos(t), is used to produce a simple related solutionC q (t) = q Cos(t/ √ q) which is an eigensolution of a MADE in the sense of the Abstract. In turn, C q (t) is then used to obtain a new q-advanced Airy function Aiq(t) satisfying a MADE analogue of the Airy differential equation. Then Aiq(t) itself is used along with convolution to generate families of functions φ q (x, t) solving a q-advanced PDE.
In Section 3, a family of MADE solutions, under convolution and auto-correlation, are seen to produce related solutions of new MADEs. Furthermore, the least-element method in Poincare asymptotics is deployed to find natural extensions to related MADE solutions on the negative real line. A theory of asymptotic extensions to t < 0 is developed to clarify the notion that solutions to MADEs behave smoothly in a neigborhood of the origin. We also give conditions that ensures a natural extension to all of R, as is needed to even consider a Fourier transform. An investigation of the inhomogeneous MADEs that these solve is begun.
In Section 4 we focus on considering solutions of MADEs as perturbations of classical solutions, and, mirroring a more direct convergence proof in Section 2, we exhibit MADE solutions which converge to a classical solution of a damped-oscillation equation-the convergence being uniform on compact subsets of [0, ∞).
In Section 5, we return to the topics of convolution and auto-correlation to observe their impact when applied to MADE solutions. In this paper, we will discuss convolutions, correlations, and Fourier transforms for MADEs. A table of Fourier transforms of global MADE solutions under study here is provided in Section 6.  These will be solutions of new MADEs, for which we obtain new elements in a table of Fourier  transforms. This new table mimics what is often done for Laplace transforms, in the study of linear  constant coefficient ODEs. In various theories of differential equations, convolutions provide a useful tool since general solutions can be determined from fundamental solutions, as demonstrated here in Equation (33). This is one motivation for obtaining solutions to homogenous equations, as appears in Proposition 2.

A Normalized Cosine Example and Extensions
From [25], consider the following Schwartz functions, for q > 1 and all t ∈ R, Next defineC There are several properties that we note. In particular, the functionC q (t) is normalized, in that the uniform bound C q ∞ = 1 holds, after some delicate work performed in [25], for each q > 1. It also solves the following eigen-MADE for all t ∈ R and each q > 1, From (6) we see thatC q (t) satisfies an eigen-MADE in the sense of the Abstract, with E = −1 independently of the advancing parameter q > 1. Note that q Cos (t) = −q q Cos(qt) (as recorded in (10) below) does not have an eigenvalue (−q) independent of q, thus we rely onC q (t) as the appropriate eigen-MADE solution.
SinceC q (t) is not only C ∞ and bounded, but in fact Schwartz, we can obtain its Fourier transform, an operation defined for any f ∈ L 1 (R), aŝ In [25] it was found that where N q was defined in Equation (4) above, and the other normalizing constant is To express the Fourier transform of linear, homogeneous MADEs, we found multiple uses of the Jacobi theta function which allows the association that N q = θ(q 2 ; −1/q), and which ensures that N q = 0 for all q > 1, due to the product formula. It will be of significance to note that the reciprocal 1/θ(q; u), for u ≥ 0, is Schwartz when extended to be identically 0 for u < 0. Critical algebraic properties that we use are θ(q; q p u) = q p(p+1)/2 u p · θ(q; u) , ∀p ∈ Z, u ∈ C * , and v · θ(q; A consequence is that the only zeros of θ(q; u) are for u = −q p for all p ∈ Z. This is obvious from the product definition of θ(q; u) in Equation (8).

Uniform Convergence
Using Taylor series methods as an approach paralleling that in [25] we show: Proposition 1. On any compact subset of R,C q (t) approaches cos(t) uniformly as q → 1 + .

Proof.
A given compact set is contained in an interval [−ρ, ρ] for ρ sufficiently large, so it suffices to prove the theorem on [−ρ, ρ].
From these, by induction on the even order derivatives of q Cos(t), we obtain the higher order derivatives and q Cos (2L+1) (t) = [(−1) L q L 2 q Cos(q L t)] = (−1) L+1 q L 2 +L q Sin(q L t) .
We infer all derivatives ofC q (t) viã Evaluating the derivatives ofC q (t) at t = 0 yields for all L ≥ 0.

Remark 1.
Note that, alternatively, one can express Proposition 1 as A similar convergence proof is given in Section 4, with details related to the novelty of the result.

Application to PDE Example
We are now in a position to obtain q-versions of various equations usingC q (t) as a building block for relaxing equations. For example, define the Airy function (see page 570 in [26]) Some properties of this C ∞ (R) function are that Ai(t) → 0 as |t| → ∞, and Ai(0) > 0. We now show: Proposition 2. The q-advanced Airy function is defined here to be for q > 1. The functions Ai(t) and Aiq(t) satisfy the homogeneous ODE and MADE respectively, for t ≥ 0. Basic properties of Aiq(t) for q > 1, are that Aiq(t) is Schwartz with Aiq(0) > 0. Furthermore, for each T > 0, > 0, and R > T sufficiently large, ∃ q( , T, R) > 1 so that In other words, Aiq(t) → Ai(t) uniformly for t in compact subsets of R, as q → 1 + . (27) may seem rather straight forward, due to the uniform convergence ofC q (t) to cos(t) on compact sets. However, we need to use a careful /3 argument, as demonstrated here.

Remark 2. Verifying convergence in Equation
Proof. That Aiq(t) is Schwartz follows from the same forC q (t), whereas the property Aiq(0) > 0 requires a manipulation of theta functions, and is shown in Appendix A. We start with the second equation in (26) since the first equation is known to hold [26]. First define the functioñ Now compute, using v = q(u 3 /3 + ut), and w = q 1/3 u, for t ≥ 0, Next, to show convergence, consider any > 0 and, without loss of generality, fix T > 1. Let t be in the interval |t| ≤ T. Then, for any R > T, using integration by parts and boundedness of the sine function, we can write Thus, for all |t| ≤ T we can easily find R > T sufficiently large so that The bound in Equation (30) also holds if Ai(t) is replaced with Aiq(t) since |C q (t)| ≤ 1 and |S q (t)| ≤ 1 for all q > 1. Now, fix R > 0 sufficiently large so that the bounds in (30), and also (30) with cos replaced byC q , are less than /3. It is essential to note that this value of R is independent of q > 1.
Finally, for each t ∈ R, define the function From the uniform convergence in Equation (24) we can choose q( , T, R) > 1 so that for |t| ≤ T , |u| ≤ R, and 1 < q < q( , T, R). This is now sufficient to verify the expression in Equation (27).

A q-Advanced PDE Example
The argument in the proof of Proposition 2 shows that knowledge of one MADE can help to generate and study other MADEs. In fact, this extends to Partial Differential Equations (PDEs). For example, consider the linear constant-coefficient Airy PDE [27] for x ∈ R, t ∈ R + 0 , and constant a > 0. To obtain an advanced-type equation, consider the kernel function, defined for each t > 0, for appropriate A 0 (q) = 0, to be determined. For any integrable f (x) and any a = 0, define, (compare with Equation (2.2) of [27]). Recall that the functional operation of convolution for integrable functions g, h ∈ L 1 (R) gives a new function g * h ∈ L 1 (R) defined by where the last equality in Equation (35) is the Convolution Theorem (see [28] Theorem IX.4). To discover the PDE that φ q solves, first compute the t-partial derivative of Equation (34), to obtain Now, taking three derivatives of Equation (34) with respect to x, gives By replacing t → q 2 t in Equation (38), one can verify that the q-advanced PDE, for q > 1 and holds. To obtain consistency with the initial data f (x), first define the constant, for q > 1, which is finite since Aiq(t) is Schwartz (thus integrable) for q > 1. Then we require, The q-Airy Hypothesis: Given q > 1, the expression in Equation (40) does not vanish, ie. A 0 (q) = 0.
In Appendix B we show that the q-Airy Hypothesis holds for all q > 1. Then where convergence in Equation (41) is pointwise, and is shown in Appendix C using a mollifier-type argument. If, in addition, we have f ∈ C 1 ∩ L 1 and f ∈ L ∞ , then convergence in Equation (41) becomes uniform.

Solutions of MADEs and Natural Extensions
Define the family of Dirichlet-type functions for t ∈ R + 0 , and q > 1, as introduced in [29], For each µ ∈ Q and λ ∈ Q + , the corresponding function solves the eigen-MADE Here λ/2 = γ/δ ∈ Q + is in reduced form with γ, δ ∈ N. The function f µ,λ (t) has eigenvalue E = (−1) γ+δ q γ(γ+µ)/λ and can be normalized so that the function g µ,λ (t) ≡ f µ,λ t/q γ(γ+µ)/(δλ) solves the q-advanced eigen equation for t > 0. In this manner the q-dependence of the eigenvalue can be removed. Note that the sign of the eigenvalue (−1) γ+δ can dramatically affect the behavior of the solution.

Flat Solutions of MADEs
In [29] we found special conditions under which f µ,λ (t) extends to all t ∈ R, so that gives a Schwartz solution to an associated MADE to all t ∈ R . The essential condition is that f (n) µ,λ (0 + ) = 0 for all n ∈ N 0 , which is a property called flatness, at t = 0. It was shown in [29] that f µ,λ (t) is flat at t = 0 ⇐⇒ µ is an odd integer and λ is an even integer .
This condition for flatness can be expressed as Then, for µ, λ as in Equation (46), F µ,λ (t) all solve first-order MADEs: for t ∈ R and q > 1. See examples in Figure 1. Furthermore, the Fourier transform has a special form: where for each j ∈ {0, 1, 2, . . . , N − 1}, the points of valuation of the theta function require,

A Non-Trivial Extension of a MADE Solution
Now consider the situation where ∃n * ∈ N 0 where f Then an extension of f µ,λ (t) to the region t < 0 is not so clear. However, by truncating the series in Equation (42) an asymptotic exponential-series is obtainable that provides, what appears to be, a smooth extension to the t < 0 region. However, extending in this manner does not lead to a homogeneous, eigen-MADE in the region t < 0. This is demonstrated with a specific example.
We begin by recalling the Airy equation as given in Proposition 2 However, taking the derivative of this equation gives a generalization where the right hand side is expected to be small for t 0. Hence a solution to the constant coefficient equation see Section 4, may be considered to be an approximate solution to the Airy equation near the origin. For example, the function solves (49) with initial conditions Now we consider a q-relaxed version of (50) in the form of a solution to the MADE with parameter q > 1. Note that (52) is a multiplicatively advanced relaxed version of the approximate Airy ODE (49) for q 1 + . From Equations (42) and (43), a particular solution of Equation (52) is To extend η(t) to all of t ∈ R in a C ∞ fashion, we find that is a Schwartz function, where f 1,2/3 (z) is analytic for (z) > 0 and bounded for (z) = 0. Although there is no unique solution to MADEs in general, the function W 1,2/3 (t) constructed in Equation (53) will be called canonical, and it solves the MADE in Equation (52) for all t ∈ R.

Asymptotic Analysis of an Extension
There is an alternate continuous way to extend η(t) to the region t * < t < 0, for t * < 0 defined below, in terms of q > 1. Define the constant C + q so that where the last equality follows from (8). Note that θ(q 3 ; −q) is non-zero for real q > 1 by (9), whence C + q is well-defined and finite. For t ≥ 0 the function η(t) is defined as Now, for t ≥ 0, η(t) solves (52) with initial conditions However, for each t < 0 the function η(t) diverges, due to the rapid growth of e −q k t = e q k |t| , in k, as compared to that of q 3k(k−1)/2 in the summands of (55), as k approaches infinity. Thus, for each t < 0 the function η(t) is not defined.
To remedy this, while keeping the same summands as in (55), we truncate the upper limit of summation in (55). Thus, for all t ∈ R we define the asymptotic extensionη(t) of η(t) bỹ where the integer upper limit of the sum, and the normalizing coefficient, are defined to be Since it will follow from the definition below that N * (q, t) → ∞ as t → 0 − , continuity forη(t) is achieved at t = 0. However, as a solution to a MADE, we have thatη(t) ∈ D , where D is the space of distributions, dual to D ≡ C ∞ 0 (R), the set of compactly supported, infinitely differentiable functions.
The choice of truncation N * (q) presented here is made based on the least-term approximation from Poincaré asymptotics, as presented on p. 94 of Bender and Osrzag [26]: "We look over the individual terms in the asymptotic series; ...For every given value of ... [t]... we locate the smallest term. We then add all the preceding terms in the asymptotic series up to but not including the smallest term." Traditionally this rule gives a good estimate of the actual function, which is often the solution of a differential equation. In our case the rule above can only be applied for t < 0 sufficiently close to the origin, which for this function turns out to be |t| < 3/(e √ q ln(q)) .
This is a consequence of the following more general result.

Proposition 3.
For µ, λ ∈ R with λ > 0, define the following function on t ∈ R for any bounded sequence {a k } ∈ ∞ . Define the exponential growth portion of the summands as Then, define two constants, for fixed q > 1, For t ∈ (t * , 0), the function N * (q, t) exists uniquely as the local minimum of T µ,λ (k, |t|).

Remark 3.
The coefficients a k in Equation (62) play no part in the following analysis. However, if they decay as |k| → ∞, or if they change sign, then the asymptotic behavior may be different than what is derived here.
Proof. Differentiating the exponent h µ,λ (k, |t|) = ln T µ,λ (k, |t|) in (63) with respect to k gives the critical point condition Taking a second derivative of h µ,λ (k, |t|) with respect to k gives the inflection point condition Interpreting the middle critical point condition in (65) as the intersection of the concave up function |t|q k with the fixed line (2k − µ)/λ reveals three possibilities: Case 1: There are two critical points k 1 < k 2 with an intervening inflection point k 3 ∈ (k 1 , k 2 ) for |t| and q sufficiently small. By the first derivative test, a local maximum occurs at k 1 while the desired local minimum then occurs at k 2 .

Case 2:
An edge case occurs, in which the two critical points coalesced to one point equaling the inflection point, k 1 = k 2 = k 3 . There is no local minimum for h µ,λ (k, |t|) in this setting.

Case 3:
There are no critical points when either |t| or q is too large, resulting in no local minimum for h µ,λ (k, |t|) in this setting.
Furthermore, observe that since y = exp(x − 1) is concave up with tangent line y = x at x = 1 then the inequality exp(x − 1) ≥ x holds for all x and equality holds if and only if x = 1. Replacing x by (k − µ/2) ln(q) in our inequality gives Multiplying the inequality on the left through by 2/(λ ln(q)) gives whence (65) also holds at the same value of k = µ/2 + 1/ ln(q). Thus, the critical points and the inflection point coalesced to the common value k = µ/2 + 1/ ln(q) and Case 2 holds. We see that Case 2 holding is equivalent to −t = t * = 2/(λ e q µ/2 ln(q)) holding. Furthermore, one sees that Case 1 holds when |t| < |t * |, and a local minimum is obtained. Thus, the asymptotic phenomena occurs for |t| < |t * | where for the upper index limit N * (q, t) we take the larger of the two solutions to the transendental equation for k * in Equation (65): Then, for |t| < |t * | = 2/(λ e q µ/2 ln(q)) sufficiently small, T µ,λ (k, |t|) has a local minimum at N * (q, t) = k * , which can be found by taking a seed point greater than the value ln (2/(|t|λ ln(q))) / ln(q) of the inflection point and utilizing Newton's method.
From Figure 2 the asymptotic expansion is valid to around t ∼ −2, using N * ∼ 10, rather than t ∼ −5, using N * ∼ 6, contrary to what was expected from Equation (74). This is due to the alternation a k = (−1) k since cancelations require a more careful analysis. This is not done here, but the next example considers a comparatively simple case, which gives a better comparison.

An even Simpler Example of MADE Asymptotics
In this section, we motivate a simpler type of asymptotic extension, distinct from Section 3.4, using two examples.
To begin, we recall a MADE that was studied in [30], for q > 1 and t ≥ 0, where for t ≥ 0 set N(q; t) = ∞. Here we consider the extension to negative values of the parameter. Then, for t * < t < 0, we will choose the constant N(q; t) = N * (q, t * ). To use the asymptotic analysis, note that a k = 1, µ = 1 and λ = 2. Thus, we obtain an approximate MADE solution extension to the region t < 0. Start by defining T ≡ e q j |t| q j(j−1)/2 = e h , where: h ≡ q j |t| − ln(q)(j(j − 1)/2) .
Combining these expressions to eliminate q j |t| gives from Equation (67) which then results in from Equation (68). For t * < t < 0, we have N(q, t) > N(q, t * ) = j * . By inspection, Figure 2 (Right) indicates that we maintain a good asymptotic expansion by letting all N(q, t) = N(q, t * ) = j * . In particular, for q = 1.2 our rule suggests j ≤ j * = N * ∼ 6, which is expected to be valid for t ∈ (−1.8, 0). The Right of Figure 2 indicates a good match for t ∈ (−1, 0), using N * ∼ 10. Finally, we return to Equation (57), and consider the slightly different series, for all t > t * (where where now the integer upper-sum limit, and the normalizing coefficient, are defined to be, respectively The functionη(t) is differentiable for t ∈ (t * , ∞), and solves an inhomogeneous MADẼ wheref * ∈ D is derived in Appendix D. Note thatf * (t) is distinct fromf (t) for t > t * in Equation (59), and the corresponding weak solutionη * (t) is much easier to compute thanη(t), with little consequence to the asymptotics.

Convergence of MADEs to Classical Solutions
In this section, we present another example where we can study convergence of a MADE solution to its classical analogue. This requires an apriori uniform bound in a fixed neighborhood of t = 0 for all q > 1 sufficiently small. Obtaining a uniform-in-q bound for general f µ,λ (t) is rather deep, and complicated by the presence of the alternation (−1) m in Equation (42). Here we study a series without this alternating factor, which defines a function that behaves like a damped oscillation. The details are more challenging than what appears in the proof of Proposition 1, so a full analysis is provided.
Consider the following linear third-order MADE for q > 1, on the interval t ∈ [0, ∞), satisfying the initial conditions For small q > 1, as q → 1 + , Equations (80) and (81) can be considered to be a perturbation of the classical analogue, which is the ODE with initial conditions g(0) = 0, g (0) = 1 g (0) = −1 (83) obtained by setting q = 1 in (80) and (81). One can check directly that (82) and (83) is solved uniquely by Now, using techniques mirroring those of Theorem 3.2 of [29], a particular solution to (80) is for t ≥ 0. Note that the expression in Equation (85) does not have the alternation (−1) k , unlike the expression in Equation (55) for η(t), and this will allow a sharp bound onf (t) for all t ≥ 0, independent of q > 1. The first derivative off (t) is seen to bẽ where the fact that: = cos(2π/3) sin(x) + sin(2π/3) cos(x) = sin(x + 2π/3) , was used explicitly to obtain the last equality in (86). Using this identity implicitly, we obtain: and finally we verify: A re-indexing m = k − 1 was used to move from (88) to (89). Note that (89) gives that (80) holds. From (85)-(87), one sees that where the last equalities of (91) and (92) are obtained from (8).
one sees that f (t) now satisfies the MADE (80) along with the initial conditions (81). The last initial condition follows from the fact that where the last equality in (94) follows from the next lemma.
In addition to the last equality of (94) being proven by the first equality in (95) in Lemma 1, the second equality of (95) proves that the second derivative of W 1,2/3 (t)/W 1,2/3 (0) at t = 0 equals −q.
The following theta function bound will also be helpful.

Lemma 2.
For q > 1 the Jacobi theta function (8) satisfies Proof. Observe that while Comparing each factor in the square brackets in (98) with the corresponding factor in the square brackets of (99) one sees that for all n ≥ 0 from which one concludes that θ(q 3 ; 1) giving the left inequality in (97). The right inequality in (97) holds via the assumption that q > 1.
Next we compute all derivatives of g(t) = 2 exp(−t/2) sin( √ 3t/2) √ 3 at t = 0 and of f (t) = f (t)/f (0) at t = 0, in preparation for the computation of the Taylor series expansion at t = 0 for both g(t) and f (t). From (82) we immediately have that for k ≥ 0 and j = 0, 1, 2 From (102) and (83) one concludes that for k ≥ 0 The analogous results for f (t) =f (t)/f (0) are obtained in the following lemma.
Furthermore, at t = 0 one has Proof. We first establish (104) for the case that j = 0 by induction on k. So for j = 0 note that (104) holds as a tautology for k = 0, and for k = 1 it holds by (89). Assume that f (3k) (t) = q 3 k(k+1)/2 f (q k t) for fixed k. Then where: the inductive hypothesis gives the rightmost equality in (106), and that (89) along with the chain rule gives the first equality in (107). Thus, the j = 0 case holds for all k. Now differentiate the expression f (3k) (t) = q 3 k(k+1)/2 f (q k t) either j = 1 or j = 2 times to obtain (104) in all remaining cases. Evaluating (104) at t = 0 and relying on (90)-(94) gives (105).
Next, the 3N + 2-degree Taylor polynomials P 3N [g](t), P 3N [ f ](t) of g and f , respectively, expanded about t = 0 are given by where (108) follows from (103), and (109) follows from (105). For t ≥ 0, these have respective remainder terms for some ξ ∈ [0, t] and ζ ∈ [0, t]. The goal of uniform convergence on compact subsets is now obtained in the following proposition.

Convolutions, Correlations and Bounds
Here we briefly demonstrate that solutions of MADEs beget new solutions of different MADEs.

Distinction between Convolutions and Correlations
Let f , g ∈ L 1 (R) and recall the standard definitions: Correlation between f and g ≡ Proposition 5. Consider f , g ∈ S(R), which solve the following MADEs respectively, for q > 1, a, b ∈ N, and c f = 0, c g = 0. Then the correlation and convolution solve the following higher-order MADEs Proof. The fact that convolution and correlation preserve the Schwartz property follows from Theorem 3.3 of [31]. The MADE equations easily follow from repeated applications of integration by parts, use of Equation (129), and a change of variables.

Cross-Correlation
Let us consider an example that involves different MADE solutions, to obtain a new MADE. Knowing the Fourier transform of these functions allows us to easily derive properties of the resulting function. Compute, using Plancherel's Lemma, Now, to simplify the integrand in (130), we use the Fourier transforms from [25,32] respectively, to write: The equality in (131) follows from the fact that θ(q 2 ; ω 2 ) = θ(q; iω) θ(q; −iω) and uses the definition of the Jacobi theta function in Equation (8). The consequence is that there are simple poles when ω = −iq k for k ∈ Z, but double poles at ω = iq k . Computing the integral in (130) using residue theory, requires a careful consideration of the position of these poles off the real axis.

Expanded Table of Fourier Transforms
In this final section we establish a short table of Fourier transforms for solutions of MADEs and their relations to Jacobi theta functions. Included are well-established results, along with new functions. The positive constants K 1 and K 2 are generic, but estimates are not presented here.
The introduction of new functions are as follows: For K(t) see [32] for decay constants K 1 and K 2 in Table 1; The functions q Cos(t) and q Sin(t) are closely related toC q (t) andS q (t), respectively, introducted in [25], where constants K 3 and K 4 are obtained; The q-Bessel functions, related to J (t), were introduced in [33], along with decay constants K 5 and K 6 ; Flat wavelets F(t) have Fourier transforms that are averages of theta functions, first derived in [29], along with constants K 7 and K 8 . The functions K * C q (t) and W µ,λ (t), have Fourier transforms that involve theta functions, which can be used to obtain decay parameters K 9 and K 10 .
Note that similar tables for Laplace transforms are quite extensive, since applications only require control of function growth on R + . Here we are concerned with globally defined functions on R for which a Fourier transform can be defined.   Acknowledgments: The authors would like to thank the reviewers for very helpful comments and suggestions that aided in the completion of this project.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript:

ODE Ordinary Differential Equation PDE
Partial Differential Equation MADE Multiplicatvely Advanced Differential Equation

Appendix A. Normalization in Terms of Theta Functions
The normalization forC q (t) in Equation (4) involves a theta function, so that The last expression in Equation (A1) does not vanish for q > 1 due to the product formula in Equation (8). Similarly, we can show that Aiq(0) = 0. Indeed, from the definition, note that using the change of variables w = q (k−1/2)/3 u, ∞ 0 e −w 3 dw · θ(q 2 ; −1/q 4/3 ) θ(q 2 ; −1/q) , and the final expression clearly does not vanish for any q > 1.

Appendix B. Establishing the q-Airy Hypothesis for q > 1
To compute A 0 (q) explicitly, we will find the Fourier transform of Aiq(t) and then find its value at the origin. This requires a careful change of variables. To begin, we combine definition Equation (25) and the inverse Fourier transform of formula in Equation (7) giving To handle the double integral note that the odd power of both the k and the t variables allows the following rearrangement ∞ 0 ∞ 0 exp ik(t 3 /3 + xt) + exp i(−k)(t 3 /3 + xt) We can now obtain the Fourier transform Finally, computing at ω = 0 gives the final result which is clearly a finite, positive, non-zero quantity, for each q > 1.
for N * = N * (q, t * ) fixed. Note that C + q is defined in Equation (54), and C − q is defined in Equation (78), so thatη(0 + ) =η(0 − ) = 0, andη (0 + ) =η (0 − ) = 1. Thus, the first derivative is continuous, and the second derivative is bounded. However, the third derivative results in the appearance of a distribution, which is an inhomogeneous MADE for all t > t * , and which definesf * (t), by inspection of the quantity in the square brackets. The last three terms on the right hand side vanish as q → 1 + , where N * → ∞ in a manner described after the proof of Proposition 3.