login
Expansion of -(9520*x^9 -11504*x^8 +2840*x^7 -1040*x^6 +248*x^5 +36*x^4 -1) / ((2*x -1)*(4*x^2 +2*x -1)).
0

%I #28 Jun 13 2015 00:52:02

%S 1,4,16,56,156,248,1584,2248,18512,51856,189440,609664,2023808,

%T 6579712,21441536,69575680,225665024,731127808,2367905792,7666302976,

%U 24816189440,80321511424,259955621888,841292972032,2722599796736

%N Expansion of -(9520*x^9 -11504*x^8 +2840*x^7 -1040*x^6 +248*x^5 +36*x^4 -1) / ((2*x -1)*(4*x^2 +2*x -1)).

%H <a href="/index/Rec">Index entries for linear recurrences with constant coefficients</a>, signature (4,0,-8).

%F G.f.: -(9520*x^9 -11504*x^8 +2840*x^7 -1040*x^6 +248*x^5 +36*x^4 -1) / ((2*x -1)*(4*x^2 +2*x -1)). - _Colin Barker_, Mar 15 2013

%F a(n) = +4*a(n-1) -8*a(n-3). - _Colin Barker_, Mar 15 2013

%t M0 = {{2, 2, 0}, {1, 0, 1}, {0, 2, 2}}; M1 = {{1, 0, 1}, {2, 2, 0}, {0, 2, 2}}; M2 = {{2, 2, 0}, {0, 2, 2}, {1, 0, 1}}; M[n_] := M[n] = If[Mod[v[n][[1]], 3] == 0, M1, If[Mod[v[n][[2]], 3] == 0, M0, M2]] v[0] = {1, 1, 1}; M[0] = {{2, 2, 0}, {0, 2, 2}, {1, 0, 1}}; v[n_] := v[n] = M[n - 1].v[n - 1] a0 = Table[v[n][[1]], {n, 0, 25}]

%K nonn,easy,less

%O 0,2

%A _Roger L. Bagula_, Mar 03 2006

%E Explicit name by _Joerg Arndt_, Mar 15 2013