Showing that $\prod_{k=1}^{n} \left( 3 + 2\cos\left(\frac{2\pi}{n+1}k\right) \right)$ is the square of a Fibonacci number I was experimenting with products of the form
$$\prod_{k=1}^{n} \left( a + b\cos(ck) \right)$$
when I found that the expression
$$\prod_{k=1}^{n} \left( 3 + 2\cos\left(\frac{2\pi}{n+1}k\right) \right)$$
seems to always return a perfect square. (verified numerically).
$$\left(3 + 2\cos\left(\frac{2\pi}{2}\right)\right) = 1 = F_2^2$$
$$\left(3 + 2\cos\left(\frac{2\pi}{3}\right)\right) \left(3 + 2\cos\left(\frac{4\pi}{3}\right)\right) = 4 = F_3^2$$
$$\left(3 + 2\cos\left(\frac{2\pi}{4}\right)\right)\left(3 + 2\cos\left(\frac{4\pi}{4}\right)\right) \left(3 + 2\cos\left(\frac{6\pi}{4}\right)\right) = 9 = F_4^2$$
$$\dots$$
Furthermore, it appears these squares are the squares of Fibonacci numbers. ($F_1 = F_2 = 1$, $F_n = F_{n-1} + F_{n-2})$
I've tried to prove that this expression always gives a perfect square by considering the identity (Cassini's identity)
$$F_{n}^2 = F_{n+1}F_{n-1} + (-1)^{n-1}$$
and using induction to show that the above product satisfies this relation, and since the base cases ($F_1, F_2...$) are that of the Fibonacci sequence then this would prove that the expression is the square of Fibonacci numbers. However, this approach led nowhere.
Is there a clever way to show that the above expression is related to Fibonacci numbers?

This was a blast solving. Never in a million years would I have guessed these to be true. But the $n+1$ in the denominator really hints at what to use. I will be skipping the more well known proofs.
We start with the Chebyshev polynomials of the second kind $U_n(x)$ which satisfy the following recursion relations
$$U_0(x)=1,\ U_1(x)=2x,\ U_n(x)=2xU_{n-1}(x)-U_{n-2}(x)$$
Observe that $U_n(x)$ is a $n$-degree polynomial. Easy proof. It is true for $n=0,1$. Then use induction.
Claim 1 : $$U_n(\cos\theta)=\frac{\sin(n+1)\theta}{\sin\theta}$$
Then the roots of the polynomial $U_n(x)$ are exactly where the RHS of the above vanishes, namely $\frac{k\pi}{n+1}$, and you can produce $n$ of them.
Proposition 2 : $$U_n(x)=2^n\prod_{k=1}^n\left(x-\cos\left(\frac{k\pi}{n+1}\right)\right)$$
Proof : All that's left to prove is the fact that the leading coefficient of $U_n$ is $2^n$. It is true for $n=0,1$. By induction, $U_{n-1},\ U_{n-2}$ has leading coefficients $2^{n-1},\ 2^{n-2}$ respectively and thus $U_n$ has leading coefficient $2\cdot 2^{n-1}=2^n$. The result follows. $\ \square$

Let $F_n$ denote the $n$-th Fibonacci number, $n\ge 0$.
Theorem 3 : $$F_n=i^{-n}U_n\left(\frac{i}{2}\right)$$
Proof : Induction yet again. $$1=F_0=i^{-0}U_0\left(\frac{i}{2}\right)=1$$ $$1=F_1=i^{-1}U_1\left(\frac{i}{2}\right)=i^{-1}2\left(\frac{i}{2}\right)=1$$
Let this be true for all $0\le k\le n-1$. Then for $k=n$
\begin{align*}
F_n&=F_{n-1}+F_{n-2}\\&=i^{-n+1}U_{n-1}\left(\frac{i}{2}\right)+i^{-n+2}U_{n-2}\left(\frac{i}{2}\right)\\&=i^{-n+1}2\left(\frac{i}{2}\right)U_{n-1}\left(\frac{i}{2}\right)+i^{-n+2}U_{n-2}\left(\frac{i}{2}\right)\\&=i^{-n}\left\{2\left(\frac{i}{2}\right)U_{n-1}\left(\frac{i}{2}\right)+i^2U_{n-2}\left(\frac{i}{2}\right)\right\}\\&=i^{-n}\left\{2\left(\frac{i}{2}\right)U_{n-1}\left(\frac{i}{2}\right)-U_{n-2}\left(\frac{i}{2}\right)\right\}\\&=i^{-n}U_n\left(\frac{i}{2}\right)
\end{align*}
and we are done by induction. $\ \square$
Main Result : $$F_n^2=\prod_{k=1}^n\left( 3 + 2\cos\left(\frac{2\pi}{n+1}k\right) \right)$$
Proof : We first look at $U_{n}^2\left(\frac{i}{2}\right)$.
We have
\begin{align*}
U_n^2\left(\frac{i}{2}\right)&=2^{2n}\left[\prod_{k=1}^n\left(\left(\frac{i}{2}\right)-\cos\left(\frac{k\pi}{n+1}\right)\right)\right]^2\\&=\left[2^n\prod_{k=1}^n\left(\left(\frac{i}{2}\right)-\cos\left(\frac{k\pi}{n+1}\right)\right)\right]\left[2^n\prod_{k=1}^n\left(\left(\frac{i}{2}\right)-\cos\left(\frac{k\pi}{n+1}\right)\right)\right]\\&=\left[\prod_{k=1}^n\left(i-2\cos\left(\frac{k\pi}{n+1}\right)\right)\right]\left[\prod_{k=1}^n\left(i-2\cos\left(\frac{k\pi}{n+1}\right)\right)\right]
\end{align*}
Now to calculate this, pair up the $j$-th term of the first product with the $n+1-j$-th term of the second product. This produces terms like
$$\left(i-2\cos\left(\frac{j\pi}{n+1}\right)\right)\left(i-2\cos\left(\frac{(n+1-j)\pi}{n+1}\right)\right)=\left(i-2\cos\left(\frac{j\pi}{n+1}\right)\right)\left(i+2\cos\left(\frac{j\pi}{n+1}\right)\right)=-\left(1+4\cos^2\left(\frac{j\pi}{n+1}\right)\right)$$
Plugging this back, we have
$$U_n^2\left(\frac{i}{2}\right)=(-1)^n\prod_{k=1}^n\left(1+4\cos^2\left(\frac{k\pi}{n+1}\right)\right)$$
Using Theorem 3, we get
$$F_n^2=(-1)^{-n}\prod_{k=1}^n\left(1+4\cos^2\left(\frac{k\pi}{n+1}\right)\right)$$
Now using the standard cosine double angle identity, we get $$F_n^2=\prod_{k=1}^n\left(3+2\cos\left(\frac{2k\pi}{n+1}\right)\right)$$ and we are done. $\ \square$

Some random musings:
This seems generalizable in quite a few directions. What if we use $\sin$ instead of $\cos$? Do we still get something that has some recurrence? What if we used Chebyshev polynomials of the first kind? What kind of identities will that lead to?
More algebraically, does this suggest some deeper connection between cyclotomic extensions and quadratic extensions (Since Fibonacci numbers are related to $\mathbb Z[\sqrt 5]$ and we have a bunch of products of cosines)? [Kronecker-Weber is a step in this direction, but I am not sure how it related honestly]. Or more intriguing, is there a relation between cyclotomic extensions and Chebyshev polynomials?
I am not learned enough to even try and make a guess. Maybe people in the comment can guide me?

I just now saw the comment by Thomas Andrews mentioning Chebyshev polynomials. Interesting how we came up with the same idea. The $n+1$ denominator suggested this approach. I would be interested in learning how they came up with it though.
