login
Exponential Riordan array [exp(x*exp(-x)),x].
5

%I #15 Feb 13 2017 10:49:53

%S 1,1,1,-1,2,1,-2,-3,3,1,9,-8,-6,4,1,-4,45,-20,-10,5,1,-95,-24,135,-40,

%T -15,6,1,414,-665,-84,315,-70,-21,7,1,49,3312,-2660,-224,630,-112,-28,

%U 8,1,-10088,441,14904,-7980,-504,1134,-168,-36,9,1

%N Exponential Riordan array [exp(x*exp(-x)),x].

%C For commuting lower unitriangular matrices A and B we define A raised to the matrix power B, denoted A^^B, to be the matrix Exp(B*Log(A)). Here Exp denotes the matrix exponential and the matrix logarithm Log(A) is defined as sum {n >= 1} (-1)^(n+1)*(A-1)^n/n. Call the present triangle X and let P denote Pascal's triangle A007318. Then X solves the matrix equation X^^P = P. Equivalently, the infinite tower of matrix powers X^^(X^^(X^^(....))) equals P. Note that the infinite tower of powers P^^(P^^(P^^(...))) of the Pascal triangle equals the hyperbinomial array A088956. Thus we might view the present array as the hypobinomial triangle.

%H G. Helms, <a href="http://go.helms-net.de/math/tetdocs/PascalMatrixTetrated.pdf">Pascalmatrix tetrated</a>

%F T(n,k) = binomial(n,k)*A003725(n-k).

%F The triangle equals P^^Q, where P is Pascal's triangle and Q is the inverse of P. Column 0 equals A003725.

%F E.g.f.: exp(x*t)*exp(x*exp(-x)) = 1 + (1 + t)*x + (-1 + 2*t + t^2)*x^2/2! + (-2 - 3*t + 3*t^2 + t^3)*x^3/3! + ....

%F The infinitesimal generator for this triangle is the generalized exponential Riordan array [x*exp(-x),x], which factors as [x,x]*[exp(-x),x] = A132440*A007318^(-1). The infinitesimal generator begins

%F ..0

%F ..1....0

%F .-2....2....0

%F ..3...-6....3....0

%F .-4...12..-12....4....0

%F This is a signed version of the triangle of denominators from Leibniz's harmonic triangle - see A003506.

%e Triangle begins

%e .n\k.|....0.....1.....2.....3.....4.....5.....6.....7

%e = = = = = = = = = = = = = = = = = = = = = = = = = = =

%e ..0..|....1

%e ..1..|....1.....1

%e ..2..|...-1.....2.....1

%e ..3..|...-2....-3.....3.....1

%e ..4..|....9....-8....-6.....4.....1

%e ..5..|...-4....45...-20...-10.....5.....1

%e ..6..|..-95...-24...135...-40...-15.....6.....1

%e ..7..|..414..-665...-84...315...-70...-21.....7.....1

%e ...

%t max = 9; MapIndexed[ Take[#1, #2[[1]]]&, CoefficientList[ Series[ Exp[x*t]*Exp[x*Exp[-x]], {x, 0, max}, {t, 0, max}], {x, t}]*Range[0, max]!, 1] // Flatten (* _Jean-François Alcover_, Jan 08 2013 *)

%Y Cf. A003506, A003725 (column 0), A007318, A088956.

%K sign,easy,tabl

%O 0,5

%A _Peter Bala_, Sep 11 2012