Next Article in Journal
Unsteady Electro-Hydrodynamic Stagnating Point Flow of Hybridized Nanofluid via a Convectively Heated Enlarging (Dwindling) Surface with Velocity Slippage and Heat Generation
Next Article in Special Issue
Poly-Cauchy Numbers with Higher Level
Previous Article in Journal
DualAC2NN: Revisiting and Alleviating Alert Fatigue from the Detection Perspective
Previous Article in Special Issue
Implementation of Two-Mode Gaussian States Whose Covariance Matrix Has the Standard Form
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bell’s Polynomials and Laplace Transform of Higher Order Nested Functions

by
Diego Caratelli
1,2,* and
Paolo Emilio Ricci
3
1
Department of Electrical Engineering, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
2
Department of Research and Development, The Antenna Company, High Tech Campus 29, 5656 AE Eindhoven, The Netherlands
3
Department of Mathematics, International Telematic University UniNettuno, Corso Vittorio Emanuele II, 39, 00186 Rome, Italy
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(10), 2139; https://doi.org/10.3390/sym14102139
Submission received: 25 July 2022 / Revised: 3 October 2022 / Accepted: 8 October 2022 / Published: 13 October 2022
(This article belongs to the Special Issue Theory and Applications of Special Functions II)

Abstract

:
Using Bell’s polynomials it is possible to approximate the Laplace Transform of composite functions. The same methodology can be adopted for the evaluation of the Laplace Transform of higher-order nested functions. In this case, a suitable extension of Bell’s polynomials, as previously introduced in the scientific literature, is used, namely higher order Bell’s polynomials used in the representation of the derivatives of multiple nested functions. Some worked examples are shown, and some of the polynomials used are reported in the Appendices.
MSC:
44A10; 05A10; 11B65; 11B83

1. Introduction

In this study, we illustrate a procedure for the evaluation of the Laplace Transform (LT) of multi-nested analytic functions. To this end, we make use of Bell’s polynomials [1,2,3,4,5], which constitute the essential tool for computing the subsequent derivatives of composite functions.
The Bell’s polynomials appear in many different fields, ranging from number theory [6,7,8] to operator theory [9], and from differential equations [4] to integral transforms [10,11]. It is worth noting here that Bell’s polynomials are closely related to and can be written in terms of symmetric functions in combinatorial Hopf algebras [12].
The importance of the LT [13,14] is well known and it is redundant to remind it here.
We use the classic definition of the LT:
L ( f ) : = 0 exp ( s t ) f ( t ) d t = L ( s ) .
The LT converts a function of a real variable t (usually representing the time) to a function of a complex variable s (which represents the complex frequency). The LT holds for locally integrable functions on [ 0 , + ) . It is convergent in every half-plane R e ( s ) > a , where a is the so-called convergence abscissa, depending on the growth rate at infinity of f ( t ) .
Our procedure is as follows: we use Taylor’s expansion of the considered analytic function, and express the relevant coefficients in terms of Bell’s polynomials; then, we approximate the LT of the given nested function by a series expansion, which provides an asymptotic representation of the LT when that exists.
We start from the easier case of the LT of a nested exponential function, considering the first few values of the complete Bell’s polynomials. The result is a Laurent expansion which approximates the relevant LT.
Then, we consider the case of the LT of general nested functions. The main problem is to provide a table of Bell’s polynomials. These exhibit higher complexity, but their evaluation can be easily performed through a dedicated computer code.
Our results can be compared with the LT of nested functions appearing in the literature only in a few cases [15], but the results we have obtained in these cases are completely satisfactory.
All the computations reported in this study have been performed using the computer algebra program Mathematica © .
The second-order Bell’s polynomials Y n [ 2 ] , representing the derivatives of nested functions of the type f ( g ( h ( t ) ) are then introduced, and two examples of LT of these functions are given.
In Appendix A a table of the second-order Bell’s polynomials is reported.
Lastly, we give some examples to show that the same methodology can be used even for the LT of higher-order nested functions. The first few terms of the corresponding generalized Bell’s polynomials, of order 4, Y n [ 4 ] , are shown in Appendix B The polynomials Y n [ 7 ] have been computed in the same way but are not reported here owing to the lack of space.
It is worth noting that more general extensions of Bell’s polynomials have been introduced in the past, including those appearing in the two-variable case [16], as well as the multi-variable case [17]. Since all the aforementioned extensions have been proven through the classical case, more general results could be obtained by applying the methods described in this article.

2. Definition of Bell’s Polynomials

The n-th derivative of the composite (differentiable) function Φ ( t ) : = f ( g ( t ) ) , as evaluated by the chain rule, is expressed by Bell’s polynomials as follows
Φ n : = D t n Φ ( t ) = Y n ( f 1 , g 1 ; f 2 , g 2 ; ; f n , g n ) = k = 1 n B n , k ( g 1 , g 2 , , g n k + 1 ) f k ,
where
f h : = D x h f ( x ) | x = g ( t ) , g k : = D t k g ( t ) .
The coefficients B n , k , for all k = 1 , , n , are polynomials of the variables g 1 , g 2 , , g n k + 1 , that are homogeneous of degree k and isobaric of weight n (i.e., they are a linear combination of monomials g 1 k 1 g 2 k 2 g n k n whose weight is constantly given by k 1 + 2 k 2 + + n k n = n ); in the literature, they are also referred to as partial Bell’s polynomials.
Bell’s polynomials satisfy the recursion
Y 0 : = f 1 ; Y n + 1 ( f 1 , g 1 ; ; f n , g n ; f n + 1 , g n + 1 ) = = k = 0 n ( n k ) Y n k ( f 2 , g 1 ; f 3 , g 2 ; ; f n k + 1 , g n k ) g k + 1 .
An explicit representation is given by the Faà di Bruno’s formula
Y n ( f 1 , g 1 ; f 2 , g 2 ; ; f n , g n ) = π ( n ) n ! r 1 ! r 2 ! r n ! f r g 1 1 ! r 1 g 2 2 ! r 2 g n n ! r n ,
where the sum runs over all the partitions π ( n ) of the integer n, r i denotes the number of parts of size i, and r = r 1 + r 2 + + r n denotes the number of parts of the considered partition [5].
The B n , k coefficients satisfy the recursion n
B n , 1 = g n , B n , n = g 1 n , B n , k ( g 1 , g 2 , , g n k + 1 ) = h = 0 n k n 1 h B n h 1 , k 1 ( g 1 , g 2 , , g n k h + 1 ) g h + 1 .

3. LT of Composite Functions

Let f ( g ( t ) ) be a composite function that is analytic in a neighborhood of the origin, and whose Taylor’s expansion is given by
f ( g ( t ) ) = n = 0 a n t n n ! , a n = D t n [ f ( g ( t ) ) ] t = 0 .
According to the preceding equations, it results in
a 0 = f ( g 0 ) , a n = D t n [ f ( g ( t ) ) ] t = 0 = k = 1 n B n , k ( g 1 , g 2 , , g n k + 1 ) f k , ( n 1 ) ,
where
f k : = D x k f ( x ) | x = g ( 0 ) , g h : = D t h g ( t ) | t = 0 .
Then, the following result easily follows.
Theorem 1.
Consider a composite function f ( g ( t ) ) that is analytic in a neighborhood of the origin, and can be expressed by Taylor’s expansion in (6). For its LT the following asymptotic representation holds
0 + f ( g ( t ) ) e t s d t f ( g 0 ) s + n = 1 N 0 + k = 1 n B n , k ( g 1 , g 2 , , g n k + 1 ) f k t n n ! e t s d t = = f ( g 0 ) s + n = 1 N ( k = 1 n B n , k ( g 1 , g 2 , , g n k + 1 ) f k ) 0 + t n n ! e t s d t = = f ( g 0 ) s + n = 1 N k = 1 n B n , k ( g 1 , g 2 , , g n k + 1 ) f k 1 s n + 1 ,
where N denotes a finite expansion order.

3.1. The Particular Case of the Exponential Function

In the particular case when f ( x ) = e x , that is considering the function e g ( t ) , and assuming g ( 0 ) = 0 , we have the simple form
k = 1 n B n , k ( g 1 , g 2 , , g n k + 1 ) f k = k = 1 n B n , k ( g 1 , g 2 , , g n k + 1 ) = B n ( g 1 , g 2 , , g n ) ,
where the B n are the complete Bell’s polynomials. It results B 0 ( g 0 ) : = f ( g 0 ) , and the first few values of B n , for n = 1 , 2 , , 5 , are given by
B 1 = g 1 , B 2 = g 1 2 + g 2 , B 3 = g 1 3 + 3 g 1 g 2 + g 3 , B 4 = g 1 4 + 6 g 1 2 g 2 + 4 g 1 g 3 + 3 g 2 2 + g 4 , B 5 = g 1 5 + 10 g 1 3 g 2 + 15 g 1 g 2 2 + 10 g 1 2 g 3 + 10 g 2 g 3 + 5 g 1 g 4 + g 5 , B 6 = g 1 6 + 15 g 1 4 g 2 + 45 g 1 2 g 2 2 + 15 g 2 3 + 20 g 1 3 g 3 + 60 g 1 g 2 g 3 + 10 g 3 2 + 15 g 1 2 g 4 + 15 g 2 g 4 + 6 g 1 g 5 + g 6 , B 7 = g 1 7 + 21 g 1 5 g 2 + 105 g 1 3 g 2 2 + 105 g 1 g 2 3 + 35 g 1 4 g 3 + 210 g 1 2 g 2 g 3 + 105 g 2 2 g 3 + 70 g 1 g 3 2 + 35 g 1 3 g 4 + 105 g 1 g 2 g 4 + 35 g 3 g 4 + 21 g 1 2 g 5 + 21 g 2 g 5 + 7 g 1 g 6 + g 7 , B 8 = g 1 8 + 28 g 1 6 g 2 + 210 g 1 4 g 2 2 + 420 g 1 2 g 2 3 + 105 g 2 4 + 56 g 1 5 g 3 + 560 g 1 3 g 2 g 3 + 840 g 1 g 2 2 g 3 + 280 g 1 2 g 3 2 + 280 g 2 g 3 2 + 70 g 1 4 g 4 + 420 g 1 2 g 2 g 4 + 210 g 2 2 g 4 + 280 g 1 g 3 g 4 + 35 g 4 2 + 56 g 1 3 g 5 + 168 g 1 g 2 g 5 + 56 g 3 g 5 + 28 g 1 2 g 6 + 28 g 2 g 6 + 8 g 1 g 7 + g 8 , B 9 = g 1 9 + 36 g 1 7 g 2 + 378 g 1 5 g 2 2 + 1260 g 1 3 g 2 3 + 945 g 1 g 2 4 + 84 g 1 6 g 3 + 1260 g 1 4 g 2 g 3 + 3780 g 1 2 g 2 2 g 3 + 1260 g 2 3 g 3 + 840 g 1 3 g 3 2 + 2520 g 1 g 2 g 3 2 + 280 g 3 3 + 126 g 1 5 g 4 + 1260 g 1 3 g 2 g 4 + 1890 g 1 g 2 2 g 4 + 1260 g 1 2 g 3 g 4 + 1260 g 2 g 3 g 4 + 315 g 1 g 4 2 + 126 g 1 4 g 5 + 756 g 1 2 g 2 g 5 + 378 g 2 2 g 5 + 504 g 1 g 3 g 5 + 126 g 4 g 5 + 84 g 1 3 g 6 + 252 g 1 g 2 g 6 + 84 g 3 g 6 + 36 g 1 2 g 7 + 36 g 2 g 7 + 9 g 1 g 8 + g 9 , B 10 = g 1 10 + 45 g 1 8 g 2 + 630 g 1 6 g 2 2 + 3150 g 1 4 g 2 3 + 4725 g 1 2 g 2 4 + 945 g 2 5 + 120 g 1 7 g 3 + 2520 g 1 5 g 2 g 3 + 12600 g 1 3 g 2 2 g 3 + 12600 g 1 g 2 3 g 3 + 2100 g 1 4 g 3 2 + 12600 g 1 2 g 2 g 3 2 + 6300 g 2 2 g 3 2 + 2800 g 1 g 3 3 + 210 g 1 6 g 4 + 3150 g 1 4 g 2 g 4 + 9450 g 1 2 g 2 2 g 4 + 3150 g 2 3 g 4 + 4200 g 1 3 g 3 g 4 + 12600 g 1 g 2 g 3 g 4 + 2100 g 3 2 g 4 + 1575 g 1 2 g 4 2 + 1575 g 2 g 4 2 + 252 g 1 5 g 5 + 2520 g 1 3 g 2 g 5 + 3780 g 1 g 2 2 g 5 + 2520 g 1 2 g 3 g 5 + 2520 g 2 g 3 g 5 + 1260 g 1 g 4 g 5 + 126 g 5 2 + 210 g 1 4 g 6 + 1260 g 1 2 g 2 g 6 + 630 g 2 2 g 6 + 840 g 1 g 3 g 6 + 210 g 4 g 6 + 120 g 1 3 g 7 + 360 g 1 g 2 g 7 + 120 g 3 g 7 + 45 g 1 2 g 8 + 45 g 2 g 8 + 10 g 1 g 9 + g 10 .
The values of the complete Bell’s polynomials for particular choices of the relevant parameters can be found in [6].
The complete Bell’s polynomials satisfy the identity (see, e.g., [4])
B n + 1 ( g 1 , , g n + 1 ) = k = 0 n n k B n k ( g 1 , , g n k ) g k + 1 .
In this case Equation (9) reduces to
0 + exp ( g ( t ) ) e t s d t exp ( g 0 ) s + n = 1 N B n ( g 1 , g 2 , , g n ) 1 s n + 1 .
In what follows, we evaluate the approximation of the LT of nested functions. The reported results have been obtained using the computer algebra program Mathematica © .

Examples

We first recall the case of the LT of nested exponential functions, showing two particular examples.
  • Consider the Bessel function g ( t ) : = J 1 ( t ) and the LT of the corresponding exponential function. We find
0 + exp ( J 1 ( t ) ) e t s d t = 1 s + 1 2 s 2 + 1 4 s 3 3 4 s 4 11 16 s 5 19 32 s 6 + 91 64 s 7 + + 701 128 s 8 + 953 256 s 9 15245 512 s 10 + O 1 s 11 ,
for s .
  • Consider the function g ( t ) : = arctan ( t ) and the LT of the corresponding exponential function. We find
0 + exp ( arctan ( t ) ) e t s d t = 1 s 1 s 2 + 1 s 3 + 1 s 4 7 s 5 5 s 6 + 145 s 7 + + 5 s 8 6095 s 9 + 5815 s 10 + O 1 s 11 ,
for s .

4. LT in Two Known Cases

We considered two cases concerning composite functions whose transform and anti-transform are known (see [15]). By using the computer algebra program Mathematica © , we have been able to prove the correctness of the methodology used.

4.1. Case #1

Consider the function l ( t ) = log [ cosh ( t ) ] . The LT of l ( t ) is found to be [15]:
L ( s ) = 1 2 s ψ 1 2 + s 4 ψ s 4 1 s 2 ,
for s > 0 , and where ψ ( z ) is the logarithmic derivative of the gamma function, given by
ψ ( z ) d d z ln Γ ( z ) = Γ ( z ) Γ ( z ) .
Using our methodology, we find that
L ( s ) L ˜ ( s ) = 1 s 3 2 s 5 + 16 s 7 272 s 9 + 7936 s 11 ,
so that, the inverse Laplace transformation is given by
l ˜ ( t ) t 2 2 t 4 12 + t 6 45 17 t 8 2520 + 31 t 10 14175 H t ,
with H · denoting the Heaviside distribution which can be defined as follows:
H x = x δ u d u ,
in terms of the Dirac delta distribution δ · .

4.2. Case #2

Let us consider the function l ( t ) = J 0 ( t 2 ) . The LT of l ( t ) is found to be [15]:
L ( s ) = π s 16 J 1 / 4 ( s 2 / 8 ) 2 + Y 1 / 4 ( s 2 / 8 ) 2 ,
for s > 0 .
Using our methodology, we find that
L ( s ) L ˜ ( s ) = 1 s 6 s 5 + 630 s 9 207900 s 13 + 141891750 s 17 164991726900 s 21 ,
so that, the inverse Laplace transformation is given by
l ˜ ( t ) 1 t 4 4 + t 8 64 t 12 2304 + t 16 147456 t 20 14745600 H t .

5. A First Extension of Bell’s Polynomials

We consider the second-order Bell’s polynomials, Y n [ 2 ] ( f 1 , g 1 , h 1 ; f 2 , g 2 , h 2 ; ; f n , g n , h n ) , defined by the n-th derivative of the composite function Φ ( t ) : = f ( g ( h ( t ) ) ) .
Consider the functions x = h ( t ) , z = g ( x ) , and y = f ( z ) , and suppose that h ( t ) , g ( x ) , and f ( z ) are n times differentiable with respect to their variables, so that the composite function Φ ( t ) : = f ( g ( h ( t ) ) ) can be differentiated n times with respect to t, by using the chain rule.
We use, as before, the following notation:
Φ j : = D t j Φ ( t ) , f h : = D y h f ( y ) | y = g ( x ) , g k : = D x k g ( x ) | x = h ( t ) , h r : = D t r h ( t ) .
Then, the n-th derivative can be represented by the compact symbol:
Φ n = Y n [ 2 ] ( f 1 , g 1 , h 1 ; f 2 , g 2 , h 2 ; ; f n , g n , h n ) = Y n [ 2 ] ( [ f , g , h ] n ) ,
where the Y n [ 2 ] are defined as the second order Bell’s polynomials.
The first few terms are as follows.
Y 1 [ 2 ] ( [ f , g , h ] 1 ) = f 1 g 1 h 1 ; Y 2 [ 2 ] ( [ f , g , h ] 2 ) = f 1 g 1 h 2 + f 1 g 2 h 1 2 + f 2 g 1 2 h 1 2 ; Y 3 [ 2 ] ( [ f , g , h ] 3 ) = f 1 g 1 h 3 + f 1 g 3 h 1 3 + 3 f 1 g 2 h 1 h 2 + 3 f 2 g 1 g 2 h 1 3 + f 3 g 1 3 h 1 3 ; Y 4 [ 2 ] ( [ f , g , h ] 4 ) = f 4 g 1 4 h 1 4 + 6 f 3 g 1 2 g 2 h 1 4 + 3 f 2 g 2 2 h 1 4 + 4 f 2 g 1 g 3 h 1 4 + f 1 g 4 h 1 4 + 6 f 3 g 1 3 h 1 2 h 2 + + 18 f 2 g 1 g 2 h 1 2 h 2 + 6 f 1 g 3 h 1 2 h 2 + 3 f 2 g 1 2 h 2 2 + 3 f 1 g 2 h 2 2 + 4 f 2 g 1 2 h 1 h 3 + 4 f 1 g 2 h 1 h 3 + f 1 g 1 h 4 ; Y 5 [ 2 ] ( [ f , g , h ] 5 ) = f 5 g 1 5 h 1 5 + 10 f 4 g 1 3 g 2 h 1 5 + 15 f 3 g 1 g 2 2 h 1 5 + 10 f 3 g 1 2 g 3 h 1 5 + 10 f 2 g 2 g 3 h 1 5 + + 5 f 2 g 1 g 4 h 1 5 + f 1 g 5 h 1 5 + 10 f 4 g 1 4 h 1 3 h 2 + 60 f 3 g 1 2 g 2 h 1 3 h 2 + 30 f 2 g 2 2 h 1 3 h 2 + 40 f 2 g 1 g 3 h 1 3 h 2 + + 10 f 1 g 4 h 1 3 h 2 + 15 f 3 g 1 3 h 1 h 2 2 + 45 f 2 g 1 g 2 h 1 h 2 2 + 15 f 1 g 3 h 1 h 2 2 + 10 f 3 g 1 3 h 1 2 h 3 + 30 f 2 g 1 g 2 h 1 2 h 3 + + 10 f 1 g 3 h 1 2 h 3 + 10 f 2 g 1 2 h 2 h 3 + 10 f 1 g 2 h 2 h 3 + 5 f 2 g 1 2 h 1 h 4 + 5 f 1 g 2 h 1 h 4 + f 1 g 1 h 5 .
A more extended table is given in Appendix A.
The connections to the ordinary Bell’s polynomials are highlighted below.
Theorem 2.
For every integer n, the polynomials Y n [ 2 ] are represented in terms of the ordinary Bell’s polynomials by the following equation, where a compact notation similar to the one in (23) is used:
Y n [ 2 ] ( [ f , g , h ] n ) = = Y n f 1 , Y 1 ( [ g , h ] 1 ) ; f 2 , Y 2 ( [ g , h ] 2 ) ; ; f n , Y n ( [ g , h ] n )
Proof. 
Using induction, we can conclude that (24) is true for n = 1 , since
Y 1 [ 2 ] ( [ f , g , h ] 1 ) = f 1 g 1 h 1 = f 1 Y 1 ( [ g , h ] 1 ) = Y 1 f 1 , Y 1 ( [ g , h ] 1 ) .
Then, assuming that (24) is true for every n, it follows that
Y n + 1 [ 2 ] ( [ f , g , h ] n + 1 ) = D t Y n [ 2 ] ( [ f , g , h ] n ) = D t Y n f 1 , Y 1 ( [ g , h ] 1 ) ; ; f n , Y n ( [ g , h ] n ) = = Y n + 1 f 1 , Y 1 ( [ g , h ] 1 ) ; f 2 , Y 2 ( [ g , h ] 2 ) ; ; f n + 1 , Y n + 1 ( [ g , h ] n + 1 ) .
Consequently, we have the theorem:
Theorem 3.
The second-order Bell’s polynomials verify the recursion
Y 0 [ 2 ] = f 1 ; Y n + 1 [ 2 ] ( [ f , g , h ] n + 1 ) = = k = 0 n ( n k ) Y n k [ 2 ] f 2 , g 1 , h 1 ; f 3 , g 2 , h 2 ; ; f n k + 1 , g n k , h n k Y k + 1 ( [ g , h ] k + 1 ) .
Proof. 
By means of (24) we express Y n + 1 [ 2 ] ( [ f , g , h ] n + 1 ) in terms of
Y n + 1 f 1 , Y 1 ( [ g , h ] 1 ) ; ; f n + 1 , Y n + 1 ( [ g , h ] n + 1 ) .
Then, by using the recursion (9) and again Equation (24), the expression (26) follows. □

6. LT of Second-Order Nested Functions

Let be f ( g ( h ( t ) ) ) be a composite function that is analytic in a neighborhood of the origin and, therefore, can be expressed by the Taylor’s expansion
f ( g ( ( h ( t ) ) ) = n = 0 a n t n n ! , a n = D t n [ f ( g ( ( h ( t ) ) ) ] t = 0 .
According to the preceding equations, it results
a 0 = f 0 = f ( g ( h ( 0 ) ) , a n = D t n [ f ( g ( ( h ( t ) ) ) ] t = 0 = Y n [ 2 ] ( [ f , g , h ] n ) , ( n 1 ) ,
where
f h : = D x h f ( y ) | y = g ( 0 ) , g k : = D t k g ( x ) | x = h ( 0 ) , h r : = D t r h ( t ) | t = 0 .
This expansion can be used to evaluate the LT of analytic nested functions.
Theorem 4.
Consider a nested function f ( g ( ( h ( t ) ) ) that is analytic in a neighborhood of the origin, and whose Taylor’s expansion is given by (27). For its LT, the following asymptotic representation holds
0 + f ( g ( ( h ( t ) ) ) e t s d t f 0 s + n = 1 N Y n [ 2 ] ( [ f , g , h ] n ) 0 + t n n ! e t s d t = = f 0 s + n = 1 N Y n [ 2 ] ( [ f , g , h ] n ) 1 s n + 1 ,
where N denotes a finite expansion order.
Proof. 
It is a straightforward application of the definition of the second-order Bell’s polynomials. □
Example 1.
• Assuming f ( x ) = e x 1 , g ( y ) = cos ( y ) , h ( t ) = sin ( t ) , it results in (see Figure 1)
0 + exp [ cos ( sin ( t ) ) 1 ] e t s d t = 1 s 1 s 3 + 8 s 5 127 s 7 + 3523 s 9 146964 s 11 + O 1 s 13 ,
for s . The corresponding inverse LT is approximated by (see Figure 2)
l ˜ ( t ) 1 1 2 t 2 + 1 3 t 4 127 720 t 6 + 3523 40320 t 8 12247 302400 t 10 H t .
Example 2.
• Upon assuming f ( x ) = log 1 + x 2 , g ( y ) = cosh ( y ) 1 , h ( t ) = sin ( t ) , it results in (see Figure 3)
0 + log 1 + cosh ( sin ( t ) ) 1 2 e t s d t = 1 2 s 3 9 4 s 5 27 2 s 7 + 1169 8 s 9 5869 2 s 11 + O 1 s 13 ,
for s . The corresponding inverse LT can be approximated as (see Figure 4)
l ˜ ( t ) 1 4 t 2 3 32 t 4 + 3 160 t 6 167 46080 t 8 + 5869 7257600 t 10 H t .

7. Higher Order Bell’s Polynomials

Consider the nested function Φ ( t ) : = f ( 1 ) ( f ( 2 ) ( ( f ( M ) ( t ) ) ) ) , i.e., the composition of the functions x M 1 = f ( M ) ( t ) , …, x 1 = f ( 2 ) ( x 2 ) , y = f ( 1 ) ( x 1 ) , and suppose that f ( M ) , …, f ( 2 ) , f ( 1 ) are n times differentiable with respect to their independent variables. Then, Φ ( t ) can be differentiated n times with respect to t using the chain rule. By definition we put x M : = t , so that y = Φ ( t ) .
We use the following notation:
Φ h : = D t h Φ ( t ) , f ( 1 ) , h : = D x 1 h f ( 1 ) | x 1 = f ( 2 ) ( f ( 3 ) ( ( f ( M ) ( t ) ) ) ) , f ( 2 ) , k : = D x 2 k f ( 2 ) | x 2 = f ( 3 ) ( f ( 4 ) ( ( f ( M ) ( t ) ) ) ) , . . . . . . . . . . . . . . f ( M ) , j : = D x M j f ( M ) | x M = t .
Then, the n-th derivative can be represented as
Φ n = Y n [ M 1 ] ( f ( 1 ) , 1 , , f ( M ) , 1 ; f ( 1 ) , 2 , , f ( M ) , 2 ; ; f ( 1 ) , n , , f ( M ) , n ) ,
where the Y n [ M 1 ] are, by definition, Bell’s polynomials of order M 1 .
The above Theorems 2 and 3 can be generalized as follows.
Theorem 5.
For every integer n, the polynomials Y n [ M 1 ] are expressed in terms of Bell’s polynomials of a lower order, through the following equation:
Y n [ M 1 ] ( f ( 1 ) , 1 , , f ( M ) , 1 ; ; f ( 1 ) , n , , f ( M ) , n ) = = Y n ( f ( 1 ) , 1 , Y 1 [ M 2 ] ( f ( 2 ) , 1 , , f ( M ) , 1 ) ; f ( 1 ) , 2 , Y 2 [ M 2 ] ( f ( 2 ) , 1 , , f ( M ) , 1 ; f ( 2 ) , 2 , , f ( M ) , 2 ) ; ; f ( 1 ) , n , Y n [ M 2 ] ( f ( 2 ) , 1 , , f ( M ) , 1 ; ; f ( 2 ) , n , , f ( M ) , n ) .
Theorem 6.
The following recurrence relation for the Bell’s polynomials Y n [ M 1 ] of order M 1 holds true:
Y 0 [ M 1 ] = f ( 1 ) , 1 ; Y n + 1 [ M 1 ] ( f ( 1 ) , 1 , , f ( M ) , 1 ; ; f ( 1 ) , n + 1 , , f ( M ) , n + 1 ) = = k = 0 n n k Y n k [ M 1 ] f ( 1 ) , 2 , f ( 2 ) , 1 , , f ( M ) , 1 ; f ( 1 ) , 3 , f ( 2 ) , 2 , , f ( M ) , 2 ; ; f ( 1 ) , n k + 1 , f ( 2 ) , n k , , f ( M ) , n k × × Y k + 1 [ M 2 ] f ( 2 ) , 1 , , f ( M ) , 1 ; ; f ( 2 ) , k + 1 , , f ( M ) , k + 1 .
Example 3.
We apply the above results to the case of the LT of nested sine functions, assuming M = 4 and M = 7 .
• Let be M = 4 . We have (see Figure 5):
f 4 ( t ) = f 3 ( t ) = f 2 ( t ) = f 1 ( t ) = sin ( t ) , f ( t ) = sin ( sin ( sin ( sin ( t ) ) ) ) ,
0 exp ( s t ) sin ( sin ( sin ( sin ( t ) ) ) ) d t = 1 s 2 4 s 4 + 64 s 6 2160 s 8 + 121600 s 10 + O 1 s 12 ,
for s .
The corresponding inverse LT is approximated by (see Figure 6).
l ˜ ( t ) t 2 3 t 3 + 8 15 t 5 3 7 t 7 + 190 567 t 9 H t .
• Let be M = 7 . We have (see Figure 7):
f 7 ( t ) = f 6 ( t ) = = f 1 ( t ) = sin ( t ) , f ( t ) = sin ( sin ( sin ( sin ( t ) ) ) ) ,
0 exp ( s t ) f ( t ) d t = 1 s 2 7 s 4 + 217 s 6 14903 s 8 + 1776817 s 10 + O 1 s 12 ,
for s .
The corresponding inverse LT is approximated by (see Figure 8).
l ˜ ( t ) t 7 6 t 3 + 217 120 t 5 2129 720 t 7 + 253831 51840 t 9 H t .

8. Conclusions

We have presented a method for approximating the integral of analytic composite functions. We started from the Taylor expansion of the considered function in a neighborhood of the origin. Since the coefficients can be expressed in terms of Bell’s polynomials, the integral is reduced to the computation of an approximating series, which obviously converges if the integral is convergent. Then, this methodology has been applied to the case of the LT of an analytic composite function, starting from the case of analytic nested exponential functions. Furthermore, the evaluation of the LT of analytic nested functions is discussed, and the first few second-order Bell’s polynomials used in the framework of the presented methodology are reported in Appendix A, whereas those of order 4 are given in Appendix B. A graphical verification of the proposed technique, performed in the case when the analytical forms of both the transform and anti-transform are known, proved the correctness of our results. In future studies, attention will be devoted to the evaluation of more complex functions, such as the basic class of symmetric orthogonal polynomials (BCSOP) introduced in [18].

Author Contributions

The authors have contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Symmetry 14 02139 i001aSymmetry 14 02139 i001b

Appendix B

Symmetry 14 02139 i002aSymmetry 14 02139 i002b

References

  1. Bell, E.T. Exponential polynomials. Ann. Math. 1934, 35, 258–277. [Google Scholar] [CrossRef]
  2. Comtet, L. Advanced Combinatorics: The Art of Finite and Infinite Expansions; D. Reidel Publishing Co.: Dordrecht, The Netherlands, 1974. [Google Scholar] [CrossRef]
  3. Faà di Bruno, F. Théorie des Formes Binaires; Brero: Turin, Italy, 1876. [Google Scholar]
  4. Orozco López, R. Solution of the Differential Equation y(k) = eay, Special Values of Bell’s Polynomials, and (k, a)-Autonomous Coefficients. J. Integer Seq. 2021, 24, 21.8.6. [Google Scholar]
  5. Riordan, J. An Introduction to Combinatorial Analysis; J. Wiley & Sons: Chichester, UK, 1958. [Google Scholar]
  6. Qi, F.; Niu, D.-W.; Lim, D.; Yao, Y.-H. Special values of the Bell’s polynomials of the second kind for some sequences and functions. J. Math. Anal. Appl. 2020, 491, 124382. [Google Scholar] [CrossRef]
  7. Roman, S.M. The Faà di Bruno Formula. Am. Math. Mon. 1980, 87, 805–809. [Google Scholar] [CrossRef]
  8. Roman, S.M.; Rota, G.C. The umbral calculus. Adv. Math. 1978, 27, 95–188. [Google Scholar] [CrossRef] [Green Version]
  9. Robert, D. Invariants orthogonaux pour certaines classes d’operateurs. Ann. Mathém. Pures Appl. 1973, 52, 81–114. [Google Scholar]
  10. Caratelli, D.; Cesarano, C.; Ricci, P.E. Computation of the Bell-Laplace transforms. Dolomites Res. Notes Approx. 2021, 14, 74–91. [Google Scholar]
  11. Ricci, P.E. Bell’s polynomials and generalized Laplace transforms. arXiv 2021, arXiv:2103.07267. [Google Scholar]
  12. Aboud, A.; Bultel, J.-P.; Chouria, A.; Luque, J.-G.; Mallet, O. Bell Polynomials in Combinatorial Hopf Algebras. arXiv 2014, arXiv:1402.2960. [Google Scholar]
  13. Beerends, R.J.; Ter Morsche, H.G.; Van Den Berg, J.C.; Van De Vrie, E.M. Fourier and Laplace Transforms; Cambridge University Press: Cambridge, UK, 2003. [Google Scholar]
  14. Ghizzetti, A.; Ossicini, A. Trasformate di Laplace e Calcolo Simbolico; UTET: Torino, Italy, 1971. (In Italian) [Google Scholar]
  15. Oberhettinger, F.; Badii, L. Tables of Laplace Transforms; Springer: Berlin/Heidelberg, Germany; New York, NY, USA, 1973. [Google Scholar]
  16. Noschese, S.; Ricci, P.E. Differentiation of multivariable composite functions and Bell’s polynomials. J. Comput. Anal. Appl. 2003, 5, 333–340. [Google Scholar]
  17. Bernardini, A.; Natalini, P.; Ricci, P.E. Multi-dimensional Bell’s polynomials of higher order. Comput. Math. Appl. 2005, 50, 1697–1708. [Google Scholar] [CrossRef]
  18. Masjed-Jamei, M. A basic class of symmetric orthogonal polynomials using the extended Sturm–Liouville theorem for symmetric functions. J. Math. Anal. Appl. 2007, 325, 753–775. [Google Scholar] [CrossRef]
Figure 1. Magnitude (a) and argument (b) of the Laplace transform of exp [ cos ( sin ( t ) ) 1 ] as evaluated through the approximant L ˜ ( s ) and the rigorous integral expression L ( s ) for s = 5 + i ω .
Figure 1. Magnitude (a) and argument (b) of the Laplace transform of exp [ cos ( sin ( t ) ) 1 ] as evaluated through the approximant L ˜ ( s ) and the rigorous integral expression L ( s ) for s = 5 + i ω .
Symmetry 14 02139 g001
Figure 2. Distribution of l ( t ) = exp [ cos ( sin ( t ) ) 1 ] and the relevant approximant l ˜ ( t ) .
Figure 2. Distribution of l ( t ) = exp [ cos ( sin ( t ) ) 1 ] and the relevant approximant l ˜ ( t ) .
Symmetry 14 02139 g002
Figure 3. Magnitude (a) and argument (b) of the Laplace transform of log 1 + cosh ( sin ( t ) ) 1 2 as evaluated through the approximant L ˜ ( s ) and the rigorous integral expression L ( s ) for s = 5 + i ω .
Figure 3. Magnitude (a) and argument (b) of the Laplace transform of log 1 + cosh ( sin ( t ) ) 1 2 as evaluated through the approximant L ˜ ( s ) and the rigorous integral expression L ( s ) for s = 5 + i ω .
Symmetry 14 02139 g003
Figure 4. Distribution of l ( t ) = log 1 + cosh ( sin ( t ) ) 1 2 and the relevant approximant l ˜ ( t ) .
Figure 4. Distribution of l ( t ) = log 1 + cosh ( sin ( t ) ) 1 2 and the relevant approximant l ˜ ( t ) .
Symmetry 14 02139 g004
Figure 5. Magnitude (a) and argument (b) of the Laplace transform of sin ( sin ( sin ( sin ( t ) ) ) ) as evaluated through the approximant L ˜ ( s ) and the rigorous integral expression L ( s ) for s = 10 + i ω .
Figure 5. Magnitude (a) and argument (b) of the Laplace transform of sin ( sin ( sin ( sin ( t ) ) ) ) as evaluated through the approximant L ˜ ( s ) and the rigorous integral expression L ( s ) for s = 10 + i ω .
Symmetry 14 02139 g005
Figure 6. Distribution of l ( t ) = sin ( sin ( sin ( sin ( t ) ) ) ) and the relevant approximant l ˜ ( t ) .
Figure 6. Distribution of l ( t ) = sin ( sin ( sin ( sin ( t ) ) ) ) and the relevant approximant l ˜ ( t ) .
Symmetry 14 02139 g006
Figure 7. Magnitude (a) and argument (b) of the Laplace transform of sin ( sin ( ( sin ( t ) ) ) ) as evaluated through the approximant L ˜ ( s ) and the rigorous integral expression L ( s ) for s = 10 + i ω .
Figure 7. Magnitude (a) and argument (b) of the Laplace transform of sin ( sin ( ( sin ( t ) ) ) ) as evaluated through the approximant L ˜ ( s ) and the rigorous integral expression L ( s ) for s = 10 + i ω .
Symmetry 14 02139 g007
Figure 8. Distribution of l ( t ) = sin ( sin ( ( sin ( t ) ) ) ) and the relevant approximant l ˜ ( t ) .
Figure 8. Distribution of l ( t ) = sin ( sin ( ( sin ( t ) ) ) ) and the relevant approximant l ˜ ( t ) .
Symmetry 14 02139 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Caratelli, D.; Ricci, P.E. Bell’s Polynomials and Laplace Transform of Higher Order Nested Functions. Symmetry 2022, 14, 2139. https://doi.org/10.3390/sym14102139

AMA Style

Caratelli D, Ricci PE. Bell’s Polynomials and Laplace Transform of Higher Order Nested Functions. Symmetry. 2022; 14(10):2139. https://doi.org/10.3390/sym14102139

Chicago/Turabian Style

Caratelli, Diego, and Paolo Emilio Ricci. 2022. "Bell’s Polynomials and Laplace Transform of Higher Order Nested Functions" Symmetry 14, no. 10: 2139. https://doi.org/10.3390/sym14102139

APA Style

Caratelli, D., & Ricci, P. E. (2022). Bell’s Polynomials and Laplace Transform of Higher Order Nested Functions. Symmetry, 14(10), 2139. https://doi.org/10.3390/sym14102139

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop