OFFSET
0,2
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..402
FORMULA
(6*n-5)*a(n) = (12*n^2 +8*n +3)*a(n-1) + (6*n+1)*a(n-2). - G. C. Greubel, Oct 07 2023
MATHEMATICA
a[n_]:= a[n]= If[n<2, 16^n, ((12*n^2+8*n+3)*a[n-1] +(6*n+1)*a[n-2])/(6*n-5)];
Table[a[n], {n, 0, 30}] (* G. C. Greubel, Oct 07 2023 *)
PROG
(Magma) [n le 2 select (16)^(n-1) else ((12*(n-1)^2+8*(n-1)+3)*Self(n-1) + (6*(n-1)+1)*Self(n-2) )/(6*(n-1) -5): n in [1..30]]; // G. C. Greubel, Oct 07 2023
(SageMath)
@CachedFunction
def a(n): # a = A144499
if n<2: return (16)^n
else: return ((12*n^2+8*n+3)*a(n-1) +(6*n+1)*a(n-2))/(6*n-5)
[a(n) for n in range(31)] # G. C. Greubel, Oct 07 2023
CROSSREFS
KEYWORD
nonn
AUTHOR
David Applegate and N. J. A. Sloane, Dec 13 2008
STATUS
approved