|
|
A080882
|
|
a(n)*a(n+3) - a(n+1)*a(n+2) = 2^n, given a(0)=1, a(1)=3, a(2)=7.
|
|
6
|
|
|
1, 3, 7, 22, 52, 164, 388, 1224, 2896, 9136, 21616, 68192, 161344, 508992, 1204288, 3799168, 8988928, 28357376, 67094272, 211662336, 500798464, 1579869184, 3738010624, 11792304128, 27900891136, 88018956288, 208255086592
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,2
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (-2*x^3 - x^2 + 3*x + 1)/(4*x^4 - 8*x^2 + 1).
a(n) = (7/24*3^(1/2) + 1/2)*((1 + sqrt(3)))^n + ( - 7/24*3^(1/2) + 1/2)*((1 - sqrt(3)))^n + ( - 1/24*3^(1/2))*( - (1 + sqrt(3)))^n + (1/24*3^(1/2))*( - ((1 - sqrt(3))))^n. - Richard Choulet, Dec 06 2008
|
|
MAPLE
|
a:= n-> (Matrix([[22, 7, 3, 1]]). Matrix(4, (i, j)-> if (i=j-1) then 1 elif j=1 then [0, 8, 0, -4][i] else 0 fi)^(n))[1, 4]: seq(a(n), n=0..26); # Alois P. Heinz, Aug 23 2008
|
|
MATHEMATICA
|
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|