login
A052459
a(n) = n*(2*n^2 + 1)*(n^2 + 1)/6.
4
1, 15, 95, 374, 1105, 2701, 5775, 11180, 20049, 33835, 54351, 83810, 124865, 180649, 254815, 351576, 475745, 632775, 828799, 1070670, 1366001, 1723205, 2151535, 2661124, 3263025, 3969251, 4792815, 5747770, 6849249, 8113505
OFFSET
1,2
COMMENTS
2-magic series constant.
LINKS
Eric Weisstein's World of Mathematics, Magic constant.
FORMULA
G.f.: x*(1 + 3*x + x^2)*(1 + 6*x + x^2)/(1-x)^6. - Chai Wah Wu, Dec 17 2016
E.g.f.: x*(6 +39*x +53*x^2 +20*x^3 +2*x^4)*exp(x)/6. - G. C. Greubel, Sep 23 2019
MAPLE
seq(n*(2*n^2 +1)*(n^2 +1)/6, n=1..30); # G. C. Greubel, Sep 23 2019
MATHEMATICA
Table[n*(2*n^2 +1)*(n^2 +1)/6, {n, 30}] (* G. C. Greubel, Sep 23 2019 *)
PROG
(PARI) vector(30, n, n*(2*n^2 +1)*(n^2 +1)/6) \\ G. C. Greubel, Sep 23 2019
(Magma) [n*(2*n^2 +1)*(n^2 +1)/6: n in [1..30]]; // G. C. Greubel, Sep 23 2019
(Sage) [n*(2*n^2 +1)*(n^2 +1)/6 for n in (1..30)] # G. C. Greubel, Sep 23 2019
(GAP) List([1..30], n-> n*(2*n^2 +1)*(n^2 +1)/6); # G. C. Greubel, Sep 23 2019
CROSSREFS
Sequence in context: A189657 A190274 A374808 * A336278 A044266 A044647
KEYWORD
nonn,easy
EXTENSIONS
Formula from Vladeta Jovovic, Jun 15 2002
STATUS
approved