## 24685Re: Twin primes and Pi

Expand Messages
• Nov 26, 2012
• 0 Attachment
"Alexander" <werner.sand@...> wrote:

> > [3.196178168631, 2097257]
> > [3.196178168641, 4194581]
> > [3.196178168628, 8388617]
> > [3.196178168631, 16777289]
> > [3.196178168630, 33554501]
> > [3.196178168629, 67109321]
> > [3.196178168630, 134217779]
>
> I would appreciate if you could also use a
> simple program without Brun's constant.

OK, but now it is of course a dumb method, since
the HL correction for truncation is much bigger,
when you wilfully throw away all of that hard
work by Pascal Sebah to evaluate B:

{Sdumb(t)=log(t/(t-1));}

{Rdumb(t)=2*0.66016181584686957/log(t);}

{T=134217779;default(primelimit,T);
default(realprecision,7);e=2^21;s=0;
forprime(t=3,T,if(isprime(t+2),s+=Sdumb(t);
if(t>e,print([exp(s+Rdumb(t)),t]);e*=2)));}

[3.196175, 2097257]
[3.196230, 4194581]
[3.196199, 8388617]
[3.196216, 16777289]
[3.196203, 33554501]
[3.196185, 67109321]
[3.196191, 134217779]

David
• Show all 24 messages in this topic