OFFSET
1,3
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..250
MATHEMATICA
m=42;
M= Table[If[k>n-1, 0, If[k==n-1, n, -1]], {n, 0, m+1}, {k, 0, m+1}];
T:= T= Sum[MatrixPower[M, j]/j, {j, m+1}];
Table[T[[n+1, 2]]/2, {n, 2, 30}] (* G. C. Greubel, Sep 06 2022 *)
PROG
(PARI) {a(n)=local(M=matrix(n+4, n+4, r, c, if(r>=c, if(r==c+1, -c, 1))), L=sum(m=1, n+4, (M^0-M)^m/m)); L[n+2, 2]/2}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Mar 11 2006
STATUS
approved