|
| |
|
|
A106569
|
|
First entry of the vector (M^n)v, where M is the 2 X 2 matrix [[0,4],[1,4]] and v is the column vector [0,1].
|
|
0
| |
|
|
0, 3, 15, 84, 465, 2577, 14280, 79131, 438495, 2429868, 13464825, 74613729, 413463120, 2291156787, 12696173295, 70354336836, 389860204065, 2160364030833, 11971400766360, 66338095924299, 367604681920575
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 0,2
|
|
|
COMMENTS
| Real Pisot roots (the eigenvalues of M): -0.541381, 5.54138. a(n)=3*A015536(n).
|
|
|
FORMULA
| a(n)=first entry of v[n], where v[n]=Mv[n-1], M is the 2 X 2 matrix [[0, 4], [1, 4]] and v[0] is the column vector [0,1]. a(n)=5a(n-1)+3a(n-2); a(0)=0, a(1)=3.
a(n)=(3/37)*[5/2+(1/2)*sqrt(37)]^n*sqrt(37)-(3/37)*[5/2-(1/2)*sqrt(37)]^n*sqrt(37), with n>=0 [From Paolo P. Lava (paoloplava(AT)gmail.com), Oct 07 2008]
|
|
|
MAPLE
| a[0]:=0: a[1]:=3: for n from 2 to 23 do a[n]:=5*a[n-1]+3*a[n-2] od: seq(a[n], n=0..23);
|
|
|
MATHEMATICA
| M = {{0, 3}, {1, 5}} v[1] = {0, 1} v[n_] := v[n] = M.v[n - 1] a = Table[v[n][[1]], {n, 1, 50}]
|
|
|
CROSSREFS
| Equals 3*A015536.
Sequence in context: A193658 A195885 A115910 * A026032 A005809 A067122
Adjacent sequences: A106566 A106567 A106568 * A106570 A106571 A106572
|
|
|
KEYWORD
| nonn
|
|
|
AUTHOR
| Roger L. Bagula (rlbagulatftn(AT)yahoo.com), May 30 2005
|
|
|
EXTENSIONS
| Edited by N. J. A. Sloane (njas(AT)research.att.com), Apr 30 2006
|
| |
|
|
