OFFSET
1,3
COMMENTS
LINKS
Robert Israel, Table of n, a(n) for n = 1..1663
MAPLE
R:= gfun:-rectoproc({r(n) = (r(n - 1))/(n - 1) + r(n - 2), r(1) = 0, r(2) = 1}, r(n), remember):
seq(numer(R(n)), n=2..30); # Robert Israel, May 14 2017
MATHEMATICA
Numerator[RecurrenceTable[{r[n] == (r[n - 1])/(n - 1) + r[n - 2], r[1] == 0, r[2] == 1}, r, {n, 2, 30}]] (* Terry D. Grant, May 07 2017, fixed by Vaclav Kotesovec, Aug 14 2021 *)
Table[Numerator[ContinuedFractionK[1, 1/k , {k, 1, n}]], {n, 1, 30}] (* Vaclav Kotesovec, Aug 14 2021 *)
PROG
(PARI) a(n)=local(A); CF=contfracpnqn(vector(n, k, 1/k)); A=denominator(CF[1, 1]/CF[2, 1])
CROSSREFS
KEYWORD
frac,nonn
AUTHOR
Paul D. Hanna, Feb 19 2004
STATUS
approved