Como un apéndice de la mjqxxxx excelente respuesta, presento un enfoque diferente que ofrece una menor mejora en la precisión (aunque la diferencia de $\approx 2\%$ es lo suficientemente grande para ser notable, teniendo en cuenta cómo poco a poco el producto converge en un gran $N$).
Deje $\mathcal {P} (\mathbb{P}_2) $ representan el poder conjunto de todos los
doble de los números primos, $\mathcal {P} (\mathbb{P}{_2}(N)) $ el juego de poder de la primera $N$
doble de los números primos, y $\mathcal {P}_\kappa (\mathbb{P}{_2}(N)) $ el conjunto de los subconjuntos
de cardinalidad $\kappa.$ También vamos a $\mathcal {P}_\kappa \small{\left(\prod\frac{1}{p\in \mathbb{P}_2}\right)}$ representa el subconjunto de los productos de reciptocals de dos números primos en la que se especifica la cardinalidad.
Por ejemplo, $A=\{3,5,7,11\},\ \mathcal {P}_2 {\left(\prod\frac{1}{A}\right)}$ representaría el conjunto de $\left\{\frac{1}{15},\frac{1}{21},\frac{1}{33},\frac{1}{35},\frac{1}{55},\frac{1}{77}\right\}.$
Desde
\begin{align}
&\quad \prod_{p\in \mathbb{P}_2}^{N}\left(1-\frac{1}{p}\right)&=&\quad 1-\sum^{N} \left (\large\mathcal {P} _ {\text {odd}} \small{\left(\prod\frac{1}{p\in \mathbb{P}_2}\right)} - \large\mathcal {P} _ {\text {even}} \small{\left(\prod\frac{1}{p\in \mathbb{P}_2}\right)} \right) \\
\end{align}
como puede verse fácilmente en el caso de $N=3:$
\begin{align}
&\left(\frac{1}{3}+\frac{1}{5}+\frac{1}{7}-\frac{1}{15}-\frac{1}{21}-\frac{1}{35}+\frac{1}{105}\right)= 1-\left(1-\frac{1}{3}\right)\left(1-\frac{1}{5}\right)\left(1-\frac{1}{7}\right)\\
\end{align}
de ello se sigue que
\begin{align}
&\quad \prod_{p\in \mathbb{P}_2}\left(1-\frac{1}{p}\right)&=
&\quad \quad \sum _{p\in \mathbb{P}_2} \frac{1}{p}\\
&&&-\quad \frac{1}{2} \left(\sum _{p\in \mathbb{P}_2} \frac{1}{p}\right)^2-\frac{1}{2} \sum _{p\in \mathbb{P}_2} \frac{1}{p^2}\\
&&&+\quad \frac{1}{6} \left(\sum _{p\in \mathbb{P}_2} \frac{1}{p}\right)^3-\frac{1}{2} \left(\sum _{p\in \mathbb{P}_2} \frac{1}{p^2}\right) \sum _{p\in \mathbb{P}_2} \frac{1}{p}+\frac{1}{3} \sum _{p\in \mathbb{P}_2} \frac{1}{p^3}\\
&&&-\quad \dots
\end{align}
donde los coeficientes se muestran en la Tabla 24.2 en 1 (multinomials M2) multiplicado por el $(-1)^{q},$ donde $q$ es el número de elementos en el entero correspondiente partición.
Esta representación resulta ser beneficioso computacionalmente, ya que el poder de las sumas convergen rápidamente. $\infty$ en la suma puede ser reemplazado por un razonablemente pequeño $N,$ $\sum _{p\in \mathbb{P}_2} \frac{1}{p}$ reemplazado por Brun constante (vea la nota de abajo), para dar a la mejora ligeramente la aproximación de las $\approx 83.83 \%.$
Nota: Como Erick Wong notas en los comentarios de abajo, la actual "conocido" el valor de la constante de Brun se basa en un argumento Heurístico (Hardy Y Littlewood) que $\pi_2(x) \approx 2C_2 \int_2^x \frac{dt}{\left(\log t \right)^2},$ donde $C$ es el doble prime constante. Muy bien le da aquí una estimación de $B_2$ $1.9021605823\pm 8\times10^{-10}.$
- Milton Abramowitz y Irene A. Stegun (eds.), Manual de funciones matemáticas, 9ª ed.,
Publicaciones De Dover, Nueva York, 1972.