Next Article in Journal
Tail Dependence and Risk Spillover from the US to GCC Banking Sectors
Next Article in Special Issue
Some Relationships for the Generalized Integral Transform on Function Space
Previous Article in Journal
Optimal Exploitation of a General Renewable Natural Resource under State and Delay Constraints
Previous Article in Special Issue
Transformations of the Hypergeometric 4F3 with One Unit Shift: A Group Theoretic Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Laguerre-Type Exponentials, Laguerre Derivatives and Applications. A Survey

by
Paolo Emilio Ricci
Department of Mathematics, International Telematic University UniNettuno, Corso Vittorio Emanuele II, 39, 00186 Roma, Italy
Mathematics 2020, 8(11), 2054; https://doi.org/10.3390/math8112054
Submission received: 19 October 2020 / Revised: 12 November 2020 / Accepted: 13 November 2020 / Published: 18 November 2020
(This article belongs to the Special Issue Special Functions with Applications to Mathematical Physics)

Abstract

:
Laguerrian derivatives and related autofunctions are presented that allow building new special functions determined by the action of a differential isomorphism within the space of analytical functions. Such isomorphism can be iterated every time, so that the resulting construction can be re-submitted endlessly in a cyclic way. Some applications of this theory are made in the field of population dynamics and in the solution of Cauchy’s problems for particular linear dynamical systems.

1. Introduction

This survey article is dedicated to a topic that has received little attention in the past, and therefore seems not to be very well known by the mathematical community.
Recently the role of the Laguerre derivative was considered in a few papers.
In [1], the authors introduce an interesting application of Wright functions of the first kind to solve fractional ordinary differential equations, with variable coefficients, generalizing the Bessel-type equations.
In [2], the authors use the same tool in Combinatorics, a completely different area, and in [3] an operational approach to the subject has been examined, in the framework of Clifford algebras.
Actually, in past time, the Laguerre-type exponentials and the related Laguerre derivative were introduced and studied in several articles (see [4,5,6,7,8,9,10,11,12,13,14,15]) and applications to Special functions, have been obtained. In particular, Laguerre-type functions of Bessel, Appell, Bell and multivariate functions were defined.
The operator D x D = D + x D 2 determines a linear differential isomorphism, acting onto the space of analytic functions of the x variable. By using this isomorphism, a sort of parallel structure is created within this space, in such a way that the differentiation properties have their counterpart, which can be immediately derived.
Furthermore, iterations of the Laguerre derivative can be defined, so that this parallelism with the space of analytic functions can be iterated too, in an endless way.
Therefore, a cyclic construction is created within the space that repeats the same structure at a higher level of differentiation order. It is one of the great cycles that sometimes occur within mathematical theories: for example, in Number theory the Fibonacci numbers F n with Fibonacci indexes constitute a higher sequence of Fibonacci numbers which still satisfies the same recursion, i.e., F F n + 2 = F F n + 1 + F F n , and this property can be iterated at infinity.
However the operators D L = D x D and its iterates as D n L = D x D x D x D x D are not completely new, since they can be considered to be particular cases of the hyper-Bessel differential operators when α 0 = α 1 = = α n = 1 (the special case considered in operational calculus by Ditkin and Prudnikov [16]). In general, the Bessel-type differential operators of arbitrary order n were introduced by Dimovski, in 1966 [17] and later called by Kiryakova hyper-Bessel operators, because are closely related to their eigenfunctions, called hyper-Bessel by Delerue [18], in 1953. These operators were studied in 1994 by Kiryakova in her book [19] (Ch. 3).
Since the Laguerrian exponentials on the positive semi-axis of the abscissas are convex increasing functions, with a growth lower than the exponential one, in Section 7 a natural application was made in the context of population dynamics.
Laguerre-type linear dynamical systems were also considered in Section 8.

2. The Laguerre Derivative and the Relevant Exponentials

The Laguerre derivative, is defined by
D L : = D x D = D + x D 2 ,
where D = D x = d / d x .
It is an interesting operator. In fact, as the exponential function e a x (a constat) is an eigenfunction of the derivative operator D = D x , i.e.,
D e a x = a e a x ,
equally the function
e 1 ( x ) : = k = 0 x k ( k ! ) 2 = C 0 ( x ) ,
where C 0 ( x ) is the Tricomi function of order zero, is an eigenfunction of the Laguerre derivative D L , since:
D L e 1 ( a x ) = a e 1 ( a x ) .
The proof easily follows, by noting that:
D L e 1 ( a x ) = D + x D 2 k = 0 a k x k ( k ! ) 2 = 1 p t 26 p t = k = 1 k + k ( k 1 ) a k x k 1 ( k ! ) 2 = k = 1 k 2 a k x k 1 ( k ! ) 2 = 1 p t 26 p t = a k = 0 a k x k ( k ! ) 2 = a e 1 ( a x ) .
For this reason, the function e 1 ( x ) is called the Laguerre-type exponential (of order 1).
In preceding articles, the role of the Laguerre derivative, in connection with the monomiality principle—an important technique introduced by G. Dattoli [20]—and its application to the multidimensional Hermite (Hermite-Kampé de Fériet or Gould-Hopper polynomials, see [21,22,23]) or Laguerre polynomials [14,24], has been shown.
The above technique can be iterated, producing Laguerre classes of exponential-type functions, of higher order, called L-exponentials, and the relevant L-circular, L-hyperbolic, L-Gaussian functions (see [4]).
Similar generalized hypergeometric functions, called trigonometric/Bessel type, exponential/ confluent type and Gauss/Beta-distribution, can be found in a book by Kiryakova [19] and also in [25].
Before going on, we notice that the Laguerre derivative verifies [26]:
D x D n = D n x n D n ,
an equation which can be easily proven by recursion.

2.1. L-Exponentials of Higher Order

We consider the operator:
D 2 L : = D x D x D = D x D + x 2 D 2 = D + 3 x D 2 + x 2 D 3 ,
and the function:
e 2 ( x ) : = k = 0 x k ( k ! ) 3 .
The following theorem holds:
Theorem 1.
The function e 2 ( a x ) is an eigenfunction of the operator D 2 L , i.e.,
D 2 L e 2 ( a x ) = a e 2 ( a x )
The proof (see [4]) depends on the identity: k + 3 k ( k 1 ) + k ( k 1 ) ( k 2 ) = k 3 , so that, it can be recognized that the coefficients of the combination in Equation (7) are the Stirling numbers of the second kind, S ( 3 , 1 ) , S ( 3 , 2 ) , S ( 3 , 3 ) , (see [27], and [28] (p. 835 for an extended table)).
In general, we can state the following theorem:
Theorem 2.
The function
e n ( x ) : = k = 0 x k ( k ! ) n + 1 .
is an eigenfunction of the operator
D n L : = D x D x D x D = D x D + x 2 D 2 + + x n D n = = S ( n + 1 , 1 ) D + S ( n + 1 , 2 ) x D 2 + + S ( n + 1 , n + 1 ) x n D n + 1 ,
i.e., for every constant a it results:
D n L e n ( a x ) = a e n ( a x ) .
Remark 1.
The above results show that, for every positive integer n, we can define a Laguerre-exponential function, satisfying an eigenfunction property, which is an analog of the elementary property (2) of the exponential. The function e n ( x ) reduces to the exponential function when n = 0 , so that we put by definition:
e 0 ( x ) : = e x , D 0 L : = D .
Obviously, D 1 L : = D L .
Examples of the L-exponential functions are given in Figure 1.

2.2. L-Circular and L-Hyperbolic Functions

Starting from the equation
e 1 ( i x ) = h = 0 ( 1 ) h x 2 h ( ( 2 h ) ! ) 2 + i h = 0 ( 1 ) h x 2 h + 1 ( ( 2 h + 1 ) ! ) 2 ,
we can define the 1 L -circular functions as follows
cos 1 ( x ) : = e 1 ( i x ) = h = 0 ( 1 ) h x 2 h ( ( 2 h ) ! ) 2 ,
sin 1 ( x ) : = e 1 ( i x ) = h = 0 ( 1 ) h x 2 h + 1 ( ( 2 h + 1 ) ! ) 2 ,
so that we find the Euler-type formulas
cos 1 ( x ) = e 1 ( i x ) + e 1 ( i x ) 2 , sin 1 ( x ) = e 1 ( i x ) e 1 ( i x ) 2 i ,
Recalling Equation (6), we find the result:
Theorem 3.
The 1 L -circular functions (14) and (15) are solutions of the differential equation
D L 2 v + v = D 2 x 2 D 2 v + v = 0 .
The above results hold even for the generalized case.
Write the n L -exponential in the form:
e n ( i x ) = h = 0 ( 1 ) h x 2 h ( ( 2 h ) ! ) n + 1 + i h = 0 ( 1 ) h x 2 h + 1 ( ( 2 h + 1 ) ! ) n + 1 .
Then we can define the n L -circular functions by putting
Definition 1.
cos n ( x ) : = e n ( i x ) = h = 0 ( 1 ) h x 2 h ( ( 2 h ) ! ) n + 1 ,
sin n ( x ) : = e n ( i x ) = h = 0 ( 1 ) h x 2 h + 1 ( ( 2 h + 1 ) ! ) n + 1 ,
and we find again the Euler-type formulas:
cos n ( x ) = e n ( i x ) + e n ( i x ) 2 , sin n ( x ) = e n ( i x ) e n ( i x ) 2 i .
Theorem 3 becomes, in general:
Theorem 4.
The n L -circular functions (18) and (19) are solutions of the differential equation
D n L 2 v + v = 0 .
and satisfy the conditions:
cos n ( 0 ) = 1 , sin n ( 0 ) = 0 .
Furthermore, we find:
Theorem 5.
The n L -circular functions satisfy
D n L cos n ( x ) = sin n ( x ) , D n L sin n ( x ) = cos n ( x ) .
Examples of the L-circular functions are given in Figure 2 and Figure 3.
In a similar way we can define the n L -hyperbolic functions, putting
cosh n ( x ) : = h = 0 x 2 h ( ( 2 h ) ! ) n + 1 ,
sinh n ( x ) : = h = 0 x 2 h + 1 ( ( 2 h + 1 ) ! ) n + 1 ,
and the formulas analogues of that of the circular functions are easily derived (see [4]).
All the eigenfunctions e 1 ( x ) , e 2 ( x ) , , e n ( x ) can be expressed as generalized hypergeometric functions p F q , [29], namely: e 1 ( x ) = 0 F 1 ( x ) , e 2 ( x ) = 0 F 2 ( x ) , ..., e n ( x ) = 0 F n ( x ) . In practice, starting from the Bessel function e 1 ( x ) , all these eigenfunctions are special cases of the hyper-Bessel functions of Delerue [18], which are shown to be eigenfunctions of Dimovski’s operators mentioned above.
Naturally, the cos n , sin n functions, in Equations (19) and (20), and their hyperbolic variants, are special cases of the trigonometric type generalized hypergeometric functions considered in the Kiryakova book [19].

3. The Isomorphism T x and Its Iterations

It was previously noted (see e.g., [14]) that, in the space A x of analytic functions, it is possible to define an isomorphism T x that preserves the differentiation properties, by means of correspondence:
D D L , x · D x 1 ,
where
D x 1 f ( x ) = 0 x f ( ξ ) d ξ , D x n f ( x ) = 1 ( n 1 ) ! 0 x ( x ξ ) n 1 f ( ξ ) d ξ ,
so that
T x ( x n ) = D x n ( 1 ) = 1 ( n 1 ) ! 0 x ( x ξ ) n 1 d ξ = x n n ! .
It is worth noting that this kind of isomorphism is widely used in operational calculus and differential equations also under the name of Transmutation or Similarity operator, since it transforms one operator into another, and eigenfunctions into each other.
In fact, in such an isomorphism we have the correspondences:
  • The exponential function is transformed into the function e 1 ( x ) , since
    T x ( e x ) = k = 0 T x ( x k ) k ! = k = 0 x k ( k ! ) 2 = e 1 ( x ) .
  • The Hermite polynomial H n ( 1 ) ( x , y ) : = ( x y ) n becomes the Laguerre polynomial
    L n ( x , y ) : = n ! r = 0 n ( 1 ) r y n r x r ( n r ) ! ( r ! ) 2
    and by using the monomiality principle we can prove thate all the relations valid in the polynomial space still hold after the substitutions stated in Equation (23).
Furthermore, an iterative application of Equation (23) gives in sequence the functions e 1 ( x ) , e 2 ( x ) , e 3 ( x ) , .
We have, for example:
T x 2 ( e x ) = k = 0 T x ( x k ) ( k ! ) 2 = k = 0 x k ( k ! ) 3 = e 2 ( x ) ,
and so on.
We already noticed that the isomorphism connected with the Laguerre derivative can be iterated as many times as we wish.
Correspondently, the derivative operator is transformed into
D L = D x D , D 2 L = D L D x 1 D L = D x D x D , D 3 L = D L D x 1 D L D x 1 D L = D x D x D x D , ,
and so on.
We can conclude that the L-exponentials (and the relevant L-circular and L-hyperbolic functions) are determined by an iterative application of the considered differential isomorphism.

4. Examples of Laguerre-Type Problems

4.1. L-Diffusion Equations

Theorem 6.
For any fixed integer n , consider the problem (see [4] (Theorem 5.1)):
D n L S ( x , t ) = t S ( x , t ) , in the half plane t > 0 , S ( 0 , t ) = s ( t ) ,
with analytic boundary condition s ( t ) .
The operational solution of problem (27) is given by:
S ( x , t ) = e n x t s ( t ) = k = 0 x k ( k ! ) n + 1 d k d t k s ( t )
Representing s ( t ) = k = 0 a k t k , from Equation (28) we find, in particular:
S ( x , 0 ) = k = 0 a k x k ( k ! ) n .
Please note that the operational solution becomes an effective solution whenever the series in Equation (28) is convergent. The validity of this condition depends on the growth of the coefficients a k of the boundary data s ( t ) , but it is usually satisfied in physical problems.
More general problems are shown in [4,10], where evolution problems related to an operator of the type
D p 1 x q 1 D p 2 x q 2 D p r x q r D p r + 1 ,
where p 1 , p 2 , , p r + 1 ; q 1 , q 2 , q r are fixed integers, have been considered.
An operational solution of the problem
D p 1 x q 1 D p 2 x q 2 D p r x q r D p r + 1 S ( x , t ) = D t S ( x , t ) , in the half plane t > 0 ,
with suitable initial conditions have been determined, in terms of the eigenfunctions of the same operator.
Remark 2.
Please note that the above operators generalize the subsequent Laguerre-type derivatives, since they are written as:
D n L r = ( D x D x D x D ) r ( n + 1 ) D e r i v a t i v e s = D r x r D r x r D r x r D r ,
which is an equation extending(6).
The operator (30) closely recalls the general case of hyper-Bessel B operators, in [17], since integers q 1 , q 2 , q r could be replaced by arbitrary real numbers, as are parameters α 0 , α 1 , , α n , considered in [17,30]. The solutions of the general differential equation B y ( x ) + λ y ( x ) = f ( x ) are given by Kiryakova et al. in [31].

4.2. L-Hyperbolic-Type Problems

Theorem 7.
Let Ω ^ x be a 2nd order differential operator with respect to the x variable, D n L : = D n L t the n L -derivative with respect to the t variable, and denote by ψ ( t ) and χ ( t ) two functions such that:
D n L ψ ( t ) = χ ( t ) , D n L χ ( t ) = ψ ( t ) 1 p t 20 p t ψ ( 0 ) = 1 , χ ( 0 ) = 0
then the abstract L-hyperbolic-type problem:
Ω ^ x 2 S ( x , t ) = D n L 2 S ( x , t ) , in the half plane t > 0 , S ( x , 0 ) = q ( x ) , D n L S ( x , t ) t = 0 = v ( x )
with analytic initial condition q ( x ) , v ( x ) , admits the operational solution(see [4], Theorem 5.3):
S ( x , t ) = ψ t Ω ^ x q ( x ) + χ t Ω ^ x w ( x ) ,
where w ( x ) : = Ω ^ x 1 v ( x ) .
Please note that conditions in (32) are satisfied, for any fixed integer n, assuming:
ψ ( x ) : = cosh n L ( x ) , χ ( x ) : = sinh n L ( x ) .

4.3. L-Elliptic-Type Problems

Theorem 8.
Let Ω ^ x be a 2nd order differential operator with respect to the x variable, D n L : = D n L y the n L -derivative with respect to the y variable, and denote by φ ( y ) and τ ( y ) two functions such that:
D n L φ ( y ) = τ ( y ) , D n L τ ( y ) = φ ( y ) 1 p t 20 p t φ ( 0 ) = 1 , τ ( 0 ) = 0
then the abstract L-elliptic-type problem:
Ω ^ x 2 S ( x , y ) + D n L 2 S ( x , y ) = 0 , in the half plane t > 0 , S ( x , 0 ) = q ( x ) ,
with analytic boundary condition q ( x ) , admits the operational solution (see [4], Theorem 5.4):
S ( x , y ) = φ y Ω ^ x q ( x ) .
Please note that conditions in (35) are satisfied, for any fixed integer n, assuming:
φ ( x ) : = cos n L ( x ) , τ ( x ) : = sin n L ( x ) .
Further examples of PDE’s problems involving the Laguerre derivatives can be found in [10,11].

5. Laguerre-Type Special Functions

5.1. Laguerre-Type Bessel Functions

The Laguerre-type Bessel functions, of order 1, (shortly L-Bessel functions), denoted by L J n ( x ) , are obtained substituting the exponential with the L-exponential e 1 ( x ) in the classic generating function, i.e., by putting
e 1 x 2 t 1 t = n = + L J n ( x ) t n .
We can derive the explicit expression by applying the isomorphism T x to both sides of the explicit expression of the Bessel functions, so that we find:
L J n ( x ) : = n = 0 ( 1 ) h x n + 2 h 2 n + 2 h h ! ( n + h ) ! ( n + 2 h ) ! .
We proved the results:
Theorem 9.
The L-Bessel functions L J n ( x ) satisfy the recurrence relation (see [8], Theorem 2.3):
D ^ x 1 L J n 1 ( x ) + L J n + 1 ( x ) = 2 n L J n ( x ) , L J n 1 ( x ) L J n + 1 ( x ) = 2 D ^ L L J n ( x ) .
Theorem 10.
The differential equation satisfied by the L-Bessel functions L J n ( x ) is(see [8], Theorem 2.5):
D ^ L 2 + D ^ x D ^ L n 2 D ^ x 2 + I ^ L J n ( x ) = 0 ,
where I ^ denotes the identity operator. This equation can be derived by applying the isomorphism T x to both sides of the differential equation of the ordinary first kind Bessel functions.

5.2. Laguerre-Type Hypergeometric Functions

By using the isomorphism technique it is possible to define in general Laguerre-type special functions, and in particular, the 1st order Laguerre-type hypergeometric functions.
In fact, starting from the Gauss’ hypergeometric equation:
x ( 1 x ) y + [ c ( a + b + 1 ) x ] y a b y = 0
and applying the isomorphism T x , we find the equation
x ( 1 x ) D L 2 y + [ c ( a + b + 1 ) x ] D L y a b y = 0 ,
that is:
[ x ( 1 x ) ] ( x 2 y i v + 4 x y + 2 y ) + [ c ( a + b + 1 ) x ] ( y + x y ) a b y = 0 .
The solution of Equation (38), corresponding to the Gauss’ hypergeometric equation F ( a , b , c ; x ) , is given by
L F ( a , b , c ; x ) = 1 + n = 1 a ( n ) b ( n ) c ( n ) x n ( n ! ) 2 ,
where the symbol a ( n ) denotes the rising factorial.
Of course the rth order Laguerre-type hypergeometric functions are obtained by applying to both sides of the hypergeometric equation the iterated isomorphism of order r, but the corresponding differential equation becomes more and more complicated as r increases.
The generalized hypergeometric functions have their 1st order Laguerre-type counterpart, which are given by:
L p F q ( a 1 , , a p ; b 1 , , b q , ; x ) = n = 0 a 1 ( n ) a p ( n ) b 1 ( n ) b q ( n ) x n ( n ! ) 2 ,
and those of higher order immediately follow.
Please note that the function in (41) can be viewed as a generalized hypergeometric function of the form p F q + 1 , by moving one of the n ! in the first fraction under the sum and considering Γ ( n + 1 ) / Γ ( 1 ) = 1 ( n ) = n ! = b q + 1 ( n ) as the ( q + 1 ) -th term.

5.3. Laguerre-Type Bell Polynomials

We first note that for the Laguerre derivative, the chain rule
d d t = d d x d x d t
becomes:
d d t t d d t = d d x d d t t d x d t , that is : D L t = d d x D L t x
and in general:
D n L t = d d x D n L t x .
The problem of constructing Bell polynomials can be extended in the natural way to the case of the Laguerre-type derivatives.
To this aim, we introduce the definition:
Definition 2.
The n th Laguerre-type Bell polynomial, denoted by r L Y n x ; f , g n , represents the n th r Laguerre-type derivative of the composite function f ( g ( t ) ) .
In [12] we showed that r L Y n can be expressed as a polynomial in the independent variable x, depending on f 1 , g 1 ; f 2 , g 2 ; ; f n , g n , in terms of the classical Bell polynomials.
According to Equation (6), the Leibniz rule, gives:
( D x D ) n = D n x n D n = k = 0 n n k D n k x n D n + k = 1 p t 26 p t = k = 0 n n k 2 ( n k ) ! x k D n + k = k = 0 n n ! k ! n k x k D n + k .
Therefore, the following representation formula for the Laguerre-type Bell polynomials, denoted by L Y n , holds:
Theorem 11.
The L Y n polynomials are expressed in terms of the ordinary Bell polynomials according to the equation (see [12], Theorem 4.1):
L Y n x ; f , g n = k = 0 n n ! k ! n k x k Y n + k f , g n + k .
The above results can be easily generalized, since
( D 2 L ) n = ( D x D x D ) n = D n x n D n x n D n = 1 p t 26 p t = k 1 = 0 n k 2 = 0 n n ! k 1 ! ( n + k 1 ) ! ( k 1 + k 2 ) ! n k 1 n k 2 x k 1 + k 2 D n + k 1 + k 2 .
In [12] even the general case of polynomials r L Y n Bell is considered, but we do not report here the equation which is a little more complicated.

6. The Multivariate Case

6.1. Laguerre-Type Appell Polynomials

In a preceding article [32] multivariate extensions of the Appell polynomials (including the Bernoulli and Euler cases) have been introduced, by means of the generating function [23]:
A ( t ) exp ( x t + y j ) = n = 0 R n ( j ) ( x , y ) t n n ! ,
where j is a fixed integer.
The application of the isomorphism T x , and its iterations allows defining new classes of multivariate special polynomials, the Laguerre-type Appell polynomials, and to build their main properties (recurrence relations, shift operators, differential equations, etc), in an easy and uniform way.
This has been achieved in [6] starting from generating functions of the type
A ( t ) e s ( x t ) e σ ( y t j ) = n = 0 R n ( j ) ( T x s ( x ) , T y σ ( y ) ) t n n !
where e s ( · ) and e σ ( · ) are Laguerre-type exponentials. Many properties of these functions have been derived, including recursions and differental equations.
The results obtained in this case are easily extended to the functions of r variables, since the technique works regardless of the number r.

6.2. Laguerre-Type Appell Series

We limit ourselves to the case of series in two variables, but the equations trivially extend to the general case. For | x | < 1 , | y | < 1 the double series
L F 1 ( a , b 1 , b 2 ; c ; x , y ) = m , n = 0 a ( m + n ) b 1 ( m ) b 2 ( n ) c ( m + n ) x m ( m ! ) 2 y n ( n ! ) 2
is the Laguerre-type Appell series, obtained by the classical one acting on it with the two isomorphisms T x and T y .
We avoid to consider further extension to the case of multivariate functions with several parameters, since they are trivially obtained.

7. Applications to Population Dynamics

7.1. Exponential and L-Exponential Models

In this section a possible application of the Laguerre derivative is recalled [9,13]. Since the L-exponentials for every x 0 are convex increasing functions, with a graph lower with respect to exp ( x ) , it is possible to use these function in the framework of population dynamics, as it seems that in some cases the growth of the exponential is too fast.
Consider the number N ( t ) of population individuals at time t and let N ( 0 ) = N 0 the initial number at time t = 0 .
In the Malthus model, the variation is assumed to be proportional to N ( t ) , i.e.,
d d t N ( t ) = r N ( t ) , N ( 0 ) = N 0 ,
where the growth rate r is a suitable constant.
The solution is given by the exponential function
N ( t ) = N 0 e r t .
Using the Laguerre derivative, the Laguerre-type Malthus reads:
d d t t d N d t = r N ( t ) i . e . d N d t + t d 2 N d t 2 = r N ( t ) ,
where r is a positive constant. Assuming the initial conditions
N ( 0 ) = N 0 , N ( 0 ) = N 1 = N 0 r ,
we find the solution
N ( t ) = N 0 e 1 ( r t ) = N 0 k = 0 + r k t k ( k ! ) 2 .
In this case the population growth increases according to the Laguerre exponential function e 1 ( x ) , so that the relevant increasing is slower with respect to the classical Malthus model.
In [9] it has been shown, with tables of data taken from real population dynamics, that the Laguerre-type Malthus model produces data closer to real population growth.

7.2. Logistic vs. L-logistic Model

Taking into account that the growth rate cannot be constant, since it depends on the environmental resources, Pierre Verhulst considered the so-called logistic model
d N d t = r 1 1 K N ( t ) N ( t ) , N ( 0 ) = N 0 ,
where r is called the intrinsic growth rate, and K denotes the environmental capacity.
The exact solution of this problem is given by
N ( t ) = N 0 K N 0 + ( K N 0 ) e r t ,
so that, if N 0 < K the solution is a function monotonically increasing to K, whereas, if N 0 > K , the solution is monotonically decreasing to K. In any case,
lim t N ( t ) = K ,
and the value N ( t ) = K is a stable equilibrium point for the logistic equation.
The Laguerre-logistic (shortly L-logistic) model is expressed by
N ( t ) + t N ( t ) = r N ( t ) 1 N ( t ) K , N ( 0 ) = N 0 , N ( 0 ) = N 1 .
Please note that if in the above equation N is small with respect to K, then N / K is close to 0 and consequently D t t D t N r N ( t ) .
If N K , then N / K 1 , and D t t D t N 0 .
The L-logistic equation cannot be solved explicitly, but numerically, using a Runge-Kutta method.
The behavior of the approximate solutions for the L-logistic model is shown in Figure 4. It is worth noting that the solution tends to the environmental capacity K by an oscillating behavior.
This is the main difference with respect to the ordinary logistic model, since in that case the solution was monotonically increasing or decreasing to K.
Similar results could be obtained by using the n L -derivatives, introducing suitable initial conditions which can be easily derived from the initial observations data.
Please note that as the order n increases, for x > 0 , the Laguerrian exponential attenuates its growth and for n it tends to assume the linear value 1 + x , so it can be used to model a growth as slow as it is needed.
Remark 3.
We recall that the oscillating asymptotic trend of solutions occurs in reality. For example, the classical experiment of G.F. Gause, relative to the protozoon paramecium shows such a typical behavior, represented in Figure 5. In this figure the true values, represented by a dotted line, are compared with the exponential trend of Malthus and with the logistic curve.

7.3. Modified L-Logistic Models

Many different models modifying the basic logistic model appeared in the literature: the Bernoulli, the modified logistic, the Gompertz, the Alee, and the Beverton-Holt models.
In [13] we considered the Laguerre-type version of all of them, showing that in all cases the oscillating asymptotic behavior of solutions takes the place of the monotonic one.
Instead, it was found that the model of Volterra-Lotka model is invariant under the action of the isomorphism T x , since the Laguerre derivative satisfies again the chain rule, according to Equation (42).

8. Laguerre-Type Linear Dynamical Systems

Let A be a r × r matrix and denote by u k , ( k = 1 , 2 , , r ) the invariants of A , i.e., the sum of principal minors (i.e., the elementary symmetric functions of the eigenvalues). The invariants of the matrix t A are given by u k ( t ) = t k u k , ( k = 1 , 2 , , r ) .
Consider the vectors
Z ( t ) = ( Z 1 ( t ) , , Z r ( t ) ) T Z 0 = ( Z 1 ( 0 ) , , Z r ( 0 ) ) T .
Then the solution of the linear dynamical system
Z ( t ) = A Z ( t ) , Z ( 0 ) = Z 0
writes [33]:
Z ( t ) = e t A Z 0 = h = 0 r 1 1 2 π i j = 0 r h 1 ( 1 ) j u j ( t ) γ e λ λ r h j 1 P ( λ , t ) d λ · t h Z 0 h ,
where P ( λ , t ) is the characteristic polynomial of the matrix t A and γ denotes a simple Jordan curve encircling all the eigenvalues of A . The choice of γ , without computing the eigenvalues, can be done by using the Gershgorin theorem.
In [15], a Laguerre-type version of the above classic result has been shown.
Consider the above r × r matrix, and the vectors
Z ( t ) = ( Z 1 ( t ) , , Z r ( t ) ) T Z 0 = ( Z 1 ( 0 ) , , Z r ( 0 ) ) T Z 0 = ( Z 1 ( 0 ) , , Z r ( 0 ) ) T = A · Z 0 Z 0 r 1 = ( Z 1 r 1 ( 0 ) , , Z r r 1 ( 0 ) ) T = A · Z 0 r 2 .
The following result holds (see [15] (Theorem 10)):
Theorem 12.
The Laguerre-type Cauchy problem for a homogeneous linear differential system
D L Z ( t ) = Z ( t ) + t Z ( t ) = A · Z ( t ) Z ( 0 ) = Z 0 Z 0 = A · Z 0 ,
has the solution:
Z ( t ) = e 1 ( t A ) Z 0 = h = 0 r 1 1 2 π i j = 0 r h 1 ( 1 ) j u j ( t ) γ e 1 ( λ ) λ r h j 1 P ( λ , t ) d λ · t h Z 0 h ,
where P ( λ , t ) and γ have been defined above.
The proof of this result is a straightforward application of the isomorphism T t . In [15] worked examples are reported.

9. Conclusions

The Laguerre derivative and the relevant Laguerre-type exponentials allow to associate, to any given integer n, a new class of special functions. This fact is obtained by exploiting the properties of an isomorphism, within the space of analytic functions, which acts in such a way as to preserve the differentiation properties. The successive iterations of this isomorphism produce a cyclic construction within the space that repeats the same structure at a higher level of the order of derivation.
Infinite many special functions can be defined in this way. A few of them have been presented explicitly, and the general technique to produce the others has been indicated.
This Survey has shown even possible applications of the Laguerrian derivative in the context of population dynamics and in the solution of Cauchy problems related to particular linear dynamical systems.

Funding

This research received no external funding.

Acknowledgments

The author is grateful to the referees, especially to the one who reported ties with the topics covered in the article and the hyper-Bessel operators, helping to expand the bibliography.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Garra, R.; Mainardi, F. Some applications of Wright functions in fractional differential equations. Rep. Math. Phys. 2020, in press. [Google Scholar]
  2. Penson, K.A.; Blasiak, P.; Horzela, A.; Solomon, A.I.; Duchamp, G.H.E. Laguerre-type derivatives: Dobiński relations and combinatorial identities. J. Math. Phys. 2009, 50, 3512. [Google Scholar] [CrossRef] [Green Version]
  3. Cação, I.; Falcão, M.I.; Malonek, H.R. Laguerre derivative and monogenic Laguerre polynomials: An operational approach. Math. Comput. Model. 2011, 53, 1084–1094. [Google Scholar] [CrossRef] [Green Version]
  4. Dattoli, G.; Ricci, P.E. Laguerre-type exponentials, and the relevant L-circular and L-hyperbolic functions. Georgian Math. J. 2003, 10, 481–494. [Google Scholar]
  5. Bernardini, A.; Dattoli, G.; Ricci, P.E. L-exponentials and higher order Laguerre polynomials. In Proceedings of the Fourth International Conference of the Society for Special Functions and their Applications (SSFA), Rajasthan University, Jaipur, India, 4–6 March 2003; pp. 13–26. [Google Scholar]
  6. Bretti, G.; Cesarano, C.; Ricci, P.E. Laguerre-type exponentials and generalized Appell polynomials. Comput. Math. Appl. 2004, 48, 833–839. [Google Scholar] [CrossRef] [Green Version]
  7. Bernardini, A.; Bretti, G.; Ricci, P.E. Laguerre-type exponentials, multidimensional special polynomials and applications. Lect. Notes TICMI 2004, 5, 1–28. [Google Scholar]
  8. Cesarano, C.; Germano, B.; Ricci, P.E. Laguerre-type Bessel functions. Integral Transforms Spec. Funct. 2005, 16, 315–322. [Google Scholar] [CrossRef]
  9. De Andreis, S.; Ricci, P.E. Modelling population growth via Laguerre-type exponentials. Math. Comput. Model. 2005, 42, 1421–1428. [Google Scholar] [CrossRef]
  10. Dattoli, G.; He, M.X.; Ricci, P.E. Eigenfunctions of Laguerre-type operators and generalized evolution problems. Math. Comput. Model. 2005, 42, 1263–1268. [Google Scholar] [CrossRef]
  11. Maroscia, G.; Ricci, P.E. Laguerre-type BVP and generalized Laguerre polynomials. Integral Transforms Spec. Funct. 2006, 17, 577–590. [Google Scholar] [CrossRef]
  12. Natalini, P.; Ricci, P.E. Laguerre-type Bell polynomials. Int. J. Math. Math. Sci. 2006, 2006, 45423. [Google Scholar] [CrossRef]
  13. Bretti, G.; Ricci, P.E. Laguerre-type Special functions and population dynamics. Appl. Math. Comput. 2007, 187, 89–100. [Google Scholar] [CrossRef]
  14. Ricci, P.E.; Tavkhelidze, I. An introduction to operational techniques and special polynomials. J. Math. Sci. 2009, 157, 161–189. [Google Scholar] [CrossRef]
  15. Natalini, P.; Ricci, P.E. Laguerre-type linear dynamical systems. Ilir. J. Math. 2015, 4, 24–40. [Google Scholar]
  16. Ditkin, A.P.; Prudnikov, V.A. Transforms and Operational Calculus; Pergamon Press: Oxford, UK, 1965. [Google Scholar]
  17. Dimovski, I. Operational calculus for a class of differential operators. C. R. Acad. Bulg. Sci. 1966, 19, 1111–1114. [Google Scholar]
  18. Delerue, P. Sur le calcul symbolique à n variables et fonctions hyper-besseliennes (II). Ann. Soc. Sci. Brux. 1953, 3, 229–274. [Google Scholar]
  19. Kiryakova, V. Generalized Fractional Calculus and Applications; Pitman Research Notes in Mathematics Series 301; Longman: Harlow, UK, 1994. [Google Scholar]
  20. Dattoli, G. Hermite-Bessel and Laguerre-Bessel functions: A by-product ot the monomiality principle. In Proceedings of the Melfi School on Advanced Topics in Mathematics and Physics, Melfi, Italy, 9–12 May 1999; Cocolicchio, D., Dattoli, G., Srivastava, H.M., Eds.; Advanced Special Functions and Applications. Aracne Editrice: Rome, Italy, 2000; pp. 147–164. [Google Scholar]
  21. Appell, P.; Kampé de Fériet, J. Fonctions Hypergéométriques et Hypersphériques. Polyno^mes d’Hermite; Gauthier-Villars: Paris, France, 1926. [Google Scholar]
  22. Gould, H.W.; Hopper, A.T. Operational formulas connected with two generalizations of hermite polynomials. Duke Math. J. 1962, 29, 51–62. [Google Scholar] [CrossRef]
  23. Srivastava, H.M.; Manocha, H.L. A Treatise on Generating Functions; J. Wiley & Sons: New York, NY, USA, 1984. [Google Scholar]
  24. Dattoli, G.; Ricci, P.E.; Cesarano, C. On a class of polynomials generalizing the Laguerre family. J. Comput. Appl. Math. 2005, 3, 405–412. [Google Scholar]
  25. Kiryakova, V. All the special functions are fractional differintegrals of elementary functions. J. Phys. A Math. Gen. 1997, 30, 5085–5103. [Google Scholar] [CrossRef] [Green Version]
  26. Viskov, O.V. A commutative-like noncommutation identity. Acta Sci. Math. 1994, 59, 585–590. [Google Scholar]
  27. Riordan, J. An Introduction to Combinatorial Analysis; J. Wiley & Sons: Chichester, UK, 1958. [Google Scholar]
  28. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions, with Formulas, Graphs, and Mathematical Tables; Dover: New York, NY, USA, 1965. [Google Scholar]
  29. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Higher Transcendental Functions; McGraw-Hill: New York, NY, USA, 1953; Volume 1. [Google Scholar]
  30. Kiryakova, V. From the hyper-Bessel opertors of Dimovski to the generalized fractional calculus. Fract. Calc. Appl. Anal. 2014, 17, 977–1000. [Google Scholar] [CrossRef]
  31. Kiryakova, V.S.; McBride, A.C. Explicit solution of the nonhomogeneous hyper-Bessel differential equation. C. R. Acad. Bulgare Sci. 1993, 46, 23–26. [Google Scholar]
  32. Bretti, G.; Ricci, P.E. Multidimensional extensions of the Bernoulli and Appell Polynomials. Taiwan. J. Math. 2004, 8, 415–428. [Google Scholar] [CrossRef]
  33. Natalini, P.; Ricci, P.E. Avoiding higher matrix powers in the solution of linear dynamical systems. Atlantis Trans. Geom. Model. Math. 2017, 2, 117–128. [Google Scholar]
Figure 1. e 1 ( x ) , (green) and e 2 ( x ) , (red).
Figure 1. e 1 ( x ) , (green) and e 2 ( x ) , (red).
Mathematics 08 02054 g001
Figure 2. cos 1 ( x ) , (green) and sin 1 ( x ) , (red).
Figure 2. cos 1 ( x ) , (green) and sin 1 ( x ) , (red).
Mathematics 08 02054 g002
Figure 3. cos 2 ( x ) (green) and sin 2 ( x ) (red).
Figure 3. cos 2 ( x ) (green) and sin 2 ( x ) (red).
Mathematics 08 02054 g003
Figure 4. Solutions to the L-logistic model with N ( 0 ) = N ( 0 ) < K (on the left), N ( 0 ) = N ( 0 ) > K (on the right), K = 64 , r = 0.8 , T = 100 , Δ t = 0.1 .
Figure 4. Solutions to the L-logistic model with N ( 0 ) = N ( 0 ) < K (on the left), N ( 0 ) = N ( 0 ) > K (on the right), K = 64 , r = 0.8 , T = 100 , Δ t = 0.1 .
Mathematics 08 02054 g004
Figure 5. The behavior of growth in the Gause experiment.
Figure 5. The behavior of growth in the Gause experiment.
Mathematics 08 02054 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ricci, P.E. Laguerre-Type Exponentials, Laguerre Derivatives and Applications. A Survey. Mathematics 2020, 8, 2054. https://doi.org/10.3390/math8112054

AMA Style

Ricci PE. Laguerre-Type Exponentials, Laguerre Derivatives and Applications. A Survey. Mathematics. 2020; 8(11):2054. https://doi.org/10.3390/math8112054

Chicago/Turabian Style

Ricci, Paolo Emilio. 2020. "Laguerre-Type Exponentials, Laguerre Derivatives and Applications. A Survey" Mathematics 8, no. 11: 2054. https://doi.org/10.3390/math8112054

APA Style

Ricci, P. E. (2020). Laguerre-Type Exponentials, Laguerre Derivatives and Applications. A Survey. Mathematics, 8(11), 2054. https://doi.org/10.3390/math8112054

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