OFFSET
1,1
COMMENTS
Leading zeros are not permitted, so each term is 4 digits in length.
LINKS
Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Jon Borwein, 170,000 digits of Gamma [Wayback Machine copy]
Eric Weisstein's World of Mathematics, Euler-Mascheroni Constant.
MAPLE
Digits := 420 ;
for sh from 3 do
p := floor(gamma*10^sh) mod 10000 ;
if isprime(p) and p > 999 then
printf("%d, ", p);
end if;
end do: # R. J. Mathar, Oct 31 2011
MATHEMATICA
egp[len_]:=Module[{egterms=FromDigits/@Partition[RealDigits[EulerGamma, 10, 1000][[1]], len, 1]}, Select[egterms, IntegerLength[#]==len&&PrimeQ[#]&]]; egp[4] (* Harvey P. Dale, Oct 29 2011 *)
PROG
(PARI) L=10^4; for(i=3, 999, isprime(p=Euler\.1^i%L)&p*10>L&print1(p", ")) \\ M. F. Hasler, Oct 31 2011
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Andrew G. West (WestA(AT)wlu.edu), Mar 29 2005
EXTENSIONS
Offset changed from 0 to 1 by Vincenzo Librandi, Apr 21 2013
STATUS
approved