OFFSET
1,3
FORMULA
a(n) = floor(n! * n * (n+1)! * b(n-1))) where b(0) = gamma and b(n) = b(n-1) - a(n) / (n! * n *(n+1)!). - Sean A. Irvine, May 01 2018
PROG
(Macsyma) Ouang(x, n) := block([ ans : [ ] ], for k through n do (push(floor(multthru(k!*k*(k+1)!, x)), ans), x: x-ans[ 1 ]/(k!*k*(k+1)!)), reverse(ans)); Ouang(%gamma, 25);
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
More terms from Sean A. Irvine, May 01 2018
STATUS
approved