login
Denominators of an Egyptian fraction for Euler-Mascheroni constant, using only prime numbers and allowing repetitions.
2

%I #3 Aug 24 2012 10:50:01

%S 2,13,3433,770741,9157546933,68379695298643311959,

%T 177698848753108581519320740556520784033,

%U 191429752162010655176200572455622919752517437701651231352646653127630507601

%N Denominators of an Egyptian fraction for Euler-Mascheroni constant, using only prime numbers and allowing repetitions.

%C Good approximation up to 801 decimal digits

%p P:=proc(n) local a,i; a:=evalf(gamma,100); for i from 1 by 1 to n do if 1/ithprime(i)<a then a:=a-1/ithprime(i); print(a); print(ithprime(i)); fi; od; end: P(100000);

%Y Cf. A139514-A139518, A139520-A139523.

%K easy,nonn,frac

%O 0,1

%A _Paolo P. Lava_ and _Giorgio Balzarotti_, Apr 28 2008