Operational Solution of Non-Integer Ordinary and Evolution-Type Partial Differential Equations

Konstantin V. Zhukovsky 1,* and Hari M. Srivastava 2,3 1 Faculty of Physics, Moscow State University, Leninskie Gory, Moscow 119991, Russia 2 Department of Mathematics and Statistics, University of Victoria, Victoria, BC V8W 3R4, Canada; harimsri@math.uvic.ca 3 Department of Medical Research, China Medical University Hospital, China Medical University, Taichung 40402, Taiwan * Correspondence: zhukovsk@physics.msu.ru; Tel.: +7-495-939-3177; Fax: +7-495-939-2991


Introduction
Differential equations (DE) play an important role in pure mathematics and physics.They describe a broad range of physical processes and finding their solutions is of great importance.Only a few types of DE allow explicit analytical solutions.A vast literature is dedicated to the topic, and the contribution of scientists such as A.M. Mathai can hardly be overestimated (see, for example, [1,2]).Fractional calculus has rapidly drawn increasing attention from researchers in the last decade.They study the solutions of fractional reaction-diffusion, statistical, and other equations (see, for example, [3][4][5][6].In many cases, expansion in series of orthogonal polynomials and their generalized forms with many indexes and variables as well as the usage of integral transforms are the most common tools to analytically solve DE. The method of operational solution of DE demonstrated in [7][8][9][10] is applicable to a wide spectrum of physical problems, described by linear partial differential equations (PDE), such as propagation and radiation from charged particles [11][12][13][14][15][16][17][18][19], heat diffusion [20][21][22], including processes not described by Fourier law, and others [23][24][25].In the context of the operational approach, the operational definitions for the polynomials through the operational exponent are very useful [26].The operational exponent is also applied when describing the fundamentals of structures in nature, including elementary particles and quarks [27][28][29]; such modern mathematical instruments are also used for the theoretical study of neutrino mixing [30][31][32] and for analysis of relevant experimental data [33][34][35].The obtained solutions were formulated in terms of series of generalized forms of orthogonal polynomials of Hermite, Laguerre, more general Appèl, and some other polynomials [36,37], special functions of hyperbolic, elliptic Weierstrass and Jacobi-type, cylindrical Bessel-type, and generalized Airy-type functions.
While the role of various parameters in the solutions of DE and their physical meaning is most clear in the analytical form of the solutions, this last is not always available.Modern computer methods help to solve DE.The numerical approach is widely applied nowadays due to the revolutionary breakthrough in computational technique and technical support.Advanced numerical methods for the solution of fractional differential equations, formulated, for example, in [38][39][40][41], can be effectively executed with modern computers.In this context we note also semi-analytical models and numerical simulations of relaxation of hot electrons and holes [42], the diffusion of charge carriers, and the energy relaxation and transfer with respect to the electron excited states in crystals [43,44].
Different from these numerical computations, analytical solutions, when available, give clearer insight into the underlying physical processes.In the following we will apply the operational method to obtain exact solutions for some linear ordinary DE with non-integer derivatives and for evolution-type PDE, giving examples of solutions of Schrödinger-type and Black-Scholes-type equations, and their generalized forms with the Laguerre derivative operator.
The structure of the manuscript is as follows.In the first section we will explore generalized Hermite and Laguerre polynomials, the inverse derivative operator, the Laguerre derivative, and the relations between them; we will also touch on the Appèl polynomials.In the second section we will apply the orthogonal polynomials and inverse differential operators to find the solution of some non-integer order DE.In the third section we will construct convolution forms of solutions for DE with the help of special functions and integral transforms.In the fourth section we will consider the operational solutions for some PDE; in particular, we will consider the evolution partial differential equations of Schrödinger and Black-Scholes types.In every section we will consider examples of solutions with various initial functions, such as the functions f ( . Eventually, we will provide the results and the conclusions.

Operational Approach and Orthogonal Polynomials
First of all, we note that an inverse function is one that undoes another function: For f (x) = y the inverse is g(y) = x, g( f (x)) = x.The differential operators can be treated similarly.For the DE ψ(D)F(x) = f (x), where ψ(D) is a differential operator, the inverse differential operator (ψ(D) where C is the integration constant.The inverse derivative operator of the n-th order acts according to its definition: which is complemented by the definition for its zeroth order action: and its action on the unity gives It is elementary to demonstrate that, for example, the DE ψ(D)F(x) = e αx has the following particular integral F(x) = (ψ(D)) −1 e αx = e αx (ψ(α)) −1 , and to prove the following identity: With the help of the above identity the action of the shifted inverse differential operator (ψ(D + α)) −1 on f (x) can be expressed via the inverse differential operator (ψ(D)) −1 , as follows: (5) Equation ( 5) might seem trivial, but it is particularly useful for the solution of a broad class of DE with shifted differential operators.
Traditionally, polynomial families are defined by their expansion in series.However, they can be defined operationally through the relationship with the exponential differential operators.We recall that, in general, an exponential of an operator can be viewed as the series expansion e Â = ∑ ∞ n=0 Ân /n!.
The Hermite polynomials of two variables [45], if considered in the context of the operational approach [37], can be explicitly defined by the following operational rule [36] in addition to their series expansion [46]: For the first-order polynomial we obtain simply and for the second-order polynomial we have the two-variable Hermite polynomials H n (x, y): Thus, the Hermite polynomials of two variables are defined through the action of the heat operator Ŝ: on the monomial x n .The heat operator (9) was thoroughly studied, for example, in [47].The Hermite polynomials of two variables have the following generating function: and they actually represent another form of the common Hermite polynomials of one variable: Direct application of the operational definition (8) to the Hermite polynomials yields the following identity: which consists of a shift in the y variable.
With the help of the following relation for Hermite polynomials: and with the operational identity: applied together with the operational rule (5), we obtain for the action of the heat diffusion operator Ŝ on the polynomial-exponential function the following result: x k e αx = e (αx+α 2 y) H k (x + 2αy, y).
The Hermite, Laguerre, and some other polynomials belong to a more general family of Appèl polynomials [48], if viewed in the framework of the operational approach.For example, the twovariable Hermite polynomials belong to the family of Appèl polynomials a n (x), which can be defined through the following generating function [47]: where it is assumed that a finite region of t exists, in which A(t) is expandable in Taylor series and this expansion converges.Then, with the help of the obvious identity: te xt = Dx e xt , Dx = d/dx, we can rewrite Equation ( 16) in the following operational form: Now, expanding the exponential in Equation ( 17) in series and equating the terms on the rightand left-hand sides of (17), we obtain the following definition for a n (x): where A Dx is the Appèl operator.In the case of the two-variable Hermite polynomials, the identity (18) becomes the operational definition (8) and Appèl operator for Hermite polynomials is realized by the exponential Let us assume that the inverse of the Appèl operator A Dx −1 can be defined as The main properties of the Appèl polynomials arise from the operational definition (18).For example, if the operators A Dx and Dx commute: A( Dx ), Dx = 0, then, by acting with Dx on both sides of (18), we obtain the following relation for a n (x) and a n−1 (x): Moreover, it follows from (18) that a n+1 (x) and a n (x) are related to each other as follows: This allows us to introduce the multiplicative operator M for Appèl polynomials: where M is given by the Appèl operator as follows: and on account of f ( Dx ), x = f ( Dx ), where f is the derivative of f, we write: For the Appèl polynomials a n (x) the operators M and Dx stand for the multiplicative and derivative operators and this set of operators: M, D, 1, realizes the Weyl-Heisenberg algebra.From the following relation for Appèl polynomials: it is easy to derive the following differential equation for Appèl polynomials: where A is the derivative of A. This equation is valid for all of the polynomials belonging to the Appèl family.Moreover, it is easy to recognize that Appèl polynomials satisfy the following recurrence: In the context of the Appèl polynomial family we obtain for the Hermite polynomials the following multiplicative operator M: The differential equation for the Hermite polynomials as part of the Appèl polynomial family reads as follows: 2y D2 x H n (x, y) The Laguerre polynomials of two variables can be defined in the operational way or as a finite sum: The Laguerre polynomials of two variables, as well as the Hermite polynomials of two variables, are just another way for writing proper polynomials of one variable [47]: However, there is more than just another notation behind the introduction of this form with two variables in Hermite and Laguerre polynomials.It allows us to consider proper polynomials as solutions of partial differential equations (PDE) with proper initial conditions: for Laguerre polynomials L n (x, y) and for Hermite polynomials H n (x, y).We introduce the Laguerre derivative L D x and then the two variable Laguerre polynomials can be operationally defined as follows: This operational definition is equivalent to the summation definition (30), which can be easily proved by direct execution of the action of L Dx on (−x) n n! : The differential and multiplicative operators are formed by the operators which do not commute: Moreover, in the framework of the inverse derivative (see ( 1)) the following operational relationship exists between them 10: which immediately raises associations with the relationship between the momentum and the coordinate in quantum mechanics.This relationship allows us to solve operationally the differential equations with the Laguerre derivative operator ∂ x x∂ x , as we will demonstrate in what follows.Directly from (36) and (38) we conclude that the Laguerre polynomials L n (x, y), (30) and (34), can be expressed in terms of the inverse derivative operator (1) as follows: This relation is also particularly useful for solution of some types of DE, involving the Laguerre derivative.Moreover, the operational definition (34) and the relations (38) and (39) yield the following operational rule for the Laguerre polynomials: Framing classical polynomials in the Appèl family should be done with some caution.Strictly speaking, the two-variable Laguerre polynomials can be considered members of the Appèl family with respect to the y variable only.Indeed, they are not Appèl polynomials with respect to x.However, the Laguerre polynomial family can be introduced in the context of the Appèl family in a way similar to ( 16) and ( 18) by the following substitution: In this way we obtain the following formula: With respect to the y variable, the Laguerre polynomials, as defined in (30), certainly belong the Appèl family as they can be given by the following operational rule: Moreover, a hybrid family of polynomials exists, defined by the Appèl operator C 0 (x D2 y ) or, alternatively, defined by the following sum: Further study of their properties is beyond the scope of the present paper, but they are quite interesting, being in between those of Laguerre and Hermite polynomials.Moreover, for x = 1 − y 2 /4 these polynomials reduce to the Legendre family.Studies of these and relevant polynomials were recently performed in [49][50][51].
Eventually, let us note that umbral calculus can provide a common framework for known and new identities for orthogonal polynomials.Let us recall the identity [52] which, after defining the umbral variable â 53, reads in its terms as follows: Therefore, as a consequence of the binomial theorem and of definition (46), we can write Equation (45) in the following useful form: Now with the help of identity (47), it is easy to generate new identities, such as the obvious consequence of Equation ( 47): which, together with yields the relation Moreover, from ( 50) and (47) more identities follow: We can define the operator of umbral derivative ∆a [53] by the following rule: which, together with the multiplication condition: yields the following result for the commutator of the two operators [53]: Equation ( 55) allows us to use Weyl-Heisenberg algebra when needed.We can define the associated Hermite polynomials of two variables, with the operator â as one of the variables, and thus we come to the following sum: The multiplication condition (54) does not define any new polynomial family and such Hermite polynomials H n (â, y) satisfy the following relation and the following recurrences: which are direct generalizations of the relevant terms for the two-variable Hermite polynomials H n (x, y) equations.Indeed, Equation (59) is the umbral heat equation-the direct generalization of the heat Equation (33).It can be used to define the associated polynomials (56) in terms of the following operational equation: which is the generalization of definition (8).Further study of this topic represents stand-alone research; it will be addressed elsewhere.

Operational Solution of Some Non-Integer Ordinary DE
Let us consider a differential equation where ν is not necessarily an integer, shifted by the constant Its particular integral formally reads and it can be found in the form of the integral if the well-known operational identity [26,47] is applied: For q = β 2 − D 2 we obtain the following particular solution, involving the integrated weighted action of the operator e t D 2 on the initial function f (x): Example 1.It is inviting to choose the initial function for Equation (61) in the form of the monomial f (x) = x n .The action of the heat diffusion operator Ŝ on the monomial gives the Hermite polynomials according to their operational definition (30); the action of Ŝ on the polynomial-exponential function is given in (15).With account for the generating function (10), we directly write the particular integral (62) for f (x) = x n as follows: The resulting function (65) with the Hermite polynomial of two variables is characterized by the shift of the argument x → x + 2αt .Evidently, for α = 0 we have Equation (66): whose solution in the integral form is nothing but a particular case of (64) with D → D : and involves the action of the heat diffusion operator Ŝ = e t∂ 2 x (9).

Example 2.
Let us choose the Gaussian initial function f (x) = e −x 2 .Then, by means of the operational rule 47, we immediately get the desired solution: Now, let us consider the following equation with the Laguerre derivative L D x , where ν is not necessarily an integer: Its operational solution is: The common change of variables t → e t in such cases transforms the solution of the fractional Laguerre Equation (70) into e tν e −βe t e e t L D x f (x)dt, and in the particular case of β = 1, ν = 1 we obtain the Laplace transforms 1/(1 − â) = ∞ 0 e −s(1− â) ds for the Laguerre derivative operator L D x , where the substitution a → L D x has been performed.Example 3.For Equation (71) with the initial monomial f (x) = x n the following particular integral arises: Suppose the initial function f (x) is expandable in series of the Laguerre polynomials L n (x): Then with the help of ( 40), we readily write the solution (71) of Equation (70) in the integral form: gives the solution: Now, let us consider an ordinary DE like (61), with shifted Laguerre derivative L D x instead of the common derivative d x : Let us choose the initial function for (78) in the form of the particular case of the Bessel-Wright function: where is the particular case of the Bessel-Wright function [47].In complete analogy with (63) we readily write the operational integral solution: Now we should compute the action of the heat operator with Laguerre derivative L D x on the initial function f (x) = W 0 (−x 2 , 2).With the help of the operational definition of Laguerre polynomials (30) and of the Gleisher operational rule [10], we obtain the particular integral as follows: The operational definitions of the polynomials and relevant operational rules allow writing solutions with ease for other types of equations too.For example, consider the following fractional order DE: Usage of the operational rule (63) immediately yields the integral solution for (84): which involves the operational exponent action: e tD x f (x), where we denoted the differential operator Consider the initial condition f (x) = x n .Direct application of the operational definition of the generalized Laguerre polynomials immediately gives results in the particular integral of the generalized Laguerre polynomials with the exponential power weight: Now consider the other initial condition function: f (x) = e −γx .To obtain the solution we exploit the generalized form of the Gleisher operational rule 47, which yields the solution We have demonstrated that the usage of the inverse derivative, combined with the operational formalism, provides a straightforward and easy way of solving some classes of linear DE.In what follows we will demonstrate how this technique allows solutions of partial differential equations (PDE).

Convolution Forms for Solution of DE
In what follows, we will apply the inverse differential operators in order to obtain the convolution forms of solution for Equation (61).The operational approach to the solution of Equation (61) involves the exponential operator technique, the inverse derivative formalism, and integral transforms.In general, for solution of equations with D + α operational rule (5) can be applied, where We continue with account for (66), ( 67) and ( 9), and make use of the action of the heat diffusion operator Ŝ ( 9) on e αx g(x) with the help of the following chain rule: x e αx g(x) = e α x e α 2 y e 2αy∂ x e y∂ 2 x g(x), where y and α are the parameters.This results in the following particular solution for Equation (61), expressed as the integral: where Ŝ is the heat operator ( 9) and Θ is the well-known operator of translation: The action of the operator Ŝ = e t∂ 2 x can be written in the form of the Gaussian integral transform: Therefore, apart from the phase factor, the solution (92) of Equation ( 61) consists of the integrated action of the heat operator Ŝ and in the consequent translation by Θ of the initial function f (x): where the integrand function U(x, t) is (94), shifted by Θ: Example 5.The example of the Gaussian initial condition f (x) = exp(−x 2 ) can be the illustration of the operational solution, described above.Accounting for (68), we directly write the solution of Equation ( 61), which is in turn a Gaussian: Note that from the general form of the solution (92), using the operational definition of the Hermite polynomials (8), we can directly obtain the solution (65) of the DE (61).
For the solution (64) of DE (61) with given initial function f (x) we have to calculate the action of the exponential differential operator e t D 2 in the exponential.This can be performed in a number of different ways.One of them consists in direct application of operational definitions, as we did in the case of the initial monomial x n .However, this is a rare case.The exponential operator of the second-order derivative can be reduced to the exponential of the first-order derivative if we apply the integral presentation for the exponential of a square of an operator q [54]: in our case q = √ t D. The above formula then reads as follows: Accounting for the action of the translation operator e η(∂ x +α) f (x) = e ηα f (x + η), we obtain the following particular integral (62) for the DE on non-integer order (61): Now, upon subject to the change of variables we end up with the following form of the particular solution for Equation (61): Several convolution forms are possible for the solution of (61).Indeed, for an arbitrary function f (x) in the r.h.s. of (61) and the real values of α and ν > 0 we can involve the generating function for Hermite polynomials (10) to disentangle two integrals in (102): It follows from Equation (103) that the solution of DE (61) can be written in the form of series with the power kernel: The respective coefficients in the series depend on the order of the equation, which can be a non-integer, and on the constants α, β as follows: Other convolution forms with different kernels are possible.Among them the Gaussian frequency kernel form is, perhaps, the most compact.Indeed, the integral form (102) of the solution of DE (61) can be viewed as the integral of the following convolution: where The above expression involves the convolution with the Gauss frequency function kernel.Furthermore, the remaining integral can be taken, and it gives the Bessel function of the second kind K κ (x): Note that for the integer order of the equation, ν ∈ Z, we have semi-integer index of the Bessel function of the second kind, K n−1/2 (x), the latter easily expressed in elementary functions, for example: Thus, we have obtained the particular solution for DE (61) in the form of the integral, which appears in the form of the convolution with the initial function f (x): with the kernel, containing the Bessel function of the second kind K ν−1/2 , the exponential, and the n power of x: Finally can we write the compact convolution form of the solution of DE (61) as follows: So far we have demonstrated that the usage of the inverse derivative and of the inverse differential operators constitutes a straightforward and easy way to solve some classes of linear DE.In what follows we will apply this concept to solve more complicated problems, formulated in terms of PDE.

Operational Solution for Evolution-Type Partial Differential Equations
The technique of the inverse differential and exponential operators is useful for finding solutions to a broad range of mathematical and physical problems.In what follows we shall demonstrate the solution of the evolution-type DE by the operational approach.Let us consider the Schrödinger equation for an electric charge in a constant electric field in imaginary time.It effectively corresponds to the case when the charge diffuses under a potential barrier in the electric field, so that the charge energy is lower than the height of the barrier.This process is governed by the Schrödinger equation upon the t → iτ, β → −β change: which is the common heat equation ∂ t F(x, t) = ∂ 2 x F(x, t) with the linear term βx in the r.h.s.The solution of DE (113) can be obtained operationally: and consists in the transform of the initial function F(x, 0) = f (x) by the operators Ŝ = e αt∂ 2 x and Θ = e αβt 2 ∂ x .Note that, although the solution for the Schrödinger equation in the electric field in real and in imaginary time, i.e., over and under the barrier, has the same structure (114), there is a fundamental difference between them.Indeed, the F(x, t) function for a particle in quantum mechanics is the amplitude of the probability of finding it at point x at moment t: F(x, t) = Ψ(x, t).For the charge over the barrier, the solution F(x, t) → Ψ(x, τ) of the Schrödinger equation is complex due to the complex phase Φ(x, τ; β); this does not trouble the probability |Ψ(x, τ)| 2 over the barrier for t → ∞ , which regularly converges.Example 6.Let us consider the initial polynomial f (x) = ∑ n c n x n in the context of the Fourier heat conduction of DE (113).The operational definition of the Hermite polynomials ( 8) gives e a∂ x x n = H n (x, a), and the operator Θ = e b∂ x gives the shift: F(x) ∝ H n (x + b, a).The solution immediately appears in terms of the sum of the Hermite polynomials: Example 7. Now let us choose the initial condition f (x) = ∑ k x k e γx .This function for γ < 0 represents a pulse, the shape of which depends on the values of k and γ, and varies from a sudden surge to a flat, smooth spatial wave.This choice of the initial function allows for modeling heat pulses for experimental tests (see, for example, [55]).Now applying the operational rule (15), where, in our case y → αt , and the shift by the translation operator Θ = e αβt 2 ∂ x , we obtain the solution in the form of the Hermite polynomials with the common phase Φ written in the solution (114) and ∆ 1 = γ (x + γα t +αβt 2 ) .For γ = 0 it immediately returns the result (115).For pure Fourier heat conduction β = 0 and the solution further simplifies: It is easy to follow its evolution in time: for t >> x/αγ the coordinate dependence fades out: F(x, t)| t>>x/αγ ∼ = H k (2tαγ, αt)exp γ 2 αt , and the time dependence prevails.For relatively short times of the evolution of the initial heat pulse f (x) = ∑ k x k e γx , such that t << x/αγ, the solution is approximated by F(x, t)| t<<x/αγ ∼ = e xγ 1 + tαγ 2 H k (x, αt) and for very short times αt → 0 the Hermite polynomials tend to H k (x, 0) = x k , which is in perfect agreement with our initial condition f (x) = x k e γx .Deeper consideration of the above topic is beyond the scope of the present paper.In forthcoming publications we will apply the operational method to explore and solve relativistic heat equations and other non-local extensions of the heat conduction.
Let us consider the following modification of the common Black-Scholes differential equation with the Laguerre derivative (see (30) and (36)) and the initial function g(x) = A(x, 0): where ρ, λ, and µ are constants.Equation ( 118) is in fact the general form of the equation, which unifies the Laguerre heat equation and the matter diffusion equation with Laguerre derivative, as previously explored in [10,37].In order to solve Equation (118), we employ the operational method.As usually in Black-Scholes DE, we distinguish the perfect square of the derivative, in this case of L D x = ∂ x x∂ x .
Then the solution takes the form of the exponential With the help of the operational identity (98), we reduce e (a L D x ) 2 to the first-order Laguerre derivative in the exponential and thus the following solution for A(x, t) arises: The above integral form of the solution, provided the integral converges, contains the exponential Laguerre derivative, which acts on the initial function: e −a L D x g(x).
Example 8. Let us consider the example of the polynomial initial function A(x, 0) = g(x) = ∑ n c n x n .Following the operational definition of the Laguerre polynomials (30), we directly write the solution for DE (118): Consequent integration results in the finite sum, involving gamma function Γ and hypergeometric function 1 F 1 : Now suppose the initial function A(x, 0) can be expanded in series of the Laguerre polynomials: g(x) = ∑ n a n L n (x).The operational relationships (40) and (31) in this case immediately propose the solution of DE (118) in the following form: Now let us consider the general case of the initial function A(x, 0) = g(x).Then the solution of DE (118), A(x, t) can be obtained in the following steps.With the help of the operational definitions (36) and with the inverse derivative Formula (1), we write the solution in terms of the operator of the inverse derivative D −1 x and of the function ϕ: where ϕ is the image function, determined by the integral: ϕ(x) = ∞ 0 exp(−κ)g(xκ)dκ.The exponential of the Laguerre derivative, acting on the initial function, yields the solution of the Laguerre diffusion equation 10: Hence, by applying the exponential differential operator to the function g x 1, we obtain the solution of the Laguerre diffusion equation: With account for the above relation (126) the solution of DE (118) becomes: where C 0 is the particular case of the Bessel-Tricomi function [56]: . Finally, we have obtained the solution of DE (118) with the initial condition g(x) = W 0 (−x 2 , 2) in the form of the integral (127) of the exponentially weighted function (129).

Results
We have obtained solutions for some ordinary DE of non-integer order with shifted derivatives.
In particular, we derived the integral form of the particular solution F The linear evolution-type PDEs were solved by the operational technique.In particular, the Black-Scholes equation with the Laguerre derivative L D x = ∂ x x∂ x was solved operationally.The example of the initial polynomial was considered.By using the operational definitions for Hermite polynomials we obtained explicit solutions in the form of the polynomials of x with the coefficients, given by Γ and 1 F 1 functions.The solution of the Black-Scholes type equation with Laguerre derivative L D x for the Bessel-Wright function f (x) = W 0 −x 2 , 2 is obtained in the integral form, involving Bessel-Tricomi function C n (x).We have obtained the operational solution of a Fourier-type heat equation with an additional term, describing the heat exchange with the environment, for the initial distribution f (x) = ∑ k x k e γx , which describes a heat pulse for γ < 0. We also obtained the solution of the Schrödinger equation for a charge in electric field in real and in imaginary time, i.e., over and under the potential barrier, and demonstrated that in real time, i.e. under the barrier, the solution is purely real, contrary to that over the barrier.Thus |F| 2 diverges for t → ∞ in the case of the real solution and β = 0, but converges otherwise, as a square of the amplitude |Ψ(x, t)| 2 of the probability function should behave (see also [57]).

Conclusions
In the present work we advocate the operational approach for solution of linear DE and the use of inverse differential operators, which allow direct and straightforward finding of solutions.The latter include the action of the operator of heat conduction and the operator of shift and dilatation.The operational approach involves operational definitions for Hermite and Laguerre orthogonal polynomials.In this way, we avoid cumbersome calculations and directly obtain the results of the action of proper exponential differential operators on the initial functions.If the DE contains the Laguerre derivatives, the commutation relationship between the inverse derivative operator and the Laguerre derivative operator helps in solving DE.Complemented by the usage of the integral transforms where needed, the operational technique yielded solutions of relatively complicated DE, such as Black-Scholes-type DE with Laguerre derivative, etc.Thus, our research demonstrates that the operational approach for solution of linear DEs is advantageous for its ease.The solutions are derived directly based on the operational definitions and on commutation relationships.Operational study of more complicated equations, describing heat propagation accounting for wave and ballistic heat transfer and, for equations, modeling other physical processes, is possible.It will be performed in forthcoming publications.

Example 4 .
Let us consider the initial exponential function f (x) = e −γx .Then, the usage of the generalized Gleisher operational rule 10, e −t L D x e −γx = e − γx 1−γt

Example 9 .r=0(− 1 )
Let us consider the particular case of the Bessel-Wright function[47] W n (x, m) = n ∈ 0 , for m = 2, n = 0 as the initial function: g(x) = W 0 (−x 2 , 2) = ∞ ∑ r x 2rr!(2r)! .Its image is ϕ(x) = e −x 2 .The operational identity (98) and the function (128) together yield, in accordance with the previously computed in[10], the result: g x) for real values of ν.The integrand involves the operators of heat propagation Ŝ and translation Θ, which act on the function f (x).Moreover, the convolution form of these solutions φ(x, τ) = G(x, τ) * f (η) and the integrals of other convolutions with several kernels different from each other are obtained.The comprehensive solution with the kernel, involving the Bessel function of the second kind with power-exponential weight, is obtained.Other integral forms of the solution with the convolutions with the Gaussian frequency kernel and with the monomial kernel are also obtained.We considered the examples of the Gaussian distribution f (x) = e −x 2 and of the monomial f (x) = x k and found explicit solutions for them in terms of integrals and series of Hermite polynomials.We operationally solved the DE with Laguerre derivatives:(x∂ 2 x + (α + 1)∂ x ) ν F(x) = f (x)and demonstrated the examples of solutions for the functions f (x) = exp(−γx), f (x) = x k and for the Bessel-Wright function f (x) = W 0 −x 2 , 2 .The obtained operational solutions are expressed in terms of the integrals of generalized Laguerre polynomials and Bessel functions.