OFFSET
2,1
LINKS
Robert Israel, Table of n, a(n) for n = 2..1150
MAPLE
g:= gfun:-rectoproc({(-n^2-4*n-4)*s(n+1)+(3*n^2+15*n+19)*s(n+2)+(-3*n^2-18*n-27)*s(n+3)+(n^2+7*n+12)*s(n+4), s(0) = 0, s(1) = 0, s(2) = 1/2, s(3) = 1}, s(n), remember):
f:= n -> g(n)*ilcm($1..n)^2:
map(f, [$2..100]); # Robert Israel, Aug 10 2018
MATHEMATICA
nmax = 100;
T = RecurrenceTable[{(-n^2 - 4*n - 4)*s[n+1] + (3*n^2 + 15*n + 19)*s[n+2] + (-3*n^2 - 18*n - 27)*s[n+3] + (n^2 + 7*n + 12)*s[n+4] == 0, s[0] == 0, s[1] == 0, s[2] == 1/2}, s, {n, 0, nmax}];
f[n_] := T[[n+1]]*(LCM @@ Range[n])^2;
Map[f, Range[2, nmax]] (* Jean-François Alcover, Sep 19 2022, after Robert Israel *)
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Offset corrected and more terms from Robert Israel, Aug 10 2018
STATUS
approved