Solution of Some Types of Differential Equations: Operational Calculus and Inverse Differential Operators

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.


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 [1] 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.

Inverse Derivative
For a common differential operator = / we can define an inverse derivative, such that upon the action on a function ( ) it gives another function ( ): whose derivative is ( ) = ( ). Evidently, the inverse derivative −1 is executed by an integral operator being the inverse of differential operator, acting on ( ), and its general form is ∫ ( ) = ( ) + , where is the constant of integration. The action of its th order 2 The Scientific World Journal can be complemented with the definition of the action of zero order derivative as follows: so that evidently In what follows we will appeal to various modifications of the following equation: Thus, it is important to construct the particular integral ( ) with the help of the following operational identity (see, e.g., [6]): For the operator̂= 2 − 2 we have We will also make explicit use of the generalized form of the Glaisher operational rule [9]; action of the operator̂= exp( 2 ) on the function ( ) = exp(− 2 ) yieldŝ Exponential operator mentioned above is closely related to Hermite orthogonal polynomials: as demonstrated in [4,5,10] by operational relations: Moreover, the following generating function for Hermite polynomials exists: Note also an easy to prove and useful relation [11]: Laguerre polynomials of two variables [4] ( , ) = exp (− ) are related to the following operator [10]: sometimes called Laguerre derivative . Note their noncommutative relation with the inverse derivative operator: They represent solutions of the following partial differential equation with proper initial conditions: where with proper initial conditions: 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.

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: with initial conditions The Scientific World Journal 3 We will employ operational approach, combined with integral transforms and exponential operator technique. Formal solution for our generic formulation reads as follows: We would like to underline that operatorŝand̂may not commute, so, dependently on the value of their commutator, we will obtain different sequences of operators in (23), disentanglinĝand̂ [13,14]. In the simplest case, when operatorŝand̂are multiplication and differentiation operators, respectively, they can be easily disentangled in the exponent with account for which yields Explicit expressions for the action of these operators on the initial condition function ( ) 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: Its formal solution reads where operatorŝ do not commute: Evidently, operators in the exponential disentangle: Thus, we have obtained the solution of the integrodifferential equation (24) as a sequence of exponential operators, transforming the initial condition ( ). Our further steps depend on the explicit form of this function. In the most general case of ( ) we may take advantage of the inverse derivative technique. First, consider the action of the operator exp(− /( −1 )) on ( ): where −1 is defined in (2). Equation (29) represents, in fact, the diffusion process and it is the solution of the following initial value problem: Relevant studies were performed in [10,12]. The initial condition function ( , 0) = ( ) can be written as follows: and the image ( ) is explicitly given by the following integral: which is supposed to converge. Then the result of the Laguerre diffusion (29) appears in the form of the translation of the image function : Consequently, we have to apply the exponential operator exp( −1 ), which can be expanded in series: The simplest example of the initial function ( ) = exp(− ) demonstrates the technique sketched above, resulting in and the Laguerre diffusion contribution (29) produces the following function: The inverse derivative action on the exponent reads − exp( ) = exp( )/ and eventually we obtain The Scientific World Journal 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 can be solved by Gauss transforms: In complete analogy with the above statements, we can solve the heat-type equation with differential operator (14) with, for example, the following initial condition: The Laguerre heat-type propagation problem (40) possesses the following solution (see the Appendix):

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 derivatives, can be further extended to other equation types. Consider the following example of a rather complicated differential equation: where , , and are some arbitrary constant coefficients and function ( ) = ( , = 0) is the initial condition. By introducing operator = and distinguishing the perfect square, we rewrite (42): where Thus, the following exponential solution for (42) appears: Now making use of the operational identity and applying exp( ) according to we obtain the following compact expression for ( , ): Consider two simple examples of initial condition functions. The first one is Then we immediately obtain the solution of (42) as follows: The second example is given by the following initial condition function: Trivial computations yield the following solution: ( , ) = (ln + ) exp (− ) .
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: with the initial condition: 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̂: The Scientific World Journal 5 wherê= exp (̂+̂) ,̂= 2 ,̂= .
The exponential of the evolution operator̂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̂and̂has the following nonzero value: Then we can apply the disentanglement operational identity: and the following chain rule: where , are constant parameters. With their help, we obtain the evolution operator̂for (54): which is, in fact, operator̂= exp( 2 ). The phase in̂is written as follows: The action of̂on the initial condition function ( ) yields the following solution for our problem: 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( 63) to ( ), apart from the factor Φ( , ; , ) . Now the explicit form of the solution (65) can be obtained by recalling thatΘ acts as a translation operator and that the action of̂on the function ( ) yields the solution of the ordinary heat equation, through Gauss-Weierstrass transform (A.1). Accordingly, we denote and we writeΘ Thus, (54) with initial condition (55) has the following explicit solution: ( , ) provided that the integral converges. Summarizing the above outlined procedure, we conclude that a solution for (54) consists in finding a Gauss transformed function with a shifted argument: ( , ) = Φ( , ; , ) ( + 2 , ) .
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 ( ) allows the expansion in the following series: In this case, we refer to the identity which arises from exp ( with account for the generating function of Hermite polynomials (11). Then the action (67) of operator̂on the initial function (72) produces ( , ) = ∑ ( , ) .
(75) 6 The Scientific World Journal In addition, the translation, operated byΘ, shifts the argument: ( , ) → ( + , ), = . Thus, we obtain the solution in a form of Hermite polynomial as follows: where Φ is defined by (64) and It is now evident that, for a short time, when ≪ / , the solution will be spanning in space off the initial function, modulated by ( , ) and exponentially depending on time: Note that for extended times ≫ / we have dominance of time variable: Φ ≫ Φ 1 , and the solution asymptotically behaves as exp( ), while plays minor role in ( +2 + 2 , ).
The same operational technique as employed for the treatise of (54) can be easily adopted for the solution of the Schrödinger equation: where is constant ( has the dimension of force). Indeed, rescaling variables in (79), we obtain the form of equation, similar to (54): where Following the operational methodology, developed for (54), we write the following solution of (80) in the form (56): where Φ is defined by (64). The integral form of the solution then is written as follows: Again, as well as in (70), without any assumption on the nature of the initial condition function ( ) of the Schrödinger equation, its solution where ( , ) is given by (85)  of translation operatorΘ and heat propagation operator̂, representing the evolution operator of the free particle. Now let us consider another Fokker-Plank type equation, that is, the following example: 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): where operatorŝand̂are defined as follows: Quantitieŝand̂evidently do not commute: It allows disentanglement of the operators in the exponential according to the following rule: The Scientific World Journal 7 Thus, the evolution operator̂action on ( ), given below: simply reduces to a Gauss transforms, where the parameter = ( ) reads as follows: Upon the trivial change of variables, we obtain where and, eventually, we end up with the following simple solution of (89): The same example of the initial Gaussian function ( ) = exp(− 2 ) as in the case of (54) yields (compare with (71)) where is defined in (97). Note that we can meet the following modified form of Fokker-Plank type equation (89): 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( ), as written below for a Gaussian ( ): However, differently from the solution of (54), where we had consequent transforms of the initial condition function ( ) by operators of translation and heat diffusion̂(63) (see also [15]), here we have just the action of̂alone with much more complicated dependence of the solution on time.

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.