n := 12; p:=2; PRIME 2f := x-> (1-x)*cos(Pi*x) + (1/Pi)*sin(Pi*x);bnd := b->n*(gamma+log(8*Pi))-(n/2)*int((1-f(x/b))/sinh(x/2),x=0..b)
-(n/2)*int(csch(x/2),x=b..infinity)-4*int(f(x/b)*cosh(x/2),x=0..b);lcl := b->2*sum(log(p)*f(m*log(p)/b)*p^(-m/2),m=1..6); num := 1; NUMBER OF PRIMES IN Sbb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p)); CHECK UPPER LIMIT ON SUM IN LCLevalf(exp(bnd(bb)+num*lcl(bb)));num := 2;bb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));num :=3;bb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));n := 12; p:=3; PRIME 3f := x-> (1-x)*cos(Pi*x) + (1/Pi)*sin(Pi*x);lcl := b->2*sum(log(p)*f(m*log(p)/b)*p^(-m/2),m=1..3); num := 1; NUMBER OF PRIMES IN Sbb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));num := 2;lcl := b->2*sum(log(p)*f(m*log(p)/b)*p^(-m/2),m=1..4); bb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));num :=3;bb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));n := 12; p:=5; PRIME 5f := x-> (1-x)*cos(Pi*x) + (1/Pi)*sin(Pi*x);lcl := b->2*sum(log(p)*f(m*log(p)/b)*p^(-m/2),m=1..2); num := 1; NUMBER OF PRIMES IN Sbb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));num := 2;bb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));num :=3;bb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));n := 12; p:=7; PRIME 7f := x-> (1-x)*cos(Pi*x) + (1/Pi)*sin(Pi*x);lcl := b->2*sum(log(p)*f(m*log(p)/b)*p^(-m/2),m=1..2); num := 1; NUMBER OF PRIMES IN Sbb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));num := 2;bb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));num :=3;bb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(bb/log(p));evalf(exp(bnd(bb)+num*lcl(bb)));num := 0; BOUND WITHOUT LOCAL CORRECTIONSbb := fsolve(diff(bnd(b)+num*lcl(b),b)=0,b,0..8);evalf(exp(bnd(bb)+num*lcl(bb)));