Arxiv:2103.13478V2 [Math.NT] 11 May 2021 Solution of the Differential Equation Y = E , Special Values of Bell Polynomials

Home , Tangent, Touchard polynomials

arXiv:2103.13478v2 [math.NT] 11 May 2021 ino ucin.Freape a iBuo[ derivativ Bruno to di Faa related example, closely For are polynomials functions. Bell of that tion fact known a is It Introduction 1 eiaieo h opsto ffntos lo enriiadRic and Bernardini Also, functions. of composition the of derivative fteaayi ouino h uooosdffrnilequation differential [ autonomous Orozco the of hand, solution other analytic the the of On equations. repre differential to and mathematics in tool useful very a are polynomials Bell that l [ al. h etrfield vector the ensuid r osrce sn elplnmas nteother the On h which polynomials. numbers, Bell of using class constructed special are a studied, general, In been c numb polynomials. others, tangent Bell among the reduced numbers, using as Blasius of such numbers, convergence numbers zigzag Euler known the numbers, well consider that not show will will we We function. exponential a iBuosfrua ecnte oka ieeta qain as equations differential polynomials. at Bell look of then values can special We investigating formula. Bruno’s di Faa nti ae ewl ou nfidn pca auso elplnmasw polynomials Bell of values special finding on focus will we paper this In 12 ouino h ieeta Equation Differential the of Solution il,atnmu oyoil n ( and polynomials au autonomous exponential mials, functions, autonomous exponential partial obntra bet n ioilcecet,Siln n coefficients. Stirling num autonomous various coefficients, and between binomial kind relationship and the objects show combinatorial we Finally, fined. oe eisslto fteequation the of solution series power ,Cly[ Caley ], oyoil and Polynomials nti ae pca auso elplnmasaegvnby given are polynomials Bell of values special paper this In y ( k ) 3 ,adWn [ Wang and ], f = ( x [email protected] fteatnmu ieeta equation differential autonomous the of ) eatmnod Matemáticas de Departamento e ay nvria elsAndes los de Universidad pca auso Bell of Values Special , oadOoc López Orozco Ronald 13 Coefficients Bogotá, 111711 hwdterltosi ewe elpolynomials Bell between relationship the showed ] Colombia Abstract ( 1 y ,a k, ,a k, ( k ) 5 -uoooscecet r de- are coefficients )-autonomous = ,Fis [ Foissy ], e ) ay -Autonomous nadto,cmlt and complete addition, In . 6 9 ,adRodn[ Riordan and ], tde h convergence the studied ] ooospolyno- tonomous meso second of umbers y y ( escounting bers k ( nb constructed be an ) k ) = i[ ci sn the using = ouin,but solutions, so composi- of es r,Bernoulli ers, etthe sent f ad new a hand, ( orefor source a 2 v o yet not ave f y 10 ,Yli et Yildiz ], ( yusing by ) y showed ] sthe is ) n hen -th family of numbers, called (k, a)-autonomous coefficients, is obtained for each value of k. Four conjectures about these numbers are established. This paper is divided as follows. We begin with a summary of results on complete and partial Bell polynomials, which will be used to demonstrate the main results presented here. Next, we introduce the complete and partial exponential autonomous functions, the recurrence relations of these are constructed using Bell polynomials, and some recurrence relations of solutions of various initial value prob- lems are given. In the fourth section the (k, a)-autonomous coefficients are intro- duced. From these numbers we can obtain the triangular numbers, the 8-sequence numbers of [1, n] with 2 contiguous pairs, among others. We finish this work by studying the cases k =2, 3, 4 for the autonomous differential equation y(k) = eay.

2 Bell exponential polynomials

The following basic results can be found at Comtent [4], and Riordan [11]. Expo- nential Bell polynomials are used to encode information on the ways in which a set can be partitioned, hence they are a very useful tool in combinatorial analysis. Bell polynomials are obtained from the derivatives of composite functions and are given by Faa Di Bruno’s formula [5]. Bell [1], Gould [7] and Mihoubi [8] provided important results on these polynomials. We start with the deﬁnition of partial Bell polynomials

Deﬁnition 1. The exponential partial Bell polynomials are the polynomials

Bn,k(x1, xn, ..., xn k+1) − in the variables x1, x2, ... deﬁned by the series expansion

n ∞ j ∞ n t t k exp u xj =1+ u Bn,k(x1, x2, ..., xn k+1). (1) j! n! − j=1 ! n=1 X X Xk=1 The following result gives the explicit way to calculate the partial Bell polynomials

Theorem 2. The partial or incomplete exponential Bell polynomials are given by

c − +1 c1 n k n! x1 xn k+1 Bn,k(x1,...,xn k+1)= − − c1!c2! cn k+1! 1! ··· (n k + 1)! X ··· − − where the summation takes place over all integers c1,c2, ..., cn k+1 0, such that − ≥

c1 +2c2 + +(n k + 1)cn k+1 = n ··· − − c1 + c2 + + cn k+1 = k ··· − The following are special cases of partial Bell polynomials and will be very useful

2 for proving results in this paper

Bn,1(x1, ..., xn)= xn, (2) n 1 1 − n Bn,2(x1, ..., xn 1)= xkxn k, (3) − 2 k − k=1 2Xa j! n n j x2 x2(a+1) j Bn,n a(x1, ..., xa+1)= x1− Ba,j a , , − , − a! j − 2 ··· 2(a + 1) j j=a+1 X − 1 a < n, (4) ≤ n Bn,n(x1)= x1 , (5)

n n 2 Bn,n 1(x1, x2)= x1− x2, (6) − 2 n n 3 n n 4 2 Bn,n 2(x1, x2, x3)= x1− x3 +3 x1− x2, (7) − 3 4 n n 4 n n 5 n n 6 3 Bn,n 3(x1, x2, x3, x4)= x1− x4 + 10 x1− x2x3 + 15 x1− x2, (8) − 4 5 6 n n 5 n n 6 2 Bn,n 4(x1, x2, x3, x4, x5)= x1− x5 +5 x1− [3x2x4 +2x3] − 5 6 n n 7 2 n n 8 4 + 105 x − x x + 105 x − x . (9) 7 1 2 3 8 1 2 Some values of partial Bell polynomials are n B (0!, 1!, ..., (n k)!) = (Stirling number of first kind), n,k − k n 1 n! B (1!, ..., (n k)!) = − (Lah number), n,k − k 1 k! − n B (1, 1, ..., 1) = (Stirling number of second kind), n,k k n n k B (1, 2, ..., n k +1) = k − (Idempotent number). n,k − k Then we can see the beautiful relationship that exists between Bell polynomials and numbers like the above. On the other hand, the partial Bell polynomials can be efficiently computed by means of the recurrence relation n k+1 − n 1 Bn,k(x1, ..., xn k+1)= − xiBn i,k 1(x1, ..., xn i k+2). (10) − i 1 − − − − i=1 X − The definition of complete Bell polynomials is as follows Definition 3. The sum n

Bn(x1, x2, ..., xn)= Bn,k(x1, x2, ..., xn k+1) (11) − Xk=1 3 is called n-th complete exponential Bell polynomials with exponential generating function given by to make u = 1 in (1)

∞ tm ∞ tn exp x = B (x , x , ..., x ) (12) m m! n 1 2 n n! m=1 ! n=0 X X and B0 = 1. Some complete Bell polynomials are

B1(x1)= x1, 2 B2(x1, x2)= x1 + x2, 3 B3(x1, x2, x3)= x1 +3x1x2 + x3, 4 2 2 B4(x1, x2, x3, x4)= x1 +6x1x2 +4x1x3 +3x2 + x4, 5 3 2 2 B5(x1, x2, x3, x4, x5)= x1 + 10x1x2 + 15x1x2 + 10x1x3 + 10x2x3 +5x1x4 + x5.

Theorem 4. The complete Bell polynomials Bn satisfy the identity n n Bn+1(x1, ..., xn+1)= Bn i(x1, ..., xn i)xi+1. (13) i − − i=0 X From this it follows that

B2n+1(0, x2, 0, .., 0, x2n+1)=0 (14) for all n 0. ≥ Another useful identity that Bell polynomials fulﬁll is as follows

n 1 n B ( x , x , x , ..., ( 1) − x )=( 1) B (x , x , x , ..., x ) (15) n − 1 2 − 3 − n − n 1 2 3 n 3 Exponential autonomous functions

We will study the solution of the equation y(k) = eay (16) for any a C. Making y = u/a we obtain the equivalent equation ∈ u(k) = aeu (17) Then without loss of generality we will focus on the equation (17). Now by applying derivative to (17) we obtain another equation equivalent to (16)

(k+1) u (k) u = ae u′ = u u′ (18)

(k 1) Denote (x1, x2, ..., xk) the initial value problem y(0) = x1, y′(0) = x2,...,y − (0) = xk. In this section the general solution and solutions with initial values (x, 0, 0, ..., 0), k 1 (x1 + k ln c,cx2, ..., c − xk), and (x1, x2, ..., x2k 1, x2k) of the equation (17) are − − − given. We will define the complete and partial exponential autonomous functions and the exponential autonomous polynomials, which are the coefficients of the power series solution of the equation (17). Moreover, we will find special values of these functions by using Bell polynomials. We begin with the following definition

4 Deﬁnition 5. Take a C. Suppose x = (x1, ..., xk). Let fn(x, a) denote the complete exponential autonomous∈ functions of order k, k 1, recursively deﬁned as ≥

f0(x, a)= x1, (19)

f1(x, a)= x2, (20) . .

fk 1(x, a)= xk, (21) − x1 fk(x, a)= ae , (22) f (x, a)= aex1 B (f (x, a), ..., f (x, a)), n 1, (23) n+k n 1 n ≥ where Bn(y1, ..., yn) are the complete Bell polynomials. When x1 = 0, we deﬁne the exponential autonomous polynomials as q (x , ..., x )= f (0, x , ..., x ), for n 1. n 2 k n 2 k ≥

When a = 1 in the above deﬁnition, we will write fn(x)= fn(x, 1). In this section we will restrict ourselves to exponential autonomous functions and autonomous polynomials will be dealt with in the next section. The following are complete exponential autonomous functions for k =1, 2, 3, 4. They will be very useful in the next section: When k = 1, f (x, a)=(n 1)!anenx. When k =2 n −

f0(x, y, a)= x,

f1(x, y, a)= y, x f2(x, y, a)= ae , x f3(x, y, a)= aye , x x 2 f4(x, y, a)= ae (ae + y ), x x 3 f5(x, y, a)= ae (4aye + y ), x 2 2x 2 x 4 f6(x, y, a)= ae (4a e + 11ay e + y ).

When k =3

f0(x,y,z,a)= x,

f1(x,y,z,a)= y,

f2(x,y,z,a)= z, x f3(x,y,z,a)= ae , x f4(x,y,z,a)= aye , x 2 f5(x,y,z,a)= ae (z + y ), x x 3 f6(x,y,z,a)= ae (ae +3yz + y ), x x 2 2 4 f7(x,y,z,a)= ae (5ye +3z +6y z + y ).

5 And ﬁnally, when k =4

f0(x,y,z,w,a)= x,

f1(x,y,z,w,a)= y,

f2(x,y,z,w,a)= z,

f3(x,y,z,w,a)= w, x f4(x,y,z,w,a)= ae , x f5(x,y,z,w,a)= aye , x 2 f6(x,y,z,w,a)= ae (z + y ), x 3 f7(x,y,z,w,a)= ae (w +3yz + y ), x x 2 2 4 f8(x,y,z,w,a)= ae (ae +3z +4yw +6y z + y ).

In the following result we will show that the exponential generating function of the complete exponential autonomous functions is solution of the equation (17)

Theorem 6. Let x =(x1, ..., xk). The series

∞ tn E (t, x, a)= f (x, a) (24) k n n! n=0 X is solution of the diﬀerential equation (17).

Proof. Taking derivative k times with respect to t the series Ek(t, x, a), using the deﬁnition of the autonomous functions fn(x, a) and the equation (1), then

∂kE (t, x, a) ∞ tn k = f (x, a) ∂tk n+k n! n=0 X ∞ tn = aex1 + aex1 B (f (x, a), ..., f (x, a)) n 1 n n! n=1 X ∞ tn = eax1 1+ B (f (x, a), ..., f (x, a)) n 1 n n! n=1 ! x1 E (t,x,aX) x1 = ae e k − x = aeEk(t, ,a).

Now we deﬁne the partial exponential autonomous functions

Deﬁnition 7. Let gn,i(x, a) denote the partial exponential autonomous functions as gn,i(x, a)= Bn,i(f1(x, a), ..., fn i+1(x, a)) (25) − with g (x, a) = 1, g (x, a) =0, for n 1, and g (x, a)=0, for i 1. Then 0,0 n,0 ≥ 0,i ≥ n x1 fn+k(x, a)= ae gn,i(x, a). (26) i=1 X 6 In the following result we establish recurrence relations for the functions fn(x, a) and gn,i(x, a). Many important results of this paper will be proved using this theorem.

Theorem 8. The autonomous functions fn(x, a) and gn,i(x, a) fulﬁll the recurrence relations n n fn+k+1(x, a)= fn i+k(x, a)fi+1(x, a) (27) i − i=0 X and n i+1 − n 1 gn,i(x, a)= − fj(x, a)gn j,i 1(x, a). (28) j 1 − − j=1 X − x1 Proof. Making yj = fj(x, a) in (10) and (13) and multiplying these by ae , we obtain the desired result.

Now we will study the behavior of the functions fn(x, a) evaluated at x = (x1, 0, 0, ..., 0). From previous result we can construct the ﬁrst important sequence arising from the diﬀerential equation (17)

Theorem 9. The functions fn(x, a) take the following values at x =(x1, 0, ..., 0)

1. fkn+1(x1, 0, ..., 0, a) = fkn+2(x1, 0, ..., 0, a) = = fkn+k 1(x1, 0, ..., 0, a)=0, ··· − n 0, ≥ (k) nx1 2. f (x , 0, ..., 0, a)= An (a)e , n 1 kn 1 ≥ (k) where A1 (a)=1 and

n (k) kn + k 1 (k) (k) An+2(a)= − An i+1(a)Ai+1(a), (29) ki + k 1 − i=0 X − n 0, k 1. ≥ ≥

Proof. Let x = (x1, 0, ..., 0). Clearly, f1(x, a) = 0, fk+1(x, a) = B1(f1(x, a)) = f (x, a) = 0. Now suppose it is true that f (x, a)=0 for2 i n 1. By the 1 ki+1 ≤ ≤ − theorem 8

fkn+1(x, a)= fk(n 1)+k+1(x, a) − k(n 1) − k(n 1) = − fkn i(x, a)fi+1(x, a). i − i=0 X

Since the product fkn i(x, a)fi+1(x, a) contain the functions fkj+1(x, a), then fkn+1(x, a)= − 0 for all n. Likewise it is proved for fkn+j(x, a), j = 2, ..., k 1. Now we will x1 2 2x1 − prove 2. We know that fk(x, a) = ae , f2k(x, a) = a e and suppose that

7 (k) nx1 fkn(x, a)= An (a)e . Then

fkn+k(x, a)= f(kn 1)+k+1(x, a) − kn 1 − kn 1 = − fkn 1 i+k(x, a)fi+1(x, a) i − − i=0 nX1 − kn 1 = − fk(n i)(x, a)fki+k(x, a) ki +1 − i=0 Xn 1 − kn 1 (k) (n i)x1 (k) (i+1)x1 = − An i(a)e − Ai+1(a)e ki +1 − i=0 X n 1 − (n+1)x1 kn 1 (k) (k) = e − An i(a)Ai+1(a) ki +1 − i=0 X (n+1)x1 (k) = e An+1(a).

(1) n It is easy to show that An (a)=(n 1)!a when k = 1. We will use the equation − (1) i (29) to prove this result. Suppose it is true that Ai (a)=(i 1)!a for i ranging between 1 and n + 1. We have − n (1) n (1) (1) An+2(a)= An i+1(a)Ai+1(a) i − i=0 Xn n n i+1 i+1 = (n i)!a − i!a i − i=0 X n n = an+2 (n i)!i! i − i=0 X = an+2n!(n +1)= an+2(n + 1)!. We can extend the above result to all k 1 ≥ Proposition 1. For all k 1 we have ≥ (k) n (k) An (a)= a An (1) (30)

(k) i (k) Proof. Suppose by induction that Ai (a) = a Ai (1) for all i ranging between 1 and n + 1, then use the same steps as in the previous proof. From the above proposition it follows that

Ek(t, (x, 0, ..., 0), a)= Ek(at, (x, 0, ..., 0), 1). (31) Then without loss of generality it is suﬃcient to study the solution of (17) with (k 1) initial conditions y(0) = x, y′(0) = y′(0) = = y − (0) = 0 and a = 1 to (k) ··· generate the sequence An (1). (k) The following corollary of theorem 9 shows us that the numbers An can be constructed using Bell polynomials

8 (k) Corollary 1. Numbers An (a) fulﬁll the recurrence relation

k 1 k 1 − − (k) (k) (k) An (a)= Bn k(0, ..., 0, A1 (a), ..., 0, ..., 0, An k(a)), n 1. (32) − − ≥ In the following theoremz }| we{ calculate somez }| special{ values of the functions gn.i(x, a) Theorem 10. Let x =(x , 0, ..., 0). For all a R we have 1 ∈ 1. g (x, a)=0, si k n. n,i 6| (k) lx1 2. glk,1(x, a)= Al (a)e .

lx1 l kl 1 (k) (k) 3. glk,2(x, a)= e j=1 kj−1 Aj (a)Al j(a). − −

4. gn,n(x, a) = gn,nP1(x, a) = gn,n 2(x, a) = gn,n 3(x, a) = gn,n 4(x, a)=0, − − − − k > 1.

Proof. Suppose k n and gn j,i 1(x, a) = 0 for all k such that k j. Using theorem 6| − − 6| 8 and 9 we have that fj(x, a) = 0. This proves 1. To prove 2 we have

lk lk 1 glk,1(x, a)= − fj(x, a)glk j,0(x, a) j 1 − j=1 X − lk 1 = − f (x, a)g (x, a) lk 1 lk 0,0 − (k) lx1 = Al (a)e . On the other hand,

lk 1 − lk 1 glk,2(x, a)= − fj(x, a)glk j,1(x, a) j 1 − j=1 X − l lk 1 = − fkj(x, a)glk kj,1(x, a) jk 1 − j=1 X − l lk 1 (k) jx1 (k) (l j)x1 = − Aj (a)e Al j(a)e − jk 1 − j=1 X − l lx1 lk 1 (k) (k) = e − Aj (a)Al j(a). jk 1 − j=1 X − Then this proves 3. To prove 4 we use the equations (5)-(9). We conclude this section with the following properties of the exponential autonomous functions Theorem 11. For all n 1, k 1 and for all a, c C is fulﬁlled ≥ ≥ ∈ k 1 n fn(x1 + k ln c,cx2, ..., c − xk, a)= c fn(x1, x2, ..., xk, a) (33)

9 k 1 Proof. Let y = (x1 + k ln c,cx2, ..., c − xk) and x = (x1, x2, ..., xk). Suppose that the result is true for i n. Then ≤

fn+1(y, a)= f(n+1 k)+k(y, a) − k x1 = ac e Bn+1 k(f1(y, a), ..., fn+1 k(y, a)) − − k x1 n+1 k = ac e Bn+1 k(cf1(x, a), ..., c − fn+1 k(x, a)) − − k n+1 k x1 = ac c − e Bn(f1(x, a), ..., fn(x, a)) n+1 = c fn+1(x, a).

The following is the corollary to the theorem 11 that allows to calculate the k 1 solutions of (17) when the initial values are (x1 + k ln c,cx2, ..., c − xk) Corollary 2.

k 1 Ek(t, (x1 + k ln c,cx2, ..., c − xk), a)= k ln c + Ek(ct, (x1, x2, ..., xk), a) (34)

k 1 Proof. Let y =(x1 + k ln c,cx2, ..., c − xk) and x =(x1, x2, ..., xk). From the above theorem and the deﬁnition of the function Ek(t, x, a) we have

∞ tn E (t, y, a)= x + k ln c + f (y, a) k 1 n n! n=1 X ∞ tn = x + k ln c + cnf (x, a) 1 n n! n=1 X ∞ (ct)n = x + k ln c + f (x, a) 1 n n! n=1 X = k ln c + Ek(ct, x, a).

Finally, we compute fn(x, a) when x =(x1, x2, ..., x2k 1, x2k) − − − Theorem 12. For all n 0 and for all exponential autonomous functions of order 2k it is satisﬁed that ≥

n fn((x1, x2, ..., x2k 1, x2k), a)=( 1) fn((x1, x2, ..., x2k 1, x2k), a) (35) − − − − −

Proof. Let y = (x1, x2, ..., x2k 1, x2k), and let x = (x1, x2, ..., x2k 1, x2k). Sup- − − pose it is true for all− values less than− or equal to n. Then

fn+1(y, a) = f(n+1 2k)+2k(y, a) − x1 = ae Bn+1 2k(f1(y, a), f2(y, a), ..., fn+1 2k(y, a)) − − x1 n+1 2k = ae Bn+1 2k( f1(x, a), f2(x, a), ..., ( 1) − fn+1 2k(x, a)) − − n+1 2k −x1 − = ( 1) − ae Bn+1 2k(f1(x, a), f2(x, a), ..., fn+1 2k(x, a)) − − − = ( 1)n+1f (x, a). − n

10 Finally we have the corollary to theorem 12 Corollary 3.

E2k( t, (x1, x2, ..., x2k 1, x2k), a)= E2k(t, (x1, x2, ..., x2k 1, x2k), a). (36) − − − − − Proof. From the above theorem and the deﬁnition of the function E2k(t, x, a) we have

∞ ( t)n E ( t, y, a)= f (y, a) − 2k − n n! n=0 X ∞ tn = ( 1)n( 1)nf (x, a) − − n n! n=0 X ∞ tn = f (x, a) n n! n=0 X = E2k(t, x, a) where y =(x1, x2, ..., x2k 1, x2k) and x =(x1, x2, ..., x2k 1, x2k). − − − − 4 (k,a)-autonomous coeﬃcients

y When k = 1 we obtain the equation y′ = ae , which is the easiest to solve for all k. Using the method of separation of variables we reach the solution

x y(t)= ln(e− at) − − with initial condition y(0) = x. On the other hand, by the theorem 6 the solution in power series becomes

∞ tn E (t, x, a)= x + A(1)(a) 1 n n! n=1 X ∞ tn = x + (n 1)!anenx − n! n=1 X ∞ (aex)n = x + = x ln(1 aext) n − − n=1 X Now we can use the results of the previous section to prove some results already known. By the definition of complete exponential autonomous functions n n+1 a(n+1)x x 1 x 2 2x n i+1 (n i+1)x n!a e = ae B (0!a e , 1!a e , ..., (n i)!a − e − ) n,i − i=1 X n = aexanenx B (0!, 1!, ..., (n i)!) n,i − i=1 X n n n = ea(n+1)xan+1 B (0!, 1!, ..., (n i)!) = ea(n+1)xan+1 n,i − i i=1 i=1 X X 11 from which follows the result relating factorials and Stirling number of first kind n n n!= . (37) i i=1 X n Furthermore gn,i(x, 1) = i and by the equation (28) we obtain the following finite- sum identity n i+1 n +1 − n n +1 j = (j 1)! − i +1 j 1 − i j=1 X − i n! j = j! i j=n X n n! j = j! i j=i X n n! j = . j! i j=0 X On the other hand, from the equation (27) we obtain the trivial result n n (n +1)! = (n i)!i! i − i=0 Xn = n!. i=0 X The Stirling numbers of the first kind originally arose algebraically from the expansion of the falling factorial (x) = x(x 1)(x 2) (x n + 1) n − − ··· − and in polynomial form is as follows n n k n i (x) = ( 1) − x . n − i i=0 X Analogously, we want to define and study the coefficients of the expansion of the autonomous exponential polynomials qn(x, a) with x = (0, x, x, ..., x). First we calculate the degree of qn(x, a).

Proposition 2. Let x = (0, x, x, ..., x). Then the degree qn(x, a) of gr(q (x, a)) = n k, n k. (38) n − ≥ Proof. By definition n 1 − i n qn+k(x, a)= a gn,i(x, a)+ a gn,n(x, a) i=1 Xn 1 − i n n = a gn,i(x, a)+ a x1 . i=1 X As gr(g (x, a)) i, then gr(q (x, a)) = n. n,i ≤ n+k 12 We now define the autonomous polynomials and autonomous coefficients

(k) Deﬁnition 13. Let An (x, a)= qn(0, x, ..., x, a) denote the autonomous polynomials of degree n k for all n k. − ≥ Using (23) we note that

(k) (k) (k) An+k(x, a)= aBn(A1 (x, a), ..., An (x, a)) (39) for all n 1. ≥ n Definition 14. We define the (k, a)-autonomous coefficients, denoted by i (k,a), (k) q y as the coefficients of the autonomous polynomials An+k(x, a), i.e.,

n n A(k) (x, a)= xi. (40) n+k si{ i=0 (k,a) X Now we will give some values of the (k, a)-autonomous coeﬃcients

n Theorem 15. Some values of the coeﬃcients i (k,a) are q y n 0, si k n; = 6| (41) s0{ an/kA(k) (1), si k n (k,a) ( n/k | 0 =0, si i 1, (42) si{(k,a) ≥ n n = a , k>l +1, 0 l < n, (43) sn l{ n l ≤ − (k,a) −

Proof. The equation (41) follows from theorem 9. By deﬁnition, B0,i = 0 for i 1. Then (42) is true. Finally, if k>l + 1, ≥

(k) (k) Bn,n l(A1 (x, a), ..., Al+1(x, a)) = Bn,n l(x, ..., x) − − n n l = x − n l − and from here follows (43). We will now show the relationship between the (k, 1)-autonomous coeﬃcients and the binomial coeﬃcients

Theorem 16. n +1 n n +1 k = − s 0 { k 1 s 0 { (k,1) − (k,1) n n n h h +1 k + − − , (44) h s 0 {(k,1)s 0 {(k,1) hX=k+1

13 for 1 i n k +3 ≤ ≤ − k 2 n +1 − n n h = − s i {(k,1) h s i 1 {(k,1) Xh=0 − n n +1 k n n k + − + − k 1 s i { k s i 1 { − (k,1) − (k,1) n n n h h +1 k + − − (45) h s j {(k,1)s l {(k,1) h=Xk+1 jX+l=i and for n k +4 i n +1 − ≤ ≤ n i+1 n +1 − n n h = − s i {(k,1) h s i 1 {(k,1) Xh=0 − n n +1 k n n k + − + − k 1 s i { k s i 1 { − (k,1) − (k,1) n n n h h +1 k + − − (46) h s j {(k,1)s l {(k,1) h=Xk+1 jX+l=i Proof. As

n (k) n (k) (k) An+1+k(x, 1) = An i+k(x, 1)Ai+1(x, 1) i − i=0 Xk 2 − n (k) n (k) = An i+k(x, 1)x + An+1(x, 1) i − k 1 i=0 − X n n (k) n (k) (k) + An (x, 1)x + An i+k(x, 1)Ai+1(x, 1), k i − i=Xk+1 then n+1 k 2 n i n +1 − n − n i xi = − xj+1 s i { i s j { i=0 (k,1) i=0 j=0 (k,1) X X X n+1 k n k n − n +1 k n − n k + − xj + − xj+1 k 1 s j { k s j { j=0 (k,1) j=0 (k,1) − X X n n i i+1 k n − n i − i +1 k + − xj − xl . i s j { s l { j=0 (k,1) (k,1) ! i=Xk+1 X Xl=0

14 We multiply the two autonomous polynomials within the last sum n+1 k 2 n i n +1 − n − n i xi = − xj+1 s i { i s j { i=0 (k,1) i=0 j=0 (k,1) X X X n+1 k n k n − n +1 k n − n k + − xj + − xj+1 k 1 s j { k s j { j=0 (k,1) j=0 (k,1) − X X n n+1 k n − n i i +1 k + − − xh. i s j {(k,1)s l {(k,1)! i=Xk+1 Xh=0 jX+l=h Then by rearranging the ﬁrst and fourth sums we obtain

n+1 n k+3 k 2 n +1 − − n n h xi = − xi s i { h s i 1 { i=0 (k,1) i=0 h=0 − (k,1)! X Xn+1 Xn+1 i − n n h + − xi h s i 1 { i=n k+4 h=0 (k,1)! X− X − n+1 k n − n +1 k + − xi k 1 s i { − i=0 (k,1) n k+1X n − n k + − xi k s i 1 { i=1 − (k,1) n+1 k Xn − n n h h +1 k + − − xi. h s j { s l { i=0 (k,1) (k,1)! X hX=k+1 jX+l=i For a suitable value of i the desired results are attained. Finally, we show without proof the relationship between the (k, 1)-autonomous coeﬃcients and the Stirling numbers of second kind.

Conjecture 1. Suppose that A(k)(1, 1) = = A(k)(1, 1) = 1. Then 1 ··· k n n B (A(k)(1, 1), ..., A(k)(1, 1)) = A(k)(1, 1), n 1 (47) n 1 n i i ≥ i=1 X Then, n n A(k) (1, 1) = A(k)(1, 1), n 1 (48) n+k i i ≥ i=1 X and n k n j k n n n − n = + (49) s i{ j j si { i=0 (k,1) j=1 i=0 (k,1) X X j=Xk+1 X Equation (48) corresponds to the number of shifts left k 1 places under Stirling − transform.

15 5 Sequences related to the equation (17)

We conclude this article by showing sequences related to equation (17) for values of k =2, 3, 4. Especially, we show that the numbers known as reduced tangent numbers, Bernoulli numbers, Euler zigzag numbers, Eulerian numbers, Blasius numbers, triangular numbers, number of shifts left 3 places under Stirling transform, and number of 8-sequences of [1, n] with 2 contiguous pairs can be constructed using Bell polynomials, Stirling numbers of second kind, binomial coeﬃcient and autonomous coeﬃcients.

5.1 Case k=2 The ﬁrst case to be studied is y y′′ = ae (50) (3) The equation (50) is equivalent to the equation y = y′′y′, whose solution is 1 1 y′ = √2√c tan √2√c t + √2√c c 1 2 1 2 1 2 and therefore 1 1 y = ln sec2 √2√c t + √2√c c + c (51) 2 1 2 1 2 3 where c1,c2 and c3 are constants in C. Since we want y(0) = x, y′(0) = y and x y′(0) = ae , then 1 ln sec2 √2√c c + c = x, 2 1 2 3 1 √2√c tan √2√c c = y, 1 2 1 2 1 c sec2 √2√c c = aex. 1 2 1 2 As result y2 c = aex , 1 − 2 1 y √2√c1c2 = arctan , 2 2aex y2 ! − y2 c = x ln 1+p . 3 − 2aex y2 − Thus, the function

x 2 2 2ae y t y E2(t, (x, y), a)= x + ln sec − + arctan 2 2aex y2 !! p − y2 ln 1+ p (52) − 2aex y2 −

16 is the solution of the equation (50) with initial value y(0) = x, y′(0) = y. The following is a list of particular solutions of (50) which are obtained from the equation (52)

√aex/2t E (t, (x, 0), a)= x + ln sec2 , a> 0, (53) 2 √ 2 √aex/2t E (t, (x, 0), a)= x + ln sech2 , a> 0, (54) 2 − √ 2 2 2 2a y t y E2(t, (0,y), a) = ln sec − + arctan 2 2a y2 !! p − y2 ln 1+ , a> 0, p (55) − 2a y2 − 2 2 2a + y t y E2(t, (0,y), a) = ln sech + arctanh 2 − p 2 2a + y !! y2 ln 1 , a> 0. p (56) − − 2a + y2 We now show the relationship between reduced tangent numbers (A002105 in OEIS) and Bell polynomials and binomial coeﬃcients Theorem 17. Let (Tn)n 1 = (1, 1, 4, 34, 496,...) (57) ≥ denote the sequence of reduced tangent numbers. Then

(2) n 1. An (a)= a Tn.

2. Tn = Bn(0, T1,..., 0, Tn 1), n 2. − ≥ n n 1 3. ( 1) Tn = Bn(0, T1,..., 0, ( 1) − Tn 1), n 2. − − − − ≥ n 2n+1 4. Tn+2 = Tn i+2Ti+1, n 0. i=0 2i+1 − ≥ Proof. AnotherP way to write (53) is

x/2 √ae t u E2(t, (x, 0), a)= x + √2 tan du. √2 Z0 Then

x/2 √ae t ∞ u2n 1 E (t, (x, 0), a)= x + T − 2 n (2n 1)! 0 n=1 Z X − ∞ (√aex/2t)2n = x + T . n (2n)! n=1 X n nx By comparison, f2n(x, 0, a) = a Tne . Thus follows 1, 2, and 3. The recurrence relation 4 follows from equation (29).

17 In general, the solution (53) is the generating function of the sequence

n 2 3 4 5 (a Tn)n 1 =(a, a , 4a , 34a , 496a ,...). ≥ n 2n 2 (2 1) b2n − | | On the other hand, it is known that Tn = n , where the b2n are the Bernoulli numbers (A000367, A002445 in OEIS). Then the theorem 17 provides a relation between Bell polynomials and Bernoulli numbers, that is

n 2n n 1 2n 2 2 (2 1) b2n 2 − (2 − 1) b2n 2 − | | = B 0, 6 b , 0, 30 b ,..., 0, − | − | (58) n n | 2| | 4| n 1 − We now show the relationship between Euler zigzag numbers (A000111 in OEIS) and Bell polynomials, binomial coeﬃcients, and Stirling numbers of second kind Theorem 18. Suppose

(An)n 0 = (1, 1, 1, 2, 5, 16, 61, 272,...) (59) ≥ the sequence of Euler zigzag numbers. Then

1. An+1 = Bn(A0,...,An 1), n 1. − ≥ n n 1 2. ( 1) An+1 = Bn( A0, A1,..., ( 1) − An 1), n 1. − − − − ≥ n n 3. An+2 = An i+1Ai, =0. i=0 i − 6 4. A = Pn n A , n 1. n+2 i=1 i i ≥ P n 2n 5. A2n+2 = Tn + i=1 2i (2,1). q y 6. A = n P2n+1 2n+3 i=1 2i+1 (2,1) q y Proof. From equationP (55), t π E (t, (0, 1), 1) = ln sec2 + ln(2) 2 2 4 − = ln(sec2(t) + sec(t) tan(t)). Then t E2(t, (0, 1), 1) = (sec(u) + tan(u))du Z0 t ∞ un = 1+ A du n n! 0 n=1 ! Z X ∞ t un = t + A du n n! n=1 0 X Z ∞ tn+1 = t + A n (n + 1)! n=1 X ∞ tn = An 1 . − n! n=1 X 18 We apply the equation (23) to obtain 1. By corollary 3 it follows that E (t, (0, 1), 1) = E( t, (0, 1), 1) = ln(sec2(t) sec(t) tan(t)). (60) 2 − − − From equation (15) follows 2. Formula 3 follows from equation (27). From equation (48) follows 4. The identities 5 and 6 follow because the Euler zigzag numbers are (2) obtained when x = 1 in An (x, 1). When k = 2 the exponential autonomous polynomials and the autonomous polynomials match. Some autonomous polynomials of the equation (50) are

q1(y, a) = y,

q2(y, a) = a,

q3(y, a) = ay, 2 q4(y, a) = a(a + y ), 3 q5(y, a) = a(4ay + y ), 2 2 4 q6(y, a) = a(4a + 11ay + y ), 2 3 5 q7(y, a) = a(34a y + 26ay + y ), 3 2 2 4 6 q8(y, a) = a(34a + 180a y + 57ay + y ). From the above we obtain the ﬁrst (2, a)-autonomous coeﬃcients

i 0 1 2 3 4 5 6 n 0 a 1 0 a 2 a2 0 a 3 0 4a2 0 a 4 4a3 0 11a2 0 a 5 0 34a3 0 26a2 0 a 6 34a4 0 180a3 0 57a2 0 a

Table 1: (2, a)-autonomous coeﬃcients

Theorem 19. Some values of (2, a)-autonomous coeﬃcients are 2n 2n +1 = =0 (61) s2i +1{(2,a) s 2i {(2,a) for all i. Proof. The equation (61) follows from theorem 9. Conjecture 2. n = a2(2n n 1) (62) sn 2{ − − − (2,a) The sequence 2n n 1=(0, 0, 1, 4, 11, 26, 57, 120, 247, 502, 1013, 2036, 4083, 8178, 16369, − − 32752, 65519, 131054, 262125, 524268, 1048555, 2097130,...) is known as Eulerian numbers (A000295 in OEIS).

19 5.2 Case k=3 When k = 3 we obtain the equation y(3) = aey (63) Solving (63) with initial conditions (0, 0, x) and a = 1 we get the solution of Blasius equation − (3) u + u′′u =0. (64) The Blasius equation [14] describes the velocity profile of the fluid in the boundary layer which forms when fluid flows along a flat plate. Using the theorem 9 and the corollary 1 we reach the following result on Blasius numbers (A018893 in OEIS) Theorem 20. Let (bn)n 1 = (1, 1, 11, 375, 27.897,...) (65) ≥ denote the sequence of Blasius numbers. Then

1. bn = Bn(0, 0, b1,..., 0, 0, bn 1), n 2. − ≥ n 3n+2 2. bn+2 = bn i+1bi+1, n 0. i=0 3i+2 − ≥ On the otherP hand, the autonomous polynomials for the equation (63) are (3) A3 (x, a)= a, (3) A4 (x, a)= ax, (3) 2 A5 (x, a)= a(x + x ), (3) 2 3 A6 (x, a)= a(a +3x + x ), (3) 2 3 4 A7 (x, a)= a(5ax +3x +6x + x ), (3) 2 3 4 5 A8 (x, a)= a(11ax + 16ax + 15x + 10x + x ), (3) 2 2 3 4 5 6 A9 (x, a)= a(11a + 84ax + (42a + 15)x + 45x + 15x + x ), (3) 2 2 3 4 5 6 7 A10 (x, a)= a(117a x + 129ax + 384ax + (99a + 105)x + 105x + 21x + x ) and from here we obtain the following table with the ﬁrst (3, a)-autonomous coeﬃcients i 0 1 2 3 4 5 6 7 n 0 a 1 0 a 2 0 a a 3 a2 0 3a a 4 0 5a2 3a 6a a 5 0 11a2 16a2 15a 10a a 6 11a3 0 84a2 42a2 + 15a 45a 15a a 7 0 117a3 129a2 384a2 99a2 + 105a 105a 21a a

Table 2: (3, a)-autonomous coeﬃcients

20 Theorem 21. Some values of (3, a)-autonomous coeﬃcients are

n = a, sn{(3,a) n n = a . sn 1{ 2 − (3,a) Proof. The results follow from theorem 15 with l = 0, 1 and by keeping in mind n n that n 1 = 2 . − Conjecture 3. n n = a 2 (66) sn 2{ 2 − (3,a) The numbers n are the triangular numbers n 2 (3,1) q − y (0, 0, 3, 15, 45, 105, 210, 378, 630, 990, 1485,...) (67)

(A050534 in OEIS).

Finally, by the equations (15), (39) and (48) we have

Theorem 22. Let

(3) (en)n 1 =(An (1, 1))n 1 = (1, 1, 1, 1, 2, 5, 15, 53, 213,...). (68) ≥ ≥ denote the number of shifts 3 places left under exponentiation (A007548 in OEIS). Then

1. e = B (e ,..., e ), n 1. n+3 n 1 n ≥ n n 1 2. ( 1) e = B ( e , e ,..., ( 1) − e ), n 1. − n+3 n − 1 2 − n ≥ 3. e = n n e , n 1. n+3 i=1 i i ≥ P 3n 3n 4. d3n = bn + i=1 i (3,1). q y 3Pn+j 3n+j 5. d3n+j = i=1 i (3,1), j =1, 2. q y P 5.3 Case k=4 The equation to be studied is y(4) = aey. (69) This equation is not commonly studied in the literature. Here we show the relation of this equation with the number of shifts left 3 places under Stirling transform,

21 (4) and also the relation with the numbers An (1). A list of exponential autonomous polynomials of the equation (69) is as follows:

q1(y,z,w,a)= y,

q2(y,z,w,a)= z,

q3(y,z,w,a)= w,

q4(y,z,w,a)= a,

q5(y,z,w,a)= ay, 2 q6(y,z,w,a)= a(z + y ), 3 q7(y,z,w,a)= a(w +3yz + y ), 2 2 4 q8(y,z,w,a)= a(a +3z +4yw +6y z + y ).

(4) From the equation (29) we calculate the ﬁrst numbers An (1),

(4) q4(0, 0, 0, 1) = A1 (1) = 1, 3 q (0, 0, 0, 1) = A(4)(1) = A(4)A(4) =1, 8 2 3 1 1 7 7 q (0, 0, 0, 1) = A(4)(1) = A(4)A(4) + A(4)A(4) = 35, 12 3 3 2 1 7 1 2 11 11 11 q (0, 0, 0, 1) = A(4)(1) = A(4)A(4) + A(4)A(4) + A(4)A(4) = 6140. 16 4 3 3 1 7 2 2 11 1 3 Following theorem 9, corollary 1 and equation (15) we have the following recurrence (4) relations for the numbers An (1)

(4) Theorem 23. Let (cn)n 1 =(An (1))n 1 = (1, 1, 35, 6140,...). Then ≥ ≥

1. cn = Bn(0, 0, 0, c1,..., 0, 0, 0, cn 1), n 2. − ≥ n 4n+3 2. cn+2 = cn i+1ci+1, n 0. i=0 4i+3 − ≥ n n 3. ( 1) cnP= Bn(0, 0, 0, c1,..., 0, 0, 0, ( 1) cn 1) − − − − The autonomous polynomials associated with the equation (69) are

(4) (4) (4) A1 (x, a)= A2 (x, a)= A3 (x, a)= x, (4) A4 (x, a)= a, (4) A5 (x, a)= ax, (4) 2 A6 (x, a)= a(x + x ), (4) 2 3 A7 (x, a)= a(x +3x + x ), (4) 2 3 4 A8 (x, a)= a(a +7x +6x + x ), (4) 2 3 4 5 A9 (x, a)= a(6ax + 10x + 25x + 10x + x ), (4) 2 3 4 5 6 A10 (x, a)= a(16ax + 32ax + 75x + 65x + 15x + x ), (4) 2 3 4 5 6 7 A11 (x, a)= a(36ax + 136ax + (64a + 175)x + 315x + 140x + 21x + x ).

22 (4) We now derive recurrence relations of the numbers An (1, 1) using the equations (15), (39), and (48). Theorem 24. Suppose

(4) (dn)n 1 =(An (1, 1))n 1 = (1, 1, 1, 1, 1, 2, 5, 15, 53, 222, 1115, 6698,...) (70) ≥ ≥ the number of shifts left 3 places under Stirling transform (A336020 in OEIS). Then 1. d = B (d ,..., d ), n 1. n+4 n 1 n ≥ n n 1 2. ( 1) d = B ( d , d ..., ( 1) − d ), n 1. − n+4 n − 1 2 − n ≥ 3. d = n n d , n 1. n+4 i=1 i i ≥ P 4n 4n 4. d4n = cn + i=1 i (4,1). q y 4Pn+j 4n+j 5. d4n+j = i=1 i (4,1), j =1, 2, 3. q y The followingP is a table of the ﬁrst (4, a)-autonomous coeﬃcients

i 0 1 2 3 4 5 6 7 n 0 a 1 0 a 2 0 a a 3 0 a 3a a 4 a2 0 7a 6a a 5 0 6a2 10a 25a 10a a 6 0 16a2 32a2 75a 65a 15a a 7 0 36a2 136a2 64a + 175 315a 140a 21a a

Table 3: (4, a)-autonomous coeﬃcients

Theorem 25. Some values of (4, a)-autonomous coeﬃcients are

n = a (71) sn{(4,a) n n = a (72) sn 1{ 2 − (4,a) n n +2 = a (73) sn 2{ n − (4,a) Proof. The equations (71)-(73) arise from theorem 15 with l =0, 1, 2. Conjecture 4. n 5a n = (n 1) , n 5 (74) sn 3{ 2 − 5 ≥ − (4,a) 23 The sequence n = (10, 75, 315, 980, 2520, 5670, 11550, 21780, 38610, 65065, sn 3{ − (4,1) 105105, 163800, 247520, 364140, 523260, 736440, 1017450, 1382535, 1850695, 2443980, 3187800, 4111250, 5247450,...) counts the number of 8-sequences of [1, n] with 2 contiguous pairs, (A027778 in the OEIS).

References

[1] E. T. Bell, Exponential polynomials, Ann. of Math. 35(2)(1934), 258–277.

[2] A. Bernardini and P. E. Ricci, Bell Polinomials and Diﬀerential Equations of Freud-Type Polynomials, Math. Comput. Model. 36(2002), 1115-1119.

[3] A. Caley, On the theory of analytical forms called trees, Philos. Mag. 13(19)(1857),4-9.

[4] L. Comtent, Analyse Combinatoire, Tome Premier, Presses Universitaires de France (1970).

[5] C. Fa`adi Bruno, Sullo sviluppo delle funzioni, Annali di Scienze Matematiche e Fisiche 6(1855), 479–480.

[6] L. Foissy, Fa´adi Bruno subalgebras of the Hopf of planar trees from combinatorial Dyson-Schwinger equations, Adv. Math. 218(1)(2008),136-162.

[7] H. W. Gould and J. Quaintance, A Linear Binomial Rrecurrence and the Bell Numbers and Polynomials, Appl. Anal. Discrete Math. 1(2007), 371-385.

[8] M. Mihoubi, M. S. Maamra, Touchard polynomials, partial Bell polynomials and polynomials of binomial type, J. Integer Seq., 14(2011) 11.3.1.

[9] R. Orozco, Solución Anal´ıtica de la Ecuación Diferencial Autónoma de orden k, Revista De Matemática Teor´ıa y Aplicaciones vol(23) fasc.1: 63-84.

[10] J. Riordan, Derivatives of composite functions, Bull. Amer. Math. Soc. 52(8)(1946), 664–667.

[11] J. Riordan, An Introduction to Combinatorial Analysis, John Wiley & Sons, Inc., New York (1958).

[12] G. Yıldız, G. Tınaztepe, and M. Sezer, Bell Polynomial Approach for the So- lutions of Fredholm Integro-Diﬀerential Equations with Variable Coeﬃcients, CMES-Computer Modeling in Engineering & Sciences, 123(3)(2020), 973–993.

[13] Y. Wang and Y. Chen, Bell polynomials approach for two higher-order KdV- type equations in ﬂuids, Nonlinear Anal. Real World Appl., 31(2016), 533-551.

24 [14] H. Blasius, Grenzschichten in Flussigkeiten mit kleiner Reibung, Zeitschr. Math. Phys. 56(1908), 1-37.

2020 Mathematics Subject Classiﬁcation: Primary 34A34; Secondary 11B37, 11B68, 11B73, 11B83. Keywords: Bell’s polynomial, autonomous function, autonomous polynomial, autonomous coeﬃcient, Stirling number.

(Concerned with sequences A000111, A000295, A000367, A002105, A002445, A007548, A018893, A027778, A050534, and A336020.)