Simplificar $\prod_{k=3}^{n-1}\cos\left(\frac{\pi}{k}\right)$

Question

Busco simplificar lo siguiente, sin utilizar la notación Pi mayúscula: $$\prod_{k=3}^{n-1}\cos\left(\frac{\pi}{k}\right)$$ Que está destinado a producir la secuencia: $\left[1,\ \frac{1}{2},\ \frac{1}{2}\frac{\sqrt{2}}{2},\frac{1}{2}\frac{\sqrt{2}}{2}\frac{1+\sqrt{5}}{4},\frac{1}{2}\frac{\sqrt{2}}{2}\frac{1+\sqrt{5}}{4}\frac{\sqrt{3}}{2}...\right]$ .

He visto identidades de estructura similar, como: $$\prod_{k=1}^{n-1}\sin\left(\frac{k\pi}{n}\right) = \frac{n}{2^{n-1}},\qquad or\qquad \prod_{k=1}^{n-1}\cos\left(\frac{k\pi}{n}\right) = \frac{\sin(\frac{\pi n}{2})}{2^{n-1}}$$

Pero, no estoy versado ni en las pruebas de estas identidades, ni en las propiedades del $\Pi$ por lo que he tenido muchas dificultades para simplificarlo por mi cuenta. Tratar con $k$ en el denominador en vez de en el numerador (como en las dos identidades mencionadas) es algo que evidentemente no estoy preparado para tratar por mi cuenta. Gracias a quien esté dispuesto a ayudarme.

Answer 1

Como se ha dicho en los comentarios y en la respuesta, una forma cerrada es más que improbable.

Sin embargo, podemos hacer algunas aproximaciones. Sea $a_k=\cos \left(\frac{\pi }{k}\right)$ y considerar $$P_n=\prod_{k=3}^{n-1} a_k =\frac{1}{2\sqrt 2}\prod_{k=5}^{n-1} a_k$$ Tomar logaritmos $$ \log(P_n)=\log\left(\frac{1}{2\sqrt 2}\right)+\sum_{k=5}^{n-1} \log(a_k )$$ Ahora, utilizando la composición de las series de Taylor $$\log(a_k)=-\frac{\pi ^2}{2 k^2}-\frac{\pi ^4}{12 k^4}-\frac{\pi ^6}{45 k^6}-\frac{17 \pi ^8}{2520 k^8}-\frac{31 \pi ^{10}}{14175 k^{10}}-\frac{691 \pi ^{12}}{935550 k^{12}}+O\left(\frac{1}{k^{14}}\right)$$ Sólo por curiosidad $$\log(a_k)=-\sum_{p=1}^\infty\frac{2^{2 p-3} \pi ^{2 p} (E_{2 p-1}(1)-E_{2 p-1}(0))}{p\, (2 p-1)!\,k^{2p}}=\sum_{p=1}^\infty \frac {c_p}{k^{2p}}$$ donde aparecen los polinomios de Euler pero me temo que no podemos hacer mucho por la suma sobre $k$ además de utilizar el hecho de que $$\sum_{k=7}^{n-1}k^{-2p}=H_{n-1}^{(2 p)}-1-2^{-2 p}-3^{-2 p}-4^{-2 p}-5^{-2 p}-6^{-2 p}$$

Usando la expansión como se ha escrito anteriormente, deberíamos llegar a alguna expresión fea como $$\sum_{k=5}^{n-1} \log(a_k )=K-\frac{\pi ^2}{2} H_{n-1}^{(2)}-\frac{\pi ^4 }{12} H_{n-1}^{(4)}-\frac{\pi ^6}{45} H_{n-1}^{(6)}-\frac{17 \pi ^8}{2520} H_{n-1}^{(8)}-\frac{31 \pi ^{10} }{14175} H_{n-1}^{(10)}-\frac{691 \pi ^{12}}{935550} H_{n-1}^{(12)}$$ donde aparecen números armónicos genralizados y donde $$K=\frac{205 \pi ^2}{288}+\frac{22369 \pi ^4}{248832}+\frac{607493 \pi ^6}{26873856}+\frac{7339467953 \pi ^8}{1083553873920}+\frac{76853883151 \pi ^{10}}{35107145515008}+\frac{6162541526621059 \pi ^{12}}{8341457774365900800}$$

Probando con algunos valores de $n$ los valores decimales son $$\left( \begin{array}{ccc} n & \text{approximation} & \text{exact} \\ 3 & 0.99918368 & 1.00000000 \\ 4 & 0.49999441 & 0.50000000 \\ 5 & 0.35355339 & 0.35355339 \\ 6 & 0.28603083 & 0.28603070 \\ 7 & 0.24770997 & 0.24770985 \\ 8 & 0.22317897 & 0.22317887 \\ 9 & 0.20619049 & 0.20619039 \\ 10 & 0.19375568 & 0.19375558 \\ 11 & 0.18427260 & 0.18427251 \\ 12 & 0.17680827 & 0.17680818 \\ 13 & 0.17078367 & 0.17078359 \\ 14 & 0.16582101 & 0.16582093 \\ 15 & 0.16166353 & 0.16166345 \\ 16 & 0.15813079 & 0.15813072 \\ 17 & 0.15509235 & 0.15509228 \\ 18 & 0.15245161 & 0.15245154 \\ 19 & 0.15013553 & 0.15013546 \\ 20 & 0.14808788 & 0.14808780 \end{array} \right)$$

Siguiendo con la expansión anterior y la asíntota de los números armónicos, convirtiendo a decimales, deberíamos tener para grandes $n$ $$P_n \sim \frac{1}{2 \sqrt{2}} \exp \left(c+\frac{\pi ^2}{2 n} +\frac{\pi ^2}{4 n^2} \right)\qquad \text{with} \qquad c \approx -1.12361$$

Aplicado a $n=20$ esta última fórmula daría $0.148018$ .

Editar

Para el producto infinito

$$\prod_{k=3}^{\infty}\cos\left(\frac{\pi}{k}\right)=0.1149420448532962007010401576568126847536004314847$$ mientras que la apoximación dada conduciría a $0.11494164$ .

Actualización

Recordando que conocemos los valores exactos hasta $k=6$ utilizando como aproximación $$\cos(x) \simeq\frac{\pi ^2-4x^2}{\pi ^2+x^2}\qquad (-\frac \pi 2 \leq x\leq\frac \pi 2)\implies \cos\left(\frac{\pi}{k}\right)=\frac{k^2-4}{k^2+1}$$ llegamos a $$P_n=\prod_{k=3}^{n-1}\cos\left(\frac{\pi}{k}\right)=\frac{\sqrt{6}+\sqrt{30}}{32} \prod_{k=7}^{n-1}\cos\left(\frac{\pi}{k}\right)=\frac{40885 \left(\sqrt{6}+\sqrt{30}\right) \pi \text{csch}(\pi )}{774144} \frac{ \Gamma (n-2)\, \Gamma (n+2)}{ \Gamma (n-i) \,\Gamma (n+i)}$$ dando $$P_\infty=\frac{40885 \left(\sqrt{6}+\sqrt{30}\right) \pi \text{csch}(\pi )}{774144}\approx 0.113881$$

Utilizando en su lugar el $[2,2]$ Padé aproximante $$\cos\left(\frac{\pi}{k}\right)=\frac{12 k^2-5 \pi ^2}{12 k^2+\pi ^2}$$ deberíamos conseguir $$P_n=\frac{K_1}{K_2} \frac{\Gamma \left(n-\frac{1}{2} \sqrt{\frac{5}{3}} \pi \right) \Gamma \left(n+\frac{1}{2} \sqrt{\frac{5}{3}} \pi \right)}{\Gamma \left(n-\frac{i \pi }{2 \sqrt{3}}\right) \Gamma \left(n+\frac{i \pi }{2 \sqrt{3}}\right)}$$ donde $$K_1=\sqrt{5} \left(\sqrt{6}+\sqrt{30}\right) \left(12+\pi ^2\right) \left(48+\pi ^2\right) \left(108+\pi ^2\right) \left(192+\pi ^2\right) \left(300+\pi ^2\right) \left(432+\pi ^2\right)$$ $$ \sin \left(\frac{1}{2} \sqrt{\frac{5}{3}} \pi ^2\right) \text{csch}\left(\frac{\pi ^2}{2 \sqrt{3}}\right)$$ $$K_2=800 \left(\pi ^2-60\right) \left(5 \pi ^2-432\right) \left(5 \pi ^2-192\right) \left(5 \pi ^2-108\right) \left(5 \pi ^2-48\right) \left(5 \pi ^2-12\right)$$ dando $$P_\infty=\frac{K_1}{K_2}\approx 0.114938$$

Es normal que el aproximante de Padé conduzca a mejores resultados sincd $$\int_0^{\frac \pi 7} \left(\cos (x)-\frac{\pi ^2-4 x^2}{\pi ^2+x^2}\right)^2\,dx\approx 1.028 \times 10^{-7}$$ mientras que $$\int_0^{\frac \pi 7}\left(\cos (x)-\frac{12-5 x^2}{12+x^2}\right)^2\,dx\approx 9.574 \times 10^{-12}$$

Answer 2

Bueno, podrías escribir $$ \prod_{k=3}^n \cos(\pi/k) = 2^{2-n} \sum_{signs} \cos\left(\pm \frac{\pi}{3} \pm \frac{\pi}{4} \pm \ldots \pm \frac{\pi}{n}\right) $$ donde la suma es sobre todos los $2^{n-2}$ posibles opciones del $\pm$ signos.
Si $n$ es moderadamente grande, los $\pm \pi/3 \pm \ldots \pm \pi/n$ será bastante desagradable múltiplos racionales de $\pi$ . Así que no es exactamente una "simplificación".

Por ejemplo, si $n=7$ Recibo $$ 16^{-1} \left(\cos \left( {\frac {11\,\pi}{420}} \right) +\cos \left( {\frac {13\, \pi}{140}} \right) +\cos \left( {\frac {27\,\pi}{140}} \right) +\cos \left( {\frac {109\,\pi}{420}} \right) \\+\cos \left( {\frac {43\,\pi}{ 140}} \right) +\cos \left( {\frac {179\,\pi}{420}} \right) +\cos \left( {\frac {59\,\pi}{420}} \right) +\cos \left( {\frac {83\,\pi}{ 140}} \right) \\+\cos \left( {\frac {57\,\pi}{140}} \right) +\cos \left( {\frac {199\,\pi}{420}} \right) +\cos \left( {\frac {97\,\pi}{ 140}} \right) +\cos \left( {\frac {319\,\pi}{420}} \right)\\ +\cos \left( {\frac {31\,\pi}{420}} \right) +\cos \left( {\frac {113\,\pi}{ 140}} \right) +\cos \left( {\frac {151\,\pi}{420}} \right) +\cos \left( {\frac {153\,\pi}{140}} \right) \right)$$