La exponencial de la generación de la función de $J_{m,n}$ fijos $n$ es
$$G_n(t) = \sum_{m=0}^\infty J_{m,n} t^m/m! = \int_0^\pi e^{tx} \dfrac{\sin(nx)^2}{\sin(x)^2}\ dx$$
Luego tenemos la $(\sin(nx)/\sin(x))^2 = n + 2 \sum_{j=1}^{n-1} (n-j) \cos(2 j x)$
y $$G_n(t) = (e^{\pi t} - 1) \left(\frac{n}{t} + 2 \sum_{j=1}^{n-1} (n-j) \dfrac{t}{t^2 + 4 j^2}\right)$$
EDITAR:
Por ejemplo, $G_3(t) = (e^{\pi t} - 1) \left(\dfrac{3}{t} + \dfrac{4t}{t^2+4} + \dfrac{2t}{t^2+16}\right)$. Ahora $$e^{\pi t} - 1 = \sum_{j=1}^\infty \pi^j \dfrac{t^j}{j!} = \pi t + \dfrac{\pi^2 t^2}{2} + \dfrac{\pi^3 t^3}{6} + \ldots$$
$$\dfrac{3}{t} + \dfrac{4t}{t^2+4} + \dfrac{2t}{t^2+16} = \dfrac{3}{t} + \sum_{j=0}^\infty (-1)^j(4^{-j}+ 16^{-j}/8) t^{2j+1} = \dfrac{3}{t} + \left(1+\dfrac{1}{8}\right) t - \left(\dfrac{1}{4} + \dfrac{1}{16 \times 8}\right) t^3 + \ldots$$
y por lo $J_{m,3}$ $m!$ multiplicado por el coeficiente de $t^m$ en el producto de estos.