A322504 a(n) = -4*a(n-1) - 3*a(n-2) + a(n-3), a(0) = 1, a(1) = -2, a(2) = 4.
1, -2, 4, -9, 22, -57, 153, -419, 1160, -3230, 9021, -25234, 70643, -197849, 554233, -1552742, 4350420, -12189221, 34152882, -95693445, 268125913, -751270435, 2105010556, -5898105006, 16526117921, -46305146110, 129744125671, -363534946433, 1018602262609, -2854060085466, 7996898607604
Offset: 0
Links
- Colin Barker, Table of n, a(n) for n = 0..1000
- Index entries for linear recurrences with constant coefficients, signature (-4,-3,1).
Programs
-
Magma
I:=[1,-2,4]; [n le 3 select I[n] else -4*Self(n-1) - 3*Self(n-2) + Self(n-3): n in [1..31]]; // Vincenzo Librandi, Jan 13 2019
-
Mathematica
LinearRecurrence[{-4,-3,1},{1,-2,4},50] (* Stefano Spezia, Jan 11 2019 *) RecurrenceTable[{a[0]==1, a[1]==-2, a[2]==4, a[n]==-4 a[n-1]-3 a[n-2]+a[n-3]}, a, {n, 30}] (* Vincenzo Librandi, Jan 13 2019 *)
-
PARI
Vec((1 + 2*x - x^2) / (1 + 4*x + 3*x^2 - x^3) + O(x^30)) \\ Colin Barker, Jan 11 2019
Formula
G.f.: (1 + 2*x - x^2) / (1 + 4*x + 3*x^2 - x^3). - Colin Barker, Jan 11 2019
Comments