Next Article in Journal
Nonexpansiveness and Fractal Maps in Hilbert Spaces
Next Article in Special Issue
Geometric Features of the Hurwitz–Lerch Zeta Type Function Based on Differential Subordination Method
Previous Article in Journal
Study on the Cascading Failure Robustness of the Belt and Road Land–Sea Transport Network under Emergencies
Previous Article in Special Issue
Several Symmetric Identities of the Generalized Degenerate Fubini Polynomials by the Fermionic p-Adic Integral on Zp
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Inverse of the Linearization Coefficients of Bessel Polynomials

by
Mohamed Jalel Atia
1,2
1
Department of Mathematics, College of Science, Qassim University, Buraidah 51452, Saudi Arabia
2
Laboratory LR17ES11, Gabes University, Gabès 6072, Tunisia
Symmetry 2024, 16(6), 737; https://doi.org/10.3390/sym16060737
Submission received: 13 May 2024 / Revised: 6 June 2024 / Accepted: 10 June 2024 / Published: 13 June 2024

Abstract

:
In this contribution, we first present a new recursion relation fulfilled by the linearization coefficients of Bessel polynomials (LCBPs), which is different than the one presented by Berg and Vignat in 2008. We will explain why this new recursion formula is as important as Berg and Vignat’s. We introduce the matrix linearization coefficients of Bessel polynomials (MLCBPs), and we present some new results and some conjectures on these matrices. Second, we present the inverse of the connection coefficients with an application involving the modified Bessel function of the second kind. Finally, we introduce the inverse of the matrix of the linearization coefficients of the Bessel polynomials (IMLCBPs), and we present some open problems related to these IMLCBPs.

1. Introduction

In 2008, Berg and Vignat presented a recursion relation fulfilled by the linearization coefficients of Bessel polynomials (LCBPs), which enabled them to prove the positivity of these coefficients. In 2013, they used the Student t-density to present an alternative evaluation of the Boros–Moll integral. Atia and Zeng presented, in 2012, an explicit single-sum formula for these LCBPs, which was missing in Berg and Vignat’s paper. Atia and Zeng derived some immediate consequences from this explicit single sum and generalized a formula on the integral evaluation from Boros and Moll, from which they derived the formula presented by Berg and Vignat. BenAbdallah and Atia discussed, in 2018, a more general aspect of the LCBPs and presented an evaluation, via the convolution of Student t-densities, of a double integral. This is the expository part of this paper.
The second part of this contribution is detailed in the abstract. Let us start by recalling Bessel polynomials. Bessel polynomials q n of degree n are defined by
q n ( u ) = k = 0 n ( n ) k 2 k ( 2 n ) k k ! u k ,
where we use the Pochhammer symbol ( z ) n : = z ( z + 1 ) ( z + n 1 ) for z C and n N . The first values are
q 0 ( u ) = 1 , q 1 ( u ) = 1 + u , q 2 ( u ) = 1 + u + u 2 3 .
Using hypergeometric functions, we have q n ( u ) = F 1 1 ( n ; 2 n ; 2 u ) . They are normalized according to q n ( 0 ) = 1 , and thus differ from the monic normalization θ n ( u ) in [1]:
θ n ( u ) = ( 2 n ) ! n ! 2 n q n ( u ) .
Polynomials θ n are sometimes called reverse Bessel polynomials, and y n ( u ) = u n θ n ( 1 u ) , ordinary Bessel polynomials. These Bessel polynomials are then written as
y n ( u ) = ( 2 n ) ! n ! 2 n u n q n ( 1 u ) = k = 0 n ( n + k ) ! 2 k k ! ( n k ) ! u k .
The so-called linearization coefficients of the Bessel polynomials LCBPs, β k ( n , m ) [2], are defined by
q n ( a u ) q m ( ( 1 a ) u ) = k = m i n ( n , m ) n + m β k ( n , m ) ( a ) q k ( u ) .
For example, we have β k ( 2 , 3 ) = 0 for k = 0 , 1 and
β 2 ( 2 , 3 ) ( a ) = a 5 , β 3 ( 2 , 3 ) ( a ) = ( 1 + a ) ( 15 a 4 20 a 3 + 15 a 2 6 a + 1 ) , β 4 ( 2 , 3 ) ( a ) = 7 a ( 1 + a ) 2 ( 5 a 2 4 a + 1 ) , β 5 ( 2 , 3 ) ( a ) = 21 a 2 ( 1 + a ) 3 ,
and we have
q 2 ( a u ) q 3 ( ( 1 a ) u ) = k = 2 5 β k ( 2 , 3 ) ( a ) q k ( u ) .

2. Berg and Vignat’s Results

In 2008, Berg and Vignat [2,3] established some important results about the following:
  • The connection coefficients c k ( n ) ( a ) . They proved their non-negativity for a [ 0 , 1 ] in the expansion
    q n ( a u ) = k = 0 n c k ( n ) ( a ) q k ( u )
    with
    c k ( n ) ( a ) = a k ( 1 a ) n k 2 n 2 k j = 0 m i n ( ( n j 1 ) , k ) n + 1 k j n k 1 j ( 1 a ) j
    = a k ( 1 a ) n k n + 1 k 2 n 2 k F 1 2 k , n + k + 1 n k + 2 ; 1 a .
  • The coefficients β k ( n , m ) ( a ) with applications. In particular, they proved Theorem 2.1 [2]: β n ( n , 0 ) ( a ) = c 0 ( n ) ( a ) β k ( n , 0 ) ( a ) = a n and, for 0 k n 1 ,
    β k ( n , 0 ) ( a ) = c k ( n ) ( a ) .
In fact, for n , m 1 , Berg and Vignat [2] proved that the coefficients β k ( n , m ) ( a ) satisfy the following recurrence relation, Lemma 3.6 [2]:
1 2 k + 1 β k + 1 ( n , m ) ( a ) = a 2 2 n 1 β k ( n 1 , m ) ( a ) + ( 1 a ) 2 2 m 1 β k ( n , m 1 ) ( a )
for k = 0 , 1 , , m + n 1 and β 0 ( n , m ) ( a ) = 0 (because n , m 1 ). From (3) and (5), they derived the positivity of β k ( n , m ) ( a ) when 0 a 1 , and also, they proved that β k ( n , m ) ( a ) = 0 for k < min ( m , n ) . However, an explicit single-sum formula for β k ( n , m ) ( a ) was missing in their paper.
Berg and Vignat [3] used the Student t-distribution and β k ( n , n ) ( 1 2 ) to present an alternative evaluation of the Boros–Moll integral [4]:
N 0 , 4 ( x ; m ) : = 0 + d s ( 1 + 2 x s 2 + s 4 ) m + 1 = π P m ( x ) 2 m + 3 2 ( x + 1 ) m + 1 2 ,
where
P m ( x ) = j = 0 m d j , m x j
and
d j , m = 2 2 m i = j m 2 i 2 m 2 i m i m + i m i j .
We present, at the end of this paper, a very short set of Maple instructions in Appendix A in order to assure readers about the correctness of this result.

3. Atia and Zeng’s Results

In 2012, Atia and Zeng presented this explicit single-sum formula for β k ( n , m ) ( a ) , providing the unique solution of recurrence system (5) with the boundary condition (3) [5] summarized in the following theorem.
Theorem 1.
For k = 0 , 1 , , n + m , we have
β k ( n , m ) ( a ) = a 2 n + m k ( 1 a ) n + k ( 1 / 2 ) n + m k ( 1 / 2 ) k ( 1 / 2 ) n ( 1 / 2 ) m ×   j = 0 2 ( n + m k ) ( 1 ) j n + m + 1 2 n + 2 m 2 k j m + k + j j a j .
Equivalently, the coefficients β k ( n , m ) ( a ) can be written in terms of a hypergeometric function as follows:
(i
If k ( n + m 1 ) / 2 , then
β k ( n , m ) ( a ) = a 2 n + m k ( 1 a ) n + k ( 1 / 2 ) n + m k ( 1 / 2 ) k ( 1 / 2 ) n ( 1 / 2 ) m n + m + 1 2 n + 2 m 2 k ×   F 1 2 2 m 2 n + 2 k , m + k + 1 2 k m n + 2 ; 1 a ;
(ii
If k ( n + m 1 ) / 2 , then
β k ( n , m ) ( a ) = a n + k + 1 ( 1 a ) n + k ( 1 / 2 ) n + m k ( 1 / 2 ) k ( 1 / 2 ) n ( 1 / 2 ) m n k 1 n + m 2 k 1 ×   ( 1 ) n + m 1 F 1 2 m n 1 , n k n + m 2 k ; 1 a .
Atia and Zeng derived some immediate consequences from the above Theorem 1. In fact, if we assume that n m and 0 a 1 , then the following apply:
  • If k < n , then the binomial coefficient n k 1 n + m 2 k 1 vanishes because n + m 2 k 1 > n k 1 . Thus, Formula (9) implies that β k ( n , m ) ( a ) = 0 for k < m n .
  • If k n , then each term in (7) is a polynomial in a. Indeed, it is clear that ( 1 a ) n + k Z [ a ] alongside the following:
    • If k m , then 2 n + m k j 2 n + m k 2 ( n + m k ) = k m 0 for 0 j 2 ( n + m k ) ;
    • If k < m , then the binomial coefficient m + k + j j does not vanish only if m + k + j < 0 , that is, j < m k ; therefore, 2 n + m k j > 2 n + m k ( m i ) = 2 n 0 .
  • From definition (1), the symmetry property is immediately derived as
    β k ( n , m ) ( a ) = β k ( m , n ) ( 1 a )
    which follows also from (8) (resp. (9)) and Pfaff’s transformation (see p. 94 [6]):
    F 1 2 A B C ; x = ( 1 x ) A F 1 2 A , C B C ; x x 1 ,
    with A = 2 n 2 m + 2 k , B = m + k + 1 , C = n m + 2 k + 2 , and x = 1 / a if k ( n + m 1 ) / 2 (resp. A = n m 1 , B = n k , C = n + m 2 k , and x = 1 / a if k ( n + m 1 ) / 2 ).
For applications of this Theorem 1, Atia and Zeng derived a formula from Berg and Vignat’s [2] when m = n , and they also proved the positivity of these coefficients, β k ( n , m ) ( a ) , m i n ( n , m ) k n + m . In fact, Theorem 1 implies this immediately following positivity result. Suppose that n < m . Then, the coefficient β k ( n , m ) ( a ) is positive for n k < m (for the proof, please see [5]).
As a particular case of the coefficient β k ( n , m ) ( a ) , Atia and Zeng found the coefficients β i ( n ) ( a ) defined by
q n ( a u ) q n ( ( 1 a ) u ) = i = 0 n β i ( n ) ( a ) q n + i ( u ) .
For 0 i n , we have
β i ( n ) ( a ) = ( 4 a ( 1 a ) ) i 4 n ( n ) i ( n + 1 2 ) i i ! ( n + 1 2 ) i F 1 2 n + i , n 1 2 ; 1 2 ; ( 2 a 1 ) 2 × ( 2 n 2 i ) ! ( 2 n + 2 i ) ! ( n i ) ! ( n + i ) ! j = 0 n i 2 n + 1 2 j n j i ( 2 a 1 ) 2 j ,
(for the proof, please see [5]).
Finally, Atia and Zeng generalized a formula on an integral evaluation of Boros and Moll’s integral [4]:
+ d y ( 1 + y 2 ) n + 1 ( 1 + ( x y ) 2 ) m + 1 = 1 2 k = m n m + n A k + 1 2 A n + 1 2 A m + 1 2 β k ( n , m ) ( 1 2 ) ( 1 + 1 4 x 2 ) k + 1
where
A ν = Γ ( ν + 1 2 ) Γ ( 1 2 ) Γ ( ν ) .
We present, at the end of this paper, a very short set of Maple instructions in Appendix B in order to assure readers about the correctness of the result.

4. BenAbdallah and Atia’s Results

In 2018, BenAbdallah and Atia presented a recursion formula (in a more general case) for the linearization coefficients β l ( N k ) ( v k ) for Bessel polynomials ( q n ) , mentioned by C. Berg and C. Vignat (p. 21 (14) [2]), in the expansion q n 1 ( a 1 u ) q n 2 ( a 2 u ) q n k ( a k u ) = l = min ( n 1 , n 2 , , n k ) L k β l ( N k ) ( v k ) q l ( u ) , where u R , N k = ( n 1 , n 2 , , n k ) N k , v k = ( a 1 , a 2 , , a k ) R + k with 1 i k a i = 1 , and L k = 1 i k n i . BenAbdallah and Atia proved that this recursion formula yields, again, the positivity of the coefficients. In addition, in the case when k = 3 , they presented an explicit formula for β l ( N 3 ) ( v 3 ) , and they derived a double-integral formula through the convolution of Student t-densities. As a bonus, two reductions for this double integral were provided; the first reduction uses an integral formula given by Atia and Zeng [5], and the second uses the one given by Boros and Moll [4].
In fact, BenAbdallah and Atia studied the general case of the Berg–Vignat linearization problem v k = ( a 1 , , a k ) R k , k 3 , with i = 1 k a i = 1 ; see (p. 21 (14) [2]).
q n 1 ( a 1 u ) q n 2 ( a 2 u ) q n k ( a k u ) = l = min ( n 1 , n 2 , , n k ) L k β l ( N k ) ( v k ) q l ( u ) .
BenAbdallah and Atia performed the following:
  • Provided the recurrence formula satisfied by the linearization coefficients β l ( N k ) ( v k ) in (12) and then proved that they are non-negative and min ( n 1 , n 2 , , n k ) l L k , with L k = 1 i k n i .
  • Provided a triple-sum formula of the linearization coefficients β l ( N 3 ) ( v 3 ) in the expansion
    q n 1 ( a 1 u ) q n 2 ( a 2 u ) q n 3 ( a 3 u ) = l = min ( n 1 , n 2 , n 3 ) L 3 β l ( N 3 ) ( v 3 ) q l ( u ) , u R .
  • Provided an evaluation through a convolution of the Student t-densities of the double integral
    I n , m , p ( x ) : = R R A n + 1 2 A m + 1 2 A p + 1 2 d t d y ( 1 + t 2 ) n + 1 ( 1 + y 2 ) m + 1 ( 1 + ( x t y ) 2 ) p + 1
where A ν is defined in (11).
In 2017, Atia and BenAbdallah generalized (1) in their paper published in the Asian–European Journal of Mathematics:
q n ( a u ) q m ( b u ) = l = 0 n + m β l ( n , m ) ( a , b ) q l ( u ) , a + b 1 ,
and also, they proved the positivity of these linearization coefficients β l ( n , m ) ( a , b ) as follows:
(i)
For 0 k n + m 1 2 we have
β k ( n , m ) ( a , b ) = ( 2 k ) ! ( 1 ) k + 1 a n b 2 k + 1 n k ! i = 0 k + 1 ( 1 ) i ( 2 k + 2 i ) ! b i i ! ( k + 1 i ) !
×   j = m a x ( 0 , n + i 2 k 1 ) n ( n ) n j ( m ) 2 k + 1 i n + j ( b a ) j ( 2 n ) n j ( 2 m ) 2 k + 1 i n + j ( n j ) ! ( 2 k + 1 i n + j ) ! .
(ii)
For n + m 1 2 k n + m we have
β k ( n , m ) ( a , b ) = ( 2 k ) ! ( 1 ) m + n k a n b m k ! i = 0 n + m k ( 1 ) i ( m + n i + 1 ) ! b i ( 2 k + 1 + i n m ) ! ( m + n k i ) !
×   j = 0 n ( n ) n j ( m ) m i + j ( b a ) j ( 2 n ) n j ( 2 m ) m i + j ( n j ) ! ( m i + j ) ! ,
for 0 k n + m and a 0 and b 0 and where the authors assumed that n m .
The linearization coefficients in (12) satisfy the following recursion formula.
For n i 1 , with 1 i k , we have
1 2 l + 1 β l + 1 ( N k ) ( v k ) = i = 1 k a i 2 2 n i 1 β l σ i ( N k ) ( v k )
where σ i ( N k ) = ( n 1 , , n i 1 , , n k ) , and l = 0 , 1 , , L k 1 . Furthermore, β 0 ( N k ) ( v k ) = 0 . For the proof, please see [7].
Atia and Benabdallah presented a triple-sum formula of the linearization coefficients β l ( N 3 ) ( v 3 ) in the expansion
q n 1 ( a 1 u ) q n 2 ( a 2 u ) q n 3 ( a 3 u ) = l = min ( n 1 , n 2 , n 3 ) n 1 + n 2 + n 3 β l ( N 3 ) ( v 3 ) q l ( u ) , u R
where β l ( N 3 ) ( v 3 ) is given by
β l ( N 3 ) ( v 3 ) = k = 0 min ( n 1 + n 2 + n 3 l , l + 1 ) δ l ( k + l ) j = max ( 0 , k + l n 3 ) min ( k + l , n 1 + n 2 ) ( i = max ( 0 , j n 2 ) min ( j , n 1 ) α i ( n 1 ) α j i ( n 2 ) a 1 i a 2 j i ) α k + l j ( n 3 ) a 3 k + l j
which, for l min ( n 1 , n 2 , n 3 ) , they wrote it under the form
β l ( N 3 ) ( v 3 ) = ( 2 l ) ! l ! k = 0 min ( L 3 l , l + 1 ) ( 1 ) k ( l + k + 1 ) ! k ! ( l k + 1 ) !
×   j = max ( 0 , k + l n 3 ) min ( k + l , n 1 + n 2 ) i = max ( 0 , j n 2 ) min ( j , n 1 ) ( n 1 ) i ( n 2 ) j i ( n 3 ) k + l j ( 2 n 1 ) i ( 2 n 2 ) j i ( 2 n 3 ) k + l j a 1 i a 2 j i a 3 k + l j i ! ( j i ) ! ( k + l j ) !
Finally, Atia and Benabdallah presented the evaluation of the following multiple integrals:
R R ( k 1 ) t i m e s A n k + 1 2 j = 1 k 1 A n j + 1 2 d x j ( 1 + x j 2 ) n j + 1 1 + ( x j = 1 k 1 x j ) 2 n k + 1 = 1 k l = min ( n 1 , n 2 , , n k ) L k A l + 1 2 β l ( N k ) ( 1 k ) ( 1 + ( x k ) 2 ) l + 1 .
This formula becomes, as a special case, a double integral when k = 3 and N 3 = ( n , m , p ) and v 3 = ( 1 3 , 1 3 , 1 3 ) :
R R A n + 1 2 A m + 1 2 A p + 1 2 d t d y ( 1 + t 2 ) n + 1 ( 1 + y 2 ) m + 1 ( 1 + ( x t y ) 2 ) p + 1 = 1 3 l = min ( n , m , p ) n + m + p A l + 1 2 β l ( N 3 ) ( 1 3 ) ( 1 + x 2 9 ) l + 1 .
We present, at the end of this paper, a very short set of Maple instructions in Appendix C in order to assure readers about the correctness of the result.

5. Results

5.1. Another Recursion Formula for LCBPs

In 2008, Berg and Vignat [2] proved that the coefficients β k ( n , m ) ( a ) satisfy the recurrence given by (5). In this section, we are going to present our first main result concerning another recurrence relation between β k ( n , m ) ( a ) , and at the end of this paper, we will show why this relation is as important as the relation given by Berg and Vignat.
Theorem 2.
For any integers n , m , n 1 , m 1 , and for any 1 k n + m , we have the following results:
( 2 m 1 ) β k ( n , m ) ( a ) + ( 2 n + 1 ) β k ( n + 1 , m 1 ) ( a ) = ( 2 k 1 ) β k 1 ( n , m 1 ) ( a ) + ( 2 n + 2 m 2 k 1 ) β k ( n , m 1 ) ( a ) ,
taking into account either β 1 ( n , m 1 ) ( a ) = 0 or β n + m ( n , m 1 ) ( a ) = 0 .
Proof. 
Let us prove that the difference
S u b s = L H S R H S = 0 .
S u b s = ( 2 m 1 ) β k ( n , m ) ( a ) + ( 2 n + 1 ) β k ( n + 1 , m 1 ) ( a ) ( 2 k 1 ) β k 1 ( n , m 1 ) ( a ) ( 2 n + 2 m 2 k 1 ) β k ( n , m 1 ) ( a ) = ( 1 2 m ) a 2 n + m k ( 1 a ) k n ( 1 2 ) k ( 1 2 ) n + m k ( 1 2 ) n ( 1 2 ) m j = 0 2 n + 2 m 2 k ( 1 ) j a j m + n + 1 2 n + 2 m 2 k j k m + j j ( 2 n + 1 ) a 2 n + m k + 1 ( 1 a ) k n 1 ( 1 2 ) k ( 1 2 ) n + m k ( 1 2 ) n + 1 ( 1 2 ) m 1 j = 0 2 n + 2 m 2 k ( 1 ) j a j m + n + 1 2 n + 2 m 2 k j k m + j + 1 j ( 2 k 1 ) a 2 n + m k ( 1 a ) k n 1 ( 1 2 ) k 1 ( 1 2 ) n + m k ( 1 2 ) n ( 1 2 ) m 1
j = 0 2 n + 2 m 2 k ( a ) j m + n 2 n + 2 m 2 k j k m + j j ( 2 n + 2 m 2 k 1 ) a 2 n + m k 1 ( 1 a ) k n ( 1 2 ) k ( 1 2 ) n + m k 1 ( 1 2 ) n ( 1 2 ) m 1 j = 0 2 n + 2 m 2 k 2 ( a ) j m + n 2 n + 2 m 2 k j 2 m + k + j + 1 j .
Using that for any integer p, ( 1 2 ) p + 1 = ( p + 1 2 ) ( 1 2 ) p and taking away the common factor the remaining terms of the S u b s that we denote by R s u b s , we obtain
R s u b s = a ( 1 a ) j = 0 2 n + 2 m 2 k ( a ) j m + n + 1 2 n + 2 m 2 k j m + k + j j + a 2 j = 0 2 n + 2 m 2 k ( a ) j m + n + 1 2 n + 2 m 2 k j m + k + j + 1 j
a j = 0 2 n + 2 m 2 k ( a ) j m + n 2 n + 2 m 2 k j m + k + j j ( 1 a ) j = 0 2 n + 2 m 2 k 2 ( a ) j m + n 2 n + 2 m 2 k j 2 m + k + j + 1 j .
The last summation can be expanded to 2 n + 2 m 2 k because m + n 1 = m + n 2 = 0 . Then, we obtain
R s u b s = a ( 1 a ) j = 0 2 n + 2 m 2 k ( a ) j m + n + 1 2 n + 2 m 2 k j m + k + j j + a 2 j = 0 2 n + 2 m 2 k ( a ) j m + n + 1 2 n + 2 m 2 k j m + k + j + 1 j a j = 0 2 n + 2 m 2 k ( a ) j m + n 2 n + 2 m 2 k j m + k + j j ( 1 a ) j = 0 2 n + 2 m 2 k ( a ) j m + n 2 n + 2 m 2 k j 2 m + k + j + 1 j .
Expanding the first and the last sums, we obtain
R s u b s = a j = 0 2 n + 2 m 2 k ( a ) j m + n + 1 2 n + 2 m 2 k j m + k + j j a 2 j = 0 2 n + 2 m 2 k ( a ) j m + n + 1 2 n + 2 m 2 k j m + k + j j + a 2 j = 0 2 n + 2 m 2 k ( a ) j m + n + 1 2 n + 2 m 2 k j m + k + j + 1 j a j = 0 2 n + 2 m 2 k ( a ) j m + n 2 n + 2 m 2 k j m + k + j j j = 0 2 n + 2 m 2 k ( a ) j m + n 2 n + 2 m 2 k j 2 m + k + j + 1 j + a j = 0 2 n + 2 m 2 k ( a ) j m + n 2 n + 2 m 2 k j 2 m + k + j + 1 j ,
Taking the sums that have the same power of a, we obtain
R s u b s = j = 0 2 n + 2 m 2 k ( a ) 2 j m + n + 1 2 n + 2 m 2 k j ( m + k + j + 1 j m + k + j j ) + j = 0 2 n + 2 m 2 k ( a ) 1 j ( m + n 2 n + 2 m 2 k j m + k + j j
m + n + 1 2 n + 2 m 2 k j m + k + j j m + n 2 n + 2 m 2 k j 2 m + k + j + 1 j ) j = 0 2 n + 2 m 2 k ( a ) j m + n 2 n + 2 m 2 k j 2 m + k + j + 1 j .
Let us denote the above by
B j ( 2 ) = m + n + 1 2 n + 2 m 2 k j ( m + k + j + 1 j m + k + j j ) , B j ( 1 ) = ( m + n 2 n + 2 m 2 k j m + k + j j m + n + 1 2 n + 2 m 2 k j m + k + j j m + n 2 n + 2 m 2 k j 2 m + k + j + 1 j ) ,
B j ( 0 ) = m + n 2 n + 2 m 2 k j 2 m + k + j + 1 j ,
Then, R s u b s becomes
R s u b s = j = 0 2 n + 2 m 2 k ( a ) 2 j B j ( 2 ) + j = 0 2 n + 2 m 2 k ( a ) 1 j B j ( 1 ) j = 0 2 n + 2 m 2 k ( a ) j B j ( 0 ) .
In these sums, we have the following:
  • Only one term with a 2 with the coefficient
    B 0 ( 2 ) = ( m + n + 1 2 n + 2 m 2 k ( m + k + 1 0 m + k 0 ) ) a 2 = 0 ,
  • Two terms with a with the coefficient
    B 1 ( 2 ) + B 0 ( 1 ) m + n + 1 2 n + 2 m 2 k 1 ( m + k + 1 + 1 1 m + k + 1 1 ) + ( m + n 2 n + 2 m 2 k 1 m + k + 1 1 m + n + 1 2 n + 2 m 2 k 1 m + k + 1 1 m + n 2 n + 2 m 2 k 1 2 m + k + 1 + 1 1 ) = 0 ,
  • Only one term with a ( 2 n + 2 m 2 k ) with the coefficient
    B 2 n + 2 m 2 k ( 0 ) = m + n 2 m + k + 2 n + 2 m 2 k + 1 2 n + 2 m 2 k = 0 ,
  • Only two terms with a ( 2 n + 2 m 2 k ) + 1 with the coefficient
    B 2 n + 2 m 2 k 1 ( 0 ) + B 2 n + 2 m 2 k ( 1 ) = m + n 1 2 n + m k + 1 2 n + 2 m 2 k + ( m + n 0 2 n + m k 2 n + 2 m 2 k m + n + 1 0 2 n + m k 2 n + 2 m 2 k m + n 2 2 n + m k + 1 2 n + 2 m 2 k ) = 0 ,
    Finally, the terms with a j , 0 j 2 n + 2 m 2 k 2 are given by
    B j + 2 ( 2 ) + B j + 1 ( 1 ) + B j ( 0 ) = 0 .

5.2. The Matrix Linearization Coefficients of Bessel Polynomials

From (1), we can write
q n ( a u ) q 0 ( ( 1 a ) u ) = β 0 ( n , 0 ) ( a ) q 0 ( u ) + β 1 ( n , 0 ) ( a ) q 1 ( u ) + + β n ( n , 0 ) ( a ) q n ( u ) , q n 1 ( a u ) q 1 ( ( 1 a ) u ) = 0 + β 1 ( n 1 , 1 ) ( a ) q 1 ( u ) + + β n ( n 1 , 1 ) ( a ) q n ( u ) , q n 2 ( a u ) q 2 ( ( 1 a ) u ) = 0 + 0 + β 2 ( n 2 , 2 ) ( a ) q 2 ( u ) + + β n ( n 2 , 2 ) ( a ) q n ( u ) , q 2 ( a u ) q n 2 ( ( 1 a ) u ) = 0 + 0 + β 2 ( 2 , n 2 ) ( a ) q 2 ( u ) + + β n ( 2 , n 2 ) ( a ) q n ( u ) , q 1 ( a u ) q n 1 ( ( 1 a ) u ) = 0 + β 1 ( 1 , n 1 ) ( a ) q 1 ( u ) + + β n ( 1 , n 1 ) ( a ) q n ( u ) , q 0 ( a u ) q n ( ( 1 a ) u ) = β 0 ( 0 , n ) ( a ) q 0 ( u ) + β 1 ( 0 , n ) ( a ) q 1 ( u ) + + β n ( 0 , n ) ( a ) q n ( u ) .
This can be written under the following matrix representation:
q n ( a u ) q 0 ( ( 1 a ) u ) q n 1 ( a u ) q 1 ( ( 1 a ) u ) q 1 ( a u ) q n 1 ( ( 1 a ) u ) q 0 ( a u ) q n ( ( 1 a ) u ) = = β 0 ( n , 0 ) ( a ) β 1 ( n , 0 ) ( a ) β n 1 ( n , 0 ) ( a ) β n ( n , 0 ) ( a ) 0 β 1 ( n 1 , 1 ) ( a ) β n 1 ( n 1 , 1 ) ( a ) β n ( n 1 , 1 ) ( a ) 0 β 1 ( 1 , n 1 ) ( a ) β n 1 ( 1 , n 1 ) ( a ) β n ( 1 , n 1 ) ( a ) β 0 ( 0 , n ) ( a ) β 1 ( 0 , n ) ( a ) β n 1 ( 0 , n ) ( a ) β n ( 0 , n ) ( a ) q 0 ( u ) q 1 ( u ) q n 1 ( u ) q n ( u ) ,
which we denote by
q n ( a u ) q 0 ( ( 1 a ) u ) q n 1 ( a u ) q 1 ( ( 1 a ) u ) q 1 ( a u ) q n 1 ( ( 1 a ) u ) q 0 ( a u ) q n ( ( 1 a ) u ) = M n β q 0 ( u ) q 1 ( u ) q n 1 ( u ) q n ( u ) ,
Please note that M n β has an order of ( n + 1 ) × ( n + 1 ) ; here are some examples.
Example 1.
M 1 β = 1 a a a 1 a , M 2 β = 1 a a ( 1 a ) a 2 0 3 a 2 3 a + 1 3 a ( 1 a ) a a ( 1 a ) ( 1 a ) 2
M 3 β = 1 a 1 5 a ( 1 a ) ( a 5 ) 6 5 a 2 ( 1 a ) a 3 0 ( 1 a ) 3 a ( 6 a 2 8 a + 3 ) 5 a 2 ( 1 a ) 0 a 3 ( 1 a ) ( 6 a 2 4 a + 1 ) 5 a ( 1 a ) 2 a 1 5 a ( 1 a ) ( a + 4 ) 6 5 a ( 1 a ) 2 ( 1 a ) 3
After computing the traces of these three first matrices, we can announce our first result.
Proposition 1.
For any integer n such that n 0 , the following result holds:
Trace ( M n β ) = Tr ( M n β ) = 1 + ( 1 2 a ) Tr ( M n 1 β ) = 1 ( 1 2 a ) n + 1 2 a
Tr ( M n β ) = 1 + ( 1 2 a ) + ( 1 2 a ) 2 + + ( 1 2 a ) n .
Proof. 
Let us find the first values. Using (7), we obtain
Tr ( M 0 β ) = β 0 ( 0 , 0 ) ( a ) = 1 , Tr ( M 1 β ) = β 0 ( 1 , 0 ) ( a ) + β 1 ( 0 , 1 ) ( a ) = ( 1 a ) + ( 1 a ) = 2 2 a = 1 + ( 1 2 a ) , Tr ( M 2 β ) = β 0 ( 2 , 0 ) ( a ) + β 1 ( 1 , 1 ) ( a ) + β 2 ( 0 , 2 ) ( a ) = ( 1 a ) + ( 3 a 2 3 a + 1 ) + ( 1 a ) 2 = 4 a 2 6 a + 3 = 1 + ( 1 2 a ) + ( 1 2 a ) 2 = 1 + ( 1 2 a ) ( 2 2 a ) , Tr ( M 3 β ) = β 0 ( 3 , 0 ) ( a ) + β 1 ( 2 , 1 ) ( a ) + β 2 ( 1 , 2 ) ( a ) + β 3 ( 0 , 3 ) ( a ) = ( 1 a ) + ( 1 a ) 3 + ( 1 a ) ( 6 a 2 4 a + 1 ) + ( 1 a ) 3 = 4 ( 1 a ) ( 2 a 2 2 a + 1 ) = 1 + ( 1 2 a ) ( 4 a 2 6 a + 3 ) .
Here, it is a known fact that if we have
n N , u n + 1 = α u n + β , α 1 , β 0
then when u n converges, the limit l fulfills l = α l + β . We define v n as the auxiliary sequence defined by n N , v n = u n l then we obtain: n N , v n + 1 = α v n . Thus, v n is a geometric sequence with ratio α , and v n = α n v 0 and u n = l + α n ( u 0 l ) . The sequence u n converges if and only if | a | < 1 or u 0 = l . In our case, if we have the recurrence
Tr ( M n β ) = 1 + ( 1 2 a ) T r ( M n 1 β ) ,
then when a 1 2 , the limit of Tr ( M n β ) is l = 1 2 a , and
Tr ( M n β ) = 1 ( 1 2 a ) n + 1 2 a .
We still need to prove that for any positive integer n, we have
Tr ( M n β ) = 1 + ( 1 2 a ) Tr ( M n 1 β ) .
Using Theorem 1, we can write
β k ( n k , k ) ( a ) = a 2 n 2 k ( 1 a ) n 2 k j = 0 2 n 2 k ( a ) j n + 1 2 n 2 k j
which, for n = 2 p , yields
Tr ( M n β ) = k = 1 p ( 1 a ) 2 p + 1 2 k + k = 0 p ( 1 a ) 2 k ( ( 1 a ) 1 + 2 p + l = 0 2 k 1 + 2 p 2 p + 1 l ( 1 ) l a 2 p + 1 l )
and for n = 2 p 1 , we have
Tr ( M n β ) = k = 1 p ( 1 a ) 2 p 2 k 1 + k = 0 p 1 ( 1 a ) 2 k + 1 ( ( 1 a ) 2 p + l = 0 2 k + 1 2 p 2 p l ( 1 ) l a 2 p l ) .
Let us prove that for n = 2 p , we have
Tr ( M n β ) = 1 ( 1 2 a ) n + 1 2 a = 1 ( 1 2 a ) 2 p + 1 2 a .
Using (19), we obtain
Tr ( M n β ) = k = 1 p ( 1 a ) 2 p + 1 2 k + k = 0 p ( 1 a ) 2 k ( ( 1 a ) 1 + 2 p + l = 0 2 k 1 + 2 p 2 p + 1 l ( 1 ) l a 2 p + 1 l ) ) = ( 1 a ) 2 p 1 + ( 1 a ) 2 p 3 + + ( 1 a ) 1 + ( 1 a ) 0 ( 1 a ) 2 p + 1 + 2 p + 1 2 p + 1 a 2 p + 1 + ( 1 a ) 2 ( 1 a ) 2 p + 1 + 2 p + 1 2 p + 1 a 2 p + 1 2 p + 1 2 p a 2 p + 2 p + 1 2 p 1 a 2 p 1 + . . . . + ( 1 a ) 2 p
We can write this under the form
Tr ( M n β ) = 2 p + 1 2 p a 2 p 2 p + 1 2 p 1 a 2 p 1 + k = 0 p 1 ( 1 a ) 2 k + 1 + a 2 k 1 ( 2 p + 1 2 k a 2 p + 1 2 k 1 ) i = 0 p k ( 1 a ) 2 i .
Using the following
i = 0 m ( 1 a ) 2 i = ( 1 a ) 2 m + 2 1 a ( a 2 )
we obtain
Tr ( M n β ) = 2 p + 1 2 p a 2 p 2 p + 1 2 p 1 a 2 p 1 + ( 1 a ) ( 1 a ) 2 ( p 1 ) + 2 1 a ( a 2 ) + k = 0 p 1 a 2 k 1 ( 2 p + 1 2 k a 2 p + 1 2 k 1 ) ( 1 a ) 2 ( p k ) + 2 1 a ( a 2 ) ,
Equivalently,
Tr ( M n β ) = 2 p + 1 2 p a 2 p 2 p + 1 2 p 1 a 2 p 1 + ( 1 a ) ( 1 a ) 2 ( p 1 ) + 2 1 a ( a 2 ) + ( a 1 ) k = 0 p 1 a 2 k 1 2 p + 1 2 k ( 1 a ) 2 ( p k ) + 2 1 a ( a 2 ) + k = 0 p 1 a 2 k 1 ( 2 p + 1 2 k 2 p + 1 2 k 1 ) ( 1 a ) 2 ( p k ) + 2 1 a ( a 2 )
This becomes
Tr ( M n β ) = ( 2 p + 1 ) a 2 p p ( 2 p + 1 ) a 2 p 1 + ( 1 + a ) ( 1 ( 1 + a ) 2 p ) a ( a 2 ) + ( a 1 ) k = 0 p 1 a 2 k 1 2 p + 1 2 k ( 1 a ) 2 ( p k ) + 2 1 a ( a 2 ) + k = 0 p 1 a 2 k 1 ( 2 p + 1 2 k 2 p + 1 2 k 1 ) ( 1 a ) 2 ( p k ) + 2 1 a ( a 2 ) .
Using
( a 1 ) k = 0 p 1 a 2 k 1 2 p + 1 2 k ( 1 a ) 2 ( p k ) + 2 1 a ( a 2 ) = ( a 1 ) 2 ( 1 + 2 a ) 2 p + 1 2 a 2 ( a 2 ) ( 1 + a ) ( 2 p + 1 ) a 2 p a + ( a 1 ) 2 p + 2 ( a 1 ) ( a + 1 ) 2 p + 1 2 a 2 ( a 2 ) ( a 1 ) 2 2 a 2 ( a 2 )
and
k = 0 p 1 a 2 k 1 ( 2 p + 1 2 k 2 p + 1 2 k 1 ) ( 1 a ) 2 ( p k ) + 2 1 a ( a 2 ) = ( 2 p + 1 ) ( p 1 ) a 2 p 1 + ( 1 2 a ) + ( 1 2 a ) 2 p + 1 + ( 1 + a ) 2 p + 1 ( a + 1 ) + ( a + 1 ) 2 p + 1 ( 1 + a ) ) 2 a 2 ( a 2 )
and computing 2 a Tr ( M n β ) 1 , we find ( 2 a 1 ) 2 p + 1 . We follow the same steps to prove the result for n = 2 p .
Please see the Maple instructions given in Appendix D at the end of this paper, which shortens the calculus presented here. □
Our second result deals with the determinant. Using (18), we can state the following result.
Proposition 2.
For any positive integer n 0 , the following result holds:
Determinant ( M n β ) = Det ( M n β ) = ( 1 2 a ) n ( n + 1 ) 2 .
Proof. 
Using (18), we can see that the trace of M n β differs from that of M n + 1 β in that the last term is ( 1 2 a ) n + 1 . For instance, using (15) and (16), we can see that the eigenvalues appear as follows:
  • Det ( M 1 β λ I 2 ) = ( λ 1 ) ( λ ( 1 2 a ) )
    Tr ( M 1 β ) = 1 + ( 1 2 a ) ,
  • Det ( M 2 β λ I 3 ) = ( λ 1 ) ( λ ( 1 2 a ) ) ( λ ( 1 2 a ) 2 )
    Tr ( M 1 β ) = 1 + ( 1 2 a ) + ( 1 2 a ) 2 ,
  • Det ( M 3 β λ I 4 ) = ( λ 1 ) ( λ ( 1 2 a ) ) ( λ ( 1 2 a ) 2 ) ( λ ( 1 2 a ) 3 )
    Tr ( M 1 β ) = 1 + ( 1 2 a ) + ( 1 2 a ) 2 + ( 1 2 a ) 3 ,
So, we can deduce directly that
Det ( M n + 1 β ) = ( 1 2 a ) n ( n + 1 ) 2 . ( 1 2 a ) n + 1 = ( 1 2 a ) ( n + 1 ) ( n + 2 ) 2 .
It is important to remark here that in using (18), we can see that the Spectrum of M n + 1 β contains the Spectrum of M n β . □
As an application on these matrices and using (10), we can write the following result:
Application A n + 1 2 A 1 2 + d y ( 1 + y 2 ) n + 1 ( 1 + ( x y ) 2 ) A n 1 2 A 3 2 + d y ( 1 + y 2 ) n ( 1 + ( x y ) 2 ) 2 A 3 2 A n 1 2 + d y ( 1 + y 2 ) 2 ( 1 + ( x y ) 2 ) n A 1 2 A n + 1 2 + d y ( 1 + y 2 ) 1 ( 1 + ( x y ) 2 ) n + 1 = β 0 ( n , 0 ) ( 1 2 ) β 1 ( n , 0 ) ( 1 2 ) β n 1 ( n , 0 ) ( 1 2 ) β n ( n , 0 ) ( 1 2 ) 0 β 1 ( n 1 , 1 ) ( 1 2 ) β n 1 ( n 1 , 1 ) ( 1 2 ) β n ( n 1 , 1 ) ( 1 2 ) 0 β 1 ( 1 , n 1 ) ( 1 2 ) β n 1 ( 1 , n 1 ) ( 1 2 ) β n ( 1 , n 1 ) ( 1 2 ) β 0 ( 0 , n ) ( 1 2 ) β 1 ( 0 , n ) ( 1 2 ) β n 1 ( 0 , n ) ( 1 2 ) β n ( 0 , n ) ( 1 2 ) A 1 2 ( 1 + x 2 4 ) A 3 2 ( 1 + x 2 4 ) 2 A n 1 2 ( 1 + x 2 4 ) n 1 A n + 1 2 ( 1 + x 2 4 ) n .

5.3. The Inverse of the Linearization Coefficients of Bessel Polynomials

5.3.1. The Inverse of the Connection Coefficients

From (2) and (3), we can talk about the inverse of the connection coefficients c k n ( a ) , which we denote by i c k ( n ) ( a ) , defined by the expansion
q n ( u ) = k = 0 n i c k ( n ) ( a ) q k ( a u ) , n 0
and we can announce the following theorem.
Theorem 3.
The connection coefficients i c k ( n ) ( a ) are given by c k ( n ) ( 1 a ) , i.e.,
q n ( u ) = k = 0 n i c k ( n ) ( a ) q k ( a u ) = k = 0 n c k ( n ) ( 1 a ) q k ( a u ) , n 0 .
Proof. 
From (2), we have
q n ( b u ) = k = 0 n c k ( n ) ( b ) q k ( u ) ,
Then, the substitution of u by u b , b 0 yields
q n ( u ) = k = 0 n c k ( n ) ( b ) q k ( u b )
and finally, if we denote 1 b by a, we find
q n ( u ) = k = 0 n c k ( n ) ( 1 b ) q k ( u ) .
By identification, we find the desired result. □

5.3.2. Application

With this result, we can find a very interesting relation between the modified Bessel function of the second kind, K ν , ν = i + 1 2 , where 0 i n is an integer [8].
Proposition 3.
We have the following relation involving the modified Bessel function of the second kind:
e ( 1 a ) u Γ ( n + 1 2 ) K n + 1 2 ( u ) = i = 0 n c i ( n ) ( 1 a ) ( 2 u ) n i a i + 1 2 Γ ( i + 1 2 ) K i + 1 2 ( a u ) .
Proof. 
This is an immediate consequence of (20), and
e u q ν 1 2 ( u ) = 2 1 ν Γ ( ν ) u ν K ν ( u ) , u 0 .
We provide, in Appendix E, a very short set of Maple instructions in order to assure readers about the correctness of the result.

5.4. The Inverse of Linearization Coefficients of Bessel Polynomials

In the previous section, we wrote (14) as follows:
q n ( a u ) q 0 ( ( 1 a ) u ) q n 1 ( a u ) q 1 ( ( 1 a ) u ) q 1 ( a u ) q n 1 ( ( 1 a ) u ) q 0 ( a u ) q n ( ( 1 a ) u ) = M n β q 0 ( u ) q 1 ( u ) q n 1 ( u ) q n ( u )
Then, it is natural to study the inverse of the matrix M n β , which we denote by M n γ and is defined as follows:
M n γ q n ( a u ) q 0 ( ( 1 a ) u ) q n 1 ( a u ) q 1 ( ( 1 a ) u ) q n 2 ( u ) q 2 ( ( 1 a ) u ) q 0 ( a u ) q n ( ( 1 a ) u ) = q 0 ( u ) q 1 ( u ) q 2 ( u ) q n ( u )
which we write as
γ 0 ( n , 0 ) ( a ) γ 1 ( n , 0 ) ( a ) γ n ( n , 0 ) ( a ) γ 0 ( n 1 , 1 ) ( a ) γ 1 ( n 1 , 1 ) ( a ) γ n ( n 1 , 1 ) ( a ) γ 0 ( n 2 , 2 ) ( a ) γ 1 ( n 2 , 2 ) ( a ) γ n ( n 2 , 2 ) ( a ) γ 0 ( 0 , n ) ( a ) γ 1 ( 0 , n ) ( a ) γ n ( 0 , n ) ( a ) q n ( a u ) q 0 ( ( 1 a ) u ) q n 1 ( a u ) q 1 ( ( 1 a ) u ) q n 2 ( u ) q 2 ( ( 1 a ) u ) q 0 ( a u ) q n ( ( 1 a ) u ) = q 0 ( u ) q 1 ( u ) q 2 ( u ) q n ( u ) .
It is straightforward to deduce that
γ i ( n k , k ) ( a ) = γ n i ( n k , k ) ( 1 a ) , 0 i , k n .
Moreover, if we take u = 0 , we see
i = 0 n γ i n k , k ( a ) = 1 , 0 k n .
Example 2.
M γ 1 = 1 2 a 1 a 1 a a a 1 ,
M γ 2 = 1 ( 2 a 1 ) 3 ( 1 + a ) 2 ( 1 3 a ) a ( 1 + a ) ( 1 + 2 a ) a 2 ( 2 + 3 a ) 3 ( 1 a ) a 2 ( a 2 a + 1 ) ( 1 + 2 a ) 3 a ( 1 + a ) 2 ( 3 a 2 3 a + 1 ) a a ( 1 a ) ( 2 a 1 ) ( 1 a ) ( 3 a 2 3 a + 1 )
and
M γ 3 = 1 ( 2 a 1 ) 5
( 10 a 2 6 a + 1 ) ( a 1 ) 3 a ( a 1 ) 2 ( 30 a 2 26 a + 5 ) 5 a 2 ( a 1 ) ( 30 a 2 34 a + 9 ) 5 10 a 3 ( a 1 ) 2 ( a 1 ) ( 6 a 4 8 a 3 + 10 a 2 6 a + 1 ) ( 6 a 4 16 a 3 + 22 a 2 14 a + 3 ) a 10 a 2 ( a 1 ) 3 ( 6 a 3 14 a 2 + 13 a 4 ) a 2 ( a 1 ) 2 ( 6 a 3 4 a 2 + 3 a 1 ) a ( 10 a 4 20 a 3 + 16 a 2 6 a + 1 ) 2 ( 15 a 3 20 a 2 + 10 a 2 ) a ( a 1 ) 5 2 ( 15 a 3 25 a 2 + 15 a 3 ) a ( a 1 ) 5
a 3 ( 10 a 2 14 a + 5 ) 10 a 2 ( a 1 ) 3 5 ( a 1 ) 2 ( 2 a 2 2 a + 1 ) a ( a 1 ) ( 10 a 4 20 a 3 + 16 a 2 6 a + 1 ) .
In the following, we are going to present all the results of the IMLCBPs only as conjectures.
In the abstract, we wrote that the new recursion formula that we provided is more important.
In fact, if we use the matrices M γ 2 and M γ 3 , we observe the following:
  • If we write the recurrence formula found by Berg and Vignat written in terms of γ k ( n , m ) ,
    1 2 k + 1 γ k + 1 ( n , m ) ( a ) = a 2 2 n 1 γ k ( n 1 , m ) ( a ) + ( 1 a ) 2 2 m 1 γ k ( n , m 1 ) ( a )
    we can easily check that it does not work;
  • If we use the matrices M γ 2 and M γ 3 and if we write the recurrence that we found written in terms of γ k ( n , m ) ,
    ( 2 m 1 ) γ k ( n , m ) ( a ) + ( 2 n + 1 ) γ k ( n + 1 , m 1 ) ( a ) = ( 2 k 1 ) γ k 1 ( n , m 1 ) ( a ) + ( 2 n + 2 m 2 k 1 ) γ k ( n , m 1 ) ( a ) ,
    we easily check that this works perfectly; for instance, when n = 0 , m = 3 , and k = 1 , we obtain
    γ 1 ( 1 , 2 ) ( a ) + 5 γ 1 ( 0 , 3 ) ( a ) γ 0 ( 0 , 2 ) ( a ) 3 γ 1 ( 0 , 2 ) ( a ) = ( 6 a 3 14 a 2 + 13 a 4 ) a 2 ( 2 a 1 ) 5
    + 5 2 ( 15 a 3 20 a 2 + 10 a 2 ) a ( a 1 ) 5 ( 2 a 1 ) 5 ( 3 a 2 3 a + 1 ) a ( 2 a 1 ) 3 3 a ( a 1 ) ( 2 a 1 ) 2 = 0 .
and we can announce the following conjecture.
Conjecture The coefficients γ k ( n , m ) of the matrix M γ fulfil the same recurrence (13): For any integers n , m , n 1 , m 1 and for any 1 k n + m , we have the following results:
( 2 m 1 ) γ k ( n , m ) ( a ) + ( 2 n + 1 ) γ k ( n + 1 , m 1 ) ( a ) = ( 2 k 1 ) γ k 1 ( n , m 1 ) ( a ) + ( 2 n + 2 m 2 k 1 ) γ k ( n , m 1 ) ( a ) ,
taking into account either γ 1 ( n , m 1 ) ( a ) = 0 or γ n + m ( n , m 1 ) ( a ) = 0 .
Conjecture All the coefficients of the matrix M γ are positive for a [ 0 , 1 ] .
Conjecture The expression
γ n n k , k : = ( 1 2 ) n ( 1 + a ) n k + 1 a n k F 1 2 ( [ 3 2 , 1 k ] , [ n + 3 2 ] , ( 1 + 2 a ) 2 ) Γ ( 2 n 1 ) ( 1 2 ) k ( 1 2 ) n k n ( 1 + 2 a ) 2 n 1 Γ ( n ) 2 2 2 k 2
fulfills
k = 0 n β k ( n l , l ) ( a ) . γ n n k , k = 1 2 δ l , 0 , 0 k , l n .
Please use the short set of Maple instructions presented in the Appendix F below.

6. Discussion

In this paper, many conjectures, mathematical statements that have not yet been rigorously proved, are presented. These conjectures arose when I noticed a pattern that holds true for many cases. However, just because a pattern holds true for many cases does not mean that the pattern will hold true for all cases. To notice a pattern that holds true is not always easy. Conjectures must be proved for the mathematical observation to be fully accepted. When a conjecture is rigorously proved, it becomes a theorem and will have many applications, maybe not in the near future but surely in the far future.

7. Conclusions

In 2008, Berg and Vignat were the first who worked on the linearization coefficients of Bessel polynomials, and since then, Atia and co-authors continued in this direction. This is an expository paper together with some new results and conjectures.

Funding

The researcher would like to thank the Deanship of Graduate Studies and Scientific Research at Qassim University for financial support (QU-APC-2024-9/1).

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declares no conflicts of interest.

Appendix A

> restart;
> A : = n u > G A M M A ( n u + 1 / 2 ) / ( G A M M A ( n u ) G A M M A ( 1 / 2 ) ) ;
> I n m : = ( x , n , m ) > i n t ( 1 / ( ( 1 + ( y x ) 2 ) ( n + 1 ) ( y 2 + 1 ) ( m + 1 ) ) ,   y = i n f i n i t y . . i n f i n i t y ) ;
> J n m : = ( x , n , m ) > ( 1 / 2 ) ( s u m ( A ( k + 1 / 2 ) B e t a n m k a ( n , m , k , 1 / 2 ) / ( A ( n + 1 / 2 )   A ( m + 1 / 2 ) ( 1 + ( 1 / 4 ) x 2 ) k + 1 ) , k = m i n ( n , m ) . . n + m ) ) ;
> f a c t o r ( s i m p l i f y ( I n m ( x , 1 , 2 ) ) ) ;
π ( 3 x 4 + 56 x 2 + 560 ) 8 ( x 8 + 16 x 6 + 96 x 4 + 256 x 2 + 256 )
> f a c t o r ( s i m p l i f y ( J n m ( x , 1 , 2 ) ) ) ;
3 π B e t a n m k a ( 1 , 2 , 1 , 1 / 2 ) 16 ( 1 + ( 1 / 4 ) x 2 ) 2 + π B e t a n m k a ( 1 , 2 , 1 , 1 / 2 ) 4 ( 1 + ( 1 / 4 ) x 2 ) 3 + 3 π B e t a n m k a ( 1 , 2 , 1 , 1 / 2 ) 10 ( 1 + ( 1 / 4 ) x 2 ) 4
> B e t a n m k a : = ( n , m , k , a ) > a 2 n + m k ( 1 a ) n + k p o c h h a m m e r ( 1 / 2 , k )
p o c h h a m m e r ( 1 / 2 , n + m k ) ( s u m ( ( a ) j
b i n o m i a l ( m + n + 1 , 2 n + 2 m 2 k j ) b i n o m i a l ( m + k + j , j ) , j = 0 . 2 m + 2   n 2 k ) ) / ( p o c h h a m m e r ( 1 / 2 , n ) p o c h h a m m e r ( 1 / 2 , m ) ) ;
> f a c t o r ( s i m p l i f y ( 3 π B e t a n m k a ( 1 , 2 , 1 , 1 / 2 ) 16 ( 1 + ( 1 / 4 ) x 2 ) 2 + π B e t a n m k a ( 1 , 2 , 1 , 1 / 2 ) 4 ( 1 + ( 1 / 4 ) x 2 ) 3 + 3 π B e t a n m k a ( 1 , 2 , 1 , 1 / 2 ) 10 ( 1 + ( 1 / 4 ) x 2 ) 4 ) ) ;
π ( 3 x 4 + 56 x 2 + 560 ) 8 ( x 8 + 16 x 6 + 96 x 4 + 256 x 2 + 256 )

Appendix B

> restart;
> d j m : = ( m , j ) > 2 2 m ( s u m ( 2 i b i n o m i a l ( 2 m 2 i , m i ) b i n o m i a l ( m + i , m )
b i n o m i a l ( i , j ) , i = j . . m ) ) ;
> P m : = m > s u m ( d j m ( m , j ) x j , j = 0 . . m ) : ; P m ( m )
j = 0 m ( 2 2 m 2 j b i n o m i a l ( 2 m 2 j , m j ) b i n o m i a l ( m + j , m )
h y p e r g e o m ( [ m + j , m + 1 + j ] , [ j m + 1 / 2 ] , 1 / 2 ) x j )
> e x p p : = m > π ( s u m ( 2 2 m 2 j b i n o m i a l ( 2 m 2 j , m j ) b i n o m i a l ( m + j , m )
h y p e r g e o m ( [ m + j , m + 1 + j ] , [ j + 1 / 2 m ] , 1 / 2 ) x j , j = 0 . . m ) ) / ( 2 m + 3 / 2 ( x + 1 ) m + 1 / 2 ) ;
s i m p l i f y ( e x p p ( 2 ) )
3 π ( 4 x 2 + 10 x + 7 ) 2 128 ( x + 1 ) 5 / 2
> a s s u m e ( x > 1 ) ; i n t m : = m > i n t ( 1 / ( s 4 + 2 s 2 x + 1 ) 1 + m , s = 0 . . i n f i n i t y ) ; s i m p l i f y ( i n t m ( 2 ) )
3 π ( 4 x 2 + 10 x + 7 ) 2 128 ( x + 1 ) 5 / 2 .

Appendix C

> restart;
> A : = n u > G A M M A ( n u + 1 / 2 ) / ( G A M M A ( n u ) G A M M A ( 1 / 2 ) ) ;
> B e t a N 3 t h i r d : = ( n , m , p ) > s u m ( f a c t o r i a l ( 2 l ) ( 1 ) k f a c t o r i a l ( k + l + 1 ) ( s u m ( p o c h h a m m e r ( p , k + l j ) ( s u m ( p o c h h a m m e r ( n , i ) p o c h h a m m e r ( m , j i ) / ( p o c h h a m m e r ( 2 n , i ) p o c h h a m m e r ( 2 m , j i ) f a c t o r i a l ( i ) f a c t o r i a l ( j i ) ) , i = m a x ( 0 , j m ) . . m i n ( j , n ) ) ) / ( p o c h h a m m e r ( 2 p , k + l j ) f a c t o r i a l ( k + l j ) ) , j = m a x ( 0 , k + l p ) . . m i n ( k + l , n + m ) ) ) / ( 3 k + l + 1 f a c t o r i a l ( l ) f a c t o r i a l ( k ) f a c t o r i a l ( l + 1 k ) ) , k = 0 . . m i n ( n + m + p l , l + 1 ) ) ;
> i n t d o u b l e : = ( n , m , p ) > i n t ( i n t ( P h i ( n + 1 / 2 ) P h i ( m + 1 / 2 ) P h i ( p + 1 / 2 ) / ( ( 1 + x 2 ) n + 1 ( 1 + y 2 ) m + 1 ( 1 + ( z x y ) 2 ) p + 1 ) , x = i n f i n i t y . . i n f i n i t y ) , y = i n f i n i t y . . i n f i n i t y ) ;
> s u m m a t i o n : = ( n , m , p ) > s u m ( P h i ( l + 1 / 2 ) B e t a N 3 t h i r d ( n , m , p , l ) / ( 1 + ( 1 / 9 ) z 2 ) ( l + 1 ) , l = m i n ( n , m , p ) . . n + m + p ) ;
> f a c t o r ( i n t d o u b l e ( 1 , 2 , 3 ) ) ;
2 ( 15 z 10 + 1015 z 8 + 30942 z 6 + 555174 z 4 + 6341571 z 2 + 45684243 ) 15 π ( z 2 + 9 ) 7
> f a c t o r ( s u m m a t i o n ( 1 , 2 , 3 ) ) ;
2 ( 15 z 10 + 1015 z 8 + 30942 z 6 + 555174 z 4 + 6341571 z 2 + 45684243 ) 15 π ( z 2 + 9 ) 7

Appendix D

> restart;
> tracep := p > ( 1 + 2 p ) a 2 p p ( 1 + 2 p ) a 2 p 1 + ( 1 + a ) ( 1   ( 1 + a ) 2 p ) / ( a ( a 2 ) ) + ( 1 / 2 ) ( 1 + 2 a ) 1 + 2 p ( 1 + a ) 2 / ( a a ( a 2 ) ) ( 1 + a )   ( 1 + 2 p ) a 2 p / a ( 1 / 2 ) ( 1 + a ) 2 / ( a 2 ( a 2 ) ) + ( 1 / 2 ) ( ( 1 + a ) 2 p + 2   ( a + 1 ) 1 + 2 p ( 1 + a ) ) / ( a 2 ( a 2 ) ) + ( 1 + 2 p ) ( p 1 ) a 2 p 1 + ( 1 / 2 )   ( ( 1 2 a ) 1 + 2 p + 1 2 a + ( 1 + a ) 1 + 2 p ( a + 1 ) + ( a + 1 ) 1 + 2 p ( 1 + a ) ) / ( a 2 ( a 2 ) ) :
> factor(simplify(2*a*tracep(3)-1));
( 1 + 2 a ) 7

Appendix E

> restart;
> K u K a u : = ( n , a ) e x p ( ( 1 a ) u ) B e s s e l K ( n + 1 / 2 , u ) / G A M M A ( n + 1 / 2 )   ( s u m ( s q r t ( a ) c n k ( n , i , 1 / a ) ( 2 / u ) ( n i ) a i B e s s e l K ( i + 1 / 2 , a u ) / G A M M A ( i + 1 / 2 ) , i = 0 . . n ) ) ;
K u K a u : = ( n , a ) e ( 1 a ) u B e s s e l K ( n + 1 / 2 , u ) Γ ( n + 1 / 2 )
i = 0 n a c n k ( n , i , 1 / a ) ( 2 u ) n i a i B e s s e l K ( i + 1 / 2 , a u ) Γ ( i + 1 / 2 )
> s i m p l i f y ( K u K a u ( 2 , a 2 ) ) ;
( 2 / 3 ) s q r t ( 2 ) e x p ( a 2 u ) ( c n k ( 2 , 2 , 1 / a 2 ) a 4 u 2 + 3 c n k ( 2 , 1 , 1 / a 2 ) a 2 u
+ 3 c n k ( 2 , 2 , 1 / a 2 ) a 2 u u 2 + 3 c n k ( 2 , 0 , 1 / a 2 ) + 3 c n k ( 2 , 1 , 1 / a 2 )
+ 3 c n k ( 2 , 2 , 1 / a 2 ) 3 u 3 ) / u 5 / 2
> c n k : = ( n , k , a ) > a k b i n o m i a l ( n , k ) / b i n o m i a l ( 2 n , 2 k ) s u m ( b i n o m i a l ( n + 1 , k + 1 r ) b i n o m i a l ( n k 1 , r 1 ) ( 1 a ) r , r = 0 . . m i n ( n k , k + 1 ) ) ;
c n k : = ( n , k , a ) >
a k b i n o m i a l ( n , k ) r = 0 m i n ( n k , k + 1 ) b i n o m i a l ( n + 1 , k + 1 r ) b i n o m i a l ( n k 1 , r 1 ) ( 1 a ) r ) b i n o m i a l ( 2 n , 2 k )
> s i m p l i f y ( ( 2 / 3 ) s q r t ( 2 ) e x p ( a 2 u ) ( c n k ( 2 , 2 , 1 / a 2 ) a 4 u 2 + 3 c n k ( 2 , 1 , 1 / a 2 )   a 2 u + 3 c n k ( 2 , 2 , 1 / a 2 ) a 2 u u 2 + 3 c n k ( 2 , 0 , 1 / a 2 ) + 3 c n k ( 2 , 1 , 1 / a 2 ) + 3 c n k ( 2 , 2 , 1 / a 2 ) 3 u 3 ) / u 5 / 2 )

Appendix F

> b e t a : = ( m , n , k ) > a ( 2 n + m k ) ( 1 a ) ( n + k ) p o c h h a m m e r ( 1 / 2 , k )
p o c h h a m m e r ( 1 / 2 , n + m k ) ( s u m ( ( a ) ( j ) b i n o m i a l ( m + n + 1 , 2 n + 2 m 2 k j ) b i n o m i a l ( m + k + j , j ) , j = 0 . 2 m + 2 n 2 k ) )
/ ( p o c h h a m m e r ( 1 / 2 , n ) p o c h h a m m e r ( 1 / 2 , m ) )
> g a m m a : = ( n , k ) > p o c h h a m m e r ( 1 / 2 , n ) 2 ( 2 n 2 k ) ( 1 + a ) ( n k + 1 ) a ( n k )
h y p e r g e o m ( [ 3 / 2 , 1 k ] , [ n + 3 / 2 ] , ( 1 + 2 a ) 2 ) G A M M A ( 2 n 1 )
/ ( p o c h h a m m e r ( 1 / 2 , k ) p o c h h a m m e r ( 1 / 2 , n k ) n ( 1 + 2 a ) ( 2 n 1 ) G A M M A ( n ) 2 2 ( 2 n 2 ) )
> s i m p l i f y ( b e t a ( 3 , 0 , 0 ) g a m m a ( 3 , 0 ) + b e t a ( 3 , 0 , 1 ) g a m m a ( 3 , 1 ) + b e t a ( 3 , 0 , 2 )
g a m m a ( 3 , 2 ) + b e t a ( 3 , 0 , 3 ) g a m m a ( 3 , 3 ) ) = 1 / 2 .

References

  1. Grosswald, E. Bessel Polynomials, Lecture Notes in Mathematics 698; Springer: New York, NY, USA, 1978. [Google Scholar]
  2. Berg, C.; Vignat, C. Linearization coefficients of Bessel polynomials and properties of Student t-Distributions. Constr. Approx. 2008, 27, 15–32. [Google Scholar] [CrossRef]
  3. Berg, C.; Vignat, C. Derivation of an integral of Boros and Moll via convolution of Student t-densities. Ramanujan J. 2012, 27, 147–150. [Google Scholar] [CrossRef]
  4. Boros, G.; Moll, V.H. An integral hidden in Gradshteyn and Ryzhik. J. Comput. Appl. Math. 1999, 106, 361–368. [Google Scholar] [CrossRef]
  5. Atia, M.J.; Zeng, J. An explicit formula for the linearization coefficients of Bessel polynomials. Ramanujan J. 2012, 28, 211–221. [Google Scholar] [CrossRef]
  6. Andrews, G.; Askey, R.; Roy, R. Special Functions, Encyclopedia of Mathematics and Its Applications 71; Cambridge University Press: Cambridge, UK, 1999. [Google Scholar]
  7. BenAbdallah, M.; Atia, M.J. Derivation of a double integral via convolution of Student t-densities. Ramanujan J. 2019, 49, 75–86. [Google Scholar] [CrossRef]
  8. Gradshteyn, I.S.; Ryzhik, I.M. Tables of Integrals, Series and Products, 6th ed.; Academic Press: San Diego, CA, USA, 2000. [Google Scholar]
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

Atia, M.J. On the Inverse of the Linearization Coefficients of Bessel Polynomials. Symmetry 2024, 16, 737. https://doi.org/10.3390/sym16060737

AMA Style

Atia MJ. On the Inverse of the Linearization Coefficients of Bessel Polynomials. Symmetry. 2024; 16(6):737. https://doi.org/10.3390/sym16060737

Chicago/Turabian Style

Atia, Mohamed Jalel. 2024. "On the Inverse of the Linearization Coefficients of Bessel Polynomials" Symmetry 16, no. 6: 737. https://doi.org/10.3390/sym16060737

APA Style

Atia, M. J. (2024). On the Inverse of the Linearization Coefficients of Bessel Polynomials. Symmetry, 16(6), 737. https://doi.org/10.3390/sym16060737

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