# Prime number counting function with error O(x^e)

eljose
i have discovered a formula for $$\pi(x^{a})$$ being Pi the prime number counting function in terms of a triple integral...but you wll say..this is already made what is this good for?..in fact if you knew $$\pi(x^{a})$$ with a total error O(x^d) by setting a=Ad and making A--->oo (infinite) the total error would be e=1/A O(x^e) with e the smallest positive number,or expanding in e powers to lineal order x^{e}=1+eln(x), so i have discovered a formula with the smallest possible error term...and this is completely new...:)

Staff Emeritus
Gold Member
There is no "smallest positive number", and the smallest possible error term is O(0).

eljose
then i made a mistake in it sorry...i mean we can always choose and a for $$\pi(x^{a})$$ so the error term goes like this O(x^e) with e can be chosen to be the smallest number it ocurs to us (for example e=10^{-100000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000} i know i,m exagerating but this is the sense of my proof,i have reduced the error term to O(x^e) with choosing a=dA and making A tends to infinite the error term would go like this O(x^e) with e an infinitesimal number ( i know this exist i have seen proofs with results a+e being e an infinitesimal number.

-The formula for $$\pi(x^a)$$ is a triple integral..how would we calculate it?..the solution is to use Gauss-Hermite Quadrature formula so the integral becomes:

$$\sum_{i,j,k}C_{i,j,k}W(x_i,x_j,x_k)F(x_i,x_j,x_k)$$

where the C.s are constants, W is a weight function in the form W(x,y,z)=exp(-x^{2}-y^{2}-z^{2}) and F is the integrand of the our triple integral...in fact is a complex integral so we first have to make the change of variable s=c+iu , q=d+iv (the i,j and k indexes are summed over the roots of Hermite Polynomials)...

-Of course the smallest error term would be O(0) but with my method the error goes like O(x^e) for very small e it is almost like 1+eln(x),that is if we ignore the lineal term,the error would be constant for every x, including the error lineal term the error would go like ln(x) and as far as i know the smallest error term for prime counting function is O(x^{2/3}ln(x)) mine is smaller...

Last edited:
Staff Emeritus
Gold Member
I'm suspicious, because my attempts to evaluate your previous integrals with Mathematica failed, returning with some astronomically huge value, which suggests that the integrals, in fact, do not converge.

But those worries aside, what is the error term for this numerical integration method? Can you find out how many terms you must add in order to get a sufficiently good upper bound on the error term? How large do the numbers in the summands become?

(And it would be nice to know how many bits of precision must be used to carry out the arithmetic)

eljose
i can not calculate the error it will be in general O(x^d) with d=a+b+c
a=number of time used to evaluate the integral.
b=time used to write and evaluate the $$\zeta(s), ln\zeta(s)$$
c=time to write $$\pi(x^a)$$
the error term in doing the integration is calculated as the error term in Gauss Hermite quadrature formula,it will be in general the product of the errors in calculating each integral...

Homework Helper
eljose said:
..in fact if you knew $$\pi(x^{a})$$ with a total error O(x^d) by setting a=Ad and making A--->oo (infinite) the total error would be e=1/A O(x^e) with e the smallest positive number...

This looks like nonsense to me. If you can find $$\pi(x^{a})$$ with a total error O(x^d) then your a and d are fixed. As this A changes, your error analyisis will no longer be correct. Before you tried to do some kind of change of variables to make the error look smaller, but actually did nothing at all, is this what you're doing again?

then there's the whole "smallest positive number" problem...

eljose said:
( i know this exist i have seen proofs with results a+e being e an infinitesimal number.

No you haven't, not in any number theory paper at least. You'll often see an error like $$O(x^{1+\epsilon})$$ and they say you can take any $$\epsilon>0$$. They are not saying that $$\epsilon$$ is an "infintessimal number", just that the big O bound of that form holds for any possible fixed $$\epsilon>0$$ that you like, though the big O constant may depend on epsilon (sometimes they write $$O_\epsilon (x^{1+\epsilon})$$ to make this dependence explicit, but not usually).

eljose said:
i can not calculate the error it will be in general O(x^d) with d=a+b+c
a=number of time used to evaluate the integral.

So you're saying the more time you spend on calculating the integral (presumably to get a more accurate value for it), the larger your error term? That sounds bad...

Homework Helper
eljose said:
i can not calculate the error it will be in general O(x^d) with d=a+b+c
a=number of time used to evaluate the integral.
b=time used to write and evaluate the $$\zeta(s), ln\zeta(s)$$
c=time to write $$\pi(x^a)$$
the error term in doing the integration is calculated as the error term in Gauss Hermite quadrature formula,it will be in general the product of the errors in calculating each integral...

? The time you take to do calculate something has nothing to do with "error" in calculating. Do you mean "how many times" you do something rather than "time used" or "time to"?

eljose
But you have the problem of calculating a complex triple integral...in general the error of calculating a triple integral will be of the form O(x^d) with d a constant fixed for every integral i mean..to calculate a triple integral in the form:

$$\int_{\Re}dxdydzF(t,x,y,z)$$ the error term will be like O(t^d) will depend only on the value t inside it...

-In the case a and d were related we would have a=f(d) so we could choose a acurately so d is the smallest possible quantity so the error will be the smallest possible number it ocurs to us.

Homework Helper
eljose said:
-In the case a and d were related we would have a=f(d) so we could choose a acurately so d is the smallest possible quantity so the error will be the smallest possible number it ocurs to us.

That's swell, but then a has changed as well, and your end result is that you've done nothing.

Let's suppose you know $$f(x)=O(x^{3})$$ (for some function f). You then know that $$f(x^{1/3})=O(x)$$. This doesn't mean you've found a smaller bound for your function f.

eljose
the qeustion is that i have given a formula for $$\pi(x^a)$$ with error O(x^d) so let,s suppose we want to calculate $$\pi(10^{1000})$$ then set a=1000 and x=10 so the total error would be O(10^d) that is the advantage of my method... if you want to judge the formula by yourselves it is:

$$\pi(x^a)=\int_{c-i\infty}^{c+i\infty}\int_0^{\infty}\int_{d-i\infty}^{d+i\infty}dqdnds\frac{x^{q}n^{s}}{qs\zeta(s)}D_{n}(n^{-1}Ln\zeta(bnq)$$

where b=1/a this is the triple integral we must calculate if we make the
change of variable c+iu=q d+iv=s the limits of integration turn real...

Last edited:
Homework Helper
eljose said:
the qeustion is that i have given a formula for $$\pi(x^a)$$ with error O(x^d) so let,s suppose we want to calculate $$\pi(10^{1000})$$ then set a=1000 and x=10 so the total error would be O(10^d) that is the advantage of my method... if you want to judge the formula by yourselves it is:

$$\pi(x^a)=\int_{c-i\infty}^{c+i\infty}\int_0^{\infty}\int_{d-i\infty}^{d+i\infty}dqdnds\frac{x^{q}n^{s}}{qs\zeta(s)}D_{n}(Ln\zeta(bnq))/n$$

where b=1/a this is the triple integral we must calculate if we make the
change of variable c+iu=q d+iv=s the limits of integration turn real...

:rofl:

"d=a+b+c "

You say a = 1000 here, now assuming b and c are both 0 or greater that means your error term is:

O(10^1000) or more and as your calculation is of Pi(10^1000) that means your error term is even greater than the value of Pi(x) your supposed to be calculating :rofl:

Or assuming you mean runtime, really hard to tell from your posts, such huge amounts of run time are totally useless.

Last edited:
eljose
No zurtex sorry it was my fault..when is put d=a+b+c a is not the power of the $$\pi(x^a)$$ d is the error made in calculating the Riemann Zeta function and the triple integral...

Homework Helper
eljose said:
the qeustion is that i have given a formula for $$\pi(x^a)$$ with error O(x^d) so let,s suppose we want to calculate $$\pi(10^{1000})$$ then set a=1000 and x=10 so the total error would be O(10^d) that is the advantage of my method... if you want to judge the formula by yourselves it is:

So you're saying that you are free to choose a to be whatever you like, apparently not affecting d or the error term at all? This is as ridiculous as having an error term in the first place for what you're also claiming is an exact formula.

Again you're missing the distinction between runtime and an error term. If you have some kind of exact algorithm, then there's no error term. What's interesting is how long it takes the algorithm to run (usually measured in number of some kind of operation, like the number of multiplications).

Maybe you're claiming that for a fixed run time you get a certain error. If this is the case, then there is no way that your error is independant of the size of the argument, which is what you'd be claiming.

Why do I get the feeling that this is going to be a complete repetition of your earlier threads on this?

eljose
I have the "exact" formula for Pi(x^a) in the form:

The following code was used to generate this LaTeX image:

$$\pi(x^a)=\int_{c-i\infty}^{c+i\infty}\int_0^{\infty}\int_{d-i\infty}^{d+i\infty}dqdnds\frac{x^{q}n^{s}}{qs\zeta(s)}D_{n}(n^{-1}Ln\zeta(bnq))$$

but unfortunately this integral can not be made exactly and also zeta can not be calculated exactly (only an approximation to it) so in general there will be a running time in the form O(x^d) (the bigger x is the most difficult is to calculate it) but we have calculated Pi(x^a) instead of Pi(x) so the true error term will be O(x^{d/a}) the running time needed to calculate the integral will depend on x...

the main advantage of my method is that you calculate Pi(x^a) so you need to calculate a low x a put a big a for example a triple integral would require a running time O(x^d) if you want to calculate Pi(10^1000) put x=10 and a=1000 so the error would be O(10^d) i don,t understand why math referees can not understand this method is better than other analytic method to calculate Pi(x) known..¡¡i really don,t like mathematicians and their nasty rigour¡¡   (Hurkyl don,t get angry or ban me,but sometimes with they saying rigour and rigour they get on my nerves,in fact this method is new as it calculates Pi(x^a) instead of Pi(x).

Last edited:
Homework Helper
Do you even know what you're claiming anymore? I sure don't, in one breath O(x^d) is the runtime, the next an error term. Fill in the blanks please:

Using this triple integral, you can calculate pi(x) with error _____ in _____ time. No unspecified constants please, i.e. no "d"'s allowed, no exponent a, this is a question about pi(x).

I'll say again, there is absolutely no way that your error term and run time are independant of the argument of the prime counting function. What will it take to get this nonsense kicked out of this forum?

eljose
-as you can see the power a of $$\pi(x^a)$$ appears as b=1/a in the form:
$$ln\zeta(bnq)$$ so the bigger b is the higher the error in calculating the integral will be but b=1/a so big b means small a,so the smallest error will get for a big (the bigger the a is the better,and the smaller error) so choosing a--->oo we will get a very small error

Homework Helper
This is the same integral you had before right? But with a missing constant and your old $$e^{st}$$ term, and a mysterious subscript n in the $$D_n$$ notation? Where did this b in the log come from? Some sort of change of variables that you aren't going to bother mentioning?

You're going to have something of the form $$x^{s}$$ in your integrand if you're trying to find pi(x), it comes from Perron's (this was your $$e^{st}$$ term, for some reason you wanted to use x=e^t). Maybe you're trying to hide this in your x^q term, but that doesn't seem dependant on a... No change of variables will let you escape the fact that as the argument of pi grows, so does the absolute value of this x^s term, either increasing your error or the time it takes to calculate the integral to a given precision.

A simple reality check is in order. Do you really believe that you are capable of computing pi(x) to within a constant error in a constant amount of time (that is neither depend on x)? Honestly? Actually you seem to be making an even more ludicrous claim that as x grows (in your notation-a is growing but x is remaining fixed), the error is actually decreasing.

Staff Emeritus
Gold Member
Why, oh why, aren't these threads in the IR section ?

eljose
the b comes from the fact that:

$$s\int_0^{\infty}dt\pi(e^{at})e^{-st}=sb\int_0^{\infty}du\pi(e^u)e^{-sbt}=sb\int_1^{\infty}dr\pi(r)r^{-sb-1}=bR(bs)$$

where R(bs) is equal to $$\sum_p{p^{-sb}}^$$ and the sum is extended over all primes,and the R(sb) function is related to $$ln\zeta(bns)$$ where b=1/a.

Last edited:
Homework Helper
eljose said:
I have the "exact" formula for Pi(x^a)[...] in fact this method is new as it calculates Pi(x^a) instead of Pi(x).

Please calculate the following with your $$\pi(x^a)$$ formula:
1. $$\pi(10^{22})$$ ($$x=10^{22},a=1$$)
2. $$\pi((10^{11})^2)$$ ($$x=10^{11},a=2$$)
3. $$\pi((10^2)^{11})$$ ($$x=10^2,a=11$$)
4. $$\pi(10^{22})$$ ($$x=10,a=22$$)

Thank you!

Homework Helper
This is hopeless. All you've done is a change of variables. This isn't going to reduce the error in your integral. If you are integrating over some range of t then set t=u/a, as a increases the range you're integrating u over changes as well. Nothing is saved. I thought I've explained this to you before when you were bent on evaluating pi(e^t) instead of pi(x)

eljose
but i have given an expression for $$\pi(x^a)$$ in the form:

$$\int_0^{\infty}dx\pi(x^a)e^{-st}=bR(sb)/s$$ it yields to a triple integral for the Pi(x^a) depending on b the bigger b is the bigger you get the error if we set O(x^d) then d=f(b) wiht f a increasing function so for bigger a (small b) the error is smaller...

Homework Helper
eljose, to help us better understand the accuracy of your formula and to help us see what you mean by "this method is new as it calculates Pi(x^a) instead of Pi(x)", would you calculate Pi(x) for the 4 values I listed above (post #20)?

eljose
i would do it ..if i had a computer program to evaluate integrals numerically or if i knew how to calculate it...:) :) :)

Staff Emeritus
Gold Member
$$\int_0^{\infty}dx\pi(x^a)e^{-st}=bR(sb)/s$$

It should be clear that this is not a convergent integral...

Homework Helper
eljose said:
i would do it ..if i had a computer program to evaluate integrals numerically or if i knew how to calculate it...:) :) :)
You don't know

a) How to write a program to do the sort of things you claim you can do

b) How to calculate your own integral numerically

And you expect people to take you seriously?

but i have given an expression for $$\pi(x^a)$$ in the form: