OFFSET
1,2
MAPLE
a := proc(n) option remember; if n < 3 then return [0, 9][n] fi;
(n + 1)^2*a(n - 2) + a(n - 1) end: print(seq(a(n), n=1..21));
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, May 20 2020
STATUS
approved