|
|
A134474
|
|
a(n) = numerator of Sum_{k=1..n} 1/A134473(k).
|
|
5
|
|
|
|
OFFSET
|
1,2
|
|
COMMENTS
|
The denominator of Sum_{k=1..n} 1/A134473(k) is A134475(n). A134474(n)/A134475(n) approaches a constant (0.6037789...) as n approaches infinity.
|
|
LINKS
|
|
|
MAPLE
|
Digits := 220 ; A134473 := proc(n) option remember ; local su, mu ; if n =1 then 2; else su := add(1/procname(k), k=1..n-1) ; mu := mul(1/(1+1/procname(j)), j=1..n-1) ; ceil( (1+su+sqrt((su-1)^2+4*mu))/2/(mu-su) ) ; fi; end: A134474 := proc(n) add(1/A134473(k), k=1..n) ; numer(%) ; end: seq(A134474(n), n=1..9) ; # R. J. Mathar, Jul 20 2009
|
|
MATHEMATICA
|
b[n_] := b[n] = If[n == 1, 2, With[{x = Product[1/(1 + 1/b[j]), {j, 1, n-1}], y = Sum[1/b[j], {j, 1, n-1}]}, Ceiling[(1 + y + Sqrt[(y-1)^2 + 4 x])/(2 (x-y))]]];
a[n_] := Sum[1/b[k], {k, 1, n}] // Numerator;
|
|
CROSSREFS
|
|
|
KEYWORD
|
frac,nonn
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|