%I #7 Nov 28 2018 19:09:50

%S 0,1,1,2,6,29,161,928,5388,31333,182269,1060358,6168762,35887625,

%T 208781333,1214615068,7066195896,41108599753,239155126777,

%U 1391318970218,8094196026894,47089136800949,273947751857609,1593730016014072

%N a(n) = 6*a(n-1)-6*a(n-3)-a(n-4).

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

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

%t m = 6 M = {{0, 1, 0, 0}, {0, 0, 1, 0}, {0, 0, 0, 1}, {-1, -m, 0, m}} Expand[Det[M - x*IdentityMatrix[4]]] NSolve[Det[M - x*IdentityMatrix[4]] == 0, x] v[1] = {0, 1, 1, 2}; v[n_] := v[n] = M.v[n - 1]; digits = 50; a = Table[v[n][[1]], {n, 1, digits}]

%t LinearRecurrence[{6,0,-6,-1},{0,1,1,2},30] (* _Harvey P. Dale_, Nov 28 2018 *)

%K nonn,easy

%O 0,4

%A _Roger L. Bagula_, May 24 2005

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