The Laplace Transform of Composed Functions and Bivariate Bell Polynomials

: The problem of computing the Laplace transform of composed functions has not found its way into the literature because it was customarily believed that there were no suitable formula to solve it. Actually, it has been shown in previous work that by making use of Bell polynomials, efﬁcient approximations can be found. Moreover, using an extension of Bell’s polynomials to bivariate functions, it is also possible to approximate the Laplace transform of composed functions of two variables. This topic is solved in this paper and some numerical veriﬁcations, due to the ﬁrst author using the computer algebra system Mathematica © , are given proving the effectiveness of the proposed method. bivariate Bell polynomials


Introduction
In the current literature, a large number of papers falsely claimed to have generalized the classical Laplace transform and the s-multiplied (or the Laplace-Carson) transform by making some obviously trivial changes in the parameter s and the argument t. Most (if not all) of such trivialities and inconsequential variations of the classical Laplace transform and the Laplace-Carson transform, which appeared and continue to appear in the literature, were pointed out and documented by H. M. Srivastava (see [1], pp. 1508-1510, and [2], pp. 36-38).
In order not to create misunderstandings, we want to say right away that the purpose of this article is not to provide generalizations of the Laplace transform (LT), but only to extend the tables of LTs often used in applied mathematics. The well known tables of Oberhettinger and Badii cited among the references [3] does not include the LT of functions considered in this paper.
We want to underline the fact that in literature it is usually considered that there are no formulas for the computation of the LT of composed functions. This belief has been refuted in a previous article [4] and the possibility of obtaining approximations in this field is shown in an even more general situation in this paper. As far as we know, there are no other articles on this topic and our approach seems to be the first in this regard.
More precisely, we show how to approximate the Laplace transform of a 2-variable composed analytic function f [ϕ(t), ψ(t)], where ϕ(t), y = ψ(t), 0 ≤ t < +∞ are analytic functions, by using the bivariate Bell polynomials, a suitable set of Bell polynomials, introduced in a preceding paper [5].
We use the classical definition: The LT converts a function of a real variable t (representing the time) to a function of a complex variable s (representing the complex frequency).
The LT can be applied to locally integrable functions on [0, +∞). It converges in each half plane Re(s) > a, where a is a constant (called the convergence abscissa), depending on the growth behavior of f (t).
Owing to the importance of this transformation in the solution of the most diverse differential problems, a large number of LT, together with the respective anti-transforms, are reported in literature (see, e.g., [3,17]).
Given a 2-variable composed function f [ϕ(t), ψ(t)] it is natural to define the relevant LT by putting: In previous articles, we have shown how to compute the LT of higher-order nested functions (see [4] and the references therein). In this article, we apply a similar method to approximate the LT of a composed analytic function of 2 variables, taking advantage of the bivariate Bell polynomials introduced in [5]. The results obtained demonstrate the correctness of the method considered, as can be seen in the numerical verifications obtained by the first author using the Mathematica © computer algebra system.
The first bivariate Bell polynomials are given in the Appendix A at the end of this article.

Recalling the Bell's Polynomials
The Bell's polynomials [18] express the nth derivative of a composed function Φ(t) := f (g(t)) in terms of the successive derivatives of the (sufficiently smooth) component functions x = g(t) and y = f (x). More precisely, if then the nth derivative of Φ(t) is represented by where B n denotes the nth Bell polynomial. The first few Bell polynomials are given by: More general results can be found in [19], p. 49. The Bell polynomials [6] are given by the equation where the B n,k satisfy the recursion [6] B n,k (g 1 , g 2 , . . . , g n−k+1 ) = The B n,k functions for any k = 1, 2, . . . , n are polynomials in the g 1 , g 2 , . . . , g n variables homogeneous of degree k and isobaric of weight n (i.e., they are linear combinations of monomials g k 1 1 g k 2 2 · · · g k n n whose weight is constantly given by k 1 + 2k 2 + . . . + nk n = n). Therefore, we have the equations and

Bivariate Bell's Polynomials
Consider the composed function f (ϕ(t), ψ(t)), under the standard assumptions on the domains of definition and differentiability. Compute the the partial derivative of the function f , h-times with respect to x and k-times with respect to y and then put Hence, the bivariate Bell polynomials are defined as follows: Note that where B n denotes the nth ordinary Bell polynomial defined in (1). Putting we find In the above formula, for the formal multiplication symbol •, we assume the definition below . . , f x m ; ϕ 1 , . . . , ϕ m ) (and B m 1 ( f y ; ψ 1 ) = B m ( f y , . . . , f y m ; ψ 1 , . . . , ψ m )) .
Therefore, we find the result [5].

Theorem 1. Let f (x, y), ϕ(t), ψ(t) be analytical functions. The nth polynomial of the system {B
(2) n ( f ; ϕ; ψ)} n∈N , can be computed with the following rule By the Leibniz rule, we can write Then, for the bivariate Bell polynomials, the recurrence relation holds

LT of a 2-Variable Composed Function
Let f (ϕ(t), ψ(t)) be a composed function, analytic in a neighborhood of the origin, so that it is expressed by the Taylor's expansion: According to the preceding equations, it results where where f k denotes the kth derivative of the function f (x(t), y(t)) with respect to t.
Using the previous formulas we are able to approximate the calculation of the LT of the function f (ϕ(t), ψ(t))) with that of a series of elementary LT of powers. In fact, we have the following Theorem. (15), for its LT the following equation holds

Theorem 3. Considering a composed function f (ϕ(t), ψ(t)), expressed by the Taylor's expansion in Equation
where the symbolic power in the above equation must be computed using the definitions in (11).
Proof. Representing the coefficients of the Taylor expansion in (15) in terms of bivariate Bell polynomials, and using the uniform convergence of series, we find Then, the result follows from elementary calculation of the LT of powers.

The Particular Case of Exponential Functions
We start considering the case of the nested exponential function Example 1 • Let f (x, y) = e x+y and ϕ(t) = sin t; ψ(t) = −t.
Then, f (ϕ(t), ψ(t)) = exp(sin t − t), and for the relevant LT, using the above described method, we find The graphical display of our approximation is shown in Figures 1-3.

The General Case
Examples of the above method for approximating the LT of general nested functions are reported in what follows.

Conclusions
We have shown a method for approximating the LT of a 2-variable composed function f [ϕ(t), ψ(t)], where ϕ(t), ψ(t), 0 ≤ t < +∞, are analytic functions, by using the bivariate Bell polynomials. Starting from the Taylor expansion in a neighborhood of the origin of the function f [ϕ(t), ψ(t)], since the coefficients can be expressed in terms of the bivariate Bell polynomials, the integral is reduced to the computation of an approximating series, which obviously converges if the integral is convergent.
We, thus, showed that just as it was possible to approximate the LT of composed functions using Bell polynomials [4], so it can be done for the LT of composed functions in two variables, using the bivariate Bell polynomials, introduced in [5].
In the last section, the proposed technique is checked in some particular cases, when the transform and the anti-transform are known, proving the correctness of our results.
Extension could be made to higher nested functions by using the results in [4].