OFFSET
0,3
PROG
(PARI) {a(n)=local(P=matrix(n+1, n+1, r, c, if(r>=c, (r^3)^(r-c)/(r-c)!)), D=matrix(n+1, n+1, r, c, if(r==c, r))); (P^-1*D*P)[n+1, 1]/(n+1)^3}
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul D. Hanna, Jun 07 2005
STATUS
approved