|
EXAMPLE
|
E.g.f.: A(x) = x^2/2! + 3*x^3/3! + 15*x^4/4! + 90*x^5/5! + 629*x^6/6! + 5019*x^7/7! + 45030*x^8/8! + 448776*x^9/9! + 4919321*x^10/10! + ...
[1],
[1, 0, 1],
[(1), 1, 2, 1, 1],
[1,(3), 5, 6, 5, 3, 1],
[1, 6, (15), 24, 28, 24, 15, 6, 1],
[1, 10, 40, (90), 139, 160, 139, 90, 40, 10, 1],
[1, 15, 91, 300, (629), 945, 1078, 945, 629, 300, 91, 15, 1], ...
where coefficients in parenthesis form the initial terms of this sequence.
|
|
PROG
|
(PARI) {a(n) = polcoeff( prod(j=0, n-1, 1 + j*x + x^2), n-2)}
for(n=2, 30, print1(a(n), ", "))
|