A120461 Expansion of x*(4-x)/( (2x-1)*(x^2-x-1)).

0, 4, 3, 15, 16, 55, 73, 206, 315, 787, 1320, 3051, 5437, 11950, 22159, 47135, 89712, 186799, 361665, 742638, 1454035, 2958619, 5835448, 11803235, 23392341, 47131150, 93701703, 188310471, 375153280, 752681287, 1501520185, 3009257486

Offset

0,2

Links

Table of n, a(n) for n=0..31.

Index entries for linear recurrences with constant coefficients, signature (1,3,-2)

Formula

a(n) = (7*2^n -A039834(n)-7*A039834(n+1))/5. - R. J. Mathar, Nov 05 2011

Maple

A039834 := proc(n)
        if n <= 1 then
                n;
        else
                -procname(n-1)+procname(n-2) ;
        end if;
end proc:
A120461 := proc(n)
        (7*2^n -A039834(n)-7*A039834(n+1))/5 ;
end proc:
seq(A120461(n), n=0..10) ; # R. J. Mathar, Nov 05 2011

Mathematica

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}]

Crossrefs

Sequence in context: A358281 A095332 A195586 * A324013 A345013 A113204

Adjacent sequences: A120458 A120459 A120460 * A120462 A120463 A120464

Keyword

nonn,easy

Author

Roger L. Bagula and Gary W. Adamson, Jun 28 2006

Status

approved