OFFSET
1,3
COMMENTS
For odd prime numbers: Product_{odd primes p} 1/(1 - 1/p^2) = Pi^2/8 = (3/4)*zeta(2) = A111003.
For odd composite numbers: Product_{odd composite numbers c} 1/(1 - 1/c^2) = (81/80) * (225/224) * (441/440) * (625/624) * (729/728) * ... = 32/Pi^3, this constant.
EXAMPLE
1.032049101862383653901505686034038...
MATHEMATICA
RealDigits[32/Pi^3, 10, 100][[1]] (* Amiram Eldar, Jan 10 2020 *)
PROG
(PARI) p = 1.0; forstep(n = 3, 10^7, 2, if(!isprime(n), p*= (1 / (1 - 1 / n^2)))); print(p)
(PARI) 32/Pi^3
CROSSREFS
KEYWORD
nonn,cons
AUTHOR
Dimitris Valianatos, Jan 08 2020
STATUS
approved