|
%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
|
