Next Article in Journal
Ensuring Topological Data-Structure Preservation under Autoencoder Compression Due to Latent Space Regularization in Gauss–Legendre Nodes
Next Article in Special Issue
A Global Method for Approximating Caputo Fractional Derivatives—An Application to the Bagley–Torvik Equation
Previous Article in Journal
Some Classical Inequalities Associated with Generic Identity and Applications
Previous Article in Special Issue
Some New Families of Finite Orthogonal Polynomials in Two Variables
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Voigt-Type Functions Extended by Neumann Function in Kernels and Their Bounding Inequalities

1
Department of Mathematics, Pondicherry University, Puducherry 605014, India
2
Institute of Applied Mathematics, Óbuda University, Bécsi út 96/b, 1034 Budapest, Hungary
3
Faculty of Maritime Studies, University of Rijeka, Studentska 2, 51000 Rijeka, Croatia
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(8), 534; https://doi.org/10.3390/axioms13080534
Submission received: 14 July 2024 / Revised: 5 August 2024 / Accepted: 6 August 2024 / Published: 7 August 2024

Abstract

:
The principal aim of this paper is to introduce the extended Voigt-type function V μ , ν ( x , y ) and its counterpart extension W μ , ν ( x , y ) , involving the Neumann function Y ν in the kernel of the representing integral. The newly defined integral reduces to the classical Voigt functions K ( x , y ) and L ( x , y ) , and to their generalizations by Srivastava and Miller, by the unification of Klusch. Following an approach by Srivastava and Pogány, we also present the multiparameter and multivariable versions V μ , ν ( r ) ( x , y ) , W μ , ν ( r ) ( x , y ) and the r positive integer of the initial extensions V μ , ν ( x , y ) , W μ , ν ( x , y ) . Several computable series expansions are obtained for the discussed Voigt-type functions in terms of Humbert confluent hypergeometric functions Ψ 2 ( r ) . Furthermore, by transforming the input extended Voigt-type functions by the Grünwald–Letnikov fractional derivative, we establish representation formulae in terms of the associated Legendre functions of the second kind Q η ν in the two-parameter and two-variable cases. Finally, functional bounding inequalities are given for V μ , ν ( x , y ) and W μ , ν ( x , y ) . Particularly interesting results are presented for the Neumann function Y ν and for the Struve H ν function in the form of several functional bounds. The article ends with a thorough discussion and closing remarks.

1. Introduction and Preliminaries

The classical Voigt functions K ( x , y ) and L ( x , y ) occur frequently in several problems of physics such as astrophysical spectroscopy, emissions, the absorption and transfer of radiation in a heated atmosphere, plasma dispersion, and also in the theory of neutron reactions [1]. The functions K ( x , y ) and L ( x , y ) , which are due essentially to Reiche [2], read as follows:
K ( x , y ) = 1 π 0 e y t t 2 4 cos ( x t ) d t , x R ; y R + ,
and
L ( x , y ) = 1 π 0 e y t t 2 4 sin ( x t ) d t , x R ; y R + .
where C , R , R + , Z 0 , and N denote the sets of complex, real, positive real, non-positive, and positive integer numbers, respectively. The Bessel function of the first kind of the order ν has the following power series definition:
J ν ( z ) = k 0 ( 1 ) k z 2 ν + 2 k k ! Γ ( ν + k + 1 ) , z N ; ν C .
where the principal branch of J ν ( z ) should be considered (it corresponds to the principal value of z ν ) and J ν ( z ) is analytic in the z-plane cut along the interval ( , 0 ] . Moreover, for ν Z , the Bessel function of the first kind is entire in z in the whole complex plane [3] (p. 15). Based on the facts [4], (Eq. (10.39.2)) is as follows:
J 1 2 ( z ) = 2 π z sin ( z ) and J 1 2 ( z ) = 2 π z cos ( z ) ,
Srivastava and Miller [5] introduced a generalization of K ( x , y ) and L ( x , y ) in the following form:
V μ , ν ( x , y ) = x 2 0 t μ e y t t 2 4 J ν ( x t ) d t , μ , x , y R + ; ( μ + ν ) > 1 ,
so that
K ( x , y ) = V 1 2 , 1 2 ( x , y ) and L ( x , y ) = V 1 2 , 1 2 ( x , y ) .
Klusch [6] defined an unification (and generalization) of the Voigt functions K ( x , y ) and L ( x , y ) in the following slightly modified form:
Ω μ ν ( x , y , z ) = x 2 0 t μ e y t z t 2 J ν ( x t ) d t , μ , x , y , z R + ; ( μ + ν ) > 1 ,
so that
V μ ν ( x , y ) = ( 2 z ) μ 1 2 Ω μ ν ( 2 x z , 2 y z , z ) ,
or, equivalently,
Ω μ ν ( x , y , z ) = ( 2 z ) μ + 1 2 V μ ν x 2 z , y 2 z .
It is worth mentioning that the extension of the Srivastava–Miller model (2) to a multiparameter and multivariable Voigt function was considered by Pathan et al. [7] (p. 253, Eq. (2.4)) and reconsidered by Srivastava and Pogány viz. [8] (p. 195, Eq. (4)), as follows:
V μ , ν ( x , y ) = j = 1 r x j 2 0 t μ e y t t 2 4 j = 1 r J ν j ( x j t ) · d t , r N ,
where μ > 0 ; x j , y R + ; ( μ + | ν | ) > 1 , u = ( u 1 , , u r ) , u { x , ν } , and | ν | = j = 1 r ν j . In turn, when assuming that r = 1 , we arrive at (2). We point out that Srivastava and Pogány, in [8], concentrate on bounding inequalities upon V μ , ν ( x , y ) .
For various other investigations involving the Voigt functions, the interested reader may be referred to several recent papers on the subject, such as [6,7,9,10,11,12,13,14,15] and the references cited therein.
We organize this paper as follows: (i) we firstly introduce the extended Voigt-type function V μ , ν ( x , y ) and its counterpart extension W μ , ν ( x , y ) , involving the Neumann function Y ν in the kernel of the representing integral; (ii) following an approach used in [8], we define the multiparameter and multivariable versions V μ , ν ( r ) ( x , y ) and W μ , ν ( r ) ( x , y ) for a positive integer r 2 ; (iii) we obtain computable series expansions for discussed Voigt-type functions in terms of the Humbert confluent hypergeometric functions Ψ 2 ( r ) of r variables; then, (iv) by transforming the input extended Voigt-type functions by the Grünwald–Letnikov fractional derivative, we establish representation formulae in terms of the associated Legendre functions of the second kind Q η ν in the case r = 1 ; finally, (v) functional bounding inequalities are given for V μ , ν ( x , y ) and W μ , ν ( x , y ) in two ways: the application of already known upper bounds for the real argument Neumann function Y ν ( x ) (see [13,16,17]) and by applying newly derived upper bounds for the Struve function H ν ( x ) . The article ends with a discussion, closing remarks, and an exhaustive references list.

2. Extended Voigt-Type Functions V , W , V ( r ) , and W ( r )

We recall that the Bessel function of the second kind (Neumann function or Weber–Bessel function) of order ν is expressible in terms of the Bessel function of the first kind as follows [18] (p. 64, Eq. (1)):
Y ν ( z ) = cos ( ν π ) J ν ( z ) J ν ( z ) sin ( ν π ) = cot ( ν π ) J ν ( z ) csc ( ν π ) J ν ( z ) , ν Z .
Also, Bessel functions of half-integer order have a connection or recurrence formula, as in [19] (p. 925, Eq. (8.465)):
Y n + 1 2 ( z ) = ( 1 ) n 1 J n 1 2 ( z ) and Y n 1 2 ( z ) = ( 1 ) n J n + 1 2 ( z ) .
On the one hand, the following is true [4] (p. 228, Eq. (10.16.1)):
J 1 2 ( z ) = Y 1 2 ( z ) = 2 π z sin ( z ) and J 1 2 ( z ) = Y 1 2 ( z ) = 2 π z cos ( z ) .
On the other hand, the asymptotic of the Neumann function when z [4] (p. 229, Entry 10.17.4) is as follows:
Y ν ( z ) = 2 π z 1 + O ( z 1 ) , | arg ( z ) | π δ , | δ | < π
This suggests extending, in the legitimate way, the Voigt-type function V μ , ν ( x , y ) by replacing the Bessel function J ν ( x t ) in (4) with the Neumann function Y ν ( x t ) of the same order in the following way:
V μ , ν ( x , y ) = x 2 0 t μ e y t t 2 4 Y ν ( x t ) d t , μ , x , y R + ; ( ν ) > 0 , ( μ + ν ) > 1 .
The form of the multiplicative constant is chosen to agree with the pioneering introduction definitions (1) and (2).
Since the Neumann function exhibits asymptotic behavior when z approaches zero, it reads as follows [4] (p. 223, Entries 10.7.4–5):
Y ν ( z ) Γ ( ν ) π 2 z ν ; ( ν ) > 0 or ν = k + 1 2 , k N , Y ν ( z ) cos ( ν π ) Γ ( ν ) π 2 z ν ; ( ν ) > 0 and ν k 1 2 , k N .
Therefore, having in mind (9), we introduce the second variant of the extended Voigt–type function, which reads as follows:
W μ , ν ( x , y ) = π Γ ( ν ) x 2 ν 0 t μ e y t t 2 4 Y ν ( x t ) d t ,
where the parameter range consists of μ , x , y R + ; ( ν ) > 0 , ( μ ν ) > 1 . Furthermore, having in mind (9) and the value Γ ( 1 2 ) = π , the constant is adapted to the multiplicative constants in (4)–(6).
Moreover, considering the generalized multiparameter Srivastava–Pogány variant (6), we also introduce the multiparameter version of the extended Voigt-type function associated with V μ , ν ( x , y ) and W μ , ν ( x , y ) , as follows:
V μ , ν ( r ) ( x , y ) = j = 1 r x j 2 · 0 t μ e y t t 2 4 j = 1 r Y ν j ( x j t ) d t ,
equipped with the parametric space μ , x j , y > 0 ; ( ν j ) > 0 , ( μ + | ν | ) > 1 , where | ν | = j = 1 r ν j . Obviously, V μ , ν ( 1 ) ( x , y ) V μ , ν 1 ( x 1 , y ) covers (8).
Next, we define simultaneously a second multiparameter extension of the Voigt–type function, now associated with W μ , ν ( x , y ) in the following form:
W μ , ν ( r ) ( x , y ) = ( 1 ) r j = 1 r π Γ ( ν j ) x j 2 ν j · 0 t μ e y t t 2 4 j = 1 r Y ν j ( x j t ) d t ,
where μ , x j , y > 0 ; ( ν j ) > 0 , ( μ | ν | ) > 1 ; here W μ , ν ( 1 ) ( x , y ) W μ , ν 1 ( x 1 , y ) . Therefore,
K ( x , y ) = V 1 2 , 1 2 ( x , y ) and L ( x , y ) = V 1 2 , 1 2 ( x , y ) .
Using the following recurrence formula [18] (p. 66, Eq. (1))
Y ν 1 ( z ) + Y ν + 1 ( z ) = 2 ν z Y ν ( z ) ,
we obtained the following recurrence formulae for extended Voigt-type functions:
2 ν x V μ , ν ( x , y ) = V μ + 1 , ν 1 ( x , y ) + V μ + 1 , ν + 1 ( x , y ) .
For the next results, we require the definition of the Humbert confluent hypergeometric functions of n variables, the series definition of which reads as follows [20] (p. 34, Eq. (9)):
Ψ 2 ( n ) [ α ; γ 1 , , γ n ; x 1 , , x n ] = m , n 0 ( α ) m 1 + + m n ( γ 1 ) m 1 ( γ n ) m n x 1 m 1 m 1 ! x n m n m n ! , max { | x 1 | , , | x n | } < .
Clearly, the confluent Humbert function of two variables Ψ 2 ( 2 ) Ψ 2 ; please see [20] (p. 26, Eq. (22)).
Theorem 1.
Let μ , y R + , x R ; | ν | < 1 and ( μ + ν ) > 1 . Then, the following representation holds true:
V μ , ν ( x , y ) = 2 μ 1 2 x ν + 1 2 cot ( ν π ) Γ ( ν + 1 ) { Γ μ + ν + 1 2 Ψ 2 μ + ν + 1 2 ; ν + 1 , 1 2 ; x 2 , y 2 2 y Γ μ + ν 2 + 1 Ψ 2 μ + ν 2 + 1 ; ν + 1 , 3 2 ; x 2 , y 2 } 2 μ 2 ν 1 2 x 1 2 ν csc ( ν π ) Γ ( 1 ν ) { Γ μ ν + 1 2 Ψ 2 μ ν 2 + 1 ; 1 ν , 1 2 ; x 2 , y 2 2 y Γ μ ν + 2 2 Ψ 2 μ ν 2 + 1 ; 1 ν , 3 2 ; x 2 , y 2 } .
Moreover, for all μ , y R + , x R , we have the following:
W μ , ν ( x , y ) = π Γ ( ν ) x 2 ν 1 2 V μ , ν ( x , y ) ,
provided that ν ( 1 , 0 ) ( 0 , 1 ) and ( μ | ν | ) > 1 .
Proof. 
Using the representation of Neumann function (7) in the definition (8), we conclude that the following holds true:
V μ , ν ( x , y ) = x 2 0 t μ e y t t 2 4 cot ( ν π ) J ν ( x t ) csc ( ν π ) J ν ( x t ) d t = x 2 0 t μ e y t t 2 4 cot ( ν π ) m 0 ( 1 ) m x t 2 2 m + ν m ! Γ ( m + ν + 1 ) csc ( ν π ) m 0 ( 1 ) m x t 2 2 m ν m ! Γ ( m ν + 1 ) d t .
The Maclaurin series expansion of e y t enables the legitimate term-wise integration of the resulting absolutely convergent double series. Hence, the following also holds true:
V μ , ν ( x , y ) = 2 μ 1 2 x ν + 1 2 cot ( ν π ) m , n 0 ( x 2 ) m ( 2 y ) n m ! n ! Γ ( m + ν + 1 ) Γ μ + ν + 1 2 + m + n 2 2 μ 1 2 x 1 2 ν csc ( ν π ) m , n 0 ( x 2 ) m ( 2 y ) n m ! n ! Γ ( m ν + 1 ) Γ μ ν + 1 2 + m + n 2 .
Now, by separating the n series into even and odd terms, we obtain the following:
V μ , ν ( x , y ) = 2 μ 1 2 x ν + 1 2 cot ( ν π ) Γ ( ν + 1 ) { Γ μ + ν + 1 2 m , n 0 ( 1 2 μ + ν + 1 ) m + n ( x 2 ) m ( y 2 ) n ( ν + 1 ) m ( 1 2 ) n m ! n ! 2 y Γ μ + ν 2 + 1 m , n 0 ( μ + ν 2 + 1 ) m + n ( x 2 ) m ( y 2 ) n ( ν + 1 ) m ( 3 2 ) n m ! n ! } 2 μ 2 ν 1 2 x 1 2 ν csc ( ν π ) Γ ( 1 ν ) { Γ μ ν + 1 2 m , n 0 ( 1 2 μ ν + 1 ) m + n ( x 2 ) m ( y 2 ) n ( ν + 1 ) m ( 1 2 ) n m ! n ! 2 y Γ μ ν 2 + 1 m , n 0 ( μ ν 2 + 1 ) m + n ( x 2 ) m ( y 2 ) n ( ν + 1 ) m ( 3 2 ) n m ! n ! } .
By using the definition of Humbert confluent hypergeometric function of two variables (12), we obtain the asserted expression. □
Remark 1.
Clearly, if we put μ = 1 / 2 , ν = 1 / 2 and μ = 1 / 2 , ν = 1 / 2 into Theorem 1, we obtain the known corrected representation by Srivastava and Miller [5] (pp. 113–114, Eq. (13–14)) and also derived by Exton [21] (p. L76, Eqs. (8–9)), respectively:
K ( x , y ) = V 1 2 , 1 2 ( x , y ) = Ψ 2 1 2 ; 1 2 , 1 2 ; x 2 , y 2 2 y π Ψ 2 1 ; 1 2 , 3 2 ; x 2 , y 2 , L ( x , y ) = V 1 2 , 1 2 ( x , y ) = 2 x π Ψ 2 1 ; 3 2 , 1 2 ; x 2 , y 2 2 x y Ψ 2 3 2 ; 3 2 , 3 2 ; x 2 , y 2 .
The next result is devoted to the multiparameter extension for the multivariable Voigt-type functions.
Theorem 2.
For all μ > 0 , ν ( 1 , 1 ) r { 0 } , such that μ + | ν | > 1 and x , y > 0 , we have the following:
V μ , ν ( r ) ( x , y ) = 2 μ 1 2 j = 1 r x j ν j + 1 2 cot ( ν j π ) Γ ( ν j + 1 ) · { Γ μ + | ν | + 1 2 Ψ 2 ( r ) μ + | ν | + 1 2 ; ν + 1 , 1 2 ; x 2 , y 2 2 y Γ μ + | ν | 2 + 1 Ψ 2 ( r ) μ + | ν | 2 + 1 ; ν + 1 , 3 2 ; x 2 , y 2 } 2 μ 1 2 j = 1 r x j 1 2 ν j csc ( ν j π ) Γ ( 1 ν j ) · { Γ μ | ν | + 1 2 Ψ 2 ( r ) μ | ν | + 1 2 ; 1 ν , 1 2 ; x 2 , y 2 2 y Γ μ | ν | 2 + 1 Ψ 2 ( r ) μ | ν | 2 + 1 ; 1 ν , 3 2 ; x 2 , y 2 } ,
where the shorthands are as follows: 1 ± ν = ( 1 ± ν 1 , , 1 ± ν r ) and x 2 = ( x 1 2 , , x r 2 ) .
Moreover, for μ | ν | > 1 , we have the following:
W μ , ν ( r ) ( x , y ) = ( 1 ) r j = 1 r π Γ ( ν j ) x j 2 ν j 1 2 · V μ , ν ( r ) ( x , y ) .
The proof is the copy of the proving procedure of Theorem 1; therefore, it is omitted.
Next, we expose another fashion computable series representation of V μ , ν ( x , y ) , for which we need the definition of the associated Legendre function of the second kind [4] (p. 354, Entry 14.3.7), which reads as follows:
Q ν μ ( x ) = e μ π i π Γ ( ν + μ + 1 ) ( x 2 1 ) μ 2 2 ν + 1 Γ ( ν + 3 2 ) x ν + μ + 1 F 1 2 ν + μ 2 + 1 , ν + μ + 1 2 ; ν + 3 2 ; 1 x 2 ,
where ( ν + μ ) N and x > 1 . Furthermore, we recall the definition of the Grünwald–Letnikov (GL) fractional derivative D x ν [ f ] of the order ρ , with respect to x of a suitable function f [22] (p. 373, Eq. (20.7)), as follows:
D x ρ [ f ] = lim h 0 1 h ρ n 0 ( 1 ) n ρ n f x + ( ρ n ) h ,
where h 0 means that, when approaching zero, h remains positive. The GL fractional derivative of order ρ of the exponential function is as follows [23] (p. 2510, Eq. (29)):
D x ρ e α x = α ρ e α x .
Here and in what follows we frequently use Laplace transform notation, as follows:
L p [ f ( t ) ] = F ( p ) = 0 e p t f ( t ) d t ,
which is of a suitable input function f ( t ) .
Theorem 3.
For all μ , x 1 , y > 0 and ν ( 0 , 1 ) , we have the following:
V μ , ν ( x , y ) = 2 x π e μ π i Γ ( ν + 1 ) n 0 ( 1 ) n + 1 4 n n ! ( ν + 1 ) 2 n · D y μ 1 ( y 2 + x 2 ) n + 1 2 Q 2 n ν y y 2 + x 2 .
Moreover, for all μ , x 1 , y > 0 and ν ( 0 , 1 ) , the following holds true:
W μ , ν ( x , y ) = 2 ν e μ π i x 2 ν n 0 ( 1 ) n 4 n n ! ( ν + 1 ) 2 n D y μ 1 ( y 2 + x 2 ) n + 1 2 Q 2 n ν y y 2 + x 2 .
Proof. 
We point out that the extended Voigt-type function V μ , ν ( x , y ) is the Laplace transform, with respect to the parameter y > 0 , viz. as follows:
V μ , ν ( x , y ) = x 2 L y t μ e t 2 4 Y ν ( x t ) .
Applying (14) to the Laplace kernel μ times, we obtain the following:
V μ , ν ( x , y ) = x 2 0 ( 1 ) μ D y μ e y t e t 2 4 Y ν ( x t ) d t
= x 2 ( 1 ) μ D y μ 0 e y t t 2 4 Y ν ( x t ) d t
= x 2 ( 1 ) μ n 0 ( 1 ) n 4 n n ! D y μ 0 e y t t 2 n Y ν ( x t ) d t .
The Laplace transform, adapted to our setting above, reads as follows [24] (p. 298, Eq. (3)):
L p t μ Y ν ( a t ) = 2 Γ ( μ + ν + 1 ) π ( p 2 + a 2 ) μ + 1 2 Q μ ν p p 2 + a 2 , ( μ ) > | ( ν ) | 1 , ( p ) > | ( a ) | ,
We transform (16), using p = y , μ = 2 n , and a = x , into the following:
V μ , ν ( x , y ) = 2 x π ( 1 ) 1 μ n 0 ( 1 ) n Γ ( 2 n + ν + 1 ) 4 n n ! D y μ 1 ( y 2 + x 2 ) n + 1 2 Q 2 n ν y y 2 + x 2 .
According to the definition of the associated Legendre function of the second kind, the parameter space excludes ( 2 n ν ) from the positive integers, and, as the Laplace transform holds for min n 0 ( 2 n ν + 1 ) = 1 ν > 0 , the first statement is proven. As for the claim (15), it is enough to consider (13). The rest is obvious. □
Another approach entails using the Grünwald–Letnikov derivative (14) applied to the newly introduced exponential term exp { β t 2 } , which replaces exp { t 2 / 4 } in the integrand of V μ , ν ( x , y ) . The related results, which follow, contain the Kummer confluent hypergeometric function F 1 1 instead of the previously used associated Legendre function of the second kind.
Theorem 4.
For all μ , x , y > 0 and ν ( 1 , 0 ) ( 0 , 1 ) , we have the following:
V μ , ν ( x , y ) = cos ( ν π ) Γ ( ν ) 2 π e μ π 2 i x 2 ν + 1 2 n 0 ( y ) n n ! { Γ 1 + n + ν 2 D β μ 2 F 1 1 n + ν + 1 2 ; 1 + ν ; x 2 4 β β n + ν + 1 2 β = 1 4 + 4 ν Γ ( ν ) Γ ( 1 + n ν 2 ) Γ ( ν ) x 2 ν D β μ 2 F 1 1 n ν + 1 2 ; 1 ν ; x 2 4 β β n ν + 1 2 β = 1 4 } .
Also, for all μ , x , y > 0 , and | ν | < 1 , the following holds true:
W μ , ν ( x , y ) = cos ( ν π ) 2 e μ π 2 i Γ ( ν ) Γ ( ν ) x 2 2 ν n 0 ( y ) n n ! { Γ 1 + n + ν 2 D β μ 2 F 1 1 n + ν + 1 2 ; 1 + ν ; x 2 4 β β n + ν + 1 2 β = 1 4 + 4 ν Γ ( ν ) Γ ( 1 + n ν 2 ) Γ ( ν ) x 2 ν D β μ 2 F 1 1 n ν + 1 2 ; 1 ν ; x 2 4 β β n ν + 1 2 β = 1 4 } .
Proof. 
Consider the Laplace transform, which reads as follows:
L y t μ e t 2 4 Y ν ( x t ) = e μ π 2 i 0 e y t D β μ 2 e β t 2 β = 1 4 Y ν ( x t ) d t = e μ π 2 i D β μ 2 0 e y t e β t 2 Y ν ( x t ) d t β = 1 4 .
By expanding the Laplace kernel exp { y t } into a Maclaurin series, interchanging the summation and integration, and using the substitution t = s , we obtain the following:
0 t μ e y t t 2 4 Y ν ( x t ) d t = 1 2 n 0 ( y ) n n ! 0 e β s s n 1 2 Y ν ( x s ) d s = 1 2 n 0 ( y ) n n ! L β s μ Y ν ( x s ) .
In turn, the Laplace transform reads as follows [24] (p. 300, Eq. 3.13.2.4.):
L p s μ Y ν ( a s ) = a ν cos ( ν π ) Γ ( ν ) Γ ( μ + ν 2 + 1 ) π 2 ν p μ + ν 2 + 1 F 1 1 μ + ν 2 + 1 ; 1 + ν ; a 2 4 p 2 ν cos ( ν π ) Γ ( ν ) Γ ( μ ν 2 + 1 ) π a ν p μ ν 2 + 1 F 1 1 μ ν 2 + 1 ; 1 ν ; a 2 4 p ,
provided 2 | ( μ ) | > | ( ν ) | 2 and ( β ) > 0 . For 2 μ = n 1 , p = β , a = x , we conclude that the following holds true:
0 t μ e y t t 2 4 Y ν ( x t ) d t = 1 2 π n 0 ( y ) n n ! { x ν cos ( ν π ) Γ ( ν ) Γ ( n + ν + 1 2 ) 2 ν β n + ν + 1 2 F 1 1 n + ν + 1 2 ; 1 + ν ; x 2 4 β + 2 ν cos ( ν π ) Γ ( ν ) Γ ( n ν + 1 2 ) x ν β n ν + 1 2 F 1 1 n ν + 1 2 ; 1 ν ; x 2 4 β } .
By now inserting (20) into (19), we complete the proving procedure of the statement.
The second statement (18) is the immediate consequence of relation (13). □

3. Functional Bounding Inequalities

First, we recall the integral definition of the Struve function, which reads as follows [18] (p. 328, Eq. (1)):
H ν ( z ) = 2 z 2 ν π Γ ( ν + 1 2 ) 0 1 ( 1 t 2 ) ν 1 2 sin ( z t ) d t , 2 ( ν ) > 1 .
It is worth mentioning that Struve considered only H 0 and H 1 . Now, the Gubler–Weber formula reads as follows: [18] (p. 165, Eq. (5))
Y ν ( z ) = 2 z 2 ν π Γ ( ν + 1 2 ) 0 1 ( 1 t 2 ) ν 1 2 sin ( z t ) d t 0 ( 1 + t 2 ) ν 1 2 e z t d t
which implies the following:
Y ν ( z ) = H ν ( z ) 2 z 2 ν π Γ ( ν + 1 2 ) 0 ( 1 + t 2 ) ν 1 2 e z t d t ,
where both representation formulae hold for ( z ) > 0 and ν > 1 / 2 , due to the convergence control of the second integral.
By splitting the ν domain into the disjoint intervals, we obtain the following:
1 2 , = 1 2 , 1 2 1 2 , 3 2 3 2 , = : U 1 U 2 U 3 ,
Baricz et al. [16] (pp. 957–958) derive the functional bounds for Y ν ( x ) as follows (see also [25] (pp. 7–8) and the monograph [17]):
Y ν ( x ) + x ν 2 ν Γ ( ν + 1 ) x 2 ν 1 π Γ ( ν + 1 2 ) , 1 2 < ν 1 2 , x 2 ν 1 π Γ ( ν + 1 2 ) + 2 ν Γ ( ν ) π x ν , 1 2 < ν 3 2 , x ν 1 2 π Γ ( ν + 1 2 ) + 2 2 ν 3 2 Γ ( ν ) π x ν , ν > 3 2
where x > 0 . To present our recent findings, we need the definition of the parabolic (or Weber) cylinder function D a ( x ) (Whittaker’s notation), which is a principal solution of the Weber differential equation [4] (p. 303, Entry 12.2.4), which reads as follows:
d 2 w d z 2 + ν + 1 2 1 4 z 2 w = 0 ; w = D ν ( z ) .
The solution is entire function of z and ν . Furthermore, in the derivation procedure, we apply the integral representation for parabolic cylinder function, as follows [19] (p. 503, Eq. 3.953.2):
D μ γ 2 β = ( 2 β ) μ 2 Γ ( μ ) e γ 2 8 β 0 t μ 1 e γ t β t 2 d t , μ , γ , β C ; min { ( β ) , ( μ ) } > 0 .
Theorem 5.
For all μ , ν + 1 2 > 0 , with ( μ ) > | ( ν ) | > 0 and for all positive x , y , we have the following:
V μ , ν ( x , y ) κ 1 ( μ , ν ) D ( μ + ν ) ( 2 y ) + κ 2 ( μ , ν ) D ( μ + ν 1 ) ( 2 y ) , 1 2 < ν 1 2 , κ 1 ( μ , ν ) D ( μ + ν ) ( 2 y ) + κ 2 ( μ , ν ) D ( μ + ν 1 ) ( 2 y ) + κ 3 ( μ , ν ) D ( μ ν ) ( 2 y ) , 1 2 < ν 3 2 , κ 1 ( μ , ν ) D ( μ + ν ) ( 2 y ) + κ 4 ( μ , ν ) D ( μ + ν 1 ) ( 2 y ) + κ 5 ( μ , ν ) D ( μ ν ) ( 2 y ) , ν > 3 2
where
κ 1 ( μ , ν ) = x ν + 1 2 e y 2 2 Γ ( μ + ν ) 2 ν μ + 1 2 Γ ( ν + 1 ) , κ 2 ( μ , ν ) = x ν 1 2 e y 2 2 Γ ( μ + ν 1 ) 2 ν μ 2 π Γ ( ν + 1 2 ) , κ 3 ( μ , ν ) = x 1 2 ν e y 2 2 Γ ( ν ) Γ ( μ ν ) 2 1 ν μ 2 π Γ ( ν + 1 ) ,
κ 4 ( μ , ν ) = x ν 1 2 e y 2 2 Γ ( μ + ν 1 ) π 2 3 ν μ 2 Γ ( ν + 1 2 ) , κ 5 ( μ , ν ) = x 1 2 ν e y 2 2 Γ ( ν ) Γ ( μ ν ) 2 4 3 ν μ 2 π Γ ( ν + 1 ) .
Proof. 
Starting with (8), using the decomposition ( 1 2 , ) = U 1 U 2 U 3 and the estimates (24), we conclude as follows:
V μ , ν ( x , y ) x 2 0 t μ 1 e y t t 2 4 Y ν ( x t ) d t V μ , U 1 ( x , y ) , 1 2 < ν 1 2 , V μ , U 2 ( x , y ) , 1 2 < ν 3 2 , V μ , U 3 ( x , y ) , ν > 3 2
where
V μ , U 1 ( x , y ) x 2 x 2 ν Γ ( ν + 1 ) 0 t μ + ν e y t t 2 4 d t + x 2 ν 1 π Γ ( ν + 1 2 ) 0 t μ + ν 1 e y t t 2 4 d t , V μ , U 2 ( x , y ) x 2 { x 2 ν Γ ( ν + 1 ) 0 t μ + ν e y t t 2 4 d t + x 2 ν 1 π Γ ( ν + 1 2 ) 0 t μ + ν 1 e y t t 2 4 d t + 2 ν Γ ( ν ) π x ν Γ ( ν + 1 ) 0 t μ ν e y t t 2 4 d t } ,
and
V μ , U 3 ( x , y ) x 2 { x 2 ν Γ ( ν + 1 ) 0 t μ + ν e y t t 2 4 d t + x ν 1 2 π Γ ( ν + 1 2 ) 0 t μ + ν 1 e y t t 2 4 d t + 2 2 ν 3 2 Γ ( ν ) π x ν Γ ( ν + 1 ) 0 t μ ν e y t t 2 4 d t } ,
which is equivalent to the statement of the theorem. In the derivation procedure, we apply (25) with γ = y and β = 1 / 4 . □
The Fox–Wright generalized hypergeometric function with p numerator parameter couples ( a 1 , A 1 ) , , ( a p , A p ) and q denominator parameter pairs ( b 1 , B 1 ) , , ( b q , B q ) has the following power series definition [26] (pp. 286–287):
Ψ q p ( a 1 , A 1 ) , , ( a p , A p ) ( b 1 , B 1 ) , , ( b q , B q ) | z = Ψ q p ( a p , A p ) ( b q , B q ) | z = n 0 j = 1 p Γ ( a j + n A j ) j = 1 q Γ ( b j + n B j ) z n n ! ,
where A j , B k 0 , j = 1 , , p , and k = 1 , , q . The series (27) converges for all z C when the following is true:
Δ : = 1 + j = 1 q B j k = 1 p A k > 0 .
When Δ = 0 , the series in (27) converges for | z | < and | z | = under the condition ( Ξ ) > 1 / 2 , where the following holds true:
: = i = 1 p A i A i j = 1 q B j B j , Ξ = j = 1 q b j i = 1 p a i + p q 2 .
The function Ψ 0 1 [ · ] is a confluent Fox–Wright function.
Theorem 6.
For all μ , ν + 1 2 > 0 and for all positive x , y > 0 , we have the following:
V μ , ν ( x , y ) λ 1 ( μ , ν ) Ψ 0 1 ( μ + ν + 1 2 , 1 2 )   ¯ | 2 y + λ 2 ( μ , ν ) Ψ 0 1 ( μ + ν 2 , 1 2 )   ¯ | 2 y , 1 2 < ν 1 2 , λ 1 ( μ , ν ) Ψ 0 1 ( μ + ν + 1 2 , 1 2 )   ¯ | 2 y + λ 2 ( μ , ν ) Ψ 0 1 ( μ + ν 2 , 1 2 )   ¯ | 2 y + λ 3 ( μ , ν ) Ψ 0 1 ( μ ν + 1 2 , 1 2 )   ¯ | 2 y , 1 2 < ν 3 2 , λ 1 ( μ , ν ) Ψ 0 1 ( μ + ν + 1 2 , 1 2 )   ¯ | 2 y + λ 4 ( μ , ν ) Ψ 0 1 ( μ + ν 2 , 1 2 )   ¯ | 2 y + λ 5 ( μ , ν ) Ψ 0 1 ( μ ν + 1 2 , 2 2 )   ¯ | 2 y , ν > 3 2
where
λ 1 ( μ , ν ) = x ν + 1 2 2 μ 1 2 Γ ( ν + 1 ) , λ 2 ( μ , ν ) = x ν 1 2 2 μ 1 2 π Γ ( ν + 1 2 ) , λ 3 ( μ , ν ) = x 1 2 ν 2 μ + 1 2 π Γ ( ν + 1 2 ) ,
λ 4 ( μ , ν ) = x ν 1 2 2 μ + ν 2 π Γ ( ν + 1 2 ) , λ 5 ( μ , ν ) = x 1 2 ν 2 μ + ν 2 Γ ( ν ) π Γ ( ν + 1 ) .
Proof. 
In the derivation procedure of Theorem 5, we apply the following integral representation for the Fox–Wright function [8] (p. 198, Equation (15)):
I μ ( y ) = 0 t μ 1 e y t t 2 4 d t = 2 μ 1 Ψ 0 1 ( μ 2 , 1 2 ) ;   ¯ ; 2 y , ( μ ) , ( y ) > 0 .
Using (26), we firstly evaluate the input extended Voigt function V μ , ν ( x , y ) , then obtain the claim of the theorem. □
In the following, we refine the bounds (24) focusing on the representation of Formula (23), which simultaneously gives novel bounds for real parameters and positive variables for either the Neumann function Y ν ( x ) or the Struve function H ν ( x ) .
Theorem 7.
For all ν 1 2 and x > 0 , we have the following:
Y ν ( x ) H ν ( x ) 2 π x , 1 2 ν 3 2 H ν ( x ) x ν 3 ( x 2 + 2 ν 1 ) π 2 ν 1 Γ ( ν + 1 2 ) ν > 3 2 .
Proof. 
For all ν 1 2 and x > 0 , the function H ν ( x ) 0 (see [4] (p. 291, Entry 11.4.13)). Indeed, the case ν = 1 2 is obvious, since the following holds true:
H 1 2 ( x ) = 2 x π 0 1 sin ( x t ) d t = 2 2 π x sin 2 x 2 0 .
In the case ν > 1 2 , consider three consecutive zeros t k = ( π / x ) k of sin ( x t ) , t 2 k , t 2 k + 1 , t 2 k + 2 . Then, having in mind that ( 1 t 2 ) ν 1 2 for t [ 0 , 1 ] , we have the following:
t 2 k t 2 k + 1 sin ( x t ) ( 1 t 2 ) ν 1 2 d t t 2 k + 1 t 2 k + 2 sin ( x t ) ( 1 t 2 ) ν 1 2 d t 0 , k N 0 .
Therefore, a fortiori, the following holds true:
H ν ( x ) = k 0 t 2 k t 2 k + 1 t 2 k + 1 t 2 k + 2 sin ( x t ) ( 1 t 2 ) ν 1 2 d t 0 .
Now, recalling (23), we see that, for 2 ν > 3 , we have, by virtue of the Bernoulli inequality [27] (p. 317), the following:
( 1 + h ) α 1 + α h , h > 1 , α ( 0 , 1 ) ,
and the following estimate:
Y ν ( x ) H ν ( x ) 2 x 2 ν π Γ ( ν + 1 2 ) 0 1 + ν 1 2 t 2 e x t d t = H ν ( x ) x ν 3 ( x 2 + 2 ν 1 ) π 2 ν 1 Γ ( ν + 1 2 ) .
The Bernoulli inequality is reversed for α ( 0 , 1 ) ; therefore, assuming ν [ 1 2 , 3 2 ] , we conclude by the A–G inequality 1 + t 2 2 t , that the following holds true:
Y ν ( x ) H ν ( x ) 2 x ν 2 π Γ ( ν + 1 2 ) 0 t ν 1 2 e x t d t = H ν ( x ) 2 π x .
Therefore, the proof is complete. □
Now it remains to apply (28) in bounding the extended Voigt and allied functions. Before this, we recall the following power series definition:
H ν ( z ) = z 2 ν + 1 n 0 ( 1 ) n z 2 2 n Γ ( n + 3 2 ) Γ ( n + ν + 3 2 ) , 2 ( ν ) > 3 .
Theorem 8.
For all μ , x , y > 0 and ν 1 2 , 3 2 , the following holds true:
V μ , ν ( x , y ) A μ , ν ( x , y ) 2 μ 2 + 1 4 π Γ μ + 1 2 e y 2 2 D ( μ + 1 2 ) 2 y ,
where
A μ , ν ( x , y ) = x ν + 3 2 e y 2 2 2 ν μ + 1 2 n 0 ( 1 ) n x 2 n Γ ( μ + ν + 2 + 2 n ) 2 n Γ ( n + 3 2 ) Γ ( n + ν + 3 2 ) D ( μ + ν + 2 + 2 n ) 2 y .
For ν > 3 2 and μ + ν > 2 , we have the following:
V μ , ν ( x , y ) A μ , ν ( x , y ) 2 μ ν + 1 2 x ν 5 2 e y 2 2 π Γ ( ν + 1 2 ) { Γ ( μ + ν ) x 2 D ( μ + ν ) ( 2 y ) + ( ν 1 2 ) Γ ( μ + ν 2 ) D ( μ + ν 2 ) ( 2 y ) } .
Proof. 
When 1 2 ν 3 2 , the integral definition (8), under upper bound (28) applied to the Neumann function in the kernel, becomes the following:
V μ , ν ( x , y ) = x 2 0 t μ e y t t 2 4 Y ν ( x t ) d t x 2 0 t μ e y t t 2 4 H ν ( x t ) 2 π x t 1 2 d t = x 2 ν + 3 2 n 0 ( 1 ) n x 2 2 n Γ ( n + 3 2 ) Γ ( n + ν + 3 2 ) 0 t μ + ν + 1 + 2 n e y t t 2 4 d t 1 π 0 t μ 1 2 e y t t 2 4 d t = x ν + 3 2 e y 2 2 2 ν μ + 1 2 n 0 ( 1 ) n x 2 n Γ ( μ + ν + 2 + 2 n ) 2 n Γ ( n + 3 2 ) Γ ( n + ν + 3 2 ) D ( μ + ν + 2 + 2 n ) 2 y 2 μ 2 + 1 4 π Γ μ + 1 2 e y 2 2 D ( μ + 1 2 ) 2 y ,
where the parabolic cylinder function formula (25) is employed twice.
Next, for ν > 3 2 , the calculation reduces to the following:
I μ , ν ( x , y ) = x 2 0 t μ e y t t 2 4 ( x t ) ν 3 [ ( x t ) 2 + 2 ν 1 ] π 2 ν 1 Γ ( ν + 1 2 ) d t = x ν 5 2 2 ν 1 2 π Γ ( ν + 1 2 ) x 2 0 t μ + ν 1 e y t t 2 4 d t + ( 2 ν 1 ) 0 t μ + ν 3 e y t t 2 4 d t = 2 μ ν + 1 2 x ν 5 2 e y 2 2 π Γ ( ν + 1 2 ) { Γ ( μ + ν ) x 2 D ( μ + ν ) ( 2 y ) + ( ν 1 2 ) Γ ( μ + ν 2 ) D ( μ + ν 2 ) ( 2 y ) } ,
where we also use the aid of (25), provided μ + ν > 2 . □
The defining relation reads as follows (13):
W μ , ν ( x , y ) = π Γ ( ν ) x 2 ν 1 2 V μ , ν ( x , y ) ,
which, when combined with Theorem 8, readily provides the refined lower bound for the counterpart extended Voigt function.
Corollary 1.
For all μ , x , y > 0 , and ν 1 2 , 3 2 , the following holds true:
W μ , ν ( x , y ) π Γ μ + 1 2 Γ ( ν ) 2 ν μ 2 3 4 x ν 1 2 e y 2 2 D ( μ + 1 2 ) 2 y B μ , ν ( x , y ) ,
where
B μ , ν ( x , y ) = π Γ ( ν ) x 2 ν + 1 e y 2 2 2 3 ν μ 2 n 0 ( 1 ) n x 2 n Γ ( μ + ν + 2 + 2 n ) 2 n Γ ( n + 3 2 ) Γ ( n + ν + 3 2 ) D ( μ + ν + 2 + 2 n ) 2 y .
In turn, for ν > 3 2 and μ + ν > 2 , we have the following:
W μ , ν ( x , y ) π 2 μ 3 ν 2 + 1 x 2 ν 3 e y 2 2 Γ ( ν ) Γ ( ν + 1 2 ) { Γ ( μ + ν ) x 2 D ( μ + ν ) ( 2 y ) + ( ν 1 2 ) Γ ( μ + ν 2 ) D ( μ + ν 2 ) ( 2 y ) } B μ , ν ( x , y ) .
Now, we significantly refine the non-negativity result H ν ( x ) 0 listed in [4] (p. 291, Entry 11.4.13), proving a more efficient functional lower bound for the Struve H ν ( x ) . With the aid of the resulting bound, we establish a lower bound for the Neumann function Y ν ( x ) .
Theorem 9.
For all 2 ( ν ) 1 and x ( 0 , π ] , there holds the following lower bound:
H ν ( x ) x ν sin ( x ) 2 ν π Γ ( ν + 3 2 ) .
Proof. 
Let t 1 = π / x be the first positive zero of sin ( x t ) in the integrand of (21). In the case t 1 1 , the secant line containing the origin and the point T 1 ( 1 , sin ( x ) ) is completely below the sine arc; that is, sin ( x t ) t · sin ( x ) , t [ 0 , 1 ] . Accordingly, taking into account the integral representation which holds true for ( ν ) > 1 2 , we conclude the following:
H ν ( x ) = 2 x 2 ν π Γ ( ν + 1 2 ) 0 1 ( 1 t 2 ) ν 1 2 sin ( x t ) d t 2 x 2 ν sin ( x ) π Γ ( ν + 1 2 ) 0 1 t ( 1 t 2 ) ν 1 2 d t ,
which evidently gives (37). □
Theorem 10.
For all 1 2 < ν < 1 , there holds the following lower bound:
Y ν ( x ) 1 2 ν π x Γ ( ν + 3 2 ) x ν + 1 2 sin ( x ) 2 ν + 1 2 Γ ( ν + 3 2 ) , x [ x 1 , x 2 ] ( 0 , π ] .
Here, x 0 , x 1 are the first two positive solutions of the following equation:
x ν + 1 2 sin ( x ) = 2 ν + 1 2 Γ ( ν + 3 2 ) .
Proof. 
Letting ν ( 1 2 , 1 ) , we apply the lower bound (37) to the relation (23) as follows:
Y ν ( x ) = H ν ( x ) 2 x 2 ν π Γ ( ν + 1 2 ) 0 ( 1 + t 2 ) ν 1 2 e x t d t x ν sin ( x ) 2 ν π Γ ( ν + 3 2 ) 2 x 2 ν π Γ ( ν + 1 2 ) 0 ( 1 + t 2 ) ν 1 2 e x t d t .
By virtue of the A–G inequality ( 1 + t 2 ) ν 1 2 ( 2 t ) ν 1 2 , | ν | < 1 2 , we clearly conclude that the estimate reads as follows:
Y ν ( x ) 1 2 ν π x Γ ( ν + 3 2 ) x ν + 1 2 sin ( x ) 2 ν + 1 2 Γ ( ν + 3 2 ) ,
which is definitely non-negative on the interval [ x 0 , x 1 ] . □

4. Discussion and Further Comments

A . The formula applied for the half-integer ν = n + 1 2 is [3] (p. 3, Eq. (2)), as follows:
Γ n + 1 2 = ( 2 n 1 ) ! ! 2 n π , n N 0 .
Specifying all parameters ν , ν as half-integers in (10) and (11), we infer the following expressions:
W μ , ν ( x , y ) = 1 2 x n + 1 2 ( 2 n 1 ) ! ! 0 t μ e y t t 2 4 Y n + 1 2 ( x t ) d t , W μ , ν ( r ) ( x , y ) = ( 1 ) r 2 r 2 j = 1 r x j n j + 1 2 ( 2 n j 1 ) ! ! · 0 t μ e y t t 2 4 j = 1 r Y n j + 1 2 ( x j t ) d t .
B . We point out the integral expression [19] (p. 1018, Eq. 9.240), which reads as follows:
0 t p 1 e y t t 2 4 d t = 2 p 1 Γ p 2 F 1 1 p 2 ; 1 2 ; y 2 2 ( p 1 ) y Γ p + 1 2 F 1 1 p + 1 2 ; 3 2 ; y 2 .
which holds for all p, ( p ) > 0 . This formula connects the parabolic cylinder function D η with the Kummer confluent hypergeometric 1 F 1 . Applying this expression for both integrals in (34), we can infer other computational series expansion formulae for the functional upper bound in Theorem 8.
C . Theorem 9 implies the questions of sharpness and possible argument and parameter extensions in (37); therefore, so does (38). Moreover, we point out that the domain [ x 1 , x 2 ] in Theorem 10 is not empty; see, for instance, the case ν = 0 , to which the numerical values x 1 1.57079633 and x 2 2.09439510 are associated, while the third root of the related equation x sin ( x ) = 2 Γ 3 2 becomes x 3 > 6.785 .
D . Future research can be focused on exhausting the Bessel–Struve functions class, which also contains the modified Bessel functions of the first and second kind, modified Struve functions of the first and second kind, the Bessel–Struve kernel function, and their various generalizations (see, e.g., the article [28] and the relevant references therein) in unifying the multiparameter and/or multivariable Voigt functions. However, we leave these study topics to another addressees.

5. Conclusions

The classical Voigt functions K ( x , y ) and L ( x , y ) are involved in several physics problems, such as astrophysical spectroscopy, emissions, the absorption and transfer of radiation in an atmosphere, plasma dispersion, and in the theory of neutron reactions, for instance. These functions were introduced by Voigt himself in 1889 and Reiche in 1913, respectively, in [1,2]. Their structure, integral forms, and the fact that both are Laplace transformations of exp t 2 4 cos ( x t ) exp t 2 4 sin ( x t ) or Fourier cosine and/or sine transforms (up to a multiplicative constant) of input functions exp y t t 2 4 —motivated the authors either to study more general oscillatory special functions, the special cases of which are the cosine and sine, or to also consider their finite products. These studies also result in multiparameter extensions; please consult the thorough introduction and section A . in the discussion section.
Our results are grouped in the following way: in Section 2. we introduce the extended Voigt-type function V μ , ν ( x , y ) and its counterpart function W μ , ν ( x , y ) using the Neumann function Y ν multiplied with the power function, obtaining t μ Y ν ( x t ) in the kernel of the observed integral. Additionally, products from kernels of r 2 such functions are also used to define the multiparameter variants V μ , ν ( r ) ( x , y ) , and W μ , ν ( r ) ( x , y ) are introduced. In Theorem 1, we establish computable series representations of V μ , ν ( x , y ) and W μ , ν ( x , y ) in terms of Humbert’s confluent hypergeometric function of Ψ 2 of two variables, while Theorem 2 is devoted to the same fashion representations in terms of Humbert’s confluent hypergeometric Ψ 2 ( r ) of r + 1 variables. Theorem 3 expresses V μ , ν ( x , y ) and W μ , ν ( x , y ) by virtue of the Grünwald–Letnikov (GL) fractional derivative of the Legendre function of the second kind, and Theorem 4 describes the GL-derivative representation of Voigt-type functions applied to the Kummer confluent hypergeometric function 1 F 1 .
Section 3 contains functional bounding inequalities for the Voigt-type functions V , W using the Gubler–Weber integral Formula (22). Theorem 5 contains the upper bound, consisting of a linear combination of parabolic cylinder function D μ , while the linear combination of Fox–Wright generalized confluent hypergeometric 1 P s i 0 functions build another form functional upper bound; see Theorem 6. Theorems 7, 8, 10, and the Corollary 1 are devoted to obtain functional upper or lower bounds for the Neumann Y ν ; finally, Theorem 9 is nothing other than a generalization of the famous Jordan inequality for the real variable Struve function H ν . The application of already known upper bounds for the real argument Neumann function Y ν ( x ) (see [13,16,17]) and applying the here-derived upper bound for the Struve function H ν ( x ) in a new way, together with introducing the GL fractional derivative in the investigation methodology, provide constructive contributions for the study of Voigt-type and allied special functions.

Author Contributions

Conceptualization, T.K.P.; Methodology, R.K.P. and T.K.P.; Validation, U.S.; Investigation, R.K.P., T.K.P. and U.S.; Writing—review & editing, R.K.P. and T.K.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The study did not report any data.

Acknowledgments

The authors highly appreciate the corrections and constructive suggestions made by all four anonymous referees, which strongly improve the completeness of the exposition herein and the readability and clear recognition of the presented ideas and results of the study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Voigt, W. Zur Theorie der Beugung ebener inhomogener Wellen an einem geradlinig begrentzen unendlichen und absolut schwarzen Schirm. Gött. Nachr. 1889, 1, 1–33. [Google Scholar]
  2. Reiche, F. Über die Emission, Absorption und Intesitätsverteilung von Spektrallinien. Ber. Deutsch. Phys. Ges. 1913, 15, 3–21. [Google Scholar]
  3. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Higher Transcendental Functions; McGraw–Hill Book Company: New York, NY, USA; Toronto, ON, Canada; London, UK, 1953; Volume I. [Google Scholar]
  4. Olver, F.W.J.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. (Eds.) NIST Handbook of Mathematical Functions; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  5. Srivastava, H.M.; Miller, E.A. A unified presentation of the Voigt functions. Astrophys. Space Sci. 1987, 135, 111–118. [Google Scholar] [CrossRef]
  6. Klusch, D. Astrophysical spectroscopy and neutron reactions: Integral transforms and Voigt functions. Astrophys. Space Sci. 1991, 175, 229–240. [Google Scholar] [CrossRef]
  7. Pathan, M.A.; Kamarujjama, M.; Alam, M.K. On multiindices and multivariables presentation of the Voigt functions. J. Comput. Appl. Math. 2003, 160, 251–257. [Google Scholar] [CrossRef]
  8. Srivastava, H.M.; Pogány, T.K. Inequalities for a unified Voigt functions in several variables. Russian J. Math. Phys. 2007, 14, 194–200. [Google Scholar] [CrossRef]
  9. Khan, S.; Agrawal, B.; Pathan, M.A. Some connections between generalized Voigt functions with the different parameters. Appl. Math. Comput. 2006, 181, 57–64. [Google Scholar] [CrossRef]
  10. Parmar, R.K. Bounding inequalities for the generalized Voigt function. J. Anal. 2020, 28, 191–197. [Google Scholar] [CrossRef]
  11. Parmar, R.K.; Saravanan, S. Extended generalized Voigt–type functions and related bounds. J. Classical Anal. 2023, 21, 45–56. [Google Scholar] [CrossRef]
  12. Pathan, M.A.; Shahwan, M.J.S. New representations of the Voigt functions. Demonstr. Math. 2006, 39, 75–80. [Google Scholar] [CrossRef]
  13. Pogány, T.K.; Süli, E. Integral representation for Neumann series of Bessel functions. Proc. Amer. Math. Soc. 2009, 137, 2363–2368. [Google Scholar] [CrossRef]
  14. Srivastava, H.M.; Pathan, M.A.; Kamarajjuma, M. Some unified presentations of the generalized Voigt functions. Commun. Appl. Anal. 1998, 2, 49–64. [Google Scholar]
  15. Srivastava, H.M.; Chen, M.P. Some unified presentations of the Voigt functions. Astrophys. Space Sci. 1992, 192, 63–74. [Google Scholar] [CrossRef]
  16. Baricz, Á.; Jankov Maširević, D.; Pogány, T.K. Integral representations for Neumann–type series of Bessel functions Iν, Yν and Kν. Proc. Amer. Math. Soc. 2012, 140, 951–960. [Google Scholar] [CrossRef]
  17. Baricz, Á.; Jankov Maširević, D.; Pogány, T.K. Series of Bessel and Kummer-Type Functions; Lecture Notes in Math; Springer: Cham, Switzerland, 2017; p. 2207. [Google Scholar]
  18. Watson, G.N. A Treatise on the Theory of Bessel Functions; University Press: Cambridge, UK, 1992. [Google Scholar]
  19. Gradshteyn, I.S.; Ryzhik, I.M. Tables of Integrals, Series, and Products, 7th ed.; Corrected and Enlarged Edition Prepared by A. Jeffrey; Academic Press: New York, NY, USA, 2007. [Google Scholar]
  20. Srivastava, H.M.; Karlsson, P.W. Multiple Gaussian Hypergeometric Series; Halsted Press: Ultimo, Australia; Ellis Horwood Limited: Chichester, UK; John Wiley and Sons: New York, NY, USA; Chichester, UK; Brisbane, Australia; Toronto, ON, Canada, 1985. [Google Scholar]
  21. Exton, H. On the reducibility of the Voigt functions. J. Phys. A. Math. Gen. 1981, 14, L75–L77. [Google Scholar] [CrossRef]
  22. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: New York, NY, USA, 1993. [Google Scholar]
  23. Ortigueira, M.D. A coherent approach to non–integer order derivatives. Signal Process. 2006, 86, 2505–2515. [Google Scholar] [CrossRef]
  24. Prudnikov, A.P.; Brychkov, Y.A.; Marichev, O.I. Integrals and Series. Direct Laplace Trandsforms; Gordon and Breach Science Publishers: New York, NY, USA, 1992; Volume 4. [Google Scholar]
  25. Jankov Maširević, D.; Pogány, T.K. Second type Neumann series of generalized Nicholson function. Results Math. 2020, 75, 14. [Google Scholar] [CrossRef]
  26. Wright, E.M. The asymptotic expansion of the generalized hypergeometric function. J. London. Math. Soc. 1935, 10, 286–293. [Google Scholar] [CrossRef]
  27. Vîjîitu, H. 103.19 Bernoulli’s inequality for negative exponents. Math. Gaz. 2019, 103, 316–317. [Google Scholar] [CrossRef]
  28. Parmar, R.K.; Pogány, T.K. On (p, q)-extension of further members of Bessel–Struve functions class. Miskolc Math. Notes 2019, 20, 451–463. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Parmar, R.K.; Pogány, T.K.; Sabu, U. On Voigt-Type Functions Extended by Neumann Function in Kernels and Their Bounding Inequalities. Axioms 2024, 13, 534. https://doi.org/10.3390/axioms13080534

AMA Style

Parmar RK, Pogány TK, Sabu U. On Voigt-Type Functions Extended by Neumann Function in Kernels and Their Bounding Inequalities. Axioms. 2024; 13(8):534. https://doi.org/10.3390/axioms13080534

Chicago/Turabian Style

Parmar, Rakesh K., Tibor K. Pogány, and Uthara Sabu. 2024. "On Voigt-Type Functions Extended by Neumann Function in Kernels and Their Bounding Inequalities" Axioms 13, no. 8: 534. https://doi.org/10.3390/axioms13080534

APA Style

Parmar, R. K., Pogány, T. K., & Sabu, U. (2024). On Voigt-Type Functions Extended by Neumann Function in Kernels and Their Bounding Inequalities. Axioms, 13(8), 534. https://doi.org/10.3390/axioms13080534

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop