login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

A061322
a(n) = a(n-1) * (1 + a(n-1)/n^2) with a(0) = 2.
1
2, 6, 15, 40, 140, 924, 24640, 12415040, 2408343949440, 71606426901226335015040, 51274803735606705472274088614112357905277056, 21728144612603201307908899563300049012978385050783094682272184269369267136230071558272
OFFSET
0,1
COMMENTS
Only first 42 terms are integers (see A003504).
LINKS
Eric Weisstein's World of Mathematics, Goebel's Sequence.
FORMULA
a(n) = a(n-1) + A003504(n+1)^2, a(n-1) = n * A003504(n+1). a(n) = A061314(2, n).
EXAMPLE
a(2) = 6 * (1 + 6/2^2) = 15.
MATHEMATICA
Block[{n = 0}, NestList[#*(1 + #/++n^2) &, 2, 11]] (* Paolo Xausa, Apr 17 2024 *)
PROG
(PARI) {a(n) = local(x); if( n<1, 2 * (n==0), (x = a(n-1)) + (x/n)^2)} /* Michael Somos, Apr 02 2006 */
CROSSREFS
Sequence in context: A206000 A316981 A319560 * A004664 A074446 A303551
KEYWORD
nonn
AUTHOR
Henry Bottomley, Apr 24 2001
STATUS
approved