OFFSET
1,1
COMMENTS
A generalization of Somos quadratic recurrence constant A112302.
LINKS
G. C. Greubel, Table of n, a(n) for n = 1..10000
FORMULA
Product_{n>0} ((2*n + 1) / (2*n - 1)) ^ (2^-n). - Michael Somos, Feb 24 2014
EXAMPLE
2.10784090263785493...
MAPLE
with(numtheory);
P:=proc(i)
local a, n;
a:=1;
for n from i by -1 to 1 do a:=(2*n-1)*sqrt(a); od;
print(evalf(sqrt(a), 1000));
end:
P(5000);
MATHEMATICA
digits = 100; Clear[p]; p[m_] := p[m] = Fold[N[Sqrt[#2*#1], digits] &, 1, Range[2*m + 1, 1, -2]] // RealDigits[#, 10, digits] & // First; p[digits]; p[m = 2*digits]; While[p[m] != p[m/2], m = 2*m]; p[m] (* Jean-François Alcover, Feb 24 2014 *)
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Paolo P. Lava, Apr 12 2011
STATUS
approved