cp's OEIS Frontend

The sequence (A228113) yields an average relative difference in absolute value, i.e. Average(Abs(A228114(n))/ (A006879(n)) = 1.03936…x10^-2 for 1<=n<=25.

Note that A057793(n) = Riemann(10^n) is not defined for n=0. Its value is set to 0.

Coefficients are A=16103485019141/2900449771918928, B=5262046568827901/29305205016290, C(0)=296261685121849/265642652464758, C(1)=38398556529727/750568780742436, C(2)=0, C(3)=-11594434149768/8254020049890781.

This sequence gives a good approximation of the number of primes with n digits (A006879); see (A228112).

As the squares of odd-indexed Fibonacci numbers F(2n+1)(1) (see A227693) are equal or close to the first values of pi(10^n) (A006880), and as F(4n)(1)=(F(2n+1)(1))^2- (F(2n-1)(1))^2, it is legitimate to ask whether the first values of the differences pi(10^n)- pi(10^(n-1)) (A006879) are also close or equal to multiple of 4 index Fibonacci numbers F(4n)(1); e.g., for n=2, F(8)(1)=21. However, when using Fibonacci polynomials, the exact relation is xF(4n)(x)=(F(2n+1)(x))^2- (F(2n-1)(x))^2.

To obtain this sequence, one switches to the product of x and multiple of 4 index Fibonacci polynomials F(4n)(x), and one obtains the sequence a(n) by computing x as a function of n such that (xF(4n)(x)) fit the values of pi(10^n)- pi(10^(n-1)) for 1 <= n <= 25, with pi(1)=0.

A228111 provides exact values of pi(10^n) - pi(10^(n-1)) for n = 1 to 3 and yields an average relative difference in absolute value, i.e. average(abs(A228112(n))/A006879(n) = 0.00375341... for 1 <= n <= 25, better than when using the 10^n/log(10^n) function, which yields 0.0469094... (see A228066) or the logarithmic integral (Li(10^n) - Li(2)) function, which yields 0.0175492... (see A228068) or the Riemann (Riemann(10^n)) function, which yields 0.0103936... (see A228114) or the Fibonacci polynomials of multiple of 4 indices, which yields 0.00473860... (see A228064) for 1 <= n <= 25.