How can I solve the following equation for $q$? I'm totally stuck. I have done everything up to this point though.
$$\left(q + \sqrt{q^2-1}\right)^{2(N+1)} = 1,$$ where N is a natural number.
Attempt: Using DeMoivre's formula for the roots of unity, we have that for $$z^M = 1,$$the $M^{th}$ roots of unity are given by $$z=\exp\left(\frac{2ik\pi}{M}\right).$$ So $$q+\sqrt{q^2-1}=\cos\left(\frac{\pi k}{N+1}\right) + i\sin\left(\frac{\pi k}{N+1}\right), \ k=1,2,3, \ldots, 2N+1$$
What next?
If $q>1$, then we have clearly no solution.
If $q\in [-1,1]$, then there is a $\vartheta\in (0,2\pi)$, such that $$ \cos\vartheta=q $$ and $$ \sqrt{q^2-1}=\pm i\sqrt{1-q^2}= i\sin \vartheta. $$ (The minus sign is absorbed in the choice of $\vartheta$.) Then your equation looks like $$ (\cos\vartheta+i\sin\vartheta)^{2N+2}=1. $$ So in this case $$ q=\cos\left(\frac{2k\pi}{2N+2}\right),\,\,\,k=0,1,\ldots,2N+1. $$ Finally, if $q<-1$, then we simply want to study the cases $$ q+\sqrt{q^2-1}=1\quad \text{and}\quad q+\sqrt{q^2-1}=-1, $$ which are dealt with easily.