login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A107382 a(n) = 4*a(n-1)-4*a(n-3)-a(n-4). 0

%I

%S 4,4,16,52,188,684,2512,9244,34052,125476,462416,1704212,6280892,

%T 23148428,85314448,314430012,1158845444,4270975556,15740867728,

%U 58013659124,213811888828,788013108844,2904256931152,10703766510172

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

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

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

%t n = 4; b4 = x /. NSolve[x^4 - n*x^3 + n*x + 1 == 0, x][[4]] b3 = x /. NSolve[x^4 - n*x^3 + n*x + 1 == 0, x][[3]] b2 = x /. NSolve[x^4 - n*x^3 + n*x + 1 == 0, x][[2]] b1 = x /. NSolve[x^4 - n*x^3 + n*x + 1 == 0, x][[1]] digits = 25 a = Table[n*(b4^m + b3^m + b1^m + b2^m)/ (b4 + b3 + b2 + b1), {n, 0, digits}]

%K nonn,easy

%O 0,1

%A _Roger L. Bagula_, May 24 2005

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

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 21 06:58 EDT 2021. Contains 348145 sequences. (Running on oeis4.)