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

 

Logo

Annual Appeal: Today, Nov 11 2014, is the 4th anniversary of the launch of the new OEIS web site. 70,000 sequences have been added in these four years, all edited by volunteers. Please make a donation (tax deductible in the US) to help keep the OEIS running.

Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A108146 a(n)= 4*a(n-1) -a(n-2) -a(n-3). 0
1, 1, 1, 2, 6, 21, 76, 277, 1011, 3691, 13476, 49202, 179641, 655886, 2394701, 8743277, 31922521, 116552106, 425542626, 1553695877, 5672688776, 20711516601, 75619681751, 276094521627, 1008046888156, 3680473349246, 13437751987201 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,4

LINKS

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

FORMULA

G.f.: (1-3*x-2*x^2)/(1-4*x+x^2+x^3). [Sep 28 2009]

MATHEMATICA

p = Expand[(x^3 - x^2 - x - 1)*(x^3 - 4*x^2 + x + 1)] v[0] = {1, 1, 1}; M = {{0, 1, 0}, {0, 0, 1}, {-1, -1, 4}}; Det[M - x*IdentityMatrix[3]] NSolve[Det[M - x*IdentityMatrix[3]] == 0, x] v[n_] := v[n] = M.v[n - 1] a = Table[v[n][[1]], {n, 0, 50}]

CROSSREFS

Sequence in context: A116819 A116782 A112091 * A116798 A116821 A116772

Adjacent sequences:  A108143 A108144 A108145 * A108147 A108148 A108149

KEYWORD

nonn,easy

AUTHOR

Roger L. Bagula, Jun 05 2005

EXTENSIONS

Definition replaced by recurrence by the Associate Editors of the OEIS, Sep 28 2009

STATUS

approved

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

Content is available under The OEIS End-User License Agreement .

Last modified December 20 11:26 EST 2014. Contains 252241 sequences.