OFFSET
0,2
PROG
(PARI) {a(n)=local(L=matrix(n+3, n+3, r, c, if(r>=c, binomial(2*r-2, 2*c-1)/2)), A);
A=sum(m=0, n, L^m/m!); A[n+3, 3]}
for(n=0, 25, print1(a(n), ", "))
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Aug 24 2014
STATUS
approved