A transform of
Gottfried Helms
Annette.Warlich at t-online.de
Fri Jan 26 14:54:29 CET 2007
Dear list -
I am currently analyzing the matrix G =
1 -1 0 0 0 0 0 0 0 0
. -1 2 -1 0 0 0 0 0 0
. . 1 -3 3 -1 0 0 0 0
. . . -1 4 -6 4 -1 0 0
. . . . 1 -5 10 -10 5 -1
. . . . . -1 6 -15 20 -15
. . . . . . 1 -7 21 -35
. . . . . . . -1 8 -28
. . . . . . . . 1 -9
. . . . . . . . . -1
which plays seemingly a multitalent role as generator
for basic number-theoretical sequences A() as
- bernoulli-numbers
- "eta"-numbers
- and all sequences like 1/binomial(k,n) (with varying k and fixed n)
if the sequence A() in question is taken as column-vector
scaled with its (column-index+1). OEIS: [1]
For instance , taking the column-vector of scaled
bernoulli-numbers
B = [ b0*1, b1*2, b2*3,... ] (where b1 = +1/2)
= [ 1*1, 1/2*2, 1/6*3,... ]
then
G * B = 0
as well as if I use
B = [ a0*1, a1*2, a2*3, ... ] (where a_k= "eta"-numbers)
= [ 1*1, 1/2*2, 0*3, -1/4*4, ...]
G * B = 0
as well as any of the reciprocals of the binomials
B = [ a0*1, a1*2, a2*3, ... ] (where a_k= binomial(k+1,1))
= [ 1*1, 1/2*2, 1/3*3, 1/4*4, ...]
B = [ a0*1, a1*2, a2*3, ... ] (where a_k= binomial(k+3,2))
= [ 1/3*1, 1/6*2, 1/10*3, 1/15*4, ...]
and so on
with the limit of
B = [ a0*1, a1*2, a2*3, ... ] (where a_k= 1/2^k)
= [ 1*1, 1/2*2, 1/4*3, 1/8*4, ...]
all
G * B = 0
means: all these very different sequences are eigenvectors of G.
(and whose generation matrix G allows to recursively compute
Bernoulli-, eta and the binomial-sequences by the same recursion
formula only using different parameters)
----------------------
I think, this alone is a remarkable property.
An additional surprise is to look at the inverse; this is
H = G^-1 =
1 -1 2 -5 14 -42 132 -429 1430 -4862
. -1 2 -5 14 -42 132 -429 1430 -4862
. . 1 -3 9 -28 90 -297 1001 -3432
. . . -1 4 -14 48 -165 572 -2002
. . . . 1 -5 20 -75 275 -1001
. . . . . -1 6 -27 110 -429
. . . . . . 1 -7 35 -154
. . . . . . . -1 8 -44
. . . . . . . . 1 -9
. . . . . . . . . -1
which is also in OEIS as the (transposed) triangle of
Catalan-numbers [2].
I tried to see, what the rowsums are; these are divergent,
but not too bad, so they can be Euler-summed. I approximated
it to
[r,r^2,r^3,..] where r=0.61..., the golden ratio
Having this amazing result I tried to find a general
expression for the transformation of the rows of this
triangle on a powerseries and I seem to have arrived at
a true powerseries-to-powerseries-transformation.
Assume a powerseries-vector V(x) = [1,x,x^2,x^3,...]~
then
H * x * V(x) = y * V(y)
where
y =sqrt( x + 1/4 ) - 1/2
for instance, if x = 3/4 then y = 1/2 and
H * [3/4, 9/16, 27/64... ] = [ 1/2 , 1/4, 1/8, ....]
Don't know, whether this is worth to be mentioned
in the database...
Gottfried Helms
-----------------------------
[1] A030528 Triangle read by rows: a(n,k)=binom(k,n-k).
The signed triangular matrix a(n,m)*(-1)^(n-m) is
the inverse matrix of the triangular Catalan convolution
matrix A033184(n+1,m+1),
n >= m >= 0, with A033184(n,m) := 0 if n<m.
Row sums A000045(n+1) (Fibonacci).
[2] A033184: Catalan triangle A009766 transposed.
http://www.research.att.com/~njas/sequences/A033184
[3]
Concerning the columnsum-/alternating columnsum-property of G and the
Fibonacci-relation:
- see also A030528, crossrefs 1
- [A010892] http://www.research.att.com/~njas/sequences/A010892
A010892:Inverse of 6th cyclotomic polynomial. A period 6 sequence.
:1, 1, 0, -1, -1, 0, 1, 1, 0, -1, -1, 0, 1
COMMENT :Any sequence b(n) satisfying the recurrence b(n)=b(n-1)-b(n-2) can be written as
b(n) = b(0) * a(n) + (b(1)-b(0)) * a(n-1) . (...)
REFERENCES :Paul Barry, A Catalan Transform and Related Transformations on Integer Sequences,
Journal of Integer Sequences, Vol. 8 (2005), Article 05.4.5.
More information about the SeqFan
mailing list