Primero de todo, al darse cuenta de $\cot^2(\pi - x) = \cot^2(x)$, podemos escribir esta identidad como
$$\sum_{k=1,3,5} \cot^2(\frac{2\pi k}{14}) = 5$$
Por escritura $\cot^2(x) = \frac{1}{\sin^2(x)} - 1$, y el uso de simetrías de $\cos$ y $\sin$ ($\cos(x)=\cos(-x)$, $\sin(\frac{\pi}{2}-x)=\cos x$), podemos escribir esta suma como sigue:
$$\sum_{k=1,3,5} \frac{1}{\cos^2(\frac{\pi k}{14})} = 8$$
Si dejamos $a_i = \cos(\frac{\pi (2i-1)}{14}), i=1,2,3$, podemos escribir esta expresión como
$$(*) \frac{(\sum_{i<j} a_i a_j)^2 - 2\prod a_i \sum a_i}{(\prod a_i)^2}$$
El 7-esima Polinomio de Chebyshev (de primera clase) se desvanece exactamente en $\cos(\frac{2k-1}{14}\pi)$, $1\le k \le 7$. Esas raíces son en realidad $\pm (a_1, a_2, a_3)$$0$, cada una es una raíz simple.
Podemos calcular el polinomio de forma recursiva y encontrar que es igual a
$$T_7(x) = 64x^7-112x^5+56x^3-7x=x(64x^6-112x^4+56x^2-7)$$
Vamos a trabajar con la $P_7(x)=\frac{T_7(x)}{64x}$, un monic polinomio con raíces $\pm(a_1,a_2,a_3)$.
Esta muestra, mediante el uso de Vieta y la simetría de las raíces (requiere de la manipulación de polinomios simétricos):
- $(\prod a_i)^2=\frac{7}{64}$ (considerando el coeficiente de $x^0$)
- $(\sum_{i<j} a_i a_j)^2 - 2\prod a_i \sum a_i = \frac{56}{64}$ (considerando el coeficiente de $x^2$ - esto requiere algunos cálculos)
De modo que la suma de $(*)$ es igual a $\frac{56}{64} / \frac{7}{64} = 8$, lo que implica su identidad. $\blacksquare$
EDIT: voy a describir algunos de la filosofía detrás de la respuesta.
La primera mitad - yo sabía que quería que los polinomios de Chebyshev de alguna manera (porque sus raíces están relacionados con la expresión), así que hice básica manipulaciones que me ayudó a utilizar los coeficientes del polinomio de Chebyshev. Yo no sabía que apriori que no son 'buenas' manipulaciones, pero yo esperaba y de hecho funcionó.
La segunda mitad - Lo que realmente quería es un polinomio $Q(x)$ cuyas raíces se $a_1,a_2,a_3$. Por desgracia, yo había conseguido sólo para construir el polinomio $P_7(x)$ lo que equivale a $-Q(x)Q(-x)$. Afortunadamente, los coeficientes de $P_7$ codificar la información suficiente acerca de los coeficientes de $Q$. De forma explícita, mediante la comparación de los coeficientes:
$$P_7[X^k] = \sum_{i+j=k} (-1)^{1+j} Q[X^i]Q[X^j]$$
He utilizado este para $k=0,2$ y fue suficiente. $k=0$ dio $P_7(0)=-Q(0)^2$, es decir, tenemos el producto de la $a_i$! (hasta firmar, pero que no es necesaria).
$k=2$ dio $P_7[X^2] = Q[X^1]^2-2Q[X^2]Q[X^0]$, que por suerte era exactamente el ingrediente que falta en el cálculo de la expresión racional $(*)$, por lo que es.
EDIT 2: me siento que necesito para ampliar la "teoría" de polinomio de Chebyshev, porque el uso que podría asustar a la gente.
El $n$'th polinomio de Chebyshev de la primera clase es el único polinomio de satisfacciones $T_n(\cos (\theta)) = \cos(n\theta)$, para cualquier $\theta$. Evidentemente, $\cos(\frac{\pi}{2n}(2k+1))$ es una raíz para cualquier $k$ - sólo tiene que conectar $\theta = \frac{\pi}{2n}(2k+1))$. Como $\deg T_n = n$ (véase el párrafo siguiente), no puede ser de otra raíces.
¿Por qué es $T_n$ necesariamente un polinomio? Así, por $n=0$ tenemos $T_0 = 1$, y para $n=1$ tenemos $T_1(x)=x$. Para $n=2$ que ya necesita algo de trigonometría: $\cos(2\theta)=2\cos^2(\theta)-1$, lo $T_2(x)=2x^2-1$. Podemos definir a la $T_n$ recursivamente por trigonométricas ideas:
$$\cos(\alpha)+\cos(\beta)=2\cos(\frac{\alpha+\beta}{2})\cos(\frac{\alpha-\beta}{2})$$
$$\implies \cos((n+1)\theta) + \cos((n-1)\theta) = 2\cos(n\theta)\cos(\theta)$$
$$\implies T_{n+1}(x) + T_{n-1}(x) = 2T_{n}(x)x$$
Esta es la forma en que calcula $T_7$. En la práctica, que sólo utiliza la relación de recurrencia $T_{n+1}(x) = 2T_{n}(x)x-T_{n-1}$ y la tabla de aquí.
Hay algunos accesos directos, ya que el coeficiente inicial de $T_n$ $2^{n-1}$ y el último coeficiente es $0$ al $n$ es impar.