Aquí hay otra manera de demostrarlo. Sorprendentemente, para $n$ integral y $k$ real, la integral en cuestión puede ser escrito como una integral indefinida. Esto da una prueba directa para los no-entero $k$, aunque, obviamente, menos clara que el contorno de método. (De hecho, es conveniente evitar entero $k$ en este método, y que se extienden a entero $k$ por la continuidad.)
Escrito $y=x/n$ e $l=n-k$, tenemos por ejemplo para $n=2$:
$$\int\frac{\sin^2(2y)}{\sin^k(ky)\sin^l(ly)}dy=
\frac{\frac{2}{l-k}\sin((l-k)y)+\sin(2y)}
{kl\sin^{k-1}(ky)\sin^{l-1}(ly)}$$
En general, para $n$ a ($n$ impar es similar con los cosenos):
$$I_{n,k}(y)=\int\frac{\sin^n(ny)}{\sin^k(ky)\sin^l(ly)}dy=
\frac{\sum_{i=0}^{n-1}\sum_{s=0}^{n-1}\lambda_{r,s}\sin(((n-1-2r)k+(n-1-2s)l)y)}
{kl\sin^{k-1}(ky)\sin^{l-1}(ly)}$$
donde
$$\lambda_{r,s}=\begin{cases}
\displaystyle
\frac{(-1)^r(n-1)^{\underline{r}}}{(r-l)^{\underline{\smash{r-s}}}s!}\lambda_{0,0},\;\;r\ge s\\
\displaystyle
\frac{(-1)^s(n-1)^{\underline{s}}}{(s-k)^{\underline{\smash{s-r}}}r!}\lambda_{0,0},\;\;s\ge r
\end{casos}
$$
$$\lambda_{0,0}=(-1)^{n/2+1}2^{1-n},$$
y $x^{\underline{r}}$ denota el poder caer $x(x-1)\ldots(x-r+1)$.
Es fácil comprobar la derivada, $I_{n,k}'(y)$ es la correcta, teniendo en cuenta el coeficiente de
de $\cos((ak+bl)y)/(\sin^k(ky)\sin^l(ly))$ por cada $a$, $b$. Si $a\neq b$, a continuación, obtendrá cero,
de lo contrario, por $a=b=n-2r$, $(ak+bl)y=(n-2r)ny$ y usted consigue $\frac{1}{2}(-1)^r\binom{n}{r}\lambda_{0,0}\cos((n-2r)ny)$. A continuación, $\sin^n(ny)$ surge de la expansión binomial:
$$(-1)^{n/2}2^{-n}\sum_{r=0}^n (-1)^r\binom{n}{r}\cos((n-2r)ny)=\sin^n(ny).$$
Tenga en cuenta que $I_{n,k}(0)=0$ ya que, al ser la integral de algo de buen comportamiento en $0$,
$I_{n,k}(y)$ debe ser continua en $0$, por lo que su numerador debe desaparecer a fin de $n-2$
al igual que su denominador. El uso de L'H\^{o}pital, teniendo en $n-2$ derivados del numerador da sólo sines, que de por sí se desvanecen en 0.
Para evaluar $I_{n,k}(\pi/n)$, tenga en cuenta que $\sin(k\pi/n)=\sin(l\pi/n)$ y
$\sin(((n-1-2r)k+(n-1-2s)l)y)=\sin(2(r-s)k\pi/n)$.
Acondicionado en $r-s=d>0$, se obtiene (por ejemplo, mediante la consideración parcial de las fracciones en $k$)
$$\sum_{s=0}^{n-1-d}\lambda_{s+d,s}=
\frac{(-1)^d(n-1)!\binom{n-2}{d-1}}{(k-1)^{\underline{\smash{n-1}}}}\lambda_{0,0},$$
y lo mismo para $d<0$ con el signo opuesto, y el uso de $-d$ en lugar de $d$.
Utilizando el binomio de expansión de $(1-e^{2\pi ik/n})^{n-2}$, se obtiene
$$\sum_{d=1}^{n-1}(-1)^d\binom{n-2}{d-1}\sin\left(\frac{2dk\pi}{n}\right)=
(-1)^{n/2}2^{n-2}\sin^{n-2}\left(\frac{k\pi}{n}\right)\sin(k\pi)$$
Así que, finalmente, poner las piezas juntas,
$$\frac{n}{\pi}I_{n,k}\left(\frac{\pi}{n}\right)=\frac{n\sum_{r=0}^{n-1}\sum_{s=0}^{n-1}\lambda_{r,s}\sin\left(\frac{2(r-s)k\pi}{n}\right)}
{\pi kl\sin^{k-1}\left(\frac{k\pi}{n}\right)\sin^{l-1}\left(\frac{l\pi}{n}\right)}=\frac{n!\el pecado(k\pi}{\pi k^{\underline{\smash{n+1}}}}$$
que (para $n$) reconocemos como $\binom{n}{k}$ por el reflejo de la fórmula para factoriales.