A Quadruple Integral Involving Chebyshev Polynomials T n ( x ) : Derivation and Evaluation

: The aim of the current document is to evaluate a quadruple integral involving the Chebyshev polynomial of the ﬁrst kind T n ( x ) and derive in terms of the Hurwitz-Lerch zeta function. Special cases are evaluated in terms of fundamental constants. The zero distribution of almost all Hurwitz-Lerch zeta functions is asymmetrical. All the results in this work are new.


Significance Statement
Named for the Russian mathematician Pafnuty Lvovich Chebyshev (1821-1894), this type of polynomial function is detailed in [1]. Chebyshev polynomials appear in almost every branch of numerical analysis, and they are especially important in recent breakthroughs in orthogonal polynomials, polynomial approximation, numerical integration, and spectral approaches [2]. Definite integrals involving Chebyshev polynomials are investigated in the works of Wang and Hu [3], where the authors derived an identity containing the integral Chebyshev Polynomials of the First-Kind. The Mellin transform of Chebyshev Polynomials is tabled in the book of Oberhettinger [4]. The definite integral of special functions is investigated in the work by McClure and Wong [5] where the authors derived an expansion for a multidimensional integral containing the Bessel function, which is used in the study of crystallography. Watson [6], investigated a quadruple integral where he expressed it in terms of the ratio of Gamma functions.
Based on current literature, there exists work conducted on definite integrals of the Chebyshev Polynomial, and also present in current literature is work conducted on quadruple integrals involving special functions though sparse. In this work we propose a systematic approach to deriving a quadruple integral involving a special function and derived it in terms of the Hurwitz-Lerch zeta function. This approach is based on the contour integral method in [7]. The goal is to expand upon work on multiple integrals involving special functions with the aim of assisting researchers where this work is useful. In this present work, we investigate the quadruple integral of the Chebyshev polynomial T n (x) and derive an integral transform which is invariant under the index n with respect to the Hurwitz-Lerch zeta function.

Introduction
The Chebyshev polynomial T n (x) see Section (22:3) in [1] has an algebraic definition given by: the polynomial has a Gauss hypergeometric function representation given by the polynomial can be expressed by the Rodrigue's formula and has a contour integral representation given in [8] by where the contour encircles the origin and no other singular points. The quadruple definite integral derived in this work is 1 0 where the parameters k, a, n, m are general complex numbers and (n) < (m). This definite integral will be used to derive special cases in terms of special functions and fundamental constants.

The Generalized Hankel Contour
An alternate derivation of Equation (2) in [9] when Re(m + w) ≤ 0 can be achieved by recalling a variant of Hankel's formula involving the Gamma function: where the contour starts from −∞ along the lower side of the negative real axis, encircles the origin once in the positive (counter-clockwise) direction, and then goes back to −∞ along the upper side of the negative real axis. Replacing t by tu in (2) gives Here, u ∈ C with suitable restriction on C u is another contour of the Hankel type depending on the arg u = θ of the polar form u = |u|e iθ : For instance, (i) If θ = 0, then

The Present Case of the Contour Integral
The derivations follow the method used in [7,9]. In the present case the cut approaches the origin from the interior of the first quadrant and the cut lies on opposite sides of the cut going round the origin with zero radius. The variable of integration in the contour integral is α = w + m. Using a generalization of Cauchy's integral formula we form the quadruple integral by replacing y by log then taking the definite integral with respect from Equation (1.9.9.5) in [4] and Equation (4.215.1) in [10] where (m + w) > 0, (n) < (m) and using the reflection Formula (8.334.3) in [10] for the Gamma function. We are able to switch the order of integration over x, y, z and t using Fubini's theorem for multiple integrals see (9.112) in [11] since the integrand is of bounded measure over the space

The Hurwitz-Lerch Zeta Function and Infinite Sum of the Contour Integral
In this section, we use Equation (2) in [9] to derive the contour integral representations for the Hurwitz-Lerch zeta function. The significance of this section is to derive a special function equivalent to the definite integral of the contour integral derived in Section 3 in terms of the same contour integral.

Definite Integrals and Invariant Index Forms
In this section we evaluate Equation (8) such that the index of the Chebyshev polynomials T n (x) is independent of the right-hand side. These types of integrals could involve orthogonal functions. We also derive special cases of Equation (8) Proof. Use Equation (8) and set m = 1/2, a = −2 and simplify using entry (4) in Table  below (64:12:7) and entry (2) in Table below (64:7) in [1].

Discussion
We provide a new approach for deriving a new integral containing the Chebyshev polynomial Tn(x) in this study. Several intriguing definite integrals can be obtained using this contour integral method. For both real and imaginary and complex values of the parameters in the integrals, the findings were mathematically validated using Wolfram Mathematica.
Author Contributions: Conceptualization, R.R.; methodology, R.R.; writing-original draft preparation, R.R.; writing-review and editing, R.R. and A.S.; funding acquisition, A.S. All authors have read and agreed to the published version of the manuscript.