login
Triangle of arctanh numbers.
4

%I #32 Mar 28 2020 14:01:56

%S 1,0,1,0,0,1,0,2,0,1,0,0,8,0,1,0,24,0,20,0,1,0,0,184,0,40,0,1,0,720,0,

%T 784,0,70,0,1,0,0,8448,0,2464,0,112,0,1,0,40320,0,52352,0,6384,0,168,

%U 0,1,0,0,648576,0,229760,0,14448,0,240,0,1

%N Triangle of arctanh numbers.

%C Sheffer triangle associated to Sheffer triangle A060524.

%C For Sheffer triangles (matrices) see the explanation and S. Roman reference given under A048854.

%C The inverse matrix of A with elements a(n,m), n,m>=0, is given in A111593.

%C In the umbral calculus notation (see the S. Roman reference) this triangle would be called associated to (1,tanh(y)).

%C The row polynomials p(n,x):=sum(a(n,m)*x^m,m=0..n), together with the row polynomials s(n,x) of A060524 satisfy the exponential (or binomial) convolution identity s(n,x+y) = sum(binomial(n,k)*s(k,x)*p(n-k,y),k=0..n), n>=0.

%C Without the n=0 row and m=0 column and signed, this will become the Jabotinsky triangle A049218 (arctan numbers). For Jabotinsky matrices see the Knuth reference under A039692.

%C The row polynomials p(n,x) (defined above) have e.g.f. exp(x*arctanh(y)).

%C Exponential Riordan array [1, arctanh(x)] = [1, log(sqrt((1+x)/(1-x)))]. - _Paul Barry_, Apr 17 2008

%C Also the Bell transform of A005359. For the definition of the Bell transform see A264428. - _Peter Luschny_, Jan 27 2016

%H Wolfdieter Lang, <a href="/A111594/a111594.txt">First 10 rows.</a>

%F E.g.f. for column m>=0: ((arctanh(x))^m)/m!.

%F a(n, m) = coefficient of x^n of ((arctanh(x))^m)/m!, n>=m>=0, else 0.

%F a(n, m) = a(n-1, m-1) + (n-2)*(n-1)*a(n-2, m), a(n, -1):=0, a(0, 0)=1, a(n, m)=0 for n<m.

%e Binomial convolution of row polynomials:

%e p(3,x)= 2*x+x^3; p(2,x)=x^2, p(1,x)= x, p(0,x)= 1,

%e together with those from A060524:

%e s(3,x)= 5*x+x^3; s(2,x)= 1+x^2, s(1,x)= x, s(0,x)= 1; therefore:

%e 5*(x+y)+(x+y)^3 = s(3,x+y) = 1*s(0,x)*p(3,y) + 3*s(1,x)*p(2,y) + 3*s(2,x)*p(1,y) +1*s(3,x)*p(0,y) = 2*y+y^3 + 3*x*y^2 + 3*(1+x^2)*y + (5*x+x^3).

%p # The function BellMatrix is defined in A264428.

%p BellMatrix(n -> `if`(n::even, n!, 0), 10); # _Peter Luschny_, Jan 27 2016

%t rows = 10;

%t t = Table[If[EvenQ[n], n!, 0], {n, 0, rows}];

%t T[n_, k_] := BellY[n, k, t];

%t Table[T[n, k], {n, 0, rows}, {k, 0, n}] // Flatten (* _Jean-François Alcover_, Jun 22 2018, after _Peter Luschny_ *)

%o (Sage) # uses[riordan_array from A256893]

%o riordan_array(1, atanh(x), 9, exp=true) # _Peter Luschny_, Apr 19 2015

%Y Row sums: A000246.

%Y Cf. A005359, A049218, A060524, A111593.

%K nonn,easy,tabl

%O 0,8

%A _Wolfdieter Lang_, Aug 23 2005