OFFSET
0,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..401
FORMULA
a(n) = ((24*n^3 + 60*n^2 + 50*n + 5)*a(n-1) + (12*n^2 + 24*n + 1)*a(n-2))/(12*n^2 - 11), with a(0) = 1, a(1) = 65. - G. C. Greubel, Oct 08 2023
MATHEMATICA
a[n_]:= a[n]= If[n<2, (65)^n, ((24*n^3+60*n^2+50*n+5)*a[n-1] +(12*n^2 + 24*n+1)*a[n-2])/(12*n^2-11)];
Table[a[n], {n, 0, 40}] (* G. C. Greubel, Oct 08 2023 *)
PROG
(Magma) [n le 2 select (65)^(n-1) else ((24*n^3-12*n^2+2*n-9)*Self(n-1) + (12*n^2-11)*Self(n-2))/(12*(n-1)^2 -11): n in [1..40]]; // G. C. Greubel, Oct 08 2023
(SageMath)
@CachedFunction
def a(n): # a = A144500
if (n<2): return (65)^n
else: return ((24*n^3 + 60*n^2 + 50*n + 5)*a(n-1) + (12*n^2 + 24*n + 1)*a(n-2))/(12*n^2 - 11)
[a(n) for n in range(41)] # G. C. Greubel, Oct 08 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
David Applegate and N. J. A. Sloane, Dec 13 2008
STATUS
approved