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 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 with Fibonacci indexes constitute a higher sequence of Fibonacci numbers which still satisfies the same recursion, i.e., , and this property can be iterated at infinity.
However the operators
and its iterates as
are not completely new, since they can be considered to be particular cases of the hyper-Bessel differential operators when
(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
where
.
It is an interesting operator. In fact, as the exponential function
(
a constat) is an eigenfunction of the derivative operator
, i.e.,
equally the function
where
is the Tricomi function of order zero, is an eigenfunction of the Laguerre derivative
, since:
The proof easily follows, by noting that:
For this reason, the function 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]:
an equation which can be easily proven by recursion.
2.1. L-Exponentials of Higher Order
We consider the operator:
and the function:
The following theorem holds:
Theorem 1. The function is an eigenfunction of the operator , i.e., The proof (see [
4]) depends on the identity:
, so that, it can be recognized that the coefficients of the combination in Equation (
7) are the
Stirling numbers of the second kind,
, (see [
27], and [
28] (p. 835 for an extended table)).
In general, we can state the following theorem:
Theorem 2. is an eigenfunction of the operator i.e., for every constant a it results: 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 reduces to the exponential function when , so that we put by definition: Obviously, .
Examples of the
L-exponential functions are given in
Figure 1.
2.2. L-Circular and L-Hyperbolic Functions
Starting from the equation
we can define the
-
circular functions as follows
so that we find the Euler-type formulas
Recalling Equation (
6), we find the result:
Theorem 3. The -circular functions (14)
and (15)
are solutions of the differential equation The above results hold even for the generalized case.
Write the
-exponential in the form:
Then we can define the
-
circular functions by putting
and we find again the Euler-type formulas:
Theorem 3 becomes, in general:
Theorem 4. The -circular functions (18)
and (19)
are solutions of the differential equation and satisfy the conditions: Furthermore, we find:
Theorem 5. The -circular functions satisfy In a similar way we can define the
-hyperbolic functions, putting
and the formulas analogues of that of the circular functions are easily derived (see [
4]).
All the eigenfunctions
can be expressed as generalized hypergeometric functions
, [
29], namely:
,
, ...,
. In practice, starting from the Bessel function
, 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
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 and Its Iterations
It was previously noted (see e.g., [
14]) that, in the space
of analytic functions, it is possible to define an isomorphism
that preserves the differentiation properties, by means of correspondence:
where
so that
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
, since
The Hermite polynomial
becomes the Laguerre polynomial
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
.
We have, for example:
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
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 , consider the problem (see [
4] (Theorem 5.1)):
with analytic boundary condition .
The operational solution of problem (
27)
is given by: Representing
, from Equation (
28) we find, in particular:
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
of the boundary data
, 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
where
are fixed integers, have been considered.
An operational solution of the problem
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:which is an equation extending(
6).
The operator (
30) closely recalls the general case of hyper-Bessel
B operators, in [
17], since integers
could be replaced by arbitrary real numbers, as are parameters
, considered in [
17,
30]. The solutions of the general differential equation
are given by Kiryakova et al. in [
31].
4.2. L-Hyperbolic-Type Problems
Theorem 7. Let be a 2nd order differential operator with respect to the x variable, the -derivative with respect to the t variable, and denote by and two functions such that:then the abstract L-hyperbolic-type problem:with analytic initial condition , admits the operational solution(see [
4], Theorem 5.3):
where .
Please note that conditions in (
32) are satisfied, for any fixed integer
n, assuming:
4.3. L-Elliptic-Type Problems
Theorem 8. Let be a 2nd order differential operator with respect to the x variable, the -derivative with respect to the y variable, and denote by and two functions such that:then the abstract L-elliptic-type problem:with analytic boundary condition , admits the operational solution (see [
4], Theorem 5.4):
Please note that conditions in (
35) are satisfied, for any fixed integer
n, assuming:
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
, are obtained substituting the exponential with the
L-exponential
in the classic generating function, i.e., by putting
We can derive the explicit expression by applying the isomorphism
to both sides of the explicit expression of the Bessel functions, so that we find:
We proved the results:
Theorem 9. The L-Bessel functions satisfy the recurrence relation (see [
8], Theorem 2.3):
Theorem 10. The differential equation satisfied by the L-Bessel functions is(see [
8], Theorem 2.5):
where denotes the identity operator. This equation can be derived by applying the isomorphism 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:
and applying the isomorphism
, we find the equation
that is:
The solution of Equation (
38), corresponding to the Gauss’ hypergeometric equation
, is given by
where the symbol
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:
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
, by moving one of the
in the first fraction under the sum and considering
as the
-th term.
5.3. Laguerre-Type Bell Polynomials
We first note that for the Laguerre derivative, the chain rule
becomes:
and in general:
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 th Laguerre-type Bell polynomial, denoted by , represents the th Laguerre-type derivative of the composite function .
In [
12] we showed that
can be expressed as a polynomial in the independent variable
x, depending on
, in terms of the classical Bell polynomials.
According to Equation (
6), the Leibniz rule, gives:
Therefore, the following representation formula for the Laguerre-type Bell polynomials, denoted by , holds:
Theorem 11. The polynomials are expressed in terms of the ordinary Bell polynomials according to the equation (see [
12], Theorem 4.1):
The above results can be easily generalized, since
In [
12] even the general case of polynomials
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]:
where
j is a fixed integer.
The application of the isomorphism , 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
where
and
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
,
the double series
is the Laguerre-type Appell series, obtained by the classical one acting on it with the two isomorphisms
and
.
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
are convex increasing functions, with a graph lower with respect to
, 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 of population individuals at time t and let the initial number at time .
In the Malthus model, the variation is assumed to be proportional to
, i.e.,
where the
growth rate r is a suitable constant.
The solution is given by the exponential function
Using the Laguerre derivative, the Laguerre-type Malthus reads:
where
r is a positive constant. Assuming the initial conditions
we find the solution
In this case the population growth increases according to the Laguerre exponential function , 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
where
r is called the
intrinsic growth rate, and
K denotes the
environmental capacity.
The exact solution of this problem is given by
so that, if
the solution is a function monotonically increasing to
K, whereas, if
, the solution is monotonically decreasing to
K. In any case,
and the value
is a stable equilibrium point for the logistic equation.
The Laguerre-logistic (shortly
L-logistic) model is expressed by
Please note that if in the above equation N is small with respect to K, then is close to 0 and consequently .
If , then , and .
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 -derivatives, introducing suitable initial conditions which can be easily derived from the initial observations data.
Please note that as the order n increases, for , the Laguerrian exponential attenuates its growth and for it tends to assume the linear value , 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
, since the Laguerre derivative satisfies again the chain rule, according to Equation (
42).
8. Laguerre-Type Linear Dynamical Systems
Let be a matrix and denote by , the invariants of , i.e., the sum of principal minors (i.e., the elementary symmetric functions of the eigenvalues). The invariants of the matrix are given by , .
Then the solution of the linear dynamical system
writes [
33]:
where
is the characteristic polynomial of the matrix
and
denotes a simple Jordan curve encircling all the eigenvalues of
. 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
matrix, and the vectors
The following result holds (see [
15] (Theorem 10)):
Theorem 12. The Laguerre-type Cauchy problem for a homogeneous linear differential system
where
and
have been defined above.
The proof of this result is a straightforward application of the isomorphism
. 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.