OFFSET
0,2
COMMENTS
These matrices resemble Hankel matrices.
LINKS
P. Steinbach, Golden fields: a case for the heptagon, Math. Mag. 70 (1997), no. 1, 22-31.
Index entries for linear recurrences with constant coefficients, signature (23,21,-4,-1).
FORMULA
a(n) = 23a(n-1)+21a(n-2)-4a(n-3)-a(n-4) for n>=4; a(0)=1, a(1)=30, a(2)=707, a(3)=16886 (follows from the minimal polynomial x^4-23x^3-21x^2+4x+1 of the matrix M).
G.f.: -(x^3+4*x^2-7*x-1) / ((x+1)*(x^3+3*x^2-24*x+1)). [Colin Barker, Dec 07 2012]
MAPLE
with(linalg): M[1]:=matrix(4, 4, [10, 9, 7, 4, 9, 8, 6, 3, 7, 6, 4, 2, 4, 3, 2, 1]): for n from 2 to 15 do M[n]:=multiply(M[1], M[n-1]) od: 1, seq(M[n][1, 1]+M[n][1, 2]+M[n][1, 3]+M[n][1, 4], n=1..15);
a[0]:=1: a[1]:=30: a[2]:=707: a[3]:=16886: for n from 4 to 15 do a[n]:=23*a[n-1]+21*a[n-2]-4*a[n-3]-a[n-4] od: seq(a[n], n=0..15);
MATHEMATICA
M = {{10, 9, 7, 4}, {9, 8, 6, 3}, {7, 6, 4, 2}, {4, 3, 2, 1}}; v[1] = {1, 1, 1, 1}; v[n_] := v[n] = M.v[n - 1]; a1 = Table[v[n][[1]], {n, 1, 50}]
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Gary W. Adamson and Roger L. Bagula, Oct 17 2006
EXTENSIONS
Edited by N. J. A. Sloane, Nov 07 2006
STATUS
approved