OFFSET
1,1
LINKS
Index entries for linear recurrences with constant coefficients, signature (2,-2,2,-1).
FORMULA
a(n) = (1/2)*(10*n - (1+2*i)*(-i)^n - (1-2*i)*i^n), where i = sqrt(-1).
a(n) = 5*n - 2*sin(Pi*n/2) - cos(Pi*n/2).
a(n) = 5*n - A228826(n-1). - Andrew Howroyd, Nov 29 2018
G.f.: x*(x^3 + x^2 + 5*x + 3) / ((x - 1)^2 *(x^2 + 1)). - Vincenzo Librandi, Dec 06 2018
MAPLE
seq(coeff(series(x*(x^3+x^2+5*x+3)/((1-x)^2*(1+x^2)), x, n+1), x, n), n = 1 .. 60); # Muniru A Asiru, Dec 06 2018
MATHEMATICA
CoefficientList[Series[(x^3 + x^2 + 5 x + 3)/((x - 1)^2 (x^2 + 1)), {x, 0, 50}], x] (* or *)
a[n_]:= (1/2) (10 n - (1 + 2 * I) (-I)^n - (1 - 2 I) I^n); Simplify[Array[a, 50]] (* Stefano Spezia, Nov 29 2018 *)
LinearRecurrence[{2, -2, 2, -1}, {3, 11, 17, 19}, 60] (* Vincenzo Librandi, Dec 06 2018 *)
PROG
(PARI) Vec((3 + 5*x + x^2 + x^3)/((1 - x)^2*(1 + x^2)) + O(x^60)) \\ Andrew Howroyd, Nov 29 2018
(Magma) I:=[3, 11, 17, 19]; [n le 4 select I[n] else 2*Self(n-1)-2*Self(n-2)+2*Self(n-3)-Self(n-4): n in [1..60]]; // Vincenzo Librandi, Dec 06 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Mark A. Thomas, Nov 29 2018
STATUS
approved