Es bien sabido que $li(x)$, la integral logaritmo es una muy buena aproximación de la $\pi(x)$, el nunmber de números primos que no exceda $x$.
Así, una muy buena aproximación de la probabilidad de que un número aleatorio con $x$ dígitos es un número primo, está dada por
$$f(x)=\frac{li(10^x)-li(10^{x-1})}{10^x-10^{x-1}}$$
Me di cuenta de que la función
$$g(x)=\frac{2000x}{(2000x^2-646x-141)\ln(10)}$$
es una excelente aproximación a $f(x)$ , utilizando una interpolación de la herramienta en un muy buen matemáticas-sitio de Arndt Brünner.
Definamos $h(x):=f(x)-g(x)$
Parece que para $x\ge 11$, la función de $h(x)$ mínimo sobre $-4.334\cdot 10^{-8}\ $ y máxima alrededor de $2.1749\cdot 10^{-8}\ $.
Son estas extremas global extremas ? Podemos demostrar rigurosamente que $|h(x)|<4.335\cdot 10^{-8}\ $ todos los $x\ge 11$ ?
Traté de verificar esto con Wolfram Alpha, pero la función parece ser demasiado complicado. Alguna idea ? El límite de $h(x)$, $x$ tiende a $\infty$, debería ser $0$.