OFFSET
0,3
COMMENTS
LINKS
Paul D. Hanna, Table of n, a(n) for n=0..50
FORMULA
Equals the main diagonal of triangle A187005.
PROG
(PARI) /* As column 0 of triangle A135080 (slower): */
{a(n)=local(F=x, M, N, P); M=matrix(n+2, n+2, r, c, F=x; for(i=1, r+c-2, F=subst(F, x, x+x^2+x*O(x^(n+2)))); polcoeff(F, c)); N=matrix(n+1, n+1, r, c, M[r, c]); P=matrix(n+1, n+1, r, c, M[r+1, c]); (P~*N~^-1)[n+1, 1]}
(PARI) /* As the main diagonal of triangle A187005 (faster): */
{a(n)=local(Rn=y); for(m=1, n+1, Rn=subst(truncate(Rn), y, y+y^2+y*O(y^m))); polcoeff(Rn/y, n, y)}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Nov 18 2007
STATUS
approved