|
|
A033595
|
|
a(n) = (n^2-1)*(2*n^2-1).
|
|
1
|
|
|
1, 0, 21, 136, 465, 1176, 2485, 4656, 8001, 12880, 19701, 28920, 41041, 56616, 76245, 100576, 130305, 166176, 208981, 259560, 318801, 387640, 467061, 558096, 661825, 779376, 911925, 1060696, 1226961
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
0,3
|
|
LINKS
|
|
|
FORMULA
|
G.f.: (1 -5*x +31*x^2 +21*x^3)/(1-x)^5. - R. J. Mathar, Feb 06 2017
E.g.f.: (1 - x + 11*x^2 + 12*x^3 + 2*x^4)*exp(x). - G. C. Greubel, Mar 05 2020
Sum_{n>=2} 1/a(n) = (Pi/sqrt(2))*cot(Pi/sqrt(2)) + 7/4.
Sum_{n>=2} (-1)^n/a(n) = (Pi/sqrt(2))*cosec(Pi/sqrt(2)) - 11/4. (End)
|
|
MAPLE
|
|
|
MATHEMATICA
|
Table[(n^2 -1)*(2*n^2 -1), {n, 0, 40}] (* G. C. Greubel, Mar 05 2020 *)
LinearRecurrence[{5, -10, 10, -5, 1}, {1, 0, 21, 136, 465}, 40] (* Harvey P. Dale, Aug 24 2020 *)
|
|
PROG
|
(PARI) vector(41, n, my(m=n-1); (m^2 -1)*(2*m^2 -1)) \\ G. C. Greubel, Mar 05 2020
(Magma) [(n^2 -1)*(2*n^2 -1): n in [0..40]]; // G. C. Greubel, Mar 05 2020
(Sage) [(n^2 -1)*(2*n^2 -1) for n in (0..40)] # G. C. Greubel, Mar 05 2020
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,easy
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|