A033574 a(n) = (2*n+1)*(10*n+1).
1, 33, 105, 217, 369, 561, 793, 1065, 1377, 1729, 2121, 2553, 3025, 3537, 4089, 4681, 5313, 5985, 6697, 7449, 8241, 9073, 9945, 10857, 11809, 12801, 13833, 14905, 16017, 17169, 18361, 19593, 20865, 22177, 23529, 24921, 26353, 27825, 29337, 30889, 32481, 34113, 35785, 37497, 39249
Offset: 0
Links
- G. C. Greubel, Table of n, a(n) for n = 0..1000
- Index entries for linear recurrences with constant coefficients, signature (3,-3,1).
Programs
-
GAP
List([0..50], n-> (2*n+1)*(10*n+1)); # G. C. Greubel, Oct 12 2019
-
Magma
[(2*n+1)*(10*n+1): n in [0..50]]; // G. C. Greubel, Oct 12 2019
-
Maple
seq((2*n+1)*(10*n+1), n=0..50); # G. C. Greubel, Oct 12 2019
-
Mathematica
Table[(2*n+1)*(10*n+1), {n, 0, 50}] (* G. C. Greubel, Oct 12 2019 *) -
PARI
a(n)=(2*n+1)*(10*n+1) \\ Charles R Greathouse IV, Jun 17 2017
-
Sage
[(2*n+1)*(10*n+1) for n in range(50)] # G. C. Greubel, Oct 12 2019
Formula
From G. C. Greubel, Oct 12 2019: (Start)
G.f.: (1 + 30*x + 9*x^2)/(1-x)^3.
E.g.f.: (1 + 32*x + 20*x^2)*exp(x). (End)
Extensions
Terms a(33) onward added by G. C. Greubel, Oct 12 2019