|
PROGRAM
|
(PARI) f(n) = for(x=0, n, print1(twinl(x!), ", ")) twinl(n) = { local(c, x); c=0; x=1; while(c<n, if(isprime(prime(x)+2), c++); x++; ); return(prime(x-1)) }
(PARI) default(primelimit, 7000000); twinl(n) = { local(c, x); c=0; x=1; while(c<n, if(isprime(prime(x)+2), c++); x++); return(prime(x-1)) } { for(n=0, 8, write("b093310.txt", n, " ", twinl(n!))); } [From Harry J. Smith (hjsmithh(AT)sbcglobal.net), Jun 20 2009]
|