TSWJ The Scientific World Journal 1537-744X 2356-6140 Hindawi Publishing Corporation 454865 10.1155/2014/454865 454865 Research Article Solution of Some Types of Differential Equations: Operational Calculus and Inverse Differential Operators Zhukovsky K. Baleanu D. Jafari H. Faculty of Physics M. V. Lomonosov Moscow State University Leninskie Gory Moscow 119899 Russia msu.ru 2014 3042014 2014 23 01 2014 23 02 2014 30 4 2014 2014 Copyright © 2014 K. Zhukovsky. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We present a general method of operational nature to analyze and obtain solutions for a variety of equations of mathematical physics and related mathematical problems. We construct inverse differential operators and produce operational identities, involving inverse derivatives and families of generalised orthogonal polynomials, such as Hermite and Laguerre polynomial families. We develop the methodology of inverse and exponential operators, employing them for the study of partial differential equations. Advantages of the operational technique, combined with the use of integral transforms, generating functions with exponentials and their integrals, for solving a wide class of partial derivative equations, related to heat, wave, and transport problems, are demonstrated.

1. Introduction

Most of physical systems can be described by appropriate sets of differential equations, which are well suited as models for systems. Hence, understanding of differential equations and finding its solutions are of primary importance for pure mathematics as for physics. With rapidly developing computer methods for the solutions of equations, the question of understanding of the obtained solutions and their application to real physical situations remains opened for analytical study. There are few types of differential equations, allowing explicit and straightforward analytical solutions. It is common knowledge that expansion into series of Hermite, Laguerre, and other relevant polynomials  is useful when solving many physical problems (see, e.g., [2, 3]). Generalised forms of these polynomials exist with many variables and indices [4, 5]. In what follows, we develop an analytical method to obtain solutions for various types of partial differential equations on the base of operational identities, employing expansions in series of Hermite, Laguerre polynomials, and their modified forms [1, 6]. The key for building these solutions will be an operational approach and development of the formalism of inverse functions and inverse differential operators, already touched in [7, 8]. We will demonstrate in what follows that when used properly and combined, in particular, with integral transforms, such an approach leads to elegant analytical solutions with transparent physical meaning without particularly cumbersome calculations.

2. Inverse Derivative

For a common differential operator D = d / d x we can define an inverse derivative, such that upon the action on a function f ( x ) it gives another function F ( x ) : (1) D - 1 f ( x ) = F ( x ) , whose derivative is F ( x ) = f ( x ) . Evidently, the inverse derivative D - 1 is executed by an integral operator being the inverse of differential operator, acting on f ( x ) , and its general form is f ( x ) = F ( x ) + C , where C is the constant of integration. The action of its n th order (2) D x - n f ( x ) = 1 ( n - 1 ) ! 0 x ( x - ξ ) n - 1 f ( ξ ) d ξ , D D D D D D D D D D D D ( n N = { 1,2 , 3 , } ) can be complemented with the definition of the action of zero order derivative as follows: (3) D x 0 f ( x ) = f ( x ) , so that evidently (4) D x - n 1 = x n n ! , ( n N 0 = N { 0 } ) . In what follows we will appeal to various modifications of the following equation: (5) ( β 2 - D 2 ) ν F ( x ) = f ( x ) . Thus, it is important to construct the particular integral F ( x ) with the help of the following operational identity (see, e.g., ): (6) q ^ - ν = 1 Γ ( ν ) 0 exp ( - q ^ t ) t ν - 1 d t , min { Re ( q ) , Re ( ν ) } > 0 . For the operator q ^ = β 2 - D 2 we have (7) ( β 2 - D 2 ) - ν f ( x ) = 1 Γ ( ν ) 0 exp ( - β 2 t ) t ν - 1 exp ( t D 2 ) f ( x ) d t . We will also make explicit use of the generalized form of the Glaisher operational rule ; action of the operator S ^ = exp ( t D x 2 ) on the function f ( x ) = exp ( - x 2 ) yields (8) S ^ f ( x ) = exp ( y 2 x 2 ) exp ( - x 2 ) = 1 1 + 4 y exp ( - x 2 1 + 4 y ) . Exponential operator mentioned above is closely related to Hermite orthogonal polynomials: (9) H n ( x , y ) = n ! r = 0 [ n / 2 ]    x n - 2 r y r ( n - 2 r ) ! r ! , as demonstrated in [4, 5, 10] by operational relations: (10) H n ( x , y ) = exp ( y 2 x 2 ) x n . Moreover, the following generating function for Hermite polynomials exists: (11) exp ( x t + y t 2 ) = n = 0 t n n ! H n ( x , y ) . Note also an easy to prove and useful relation : (12) z n H n ( x , y ) = H n ( x z , y z 2 ) .

Laguerre polynomials of two variables  (13) L n ( x , y ) = exp ( - y x x x ) ( - x ) n n ! = n ! r = 0 n ( - 1 ) r y n - r x r ( n - r ) ! ( r ! ) 2 are related to the following operator x x x : (14) D x L = x x x = D x - 1 , sometimes called Laguerre derivative D x L . Note their noncommutative relation with the inverse derivative operator: (15) [ D x L , D x - 1 ] = - 1 , ( [ A , B ] = A B - B A ) . They represent solutions of the following partial differential equation with proper initial conditions: (16) y L n ( x , y ) = - ( x x x ) L n ( x , y ) , where (17) L n ( x , y ) = n ! r = 0 n ( - 1 ) r y n - r x r ( n - r ) ! ( r ! ) 2 , with proper initial conditions: (18) L n ( x , 0 ) = ( - x ) n n ! .

In the following sections, we will investigate the possibilities to solve some partial differential equations, involving the differential operators studied above. Now we just note how the technique of inverse operator, applied for derivatives of various orders and their combinations and combined with integral transforms, allows for easy and straightforward solutions of various types of differential equations.

3. Diffusion Type, Heat Propagation Type Problems, and Inverse Derivatives

Heat propagation and diffusion type problems play a key role in the theory of partial differential equations. Combination of exponential operator technique and inverse derivative together with the operational identities of the previous section is useful for the solution of a broad spectrum of partial differential equations, related to heat and diffusion processes. Some of them have been already studied by operational method (see, e.g., [10, 12]). Below we will focus on the generalities of the solution of the following problem: (19) y F ( x , y ) = ( P ^ + M ^ ) { F ( x , y ) } , with initial conditions (20) F ( x , 0 ) = q ( x ) . We will employ operational approach, combined with integral transforms and exponential operator technique. Formal solution for our generic formulation reads as follows: (21) F ( x , y ) = exp ( y ( M ^ + P ^ ) ) q ( x ) . We would like to underline that operators P ^ and M ^ may not commute, so, dependently on the value of their commutator, we will obtain different sequences of operators in (23), disentangling P ^ and M ^ [13, 14]. In the simplest case, when operators M ^ and P ^ are multiplication and differentiation operators, respectively, they can be easily disentangled in the exponent with account for (22) [ P ^ , M ^ ] = 1 , which yields (23) F ( x , y ) = exp ( y ( M ^ + P ^ ) ) q ( x ) = exp ( - y 2 2 ) exp ( y P ^ ) exp ( y M ^ ) q ( x ) . Explicit expressions for the action of these operators on the initial condition function g ( x ) can be obtained by integral transforms, series expansions, and operational technique to be evaluated in every special case. This formulation, simple in its essence, nevertheless has wide application and allows us to frame either some of integrodifferential equations in this scheme. An elegant and interesting example is given by the following equation: (24) t F ( x , t ) = - x x x F ( x , t ) + 0 x F ( ξ , t ) d ξ , F ( x , 0 ) = g ( x ) . Its formal solution reads (25) F ( x , t ) = exp ( A ^ + B ^ ) g ( x ) , where operators (26) A ^ = - t L D x = - t D x - 1 , B ^ = t D x - 1 do not commute: (27) [ A ^ , B ^ ] = - t 2 . Evidently, operators in the exponential disentangle: (28) F ( x , t ) = exp ( A ^ + B ^ ) q ( x ) = exp ( t 2 2 ) exp ( t D x - 1 ) exp ( - t D x - 1 ) g ( x ) . Thus, we have obtained the solution of the integrodifferential equation (24) as a sequence of exponential operators, transforming the initial condition g ( x ) . Our further steps depend on the explicit form of this function. In the most general case of g ( x ) we may take advantage of the inverse derivative technique. First, consider the action of the operator exp ( - / ( D x - 1 ) ) on g ( x ) : (29) f ( x , t ) = exp ( - t D x - 1 ) g ( x ) , where D x - 1 is defined in (2). Equation (29) represents, in fact, the diffusion process and it is the solution of the following initial value problem: (30) t f ( x , t ) = - x x x f ( x , t ) , f ( x , 0 ) = g ( x ) . Relevant studies were performed in [10, 12]. The initial condition function f ( x , 0 ) = g ( x ) can be written as follows: (31) g ( x ) = φ ( D x - 1 ) 1 , and the image φ ( x ) is explicitly given by the following integral: (32) φ ( x ) = 0 exp ( - ζ ) g ( x ζ ) d ζ , which is supposed to converge. Then the result of the Laguerre diffusion (29) appears in the form of the translation of the image function φ : (33) f ( x , y ) = φ ( D x - 1 - t ) 1 . Consequently, we have to apply the exponential operator exp ( t D x - 1 ) , which can be expanded in series: (34) F ( x , t ) = exp ( t 2 2 ) n = 0 t n D x - n n ! f ( x , t ) .

The simplest example of the initial function g ( x ) = exp ( - x ) demonstrates the technique sketched above, resulting in (35) φ ( x ) = 1 ( 1 + x ) when | x | < 1 , and the Laguerre diffusion contribution (29) produces the following function: (36) f ( x , t ) = 1 1 - t exp ( - 1 1 - t x ) . The inverse derivative action on the exponent reads D x - n exp ( α x ) = exp ( α x ) / α n and eventually we obtain (37) F ( x , t ) = exp ( t 2 2 ) 1 1 - t n = 0 t n D x - n n ! exp ( - 1 1 - t x ) = 1 1 - t exp ( - x 1 - t ) exp ( - t ( 1 - t ) + t 2 2 ) .

In conclusion of the present chapter we consider the example of the solution of a heat propagation type equation by operational method, involving the inverse derivative operator and exponential operator technique. We recall that the common heat equation with initial condition problem (38) t f ( x , t ) = x 2 f ( x , t ) , f ( x , 0 ) = g ( x ) can be solved by Gauss transforms: (39) f ( x , t ) = 1 2 π t - exp ( - ( x - ξ ) 2 4 t ) g ( ξ ) d ξ . In complete analogy with the above statements, we can solve the heat-type equation with differential operator D x L (14) with, for example, the following initial condition: (40) t f ( x , t ) = D x 2 L f ( x , t ) , f ( x , 0 ) = g ( x ) . The Laguerre heat-type propagation problem (40) possesses the following solution (see the Appendix): (41) f ( x , t ) = 1 2 π t - e - ξ 2 / 4 t C 0 H ( - ξ x 2 t , - x 2 4 t ) φ ( ξ ) d ξ .

4. Operational Approach and Other Types of Differential Equations

Operational approach to solution of partial differential equations, demonstrated on the examples of diffusion-like and heat-like equations with x x x derivatives, can be further extended to other equation types. Consider the following example of a rather complicated differential equation: (42) 1 ρ t A ( x , t ) = x 2 2 x 2 A ( x , t ) + λ x x A ( x , t ) - μ A ( x , t ) , g ( x ) = A ( x , 0 ) , where ρ , λ , and μ are some arbitrary constant coefficients and function g ( x ) = A ( x , t = 0 ) is the initial condition. By introducing operator D - = x x and distinguishing the perfect square, we rewrite (42): (43) 1 ρ    t A ( x , t ) = ( ( D - + λ 2 ) 2 - ε ) A ( x , t ) , where (44) ε = μ + ( λ 2 ) 2 . Thus, the following exponential solution for (42) appears: (45) A ( x , t ) = exp { ρ t ( ( D - + λ 2 ) 2 - ε ) } g ( x ) . Now making use of the operational identity (46) exp ( p ^ 2 ) = 1 π - exp ( - ξ 2 + 2 ξ p ^ ) d ξ and applying exp ( a D - ) according to (47) exp ( a x x ) f ( x ) = f ( e a x ) , we obtain the following compact expression for A ( x , t ) : (48) A ( x , t ) = exp ( - ρ ε t ) π × - exp [ - σ 2 + σ α λ 2 ρ ] g ( x exp ( σ α ) ) d σ , where (49) α = α ( t ) = 2 ρ t . Consider two simple examples of initial condition functions. The first one is (50) g ( x ) = x n . Then we immediately obtain the solution of (42) as follows: (51) A ( x , t ) = x n exp { ρ t ( n 2 + λ n - μ ) } . The second example is given by the following initial condition function: (52) g ( x ) = ln x . Trivial computations yield the following solution: (53) A ( x , t ) = ( ln x + ρ t λ ) exp ( - ρ t μ ) .

Thus, operational technique, combined with integral transforms, operational identities, and extended forms of orthogonal polynomials, represents powerful tool for finding solutions of various classes of differential equations and initial value problems. Note that within the framework of inverse differential operators, developed and described above, the usage of the evolution operator method opens new possibilities, which we will elucidate in what follows.

Let us consider the following generalization of the heat equation: (54) t F ( x , t ) = α x 2 F ( x , t ) + β x F ( x , t ) , with the initial condition: (55) F ( x , 0 ) = f ( x ) . The evolution type equation (54) contains linear coordinate term in addition to the second order derivative. Its formal solution can be written via the evolution operator U ^ : (56) F ( x , t ) = U ^ f ( x ) , where (57) U ^ = exp ( A ^ + B ^ ) , A ^ = α t x 2 , B ^ = β t x . The exponential of the evolution operator U ^ in (56) is the sum of two noncommuting operators and it can be written as the ordered product of two exponential operators. Indeed, the commutator of A ^ and B ^ has the following nonzero value: (58) [ A ^ , B ^ ] = m A ^ 1 / 2 = 2 α β t 2 x , m = 2 β t 3 / 2 α 1 / 2 . Then we can apply the disentanglement operational identity: (59) e A ^ + B ^ = e ( m 2 / 12 ) - ( m / 2 ) A ^ 1 / 2 + A ^ e B ^ , and the following chain rule: (60) e p x 2 e q x g ( x ) = e p q 2 e q x e 2 p q x e p x 2 g ( x ) , where p , q are constant parameters. With their help, we obtain the evolution operator U ^ for (54): (61) U ^ = e Φ ( x , t ; β ) Θ ^ S ^ . Note that we have factorized two commuting operators: the operator of translation in space (62) Θ ^ = e α β t 2 x and operator (63) S ^ = e α t x 2 , which is, in fact, operator S ^ = exp ( t D x 2 ) . The phase in U ^ is written as follows: (64) Φ ( x , t ; α , β ) = 1 3 α β 2 t 3 + β t x . The action of U ^ on the initial condition function f ( x ) yields the following solution for our problem: (65) F ( x , t ) = e Φ ( x , t ; α , β ) Θ ^ S ^ f ( x ) . Thus, we conclude from the form of (65) that the problem (54) with the initial condition (55) can be solved by the consequent application of commuting operators Θ ^ (62) and S ^ (63) to f ( x ) , apart from the factor e Φ ( x , t ; α , β ) . Now the explicit form of the solution (65) can be obtained by recalling that Θ ^ acts as a translation operator (66) e s x g ( x ) = g ( x + s ) and that the action of S ^ on the function g ( x ) yields the solution of the ordinary heat equation, through Gauss-Weierstrass transform (A.1) . Accordingly, we denote (67) f ( x , t ) S ^ f ( x ) e α t x 2 f ( x ) , and we write (68) Θ ^ f ( x , t ) = f ( x + α β t 2 , t ) . Thus, (54) with initial condition (55) has the following explicit solution: (69) F ( x , t ) = e Φ ( x , t ; α , β ) 1 2 π α t -    e - ( x + α β t 2 - ξ ) 2 / 4 t α f ( ξ ) d ξ , provided that the integral converges. Summarizing the above outlined procedure, we conclude that a solution for (54) consists in finding a Gauss transformed function f with a shifted argument: (70) F ( x , t ) = e Φ ( x , t ; α , β ) f ( x + α β t 2 , t ) . Moreover, this is a general observation for this type of equation, valid for any function f ( x ) (provided the integral converges). In other words, we have obtained the solution of Fokker plank equation as a consequent action of operator S ^ of heat diffusion and operator Θ ^ of translation on the initial condition function. Note that f ( x , t ) is the solution of the heat equation, representing a natural propagation phenomenon.

The effect, produced by the translation operator Θ ^ and the operator S ^ , is best illustrated with the example of Gaussian f ( x ) = exp ( - x 2 ) evolution, when (69) becomes (71) F ( x , t ) | f ( x ) = exp ( - x 2 ) = exp ( Φ ( x , t ) ) 1 + 4 t exp ( - ( x + β t 2 ) 2 1 + 4 t ) , D D D D D D D D D D D D D D D D D D D D D D D D D D t 0 .

The above result is exactly the generalization of Gleisher rule (8), considered earlier in the context of the heat equation. Thus, (71) is the solution of the ordinary heat equation with β = 0 , when the initial function is Gaussian.

Another interesting example of solving (54) appears when the initial function f ( x ) allows the expansion in the following series: (72) f ( x ) = k c k x k e γ x . In this case, we refer to the identity (73) exp ( y D x 2 ) x k e γ    x = e ( γ    x + γ 2 y ) H k ( x + 2 γ y , y ) , which arises from (74) exp ( y m x m ) f ( x ) = f ( x + m y m - 1 x m - 1 ) 1 with account for the generating function of Hermite polynomials (11). Then the action (67) of operator S ^ on the initial function (72) produces (75) f ( x , t ) = k c k f k ( x , t ) . f k ( x , t ) = e γ ( x + γ a ) H k ( x + 2 γ a , a ) , a = α t . In addition, the translation, operated by Θ ^ , shifts the argument: f ( x , t ) f ( x + a b , t ) , b = β t . Thus, we obtain the solution in a form of Hermite polynomial as follows: (76) F ( x , t ) = k c k e Φ + Φ 1 H k ( x + α γ 2 β ( 2 t β γ + t 2 β 2 γ 2 ) , α t ) , where Φ is defined by (64) and (77) Φ 1 = γ ( x + α γ 2 β ( t β γ + ( t β γ ) 2 ) ) . It is now evident that, for a short time, when t γ / β , the solution will be spanning in space off the initial function, modulated by H k ( x , α t ) and exponentially depending on time: (78) F ( x , t ) | t γ / β k c k exp ( α γ 2 t + γ x ) H k ( x , α t ) , F ( x , t ) | t γ / β , t 1 k c k e α γ 2 t e γ x x k . Note that for extended times t γ / β we have dominance of time variable: Φ Φ 1 , and the solution asymptotically behaves as exp ( β t x ) , while x plays minor role in H k ( x + 2 γ α t + α β t 2 , α t ) .

The same operational technique as employed for the treatise of (54) can be easily adopted for the solution of the Schrödinger equation: (79) i t Ψ ( x , t ) = - 2 2 m x 2 Ψ ( x , t ) + F x Ψ ( x , t ) , Ψ ( x , 0 ) = f ( x ) , where F is constant ( F has the dimension of force). Indeed, rescaling variables in (79), we obtain the form of equation, similar to (54): (80) i τ Ψ ( x , τ ) = - x 2 Ψ ( x , τ ) + b x Ψ ( x , τ ) , where (81) τ = t 2 m , b = 2 F m 2 . Following the operational methodology, developed for (54), we write the following solution of (80) in the form (56): (82) Ψ ( x , τ ) = U - ^ f ( x ) , where (83) U - ^ = exp ( i τ H ^ ) , H ^ = x 2 - b x . Then, on account of the substitution t i τ , β - b , operators (84) Θ - ^ = exp ( b τ 2 x ) , Θ - ^ f ( x , τ ) = f ( x + b τ 2 , τ ) , (85) S - ^ = exp ( i τ x 2 ) , S - ^ f ( x ) = f ( x , i τ ) arise. Thus, the solution of the Schrödinger equation is a result of consequent action of the operator S - ^ and further action of Θ - ^ on the initial condition function: (86) Ψ ( x , t ) = exp ( - i Φ ( x , τ ; b ) ) Θ - ^ S - ^ f ( x ) , where Φ is defined by (64). The integral form of the solution then is written as follows: (87) Ψ ( x , t ) = exp ( - i Φ ( x , τ ; b ) ) 1 2 i π τ × - exp ( - ( x + b τ 2 - ξ ) 2 4 i τ ) f ( ξ ) d ξ . Again, as well as in (70), without any assumption on the nature of the initial condition function f ( x ) of the Schrödinger equation, its solution (88) Ψ ( x , τ ) = exp ( - i Φ ( x , τ ; b ) ) f ( x + b τ 2 , i τ ) , where f ( x , t ) is given by (85) and is expressed in terms of the function of two variables, obtained by the consequent application of the heat propagation and translation operators (85) and (84) to f ( x ) . So far we have demonstrated as Gauss-Weierstrass transform (A.1) describe the action of S - ^ on the initial probability amplitude and how the shift (84) finally yields the explicit form of the solution of Schrödinger equation. It means that the result of the action of evolution operator (83) on f ( x ) is the product of combined action of translation operator Θ - ^ and heat propagation operator S - ^ , representing the evolution operator of the free particle.

Now let us consider another Fokker-Plank type equation, that is, the following example: (89) t F ( x , t ) = α x 2 F ( x , t ) + β x x F ( x , t ) , with initial condition (55). Proceeding along the above outlined scheme of the solution of (54), we write the solution of (89) in general form (56): (90) F ( x , t ) = U ^ f ( x ) , U ^ = exp ( A ^ + B ^ ) , where operators A ^ and B ^ are defined as follows: (91) A ^ = t α x 2 , B ^ = t β x x . Quantities A ^ and B ^ evidently do not commute: (92) [ A ^ , B ^ ] = 2 α β t 2 x 2 = m A ^ , m = 2 t β . It allows disentanglement of the operators in the exponential according to the following rule: (93) exp ( A ^ + B ^ ) = exp ( 1 - exp ( - m ) m A ^ ) exp ( B ^ ) . Thus, the evolution operator U ^ action on f ( x ) , given below: (94) U ^ f ( x ) = exp ( σ x 2 ) exp ( β t x x ) f ( x ) = exp ( σ x 2 ) f ( e β t x ) , simply reduces to a Gauss transforms, where the parameter σ = σ ( t ) reads as follows: (95) σ ( t ) = ( 1 - exp ( - 2 β t ) ) α β . Upon the trivial change of variables, we obtain (96) U ^ f ( x ) = S ^ f ( y ) , y = x exp ( b ) ,    S ^ = exp ( ρ 2 y 2 ) , where (97) ρ ( t ) = α β ( e 2 β t - 1 ) and, eventually, we end up with the following simple solution of (89): (98) F ( x , t ) = 1 2 π ρ - e - ( exp ( β t ) x - ξ ) 2 / 2 ρ f ( ξ ) d ξ . The same example of the initial Gaussian function f ( x ) = exp ( - x 2 ) as in the case of (54) yields (compare with (71)) (99) F ( x , t ) | f ( x ) = exp ( - x 2 ) = 1 1 + 2 ρ ( t ) exp ( - e 2 β t x 2 1 + 2 ρ ( t ) ) , where ρ is defined in (97). Note that we can meet the following modified form of Fokker-Plank type equation (89): (100) t F ( x , t ) = α x 2 F ( x , t ) + β x x F ( x , t ) , in problems, related to propagation of electron beams in accelerators. Its solution arises from (99) immediately and differs from it just by a factor exp ( β t ) , as written below for a Gaussian f ( x ) : (101) F ( x , t ) | f ( x ) = exp ( - x 2 ) = e β t 1 + 2 ρ ( t ) exp ( - e 2 β t x 2 1 + 2 ρ ( t ) ) = 1 η ( t ) exp ( - x 2 η ( t ) ) , ddddddd η ( t ) = 2 α β ( 1 - e - 2 β t + β 2 α e - 2 β t ) . However, differently from the solution of (54), where we had consequent transforms of the initial condition function f ( x ) by operators of translation and heat diffusion S ^ (63) (see also ), here we have just the action of S ^ alone with much more complicated dependence of the solution on time.

5. Conclusions

Operational method is fast and universal mathematical tool for obtaining solutions of differential equations. Combination of operational method, integral transforms, and theory of special functions together with orthogonal polynomials closely related to them provides a powerful analytical instrument for solving a wide spectrum of differential equations and relevant physical problems. The technique of inverse operator, applied for derivatives of various orders and combined with integral transforms, allows for easy and straightforward solutions of various types of differential equations. With operational approach, we developed the methodology of inverse differential operators and derived a number of operational identities with them. We have demonstrated that using the technique of inverse derivatives and inverse differential operators, combined with exponential operator, integral transforms, and special functions, we can make significant progress in solution of various mathematical problems and relevant physical applications, described by differential equations.

Appendix

In complete analogy with the heat equation solution by Gauss-Weierstrass transform : (A.1) exp ( b x 2 ) g ( x ) = 1 2 π b - exp ( - ( x - ξ ) 2 4 b ) g ( ξ ) d ξ , accounting for noncommutative relation for operators of inverse derivative D x - 1 and    L D x , defined through the operational relation (14) and accounting for (31) and (32), we write the solution of (40) in the following form: (A.2) f ( x , t ) = 1 2 π t - exp ( - ( D ^ x - 1 - ξ ) 2 4 t ) 1 φ ( ξ ) d ξ , where φ is the image (32) of the initial condition function g . The kernel of the integral in the above formula can be expanded into series of two-variable Hermite polynomials H n ( x , y ) : (A.3) exp ( - ( D ^ x - 1 - ξ ) 2 4 t ) 1 = e - ξ 2 / 4 t n = 0 D ^ x - n n ! H n ( ξ 2 t , - 1 4 t ) = e - ξ 2 / 4 t n = 0 x n ( n ! ) 2 H n ( ξ 2 t , - 1 4 t ) . Taking into account formula (12) for x n H n ( ξ / 2 t , - 1 / 4 t ) in the operational identity above, which can be viewed as a generating function in terms of inverse derivative for H n ( x , y ) , we obtain the following expansion for the kernel of the integral: (A.4) exp ( - ( D ^ x - 1 - ξ ) 4 t ) 1 = e - ξ 2 / 4 t n = 0 1 ( n ! ) 2 H n ( x ξ 2 t , - x 2 4 t ) , where the series of Hermite polynomials of two arguments H n ( x , y ) can be expressed in terms of Hermite-Bessel-Tricomi functions C n H ( x , y ) —generalization of Bessel-Tricomi functions—and related to Bessel-Wright functions and to common Bessel functions : (A.5) C n H ( x , y ) = m = 0 H m ( - x , y ) r ! ( n + m ) ! , n N 0 . In particular, for n = 0 , we immediately find our series: (A.6) C 0 H ( x , y ) = m = 0 H m ( - x , y ) ( m ! ) 2 , which obviously leads to the solution of the heat-type propagation problem (40) by the following appropriate Gauss transform: (A.7) f ( x , t ) = 1 2 π t - e - ξ 2 / 4 t C 0 H ( - ξ x 2 t , - x 2 4 t ) φ ( ξ ) d ξ .

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Appèl A. Kampé de Fériet J. Fonctions Hypergéométriques et Hypersphériques: Polynômes d’Hermite 1926 Paris, France Gauthier-Villars Haimo D. T. Markett C. A representation theory for solutions of a higher order heat equation, I Journal of Mathematical Analysis and Applications 1992 168 1 89 107 2-s2.0-38249011646 Haimo D. T. Markett C. A representation theory for solutions of a higher order heat equation, II Journal of Mathematical Analysis and Applications 1992 168 2 289 305 2-s2.0-38249010207 Dattoli G. Generalized polynomials, operational identities and their applications Journal of Computational and Applied Mathematics 2000 118 1-2 111 123 2-s2.0-0040788515 Dattoli G. Srivastava H. M. Zhukovsky K. Orthogonality properties of the Hermite and related polynomials Journal of Computational and Applied Mathematics 2005 182 1 165 172 2-s2.0-19644385333 10.1016/j.cam.2004.10.021 Erdélyi A. Magnus W. Oberhettinger F. Tricomi F. G. Higher Transcendental Functions 1953 2 New York, NY, USA McGraw-Hill Book Company Dattoli G. Germano B. Martinelli M. R. Ricci P. E. Negative derivatives and special functions Applied Mathematics and Computation 2010 217 8 3924 3928 2-s2.0-78650005067 10.1016/j.amc.2010.09.056 Dattoli G. Migliorati M. Khan S. Solutions of integro-differential equations and operational methods Applied Mathematics and Computation 2007 186 1 302 308 2-s2.0-33947367482 10.1016/j.amc.2006.07.104 Srivastava H. M. Manocha H. L. A Treatise on Generating Functions 1984 Chichester, UK Halsted Press, John Wiley and Sons Mathematics and its Applications (Ellis Horwood Ltd) Dattoli G. Srivastava H. M. Zhukovsky K. A new family of integral transforms and their applications Integral Transforms and Special Functions 2006 17 1 31 37 2-s2.0-30344484568 10.1080/10652460500389081 Gould H. W. Hopper A. T. Operational formulas connected with two generalizations of Hermite polynomials Duke Mathematical Journal 1962 29 1 1 174 10.1215/S0012-7094-62-02907-1 Dattoli G. Srivastava H. M. Zhukovsky K. Operational methods and differential equations with applications to initial-value problems Applied Mathematics and Computation 2007 184 2 979 1001 2-s2.0-33846700739 10.1016/j.amc.2006.07.001 Dattoli G. Zhukovsky K. V. Appèl polynomial series expansion International Mathematical Forum 2010 5 14 649 6662 Dattoli G. Migliorati M. Zhukovsky K. Summation formulae and stirling numbers International Mathematical Forum 2009 4 41 2017 22040 Avron J. E. Herbst I. W. Spectral and scattering theory of Schrödinger operators related to the stark effect Communications in Mathematical Physics 1977 52 3 239 254 2-s2.0-0001132384 10.1007/BF01609485 Wolf K. B. Integral Transforms in Science and Engineering 1979 New York, NY, USA Plenum Press Watson G. N. A Treatise on the Theory of Bessel Functions 1944 2nd Cambridge, UK Cambridge University Press