%I #12 Aug 09 2015 15:11:25
%S 0,4,3,15,16,55,73,206,315,787,1320,3051,5437,11950,22159,47135,89712,
%T 186799,361665,742638,1454035,2958619,5835448,11803235,23392341,
%U 47131150,93701703,188310471,375153280,752681287,1501520185,3009257486
%N Expansion of x*(4-x)/( (2x-1)*(x^2-x-1)).
%H <a href="/index/Rec">Index entries for linear recurrences with constant coefficients</a>, signature (1,3,-2)
%F a(n) = (7*2^n -A039834(n)-7*A039834(n+1))/5. - _R. J. Mathar_, Nov 05 2011
%p A039834 := proc(n)
%p if n <= 1 then
%p n;
%p else
%p -procname(n-1)+procname(n-2) ;
%p end if;
%p end proc:
%p A120461 := proc(n)
%p (7*2^n -A039834(n)-7*A039834(n+1))/5 ;
%p end proc:
%p seq(A120461(n),n=0..10) ; # _R. J. Mathar_, Nov 05 2011
%t M = {{0, 1, 0, 0, 1}, {1, 0, 1, 0, 0}, {0, 1, 0, 1, 0}, {0, 0, 1, 0, 1}, {1, 0, 0, 1, 0}} v[1] = {0, 1, 1, 2, 3} v[n_] := v[n] = M.v[n - 1] a = Table[Floor[v[n][[1]]], {n, 1, 50}]
%K nonn,easy
%O 0,2
%A _Roger L. Bagula_ and _Gary W. Adamson_, Jun 28 2006
|