login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
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; text; internal format)
OFFSET

0,2

COMMENTS

Real Pisot roots (the eigenvalues of M): -0.541381, 5.54138. a(n)=3*A015536(n).

LINKS

Table of n, a(n) for n=0..20.

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, 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 * A260769 A026032 A005809

Adjacent sequences:  A106566 A106567 A106568 * A106570 A106571 A106572

KEYWORD

nonn

AUTHOR

Roger L. Bagula, May 30 2005

EXTENSIONS

Edited by N. J. A. Sloane, Apr 30 2006

STATUS

approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent | More pages
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 16 20:34 EDT 2018. Contains 316275 sequences. (Running on oeis4.)