A007761 (n+1) * a(n+1) - 2 (68*n^2+68*n+27) * a(n) + 6 * n * (772*n^2+35) * a(n-1) - 2 * (2*n-1)^2 * (68*n^2-68*n+27) * a(n-2) + (2*n-1)^2 * (n-1) * (2*n-3)^2 * a(n-3) = 0.

1, 54, 6381, 1176900, 295843509, 94263721650, 36391089828249, 16506884910849480, 8603605199199386025, 5066519768097762780270, 3326644994941284848273925, 2409605195467508091244871820

Offset

0,2

Comments

Is this always integral?

Links

G. C. Greubel, Table of n, a(n) for n = 0..250

Maple

a[0]:=1: a[1]:=54: a[2]:=6381: a[3]:=1176900: for n from 3 to 11 do a[n+1]:=(2*(68*n^2+68*n+27)*a[n]-6*n*(772*n^2+35)*a[n-1]+2*(2*n-1)^2* (68*n^2-68*n+27)*a[n-2]-(2*n-1)^2*(n-1)*(2*n-3)^2*a[n-3])/(n+1) od: seq(a[n], n=0..12); # Emeric Deutsch, Jul 20 2005

Mathematica

a[0]:=1; a[1]:=54; a[2]:=6381; a[3]:=1176900; a[n_]:= a[n]= (2*(68*(n-1)^2 + 68*(n-1) +27)*a[n-1] -6*(n-1)*(772*(n-1)^2 +35)*a[n-2] +2*(2*n-3)^2*(68*(n- 1)^2 -68*(n-1) +27)*a[n-3] -(2*n-3)^2*(n-2)*(2*n-5)^2*a[n-4])/n; Table[a[n], {n, 0, 12}] (* G. C. Greubel, Mar 04 2020 *)

Crossrefs

Sequence in context: A254620 A299862 A342893 * A178633 A369701 A299949

Adjacent sequences: A007758 A007759 A007760 * A007762 A007763 A007764

Keyword

nonn

Author

Bruno Haible

Extensions

More terms from Emeric Deutsch, Jul 20 2005

Status

approved